Research Article |
Corresponding author: Jovana M. Jasso-Martínez ( jovana.jasso@gmail.com ) Academic editor: Jose Fernandez-Triana
© 2023 Jovana M. Jasso-Martínez, Seán G. Brady, Robert R. Kula.
This is an open access article distributed under the terms of the CC0 Public Domain Dedication.
Citation:
Jasso-Martínez JM, Brady SG, Kula RR (2023) Phylogenetic affinities of the non-cyclostome subfamilies Amicrocentrinae and Dirrhopinae (Hymenoptera, Braconidae) confirmed by ultraconserved element data. Journal of Hymenoptera Research 96: 1017-1030. https://doi.org/10.3897/jhr.96.111012
|
The subfamilies Amicrocentrinae and Dirrhopinae (Hymenoptera, Braconidae) are two small, monogeneric braconid subfamilies whose species exclusively attack lepidopteran larvae. The phylogenetic placement of Amicrocentrinae as a member of the helconoid complex of subfamilies has been supported by morphological and nuclear Sanger sequence data. The subfamilial status of Dirrhopinae on the other hand has been subject to debate, although it has been suggested as closely related to the microgastroid complex based on morphology only. Here we generated for the first time genomic ultraconserved element data for members of the above subfamilies (Amicrocentrum seyrigi van Achterberg and Dirrhope americana Muesebeck) to assess their phylogenetic affinities using exhaustive taxon sampling that includes all but one of the currently valid braconid subfamilies. Our results strongly confirm the placement of both taxa within the non-cyclostome helconoid and microgastroid complexes, respectively.
Amicrocentrinae, Braconidae, Dirrhopinae, non-cyclostomes, ultraconserved elements
The parasitoid wasp family Braconidae is an extremely species-rich group within the order Hymenoptera, currently having more than 21,000 described species (
As for other megadiverse groups, the phylogenetic affinities of various taxa within Braconidae have been debated extensively due to the lack of morphological synapomorphies, limited taxon sampling and limited molecular information (
The subfamily Amicrocentrinae was established by
The placement within Braconidae of the rare, lepidopteran parasitoid genus Dirrhope Foerster has changed through time, with it suggested as a member of Microgastrinae (
In this study we included representatives of Amicrocentrinae and Dirrhopinae for the first time in a phylogenomic dataset that includes all currently extant valid braconid subfamilies except Xiphozelinae. Based on genomic-scale data obtained from UCEs, we evaluated both their proposed subfamilial status and their phylogenetic affinities within Braconidae. Additionally, we assessed previously proposed morphological apomorphic features for these taxa in relation to their recovered phylogenetic relationships.
Our taxon sampling comprises a total of 401 terminal taxa, including outgroups. Our braconid ingroup includes 233 species that belong to the cyclostome s.l. group, 156 species of the non-cyclostome group (including Amicrocentrum seyrigi and Dirrhope americana) and Apozyx penyai Mason (Apozyginae). All braconid taxa used in this study, except A. seyrigi and D. americana, are part of the dataset that was included in a recently published phylogenetic study of Braconidae (table S1 in
Locality and SRA information of Amicrocentrum seyrigi and Dirrhope americana.
Specimen voucher | Subfamily | Genus | Species | Label information | SRA accession number |
---|---|---|---|---|---|
USNMENT1322780 | Amicrocentrinae | Amicrocentrum | seyrigi | Madagascar: Antsokay. Latitude: 23°25.40'S, longitude: 43°44.51'E. Collection year: 2012. Collectors: M. Irwin, Rin’ha. | SAMN37185746 |
USNMENT1322781 | Dirrhopinae | Dirrhope | americana | USA: Virginia, Fauquier Co. Granruth Grove. Latitude: 38°49'9.60"N, longitude: 77°54'35.60"W. Collection year: 2014. Collector: Kula, R., et al. | SAMN37185747 |
Terminology for external morphology, including wing venation, follows
Genomic DNA extraction of A. seyrigi and D. americana was performed using a leg and whole specimen, respectively, with the Quiagen DNeasy Blood and Tissue kit (Qiagen Inc.,Valencia, California, U.S.A.). DNA quantitation was measured with Qubit 2.0 fluorometer (Life Technologies Inc. Carlsbad, CA). Total DNA of both samples (A. seyrigi: 41 ng; D. americana 16 ng) was used as input for shearing and library preparation. DNA was sheared for 45 seconds in a Qsonica Q800R sonicator (Qsonica LLC, Newton, CT). Genomic libraries were prepared using the Kapa Hyper Prep Kit (Roche) and the custom Tru-Seq-style dual-indexing adapter (
Raw reads were cleaned and trimmed using Illumiprocessor v 2.0.7 (Faircloth, 2013), a wrapper around Trimmomatic (
Both aligned matrices were partitioned using the Sliding-Window Site Characteristics Entropy (SWSC-EN) algorithm (
We obtained a total of 2,526 UCE loci in a complete dataset. We recovered 1,161 UCE loci with a mean length of 762.25 bp for the included species of A. seyrigi, whereas for D. americana we obtained 828 UCE loci with a mean length of 442.22 bp. A total of 1,578 and 917 UCE loci were retained for the 25% and 50% matrix, resulting in matrices with lengths of 342,464 and 170,549 bp, respectively.
The ML phylograms derived from the two different datasets analyzed (50% and 25%; Fig.
Within the helconoid complex both topologies recovered Amicrocentrinae as sister to Charmontinae with high support (BTP=100), with both subfamilies being sister to Macrocentrinae. These three subfamilies together were recovered as sister to a clade comprised by the subfamilies Homolobinae, Microtypinae and Orgilinae, followed by a clade containing the subfamilies Brachistinae, Helconinae and Acampsohelconinae (Fig.
The only topological difference at the subfamily level between both datasets analyzed was the placement of the cyclostome subfamily Pambolinae, which was recovered in a clade containing “Old World” Doryctinae with weak support (BTP = 60) or as sister to the Avgiini with moderate support (BTP = 86) in the 25% (Suppl. material
In its current status, the non-cyclostome helconoid complex contains eight subfamilies, Acampsohelconinae, Brachistinae, Charmontinae, Helconinae, Homolobinae, Macrocentrinae, Microtypinae and Orgilinae (
Amicrocentrum, the only genus of Amicrocentrinae, was previously considered a member of Macrocentrinae; however, species of Amicrocentrum lack spines on the hind trochantellus (Fig.
In this study, we confirm with strong support the phylogenetic placement of Amicrocentrinae as sister to Charmontinae, with both being sister to Macrocentrinae within the helconoid complex based on genomic UCE data, as was previously proposed based on nuclear Sanger sequence data (
Some lineages of the braconid non-cyclostome group such as Cardiochilinae and Miracinae have been treated as members of Microgastrinae (
The limits of the microgastroid complex have remained relatively stable since the publication of a phylogenetic hypothesis based on morphology, together with the recognition of the subfamilies Khoikhoiinae and Mendesellinae (
Our genomic-scale data recovered Dirrhopinae as sister to Cheloninae within the microgastroid complex of subfamilies with strong support. Dirrhope species resemble members of the chelonine tribe Adeliini in body shape; however, in addition to the spiracle synapomorphy for dirrhopines, they have forewing vein RS distinctly curved anteriorly (Fig.
Maximum likelihood phylogram resulting from the 50% complete matrix. Amicrocentrinae was recovered as sister to Charmontinae within the non-cyclostome helconoid complex (red). Dirrhopinae was recovered as sister to Cheloninae in the microgastroid complex (orange). Remaining Braconidae and outgroup taxa are in black. Numbers near nodes are BTP values that were under 100.
We thank the staff from the Smithsonian Institution Laboratories of Analytical Biology (SI-LAB) at the Smithsonian National Museum of Natural History for their general support. Thanks to Jeffrey Sosa-Calvo for his help and advice with lab work; Taina Litwak for the illustrations included in Figs
Maximum likelihood tree resulting from the 50% complete matrix without collapsed clades and with full names of the terminal taxa
Data type: jpg
Maximum likelihood tree resulting from the 25% complete matrix
Data type: jpg