Phylogenetic signal in the COI, 16S, and 28S genes for inferring relationships among genera of Microgastrinae (Hymenoptera; Braconidae): evidence of a high diversification rate in this group of parasitoids

[From GLP1 receptor agonists to triple hormone receptor activation supplemented with glucagon receptor agonism.]

Following the introduction of mono- and then dual hormone (incretin) receptor agonists into therapy, attention was turned to multiple receptor stimulation, with the additional activation of the glucagon receptor, as a new option for the pharmaceutical treatment of type 2 diabetes and obesity. In addition to its role in carbohydrate metabolism, the article reviews the other important physiological tasks of glucagon, especially its participation in intrainsular paracrine regulation, energy expenditure and the shaping of appetite and food consumption. It covers the potential benefits of the triple combination and briefly touches data on the efficacy and safety of the first triple receptor agonist drug, retatrutide, in preclinical human studies. Further confirmation of the promising results may represent progress in the treatment of these forms of disease and their accompanying conditions, such as steatosis hepatis. Orv Hetil. 2023; 164(42): 1656-1664.

Morphology and evolution of the mesopleuron in Bethylidae (Hymenoptera: Chrysidoidea) mapped on a molecular phylogeny

The mesopleuron of Bethylidae has morphostructural characters that remain poorly understood, explored, and defined. The wide range of variability of this sclerite has generated confusion both in taxonomic and cladistic studies. Aiming to solve this issue, we describe the general mesopleural anatomy of Bethylidae. Our goal is to propose primary homologies by matching external and internal structures (muscles and apodemes). We reconstruct the ancestral state of the main mesopleural structures by applying the maximum-likelihood method on ten selected character states. The phylogeny of Bethylidae is built by analyzing a dataset of COI and 28S sequences, using maximum-likelihood and Bayesian inference methods. The Bethylidae and all subfamilies are recovered as monophyletic in both resulting phylogenies, with high clade support values. Although the two analyses yielded similar results, we used the tree resulting from the Bayesian inference to map the evolution of the morphological characters, as it is better supported. The study of mesopleural anatomy allows exploration and discussion of the evolution of characters and their present states in Bethylidae and Hymenoptera in general. Reconstruction of the ancestral states shows that many characters arose independently in Bethylidae subfamilies.

A five-gene molecular phylogeny reveals Parapanteles Ashmead (Hymenoptera: Braconidae) to be polyphyletic as currently composed

Parapanteles Ashmead (Braconidae: Microgastrinae) is a medium-sized genus of microgastrine wasps that was erected over a century ago and lacks a unique synapomorphic character, and its monophyly has not been tested by any means. Parapanteles usually are parasitoids of large, unconcealed caterpillars (macrolepidoptera) and have been reared from an unusually large diversity of hosts for a relatively small microgastrine genus. We used Cytochrome Oxidase I sequences ("DNA barcodes") available for Parapanteles and other microgastrines to sample the generic diversity of described and undescribed species currently placed in Parapanteles, and then sequenced four additional genes for this subsample (wingless, elongation factor 1-alpha, ribosomal subunit 28s, and NADH dehydrogenase subunit 1). We constructed individual gene trees and concatenated Bayesian and maximum-likelihood phylogenies for this 5-gene subsample. In these phylogenies, most Parapanteles species formed a monophyletic clade within another genus, Dolichogenidea, while the remaining Parapanteles species were recovered polyphyletically within several other genera. The latter likely represent misidentified members of other morphologically similar genera. Species in the monophyletic clade containing most Parapanteles parasitized caterpillars from only five families - Erebidae (Arctiinae), Geometridae, Saturniidae, Notodontidae, and Crambidae. We do not make any formal taxonomic decisions here because we were not able to include representatives of type species for Parapanteles or other relevant genera, and because we feel such decisions should be reserved until a comprehensive morphological analysis of the boundaries of these genera is accomplished.

The evolutionary history of the cellophane bee genus Colletes Latreille (Hymenoptera: Colletidae): Molecular phylogeny, biogeography and implications for a global infrageneric classification

Colletes Latreille (Hymenoptera: Colletidae) is a diverse genus with 518 valid species distributed in all biogeographic realms, except Australasia and Antarctica. Here we provide a comprehensive dated phylogeny for Colletes based on Bayesian and maximum likelihood-based analyses of DNA sequence data of six loci: 28S rDNA, cytochrome c oxidase subunit 1, elongation factor-1α copy F2, long-wavelength rhodopsin, RNA polymerase II and wingless. In total, our multilocus matrix consists of 4824 aligned base pairs for 143 species, including 112 Colletes species plus 31 outgroups (one stenotritid and a diverse array of colletids representing all subfamilies). Overall, analyses of each of the six single-locus datasets resulted in poorly resolved consensus trees with conflicting phylogenetic signal. However, our analyses of the multilocus matrix provided strong support for the monophyly of Colletes and show that it can be subdivided into five major clades. The implications of our phylogenetic results for future attempts at infrageneric classification for the Colletes of the world are discussed. We propose species groups for the Neotropical species of Colletes, the only major biogeographic realm for which no species groups have been proposed to date. Our dating analysis indicated that Colletes diverged from its sister taxon, Hemicotelles Toro and Cabezas, in the early Oligocene and that its extant lineages began diversifying only in the late Oligocene. According to our biogeographic reconstruction, Colletes originated in the Neotropics (most likely within South America) and then spread to the Nearctic very early in its evolutionary history. Geodispersal to the Old World occurred soon after colonization of the Northern Hemisphere. Lastly, the historical biogeography of Colletes is analyzed in light of available geological and palaeoenvironmental data.

Molecular sequencing and morphological identification reveal similar patterns in native bee communities across public and private grasslands of eastern North Dakota

Bees play a key role in the functioning of human-modified and natural ecosystems by pollinating agricultural crops and wild plant communities. Global pollinator conservation efforts need large-scale and long-term monitoring to detect changes in species’ demographic patterns and shifts in bee community structure. The objective of this project was to test a molecular sequencing pipeline that would utilize a commonly used locus, produce accurate and precise identifications consistent with morphological identifications, and generate data that are both qualitative and quantitative. We applied this amplicon sequencing pipeline to native bee communities sampled across Conservation Reserve Program (CRP) lands and native grasslands in eastern North Dakota. We found the 28S LSU locus to be more capable of discriminating between species than the 18S SSU rRNA locus, and in some cases even resolved instances of cryptic species or morphologically ambiguous species complexes. Overall, we found the amplicon sequencing method to be a qualitatively accurate representation of the sampled bee community richness and species identity, especially when a well-curated database of known 28S LSU sequences is available. Both morphological identification and molecular sequencing revealed similar patterns in native bee community structure across CRP lands and native prairie. Additionally, a genetic algorithm approach to compute taxon-specific correction factors using a small subset of the most concordant samples demonstrated that a high level of quantitative accuracy could be possible if the specimens are fresh and processed soon after collection. Here we provide a first step to a molecular pipeline for identifying insect pollinator communities. This tool should prove useful for future national monitoring efforts as use of molecular tools becomes more affordable and as numbers of 28S LSU sequences for pollinator species increase in publicly-available databases.

Patterns and drivers of genetic diversity and structure in the biological control parasitoid Habrobracon hebetor in Niger

When a promising natural enemy of a key pest exists locally, it is a common practice in biological control (BC) to rear and release it for supplementary control in the targeted agroecosystem even though significant knowledge gaps concerning pre/post release may still exist. Incorporating genetic information into BC research fills some of these gaps. Habrobracon hebetor, a parasitoid of many economically important moths that infest stored and field crops worldwide is commonly used, particularly against the millet head miner (MHM), a key pest of millet in Sahelian countries. To advance our knowledge on how H. hebetor that occurs naturally in open-field cropping systems and grain stores as well as being mass-produced and released for MHM control, performs in millet agroecosystems in Niger we evaluated its population genetics using two mitochondrial and 21 microsatellite markers. The field samples were genetically more diverse and displayed heterozygote excess. Very few field samples had faced significant recent demographic bottlenecks. The mating system (i.e. nonrandom mating with complementary sex determination) of this species may be the major driver of these findings rather than bottlenecks caused by the small number of individuals released and the scarcity of hosts during the longlasting dry season in Niger. H. hebetor population structure was represented by several small patches and genetically distinct individuals. Gene flow occurred at local and regional scales through human-mediated and natural short-distance dispersal. These findings highlight the importance of the mating system in the genetic diversity and structure of H. hebetor populations, and contribute to our understanding of its reported efficacy against MHM in pearl millet fields.

Phylogeny and biogeography of the cleptoparasitic bee genus Epeolus (Hymenoptera: Apidae) and cophylogenetic analysis with its host bee genus Colletes (Hymenoptera: Colletidae)

The bee genus Epeolus Latreille (Hymenoptera: Apidae) consists of 109 species, which are known to be exclusively cleptoparasites of polyester (or cellophane) bees of the genus Colletes Latreille (Hymenoptera: Colletidae). Both genera have a nearly cosmopolitan distribution and are represented on all continents except Antarctica and Australia. We present the most comprehensive phylogeny for Epeolus to date, based on combined molecular and morphological data. In total, 59 ingroup taxa (species of Epeolus) and 7 outgroup taxa (other Epeolini) were scored for 99 morphological characters, and sequence data were obtained for seven genes (one mitochondrial and six nuclear, 5399 bp in total). Epeolus was found to be monophyletic, with a crown age estimated to be 25.0-13.4 Ma (95% HPD) and its origins traced to the Nearctic region. Epeolus was found to contain six major clades, five of which were well supported. The evolutionary history of Epeolus is explored in the context of earth history events and the evolutionary history of its host genus Colletes, for which a molecular phylogeny was constructed based on the same seven genes. A comparison of Epeolus and Colletes phylogenies limited to taxa for which there is evidence of an association suggests there was some cospeciation. However, more cladogenetic events in Epeolus were linked to instances of dispersal/vicariance. It is not yet clear the extent to which allopatric speciation contributed to diversification in Colletes, but the genus' success in having colonized and diversified across much of the globe made it possible for Epeolus to do the same.

Taxonomic integrative and phylogenetic identification of the first recorded Triatoma rubrofasciata in Zhangzhou, Fujian Province and Maoming, Guangdong Province, China

BACKGROUND

Most species of Triatominae live exclusively in Latin America. However, one species, Triatoma rubrofasciata, has been recorded in the Americas as well as in various port areas in Africa and Asia. An increasing number of T. rubrofasciata have been reported in southern China in recent years. However, the origin of this invasive insect vector in China remains unknown, therefore, accurate identification and phylogenetic analysis of the bugs are urgently needed.

METHODS

A total of seven triatomine insect specimens were found and collected from Maoming City, Guangdong Province, China (GDMM) and Zhangzhou City, Fujian Province, China (FJZZ), respectively. The obtained insect vector specimens were observed under a dissecting microscope for morphological classification and then the genomic DNA was extracted, and the 16S ribosomal RNA (rRNA), 28S rRNA as well as cytochrome oxidase subunit I (COI) genes of the species were amplified and sequenced. Subsequently, molecular phylogenetic analyses based on multiple alignments of the above genes were conducted in order to identify the species and determine the phylogenetic origin approximation accurately.

RESULTS

The triatomine insects collected from GDMM and FJZZ were identified as Triatoma rubrofasciata using morphological and genetic analyses. All of the Chinese T. rubrofasciata captured in FJZZ, GDMM and other localities in southern China, together with a Vietnamese and Brazilian strain, formed a new, cohesive clade. T. rubrofasciata in GDMM and FJZZ are likely derived from strains found in Vietnam or Brazil.

CONCLUSIONS

To the best of our knowledge, this is the first record of the invasive insect T. rubrofasciata, which is likely derived from strains native to Vietnam or Brazil, in both Maoming City, Guangdong Province and Zhangzhou City, Fujian Province of China. A comparison of the DNA sequences of the 16 s rRNA, 28 s rRNA and COI genes confirmed the specific identification of T. rubrofasciata, and its potential origin in China is based on the phylogenetic analyses undertaken in this study. More targeted interventions and improved entomological surveillance are urgently needed to control the spread of this haematophagous insect in China.

Multi-gene phylogeny and divergence estimations for Evaniidae (Hymenoptera)

Ensign wasps (Hymenoptera: Evaniidae) develop as predators of cockroach eggs (Blattodea), have a wide distribution and exhibit numerous interesting biological phenomena. The taxonomy of this lineage has been the subject of several recent, intensive efforts, but the lineage lacked a robust phylogeny. In this paper we present a new phylogeny, based on increased taxonomic sampling and data from six molecular markers (mitochondrial 16S and COI, and nuclear markers 28S, RPS23, CAD, and AM2), the latter used for the first time in phylogenetic reconstruction. Our intent is to provide a robust phylogeny that will stabilize and facilitate revision of the higher-level classification. We also show the continued utility of molecular motifs, especially the presence of an intron in the RPS23 fragments of certain taxa, to diagnose evaniid clades and assist with taxonomic classification. Furthermore, we estimate divergence times among evaniid lineages for the first time, using multiple fossil calibrations. Evaniidae radiated primarily in the Early Cretaceous (134.1-141.1 Mya), with and most extant genera diverging near the K-T boundary. The estimated phylogeny reveals a more robust topology than previous efforts, with the recovery of more monophyletic taxa and better higher-level resolution. The results facilitate a change in ensign wasp taxonomy, with Parevania, and Papatuka, syn. nov. becoming junior synonyms of Zeuxevania, and Acanthinevania, syn. nov. being designated as junior synonym of Szepligetella. We transfer 30 species to Zeuxevania, either reestablishing past combinations or as new combinations. We also transfer 20 species from Acanthinevania to Szepligetella as new combinations.

Systematics, Biology, and Evolution of Microgastrine Parasitoid Wasps

The braconid parasitoid wasp subfamily Microgastrinae is perhaps the most species-rich subfamily of animals on Earth. Despite their small size, they are familiar to agriculturalists and field ecologists alike as one of the principal groups of natural enemies of caterpillars feeding on plants. Their abundance and nearly ubiquitous terrestrial distribution, their intricate interactions with host insects, and their historical association with mutualistic polydnaviruses have all contributed to Microgastrinae becoming a key group of organisms for studying parasitism, parasitoid genomics, and mating biology. However, these rich sources of data have not yet led to a robust genus-level classification of the group, and some taxonomic confusion persists as a result. We present the current status of understanding of the general biology, taxonomic history, diversity, geographical patterns, host relationships, and phylogeny of Microgastrinae as a stimulus and foundation for further study. Current progress in elucidating the biology and taxonomy of this important group is rapid and promises a revolution in the classification of these wasps in the near future.

Infochemical use and dietary specialization in parasitoids: a meta-analysis

Many parasitoid species use olfactory cues to locate their hosts. In tritrophic systems, parasitoids of herbivores can exploit the chemical blends emitted by plants in reaction to herbivore-induced damage, known as herbivore-induced plant volatiles (HIPVs). In this study, we explored the specificity and innateness of parasitoid responses to HIPVs using a meta-analysis of data from the literature. Based on the concept of dietary specialization and infochemical use, we hypothesized that (i) specialist parasitoids (i.e., with narrow host ranges) should be attracted to specific HIPV signals, whereas generalist parasitoids (i.e., with broad host ranges) should be attracted to more generic HIPV signals and (ii) specialist parasitoids should innately respond to HIPVs, whereas generalist parasitoids should have to learn to associate HIPVs with host presence. We characterized the responses of 66 parasitoid species based on published studies of parasitoid behavior. Our meta-analysis showed that (i) as predicted, specialist parasitoids were attracted to more specific signals than were generalist parasitoids but, (ii) contrary to expectations, response innateness depended on a parasitoid's target host life stage rather than on its degree of host specialization: parasitoids of larvae were more likely to show an innate response to HIPVs than were parasitoids of adults. This result changes our understanding of dietary specialization and highlights the need for further theoretical research that will help clarify infochemical use by parasitoids.

Phylogenetic Relationships Among Malaria Vectors and Closely Related Species in Thailand Using Multilocus DNA Sequences

The evolutionary and taxonomic status is important for understanding speciation events and phylogenetic relationships between closely related vector and nonvector species. This information is useful for targeting important disease vector species groups for the development of novel genetic-based vector and pathogen control methods. In this study, different phylogenetic analyses were performed to reconstruct phylogenetic trees for the primary malaria vectors in Thailand based on sequence information of 4 DNA fragments from the nuclear and mitochondrial regions. The primary Anopheles species in the subgenus Cellia involved in malaria transmission in Thailand separate clearly into 3 distinct clades: the Leucosphyrus group, Minimus subgroup, and Maculatus group. The phylogenetic trees based on different reconstructed algorithms and different gene regions provided congruent phylogenetic status of the mosquito species studied. The phylogenetic relationships of malaria vector species examined followed similar patterns based on morphological characters. An estimate of the divergence time among the Anopheles species infers that they were present during the Eocene and Miocene periods (>41 million years ago). Congruent phylogenetic analysis of malaria vectors is presented with different algorithms and gene regions. The nuclear TOLL6 fragment appears useful for molecular phylogenetic, species DNA barcode, and Anopheles population genetic analyses.

Molecular and morphological phylogenetics of chelonine parasitoid wasps (Hymenoptera: Braconidae), with a critical assessment of divergence time estimations

Parasitoid wasps of the subfamily Cheloninae are both species rich and poorly known. Although the taxonomy of Cheloninae appears to be relatively stable, there is no clear understanding of relationships among higher-level taxa. We here applied molecular phylogenetic analyses using three markers (COI, EF1α, 28S) and 37 morphological characters to elucidate the evolution and systematics of these wasps. Analyses were based on 83 specimens representing 13 genera. All genera except Ascogaster, Phanerotoma, and Pseudophanerotoma formed monophyletic groups; Furcidentia (stat. rev.) is raised to generic rank. Neither Chelonus (Chelonus) nor Chelonus (Microchelonus) were recovered as monophyletic, but together formed a monophyletic lineage. The tribes Chelonini and Odontosphaeropygini formed monophyletic groups, but the Phanerotomini sensu Zettel and Pseudophanerotomini were retrieved as either para- or polyphyletic. The genera comprising the former subfamily Adeliinae were confirmed as being nested within the Cheloninae. To estimate the age of the subfamily, we used 16 fossil taxa. Three approaches were compared: fixed-rate dating, node dating, and total-evidence dating, with age estimates differing greatly between the three methods. Shortcomings of each approach in relation to our dataset are discussed, and none of the age estimates is deemed sufficiently reliable. Given that most dating studies use a single method only, in most cases without presenting analyses on the sensitivity to priors, it is likely that numerous age estimates in the literature suffer from a similar lack of robustness. We argue for a more rigorous approach to dating analyses and for a faithful presentation of uncertainties in divergence time estimates. Given the results of the phylogenetic analysis the following taxonomic changes are proposed: Furcidentia Zettel (stat. rev.), previously treated as a subgenus of Pseudophanerotoma Zettel is raised to generic rank; Microchelonus Szépligeti (syn. nov.), variously treated by previous authors, is proposed as a junior synonym of Chelonus Jurine; the following subgenera of Microchelonus - Baculonus Braet & van Achterberg (syn. nov.), Carinichelonus Tobias (syn. nov.) and Scabrichelonus He, Chen & van Achterberg (syn. nov.), are proposed as junior synonyms of Chelonus; a number of new species names are proposed due to homonyms resulting from the above changes and these are listed in the paper.

Patterns of Protein Evolution in Cytochrome c Oxidase 1 (COI) from the Class Arachnida

Because sequence information is now available for the 648bp barcode region of cytochrome c oxidase 1 (COI) from more than 400,000 animal species, this gene segment can be used to probe patterns of mitochondrial evolution. The present study examines levels of amino acid substitution and the frequency of indels in COI from 4177 species of arachnids, including representatives from all 16 orders and 43% of its families (267/625). It examines divergences at three taxonomic levels—among members of each order to an outgroup, among families in each order and among BINs, a species proxy, in each family. Order Distances vary fourfold (0.10–0.39), while the mean of the Family Distances for the ten orders ranges fivefold (0.07–0.35). BIN Distances show great variation, ranging from 0.01 or less in 12 families to more than 0.25 in eight families. Patterns of amino acid substitution in COI are generally congruent with previously reported variation in nucleotide substitution rates in arachnids, but provide some new insights, such as clear rate acceleration in the Opiliones. By revealing a strong association between elevated rates of nucleotide and amino acid substitution, this study builds evidence for the selective importance of the rate variation among arachnid lineages. Moreover, it establishes that groups whose COI genes have elevated levels of amino acid substitution also regularly possess indels, a dramatic form of protein reconfiguration. Overall, this study suggests that the mitochondrial genome of some arachnid groups is dynamic with high rates of amino acid substitution and frequent indels, while it is ‘locked down’ in others. Dynamic genomes are most prevalent in arachnids with short generation times, but the possible impact of breeding system deserves investigation since many of the rapidly evolving lineages reproduce by haplodiploidy, a mode of reproduction absent in ‘locked down’ taxa.

Cladistic analysis of self-grooming indicates a single origin of eusociality in corbiculate bees (Hymenoptera: Apidae)

Behavioural traits have been used extensively in recent years as an important character source for making phylogenetic inferences. The phylogenetic positions of the members of the Apini subtribe are increasingly being debated, and new characters must be examined. We analysed the presence and absence of certain behavioural patterns, as well as the sequences of some of these patterns, to generate 79 characters. Eleven species comprised the ingroup, and Xylocopini comprised the outgroup. Parsimony analysis showed that the most parsimonious tree was (Euglossina(Bombina(Apina+Meliponina))). This topology is consistent with most studies that use morphological data and the few that use behavioural data, which suggests that advanced eusociality arose only once in a common ancestor of the clade Apina plus Meliponina; however, this hypothesis is inconsistent with our molecular data. Thus we considered behavioural, molecular, and morphological data and recovered the same topology, in which eusociality has a single origin in corbiculate bees.

Molecular identification of sibling species of Sclerodermus (Hymenoptera: Bethylidae) that parasitize buprestid and cerambycid beetles by using partial sequences of mitochondrial DNA cytochrome oxidase subunit 1 and 28S ribosomal RNA gene

The species belonging to Sclerodermus (Hymenoptera: Bethylidae) are currently the most important insect natural enemies of wood borer pests, mainly buprestid and cerambycid beetles, in China. However, some sibling species of this genus are very difficult to distinguish because of their similar morphological features. To address this issue, we conducted phylogenetic and genetic analyses of cytochrome oxidase subunit I (COI) and 28S RNA gene sequences from eight species of Sclerodermus reared from different wood borer pests. The eight sibling species were as follows: S. guani Xiao et Wu, S. sichuanensis Xiao, S. pupariae Yang et Yao, and Sclerodermus spp. (Nos. 1-5). A 594-bp fragment of COI and 750-bp fragment of 28S were subsequently sequenced. For COI, the G-C content was found to be low in all the species, averaging to about 30.0%. Sequence divergences (Kimura-2-parameter distances) between congeneric species averaged to 4.5%, and intraspecific divergences averaged to about 0.09%. Further, the maximum sequence divergences between congeneric species and Sclerodermus sp. (No. 5) averaged to about 16.5%. All 136 samples analyzed were included in six reciprocally monophyletic clades in the COI neighbor-joining (NJ) tree. The NJ tree inferred from the 28S rRNA sequence yielded almost identical results, but the samples from S. guani, S. sichuanensis, S. pupariae, and Sclerodermus spp. (Nos. 1-4) clustered together and only Sclerodermus sp. (No. 5) clustered separately. Our findings indicate that the standard barcode region of COI can be efficiently used to distinguish morphologically similar Sclerodermus species. Further, we speculate that Sclerodermus sp. (No. 5) might be a new species of Sclerodermus.

Species delimitation using morphology, morphometrics, and molecules: definition of the Ophionscutellaris Thomson species group, with descriptions of six new species (Hymenoptera, Ichneumonidae)

The diverse genus Ophion is almost entirely undescribed in the Nearctic region. In this paper we define the Ophionscutellaris species group. This species group is well-supported by analysis of DNA (ITS2, COI, and 28S D2-D3) and morphology. It includes the Palearctic species Ophionscutellaris and the Nearctic species Ophionidoneus. An integrative analysis of DNA, geometric wing morphometrics, classical morphometrics and qualitative morphology indicates that this species group contains a minimum of seven species in North America, although the full diversity of the group has likely not been sampled. Ophionclave Schwarzfeld, sp. n., Ophionaureus Schwarzfeld, sp. n., Ophionbrevipunctatus Schwarzfeld, sp. n., Ophiondombroskii Schwarzfeld, sp. n., Ophionkeala Schwarzfeld, sp. n. and Ophionimportunus Schwarzfeld, sp. n. are described, and a key to the known Nearctic species of the Ophionscutellaris group is provided.

Patterns of DNA barcode variation in Canadian marine molluscs

BACKGROUND

Molluscs are the most diverse marine phylum and this high diversity has resulted in considerable taxonomic problems. Because the number of species in Canadian oceans remains uncertain, there is a need to incorporate molecular methods into species identifications. A 648 base pair segment of the cytochrome c oxidase subunit I gene has proven useful for the identification and discovery of species in many animal lineages. While the utility of DNA barcoding in molluscs has been demonstrated in other studies, this is the first effort to construct a DNA barcode registry for marine molluscs across such a large geographic area.

METHODOLOGY/PRINCIPAL FINDINGS

This study examines patterns of DNA barcode variation in 227 species of Canadian marine molluscs. Intraspecific sequence divergences ranged from 0-26.4% and a barcode gap existed for most taxa. Eleven cases of relatively deep (>2%) intraspecific divergence were detected, suggesting the possible presence of overlooked species. Structural variation was detected in COI with indels found in 37 species, mostly bivalves. Some indels were present in divergent lineages, primarily in the region of the first external loop, suggesting certain areas are hotspots for change. Lastly, mean GC content varied substantially among orders (24.5%-46.5%), and showed a significant positive correlation with nearest neighbour distances.

CONCLUSIONS/SIGNIFICANCE

DNA barcoding is an effective tool for the identification of Canadian marine molluscs and for revealing possible cases of overlooked species. Some species with deep intraspecific divergence showed a biogeographic partition between lineages on the Atlantic, Arctic and Pacific coasts, suggesting the role of Pleistocene glaciations in the subdivision of their populations. Indels were prevalent in the barcode region of the COI gene in bivalves and gastropods. This study highlights the efficacy of DNA barcoding for providing insights into sequence variation across a broad taxonomic group on a large geographic scale.

Molecular phylogeny, biogeography, and host plant shifts in the bee genus Melitta (Hymenoptera: Anthophila)

New molecular studies suggested that the family Melittidae is either a paraphyletic group from which all the other bees are derived, or the sister clade to all other existing bees. Studying the historical biogeography and evolution of each major lineage within this group is a key step to understand the origin and early radiation of bees. Melitta is the largest genus of melittid bees, for which a robust molecular phylogeny and a biogeographic analysis are still lacking. Here, we derive a phylogenetic hypothesis from the sequences of seven independent DNA fragments of mitochondrial and nuclear origin. This phylogenetic hypothesis is then used to infer the evolution of the species range and of the host-plant shifts in Melitta. Our results confirmed the monophyly of Melitta, but did not recover all previously defined clades within the genus. We propose new taxa by splitting the genus in three subgenera (including two new subgenera described in the Appendix: Afromelitta subgen. nov., Plesiomelitta subgen. nov.) and describe two new species: Melitta avontuurensis sp. n. and M. richtersveldensis sp. n. Regarding the evolution of host-plant use, our analysis suggests that all species currently specialized on one plant family originated from an ancestor that was specialized on Fabaceae plants. The inferred biogeographic history for the genus supported an African origin. In concordance with previous studies identifying Africa as the geographic origin for many clades of bees, our data bring new evidence for an African origin of melittid bees.

The evolution of host associations in the parasitic wasp genus Ichneumon (Hymenoptera: Ichneumonidae): convergent adaptations to host pupation sites

Background

The diversification of organisms with a parasitic lifestyle is often tightly linked to the evolution of their host associations. If a tight host association exists, closely related species tend to attack closely related hosts; host associations are less stable if associations are determined by more plastic traits like parasitoid searching and oviposition behaviour. The pupal-parasitoids of the genus Ichneumon attack a variety of macrolepidopteran hosts. They are either monophagous or polyphagous, and therefore offer a promissing system to investigate the evolution of host associations. Ichneumon was previously divided into two groups based on general body shape; however, a stout shape has been suggested as an adaptation to buried host pupation sites, and might thus not represent a reliable phylogenetic character.

Results

We here reconstruct the first molecular phylogeny of the genus Ichneumon using two mitochondrial (CO1 and NADH1) and one nuclear marker (28S). The resulting phylogeny only supports monophyly of Ichneumon when Ichneumon lugens Gravenhorst, 1829 (formerly in Chasmias, stat. rev.) and Ichneumon deliratorius Linnaeus, 1758 (formerly Coelichneumon) are included. Neither parasitoid species that attack hosts belonging to one family nor those attacking butterflies (Rhopalocera) form monophyletic clades. Ancestral state reconstructions suggest multiple transitions between searching for hosts above versus below ground and between a stout versus elongated body shape. A model assuming correlated evolution between the two characters was preferred over independent evolution of host-searching niche and body shape.

Conclusions

Host relations, both in terms of phylogeny and ecology, evolved at a high pace in the genus Ichneumon. Numerous switches between hosts of different lepidopteran families have occurred, a pattern that seems to be the rule among idiobiont parasitoids. A stout body and antennal shape in the parasitoid female is confirmed as an ecological adaptation to host pupation sites below ground and has evolved convergently several times. Morphological characters that might be involved in adaptation to hosts should be avoided as diagnostic characters for phylogeny and classification, as they can be expected to show high levels of homoplasy.

Utility of the DNA barcoding gene fragment for parasitic wasp phylogeny (Hymenoptera: Ichneumonoidea): data release and new measure of taxonomic congruence

The enormous cytochrome oxidase subunit I (COI) sequence database being assembled from the various DNA barcoding projects as well as from independent phylogenetic studies constitutes an almost unprecedented amount of data for molecular systematics, in addition to its role in species identification and discovery. As part of a study of the potential of this gene fragment to improve the accuracy of phylogenetic reconstructions, and in particular, exploring the effects of dense taxon sampling, we have assembled a data set for the hyperdiverse, cosmopolitan parasitic wasp superfamily Ichneumonoidea, including the release of 1793 unpublished sequences. Of approximately 84 currently recognized Ichneumonoidea subfamilies, 2500 genera and 41,000 described species, barcoding 5'-COI data were assembled for 4168 putative species-level terminals (many undescribed), representing 671 genera and all but ten of the currently recognized subfamilies. After the removal of identical and near-identical sequences, the 4174 initial sequences were reduced to 3278. We show that when subjected to phylogenetic analysis using both maximum likelihood and parsimony, there is a broad correlation between taxonomic congruence and number of included sequences. We additionally present a new measure of taxonomic congruence based upon the Simpson diversity index, the Simpson dominance index, which gives greater weight to morphologically recognized taxonomic groups (subfamilies) recovered with most representatives in one or a few contiguous groups or subclusters.

Parasites of vectors--Ixodiphagus hookeri and its Wolbachia symbionts in ticks in The Netherlands

Background

Ixodiphagus hookeri is a parasitic wasp of ixodid ticks around the world. It has been studied as a potential bio-control agent for several tick species. We suspected that the presence of Wolbachia infected I. hookeri eggs in ticks is responsible for incidental detection of Wolbachia DNA in tick samples.

Methods

The 28S rRNA and 16S rRNA genes of a specimen of I. hookeri was amplified and sequenced. PCR on part of the 28S rRNA gene was used to detect parasitic wasp DNA in 349 questing Ixodes ricinus ticks from various sampling sites. Furthermore, the wsp gene of Wolbachia was sequenced from the I. hookeri specimen and a subset of ticks was tested using this marker.

Results

Several sequences from tick specimens were identical to the Wolbachia sequence of the I. hookeri specimen. Ixodiphagus hookeri was detected in 9.5% of all tested ticks, varying between 4% and 26% depending on geographic location. Ten out of eleven sampling sites throughout the Netherlands were positive for I. hookeri. Eighty-seven percent of I. hookeri-positive but only 1.6% of I. hookeri-negative ticks were Wolbachia positive. Detection of I. hookeri DNA was strongly associated with the detection of Wolbachia in ticks.

Conclusion

This is the first reported case of I. hookeri in the Netherlands. Furthermore I. hookeri harbours Wolbachia species and is broadly distributed in the Netherlands. While detection of Wolbachia DNA in ticks might often be due to parasitism with this wasp, other sources of Wolbachia DNA in ticks might exist as well.

Mitochondrial heteroplasmy and DNA barcoding in Hawaiian Hylaeus (Nesoprosopis) bees (Hymenoptera: Colletidae)

BACKGROUND

The past several years have seen a flurry of papers seeking to clarify the utility and limits of DNA barcoding, particularly in areas such as species discovery and paralogy due to nuclear pseudogenes. Heteroplasmy, the coexistence of multiple mitochondrial haplotypes in a single organism, has been cited as a potentially serious problem for DNA barcoding but its effect on identification accuracy has not been tested. In addition, few studies of barcoding have tested a large group of closely-related species with a well-established morphological taxonomy. In this study we examine both of these issues, by densely sampling the Hawaiian Hylaeus bee radiation.

RESULTS

Individuals from 21 of the 49 a priori morphologically-defined species exhibited coding sequence heteroplasmy at levels of 1-6% or more. All homoplasmic species were successfully identified by COI using standard methods of analysis, but only 71% of heteroplasmic species. The success rate in identifying heteroplasmic species was increased to 86% by treating polymorphisms as character states rather than ambiguities. Nuclear pseudogenes (numts) were also present in four species, and were distinguishable from heteroplasmic sequences by patterns of nucleotide and amino acid change.

CONCLUSIONS

Heteroplasmy significantly decreased the reliability of species identification. In addition, the practical issue of dealing with large numbers of polymorphisms- and resulting increased time and labor required - makes the development of DNA barcode databases considerably more complex than has previously been suggested. The impact of heteroplasmy on the utility of DNA barcoding as a bulk specimen identification tool will depend upon its frequency across populations, which remains unknown. However, DNA barcoding is still likely to remain an important identification tool for those species that are difficult or impossible to identify through morphology, as is the case for the ecologically important solitary bee fauna.

An evaluation of phylogenetic informativeness profiles and the molecular phylogeny of diplazontinae (Hymenoptera, Ichneumonidae)

How to quantify the phylogenetic information content of a data set is a longstanding question in phylogenetics, influencing both the assessment of data quality in completed studies and the planning of future phylogenetic projects. Recently, a method has been developed that profiles the phylogenetic informativeness (PI) of a data set through time by linking its site-specific rates of change to its power to resolve relationships at different timescales. Here, we evaluate the performance of this method in the case of 2 standard genetic markers for phylogenetic reconstruction, 28S ribosomal RNA and cytochrome oxidase subunit 1 (CO1) mitochondrial DNA, with maximum parsimony, maximum likelihood, and Bayesian analyses of relationships within a group of parasitoid wasps (Hymenoptera: Ichneumonidae, Diplazontinae). Retrieving PI profiles of the 2 genes from our own and from 3 additional data sets, we find that the method repeatedly overestimates the performance of the more quickly evolving CO1 compared with 28S. We explore possible reasons for this bias, including phylogenetic uncertainty, violation of the molecular clock assumption, model misspecification, and nonstationary nucleotide composition. As none of these provides a sufficient explanation of the observed discrepancy, we use simulated data sets, based on an idealized setting, to show that the optimum evolutionary rate decreases with increasing number of taxa. We suggest that this relationship could explain why the formula derived from the 4-taxon case overrates the performance of higher versus lower rates of evolution in our case and that caution should be taken when the method is applied to data sets including more than 4 taxa.

Phylogeny of colletid bees (Hymenoptera: Colletidae) inferred from four nuclear genes

Colletidae comprise approximately 2500 species of bees primarily distributed in the southern continents (only two colletid genera are widely distributed: Colletes and Hylaeus). Previously published studies have failed to resolve phylogenetic relationships on a worldwide basis and this has been a major barrier to the progress of research regarding systematics and evolution of colletid bees. For this study, data from four nuclear gene loci: elongation factor-1alpha (F2 copy), opsin, wingless, and 28S rRNA were analyzed for 122 species of colletid bees, representing all subfamilies and tribes currently recognized; 22 species belonging to three other bee families were used as outgroups. Bayesian, maximum likelihood, and parsimony methods were employed to investigate the phylogenetic relationships within Colletidae and resulted in highly congruent and well-resolved trees. The phylogenetic results show that Colletidae are monophyletic and that all traditionally recognized subfamilies (except Paracolletinae) are also strongly supported as monophyletic. Our phylogenetic hypothesis provides a framework within which broad questions related to the taxonomy, biogeography, morphology, evolution, and ecology of colletid bees can be addressed.

Phylogeny of the basal angiosperm genus Pseuduvaria (Annonaceae) inferred from five chloroplast DNA regions, with interpretation of morphological character evolution

Phylogenetic relationships within the magnoliid basal angiosperm genus Pseuduvaria (Annonaceae) are investigated using chloroplast DNA sequences from five regions: psbA-trnH spacer, trnL-F, matK, rbcL, and atpB-rbcL spacer. Over 4000 nucleotides from 51 species (of the total 53) were sequenced. The five cpDNA datasets were analyzed separately and in combination using maximum parsimony (MP), maximum likelihood (ML), and Bayesian methods. The phylogenetic trees constructed using all three phylogenetic methods, based on the combined data, strongly support the monophyly of Pseuduvaria following the inclusion of Craibella phuyensis. The trees generated using MP were less well resolved, but relationships are similar to those obtained using the other methods. ML and Bayesian analyses recovered trees with short branch lengths, showing five main clades. This study highlights the evolutionary changes in seven selected morphological characters (floral sex, stamen and carpel numbers, inner petal color, presence of inner petal glands, flowering peduncle length, and monocarp size). Although floral unisexuality is ancestral within the genus, several evolutionary lineages reveal reversal to bisexuality. Other phylogenetic transitions include the evolution of sapromyophily, and fruit-bat frugivory and seed dispersal, thus allowing a wide range of adaptations for species survival.

Ancient rapid radiations of insects: challenges for phylogenetic analysis

Phylogenies of major groups of insects based on both morphological and molecular data have sometimes been contentious, often lacking the data to distinguish between alternative views of relationships. This paucity of data is often due to real biological and historical causes, such as shortness of time spans between divergences for evolution to occur and long time spans after divergences for subsequent evolutionary changes to obscure the earlier ones. Another reason for difficulty in resolving some of the relationships using molecular data is the limited spectrum of genes so far developed for phylogeny estimation. For this latter issue, there is cause for current optimism owing to rapid increases in our knowledge of comparative genomics. At least some historical patterns of divergence may, however, continue to defy our attempts to completely reconstruct them with confidence, at least using current strategies.

Molecular phylogeny and historical biogeography of the cosmopolitan parasitic wasp subfamily Doryctinae (Hymenoptera:Braconidae)

The phylogenetic relationships among representatives of 64 genera of the cosmopolitan parasitic wasps of the subfamily Doryctinae were investigated based on nuclear 28S ribosomal (r) DNA (~650 bp of the D2–3 region) and cytochrome c oxidase I (COI) mitochondrial (mt) DNA (603 bp) sequence data. The molecular dating of selected clades and the biogeography of the subfamily were also inferred. The partitioned Bayesian analyses did not recover a monophyletic Doryctinae, though the relationships involved were only weakly supported. Strong evidence was found for rejecting the monophylies of both Doryctes Haliday, 1836 and Spathius Nees, 1818. Our results also support the recognition of the Rhaconotini as a valid tribe. A dispersal–vicariance analysis showed a strong geographical signal for the taxa included, with molecular dating estimates for the origin of Doryctinae and its subsequent radiation both occurring during the late Paleocene–early Eocene. The divergence time estimates suggest that diversification in the subfamily could have in part occurred as a result of continental break-up events that took place in the southern hemisphere, though more recent dispersal events account for the current distribution of several widespread taxa.

Barcoding ciliates: a comprehensive study of 75 isolates of the genus Tetrahymena

The mitochondrial cytochrome-coxidase subunit 1 (cox1) gene has been proposed as a DNA barcode to identify animal species. To test the applicability of thecox1gene in identifying ciliates, 75 isolates of the genusTetrahymenaand three non-Tetrahymenaciliates that are close relatives ofTetrahymena,Colpidium campylum,Colpidium colpodaandGlaucoma chattoni, were selected. All tetrahymenines of unproblematic species could be identified to the species level using 689 bp of thecox1sequence, with about 11 % interspecific sequence divergence. Intraspecific isolates ofTetrahymena borealis,Tetrahymena lwoffi,Tetrahymena patulaandTetrahymena thermophilacould be identified by theircox1sequences, showing <0.65 % intraspecific sequence divergence. In addition, isolates of these species were clustered together on acox1neighbour-joining (NJ) tree. However, strains identified asTetrahymena pyriformisandTetrahymena tropicalisshowed high intraspecific sequence divergence values of 5.01 and 9.07 %, respectively, and did not cluster together on acox1NJ tree. This may indicate the presence of cryptic species. The mean interspecific sequence divergence ofTetrahymenawas about 11 times greater than the mean intraspecific sequence divergence, and this increased to 58 times when all isolates of species with high intraspecific sequence divergence were excluded. This result is similar to DNA barcoding studies on animals, indicating that congeneric sequence divergences are an order of magnitude greater than conspecific sequence divergences. Our analysis also demonstrated low sequence divergences of <1.0 % between some isolates ofT. pyriformisandTetrahymena setosaon the one hand and some isolates ofTetrahymena furgasoniandT. lwoffion the other, suggesting that the latter species in each pair is a junior synonym of the former. Overall, our study demonstrates the feasibility of using the mitochondrialcox1gene as a taxonomic marker for ‘barcoding’ and identifyingTetrahymenaspecies and some other ciliated protists.

Deciphering ancient rapid radiations

A deeper phylogenetic understanding of ancient patterns of diversification would contribute to solving many problems in evolutionary biology, yet many of these phylogenies remain poorly resolved. Ancient rapid radiations pose a major challenge for phylogenetic analysis for two main reasons. First, the pattern to be deciphered, the order of divergence among lineages, tends to be supported by small amounts of data. Second, the time since divergence is large and, thus, the potential for misinterpreting phylogenetic information is great. Here, we review the underlying causes of difficulty in determining the branching patterns of diversification in ancient rapid radiations, and review novel data exploration tools that can facilitate understanding of these radiations.

Monophyly and phylogenetic relationships of Thereva and therevine genus-groups (Insecta:Diptera:Therevidae) based on EF-1α, 28S rDNA and mitochondrial 16S rDNA sequences

Phylogenetic analyses using 28S rDNA, elongation factor (EF)-1α, and mt 16S rDNA sequences were performed to test the monophyly of Thereva Latreille. Two of the three Afrotropical Thereva species groups lack the genitalia characters that unambiguously diagnose Thereva in the Holarctic Region, but phylogenetic relationships among Thereva species groups and therevine genera are poorly understood. Using an extensive taxonomic sample (39 of the 62 therevine genera) and Thereva, sensu lato (15 spp.), simultaneous analyses of all three gene partitions recovered Nearctic and Palaearctic Thereva species in a well supported clade that includes the Afrotropical seminitida-group but excludes the Afrotropical analis- and turneri-groups. Stronger phylogenetic signal from the EF-1α partition, measured by the skewness statistic and proportion of total parsimony informative characters, dominated conflicting signal from the 16S partition and weaker, but more congruent, signal from 28S. Reducing the taxonomic sample in analyses of Therevinae reduced homoplasy, increased phylogenetic structure and partitioned Bremer support values and reduced incongruence with 28S for the 16S partition. Although molecular analyses yielded partial recovery of informal therevine genus-groups, morphological diagnoses of higher-level groups are poorly supported with the exception of Cyclotelini. The ‘Holarctic radiation’ refers to a diverse clade of genera closely related to Pandivirilia Irwin & Lyneborg and Acrosathe Irwin & Lyneborg widely distributed throughout the Holarctic Region that is the sister-group to Thereva, sensu stricto. Results from these analyses underscore the importance of male and female genitalia characters in recognising monophyletic groups and regional endemism in therevine diversification.

Dissecting the ancient rapid radiation of microgastrine wasp genera using additional nuclear genes

Previous estimates of a generic level phylogeny for the ubiquitous parasitoid wasp subfamily Microgastrinae (Hymenoptera) have been problematic due to short internal branches deep in the phylogeny. These short branches might be attributed to a rapid radiation among the taxa, the use of genes that are unsuitable for the levels of divergence being examined, or insufficient quantity of data. We added over 1200 nucleotides from four nuclear genes to a dataset derived from three genes to produce a dataset of over 3000 nucleotides per taxon. While the number of well-supported short branches in the phylogeny increased, we still did not obtain strong bootstrap support for every node. Parametric and nonparametric bootstrap simulations projected that an enormous, and likely unobtainable, amount of data would be required to get bootstrap support greater than 50% for every node. However, a marked increase in the number of well-supported nodes was seen when we conducted a Bayesian analysis of a combined dataset generated from morphological characters added to the seven gene dataset. Our results suggest that, in some cases, combining morphological and genetic characters may be the most practical way to increase support for short branches deep in a phylogeny.

Evolutionary relationships among the Scarabaeini (Coleoptera: Scarabaeidae) based on combined molecular and morphological data

The Scarabaeini is an old world tribe of ball-rolling dung beetles that have origins dating back to at least the mid-upper Miocene (19-8 million years ago). The tribe has received little to no attention in morphological or molecular phylogenetics. We obtained sequence data from the mitochondrial cytochrome oxidase subunit I (1,197 bp) and 16S ribosomal RNA (461 bp) genes for 25 species of the Scarabaeini in an attempt to further resolve broad phylogenetic relationships within this tribe. Sequence data from both markers along with 216 morphological and 3 biological characters were analysed separately and combined. Independent analyses showed poorly resolved trees with many of the intermediate and basal nodes collapsed by low bootstrap values. Many sites in both genes exhibited strong A+T nucleotide bias and high interlineage divergences. The combined analysis revealed a number of well supported relationships such as the monophyly of the nocturnal species Scarabaeus satyrus, S. [Neateuchus] proboscideus, and S. zambesianus. Furthermore, the total evidence tree suggested to elevate S. (Pachysoma) to the status of an independent genus, Pachysoma, as a sister taxon to a clade containing Pachylomerus femoralis and Scarabaeus sensu lato. Within the latter, the following subgenera were maintained by the combination of data sets: S. (Scarabaeolus), S. (Sceliages), and S. (Kheper). Both, feeding specialisation and food relocation behaviour, were inferred to be polyphyletic in the Scarabaeini. Total evidence analysis found no support for common ancestry of Scarabaeini and Eucraniini.

Mosaics of convergences and noise in morphological phylogenies: what's in a viverrid-like carnivoran?

Adaptive convergence in morphological characters has not been thoroughly investigated, and the processes by which phylogenetic relationships may be misled by morphological convergence remains unclear. We undertook a case study on the morphological evolution of viverrid-like feliformians (Nandinia, Cryptoprocta, Fossa, Eupleres, Prionodon) and built the largest morphological matrix concerning the suborder Feliformia to date. A total of 349 characters grouped into four anatomical partitions were used for all species of Viverridae and viverrid-like taxa plus representatives of the Felidae, Hyaenidae, Herpestidae, and one Malagasy mongoose. Recent molecular phylogenetic analyses suggest that viverrid-like morphotypes appeared independently at least three times during feliformian evolution. We thus used a synthetic molecular tree to assess morphological evolutionary patterns characterizing the viverrid-like taxa. We examined phylogenetic signal, convergence and noise in morphological characters using (a) tree-length distribution (g1), (b) partitioned Bremer support, (c) RI values and their distribution, (d) respective contributions of diagnostic synapomorphies at the nodes for each partition, (e) patterns of shared convergences among viverrid-like taxa and other feliformian lineages, (f) tree-length differences among alternative hypotheses, and (g) the successive removal of convergent character states from the original matrix. In addition, the lability of complex morphological structures was assessed by mapping them onto the synthetic molecular tree. The unconstrained morphological analysis yielded phylogenetic groupings that closely reflected traditional classification. The use of a synthetic molecular tree (constraint) combined with our thorough morphological investigations revealed the mosaics of convergences likely to have contributed to part of the historical uncertainty over viverrid classification. It also showed that complex morphological structures could be subjected to reversible evolutionary trends. The morphological matrix proved useful in characterizing several feliformian clades with diagnostic synapomorphies. These results support the removal from the traditionally held Viverridae of several viverrid-like taxa into three distinct families: Nandiniidae (Nandinia), Prionodontidae (Prionodon), and the newly defined Eupleridae (including Cryptoprocta, Fossa, Eupleres plus all "mongoose-like" Malagasy taxa). No clearly "phylogenetically misleading" data subsets could be identified, and the great majority of morphological convergences appeared to be nonadaptive. The multiple approaches used in this study revealed that the most disruptive element with regards to morphological phylogenetic reconstruction was noise, which blured the expression of phylogenetic signal. This study demonstrates the crucial need to consider independent (molecular) phylogenies in order to produce reliable evolutionary hypotheses and should promote a new approach to the definition of morphological characters in mammals. [Constrained analysis; convergence; evolutionary scenario; Feliformia; morphology; noise; phylogenetic signal; phylogeny; Viverridae.].

Systematics of the cyclostome subfamilies of braconid parasitic wasps (Hymenoptera: Ichneumonoidea): A simultaneous molecular and morphological Bayesian approach

Phylogenetic relationships among 95 genera collectively representing 17 of the 18 currently recognized cyclostome braconid wasp subfamilies were investigated based on DNA sequence fragments of the mitochondrial COI and the nuclear 28S rDNA genes, in addition to morphological data. The treatment of sequence length variation of the 28S partition was explored by either excluding ambiguously aligned regions and indel information (28SN) or recoding them (28SA) using the 'fragment-level' alignment method with a modified coding approach. Bayesian MCMC analyses were performed for the separate and combined data sets and a maximum parsimony analysis was also carried out for the simultaneous molecular and morphological data sets. There was a significant incongruence between the two genes and between 28S and morphology, but not for morphology and COI. Different analyses with the 28SA data matrix resulted in topologies that were generally similar to the ones from the 28SN matrix; however, the former topologies recovered a higher number of significantly supported clades and had a higher mean Bayesian posterior probability, thus supporting the inclusion of information from ambiguously aligned regions and indel events in phylogenetic analyses where possible. Based on the significantly supported clades obtained from the simultaneous molecular and morphological analyses, we propose that a total of 17 subfamilies should be recognized within the cyclostome group. The subfamilial placements of several problematic cyclostome genera were also established.

Genetic variability in the three genomes of Nasonia: nuclear, mitochondrial and Wolbachia

Nasonia consists of three closely related species of parasitoid wasps that are all infected with the endosymbiotic bacteria Wolbachia, a reproductive parasite common in arthropods. This situation presents the opportunity to compare patterns of variation in three associated genomes, Wolbachia and the nuclear and mitochondrial genomes of its host. Furthermore, although Nasonia wasps are emerging as a model for evolutionary and genetic studies, little is known about their genetic variability. Using amplified fragment length polymorphisms (AFLPs), all three species present a relatively high level of nuclear polymorphism and have different patterns of variation, with one of the species, Nasonia giraulti, being divided into two divergent subgroups. In each species, the mitochondrial pattern of variation is different from the nuclear pattern, possibly due to genetic hitchhiking of the mitochondria during (cytoplasmically inherited) Wolbachia sweeps. Mitochondria in Nasonia show a synonymous substitution rate approximately 10-15-fold higher than nuclear genes, probably reflecting an elevated mitochondrial mutation rate that is among the highest found in insects. Finally, all three species are doubly infected with their own strains of Wolbachia, one each from the two major supergroups (A and B). Sequence analysis reveals that each of the three Nasonia species acquired their A and B bacteria independently by horizontal transfer events from other insects with the exception of B type Wolbachia in N. longicornis and N. giraulti, which were acquired prior to speciation and then codiverged with the host. This represents one of the few clear-cut examples of codivergence of Wolbachia during host speciation.