Mixed Yule Coalescent Model Research Articles

Species identification in the taxonomically neglected, highly diverse, neotropical parasitoid wasp genus Notiospathius (Braconidae: Doryctinae) based on an integrative molecular and morphological approach

Various DNA sequence-based methods for species delineation have recently been developed to assess the species-richness of highly diverse, neglected invertebrate taxa. These methods, however, need to be tested under a variety of conditions, including the use of different markers and parameters. Here, we explored the species diversity of a species-rich group of braconid parasitoid wasps, the Neotropical genus Notiospathius, including 233 specimens from 10 different countries. We examined sequences of two mitochondrial (mt) (COI, cyt b) and one nuclear (wg) gene fragments. We analysed them separately as well as concatenating the mt data with the general mixed Yule-coalescent (GMYC) model for species delineation using different tree-building methods and parameters for reconstructing ultrametric trees. We evaluated the performance of GMYC analyses by comparing their species delineations with our morphospecies identifications. Reconstructing ultrametric trees with a relaxed lognormal clock rate using the program BEAST gave the most congruent results with morphology for the two mt markers. A tree obtained with wg using the programs MrBayes+Pathd8 had the fewest cases of incongruence with morphology, though the performance of this nuclear marker was considerably lower than that of COI and cyt b. Species delimitation using the coalescent prior to obtain ultrametric trees was morphologically more congruent with COI, whereas the Yule prior was more congruent with cyt b. The analyses concatenating the mt datasets failed to recover some species supported both by morphology and the separate analyses of the mt markers. The highest morphological congruence was obtained with the GMYC analysis on an ultrametric tree reconstructed with cyt b using the relaxed lognormal clock rate and the Yule prior, thus supporting the importance of using alternative markers when the information of the barcoding locus (COI) is not concordant with morphological evidence. Seventy-one species were delimited based on the congruence found among COI, cyt b and morphology. Both mt markers also revealed the existence of seven potential cryptic species. This high species richness from a scattered geographical sampling indicates that there is a remarkable number of Notiospathius species that remains undiscovered.

Toward a DNA Taxonomy of Alpine Rhithrogena (Ephemeroptera: Heptageniidae) Using a Mixed Yule-Coalescent Analysis of Mitochondrial and Nuclear DNA

Aquatic larvae of many Rhithrogena mayflies (Ephemeroptera) inhabit sensitive Alpine environments. A number of species are on the IUCN Red List and many recognized species have restricted distributions and are of conservation interest. Despite their ecological and conservation importance, ambiguous morphological differences among closely related species suggest that the current taxonomy may not accurately reflect the evolutionary diversity of the group. Here we examined the species status of nearly 50% of European Rhithrogena diversity using a widespread sampling scheme of Alpine species that included 22 type localities, general mixed Yule-coalescent (GMYC) model analysis of one standard mtDNA marker and one newly developed nDNA marker, and morphological identification where possible. Using sequences from 533 individuals from 144 sampling localities, we observed significant clustering of the mitochondrial (cox1) marker into 31 GMYC species. Twenty-one of these could be identified based on the presence of topotypes (expertly identified specimens from the species' type locality) or unambiguous morphology. These results strongly suggest the presence of both cryptic diversity and taxonomic oversplitting in Rhithrogena. Significant clustering was not detected with protein-coding nuclear PEPCK, although nine GMYC species were congruent with well supported terminal clusters of nDNA. Lack of greater congruence in the two data sets may be the result of incomplete sorting of ancestral polymorphism. Bayesian phylogenetic analyses of both gene regions recovered four of the six recognized Rhithrogena species groups in our samples as monophyletic. Future development of more nuclear markers would facilitate multi-locus analysis of unresolved, closely related species pairs. The DNA taxonomy developed here lays the groundwork for a future revision of the important but cryptic Rhithrogena genus in Europe.

SPECIATION AND DISPERSAL ALONG CONTINENTAL COASTLINES AND ISLAND ARCS IN THE INDO-WEST PACIFIC TURBINID GASTROPOD GENUSLUNELLA

Species trees were produced for the Indo-West Pacific (IWP) gastropod genus Lunella using MrBayes, BEAST, and *BEAST with sequence data from four genes. Three fossil records were used to calibrate a molecular clock. Eight cryptic species were recognized using statistical methods for species delimitation in combination with morphological differences. However, our results suggest caution in interpreting ESUs defined solely by the general mixed Yule Coalescent model in genera like Lunella, with lower dispersal abilities. Four almost entirely allopatric species groups were recovered that differ in ecology and distribution. Three groups occur predominantly along continental coastlines and one occurs on island arrays. Sympatric species occur only in the torquata and coronata groups along coastlines, whereas species in the cinerea group, distributed in two-dimensional island arrays, occur in complete allopatry. Dispersal along island arcs has been important in the maintenance of species distributions and gene flow among populations in the cinerea group. The emergence of new islands and their eventual subsidence over geological time has had important consequences for the isolation of populations and the eventual rise of new species in Lunella.

Evolutionary criteria outperform operational approaches in producing ecologically relevant fungal species inventories

Analyses of the structure and function of microbial communities are highly constrained by the diversity of organisms present within most environmental samples. A common approach is to rely almost entirely on DNA sequence data for estimates of microbial diversity, but to date there is no objective method of clustering sequences into groups that is grounded in evolutionary theory of what constitutes a biological lineage. The general mixed Yule-coalescent (GMYC) model uses a likelihood-based approach to distinguish population-level processes within lineages from processes associated with speciation and extinction, thus identifying a distinct point where extant lineages became independent. Using two independent surveys of DNA sequences associated with a group of ubiquitous plant-symbiotic fungi, we compared estimates of species richness derived using the GMYC model to those based on operational taxonomic units (OTUs) defined by fixed levels of sequence similarity. The model predicted lower species richness in these surveys than did traditional methods of sequence similarity. Here, we show for the first time that groups delineated by the GMYC model better explained variation in the distribution of fungi in relation to putative niche-based variables associated with host species identity, edaphic factors, and aspects of how the sampled ecosystems were managed. Our results suggest the coalescent-based GMYC model successfully groups environmental sequences of fungi into clusters that are ecologically more meaningful than more arbitrary approaches for estimating species richness.

DNA barcoding a highly diverse group of parasitoid wasps (Braconidae: Doryctinae) from a Mexican nature reserve

Background and aims. The preliminary results of a DNA barcoding study of the doryctine fauna of parasitoid wasps from the Chamela–Cuixmala Biosphere Reserve in Mexico, a region dominated by tropical dry forest, are presented. So far, three field trips have been carried out to the reserve and 468 specimens have been collected, of which 407 cox1 sequences were obtained.Materials and methods. The general mixed Yule-coalescent model was applied to a phylogram to investigate the number of evolutionary units that can be detected from the DNA sequence data examined.Results. A total of 185 barcoding species assigned to 20 identified doryctine genera were discriminated using the above model, 115 of which belong to the speciose genus Heterospilus, pointing out the extraordinary species richness of this subfamily of insects in a Mexican tropical dry forest.Conclusion. On the basis of the DNA barcodes generated, Ptesimogastroides Braet & van Achterberg is proposed to be a junior synonym of Ptesimogaster Marsh syn. nov. Neoheterospilus was also found deeply nested within a large Heterospilus clade, suggesting the paraphyly of the latter genus.

DNA-based association of adults and larvae in Baetidae (Ephemeroptera) with the description of a new genusAdnoptilumin Madagascar

The mayfly fauna of Madagascar is highly diverse and largely endemic. Many species remain undescribed, and many species are known from only the larval or adult life stage. The high biodiversity in Madagascar and in other areas has led to an increasing reliance on DNA-based approaches to taxonomy, i.e., to define species boundaries and to associate different developmental stages. We used the general mixed Yule-coalescent (GMYC) model to combine population- and species-level sequence variation of mitochondrial deoxyribonucleic acid (mtDNA) to detect species boundaries in Baetidae mayflies (Ephemeroptera). Starting with a database of 240 sequences (57 species), significant clustering of newly sequenced larvae allowed us to establish 1 new species and 1 new combination and to associate adult and larval stages for both. A molecular phylogeny using additional nuclear (18S) and mtDNA (rrnS, rrnL) gene regions recovered the new species and new combination as a monophyletic group, distinct from other Afrotropical lineages. Therefore, we established a new genus, Adnoptilum, endemic to Madagascar. Adnoptilum gen. n. can be distinguished from other species in both the imaginal and adult stages and appears to belong to the Bugilliesia complex. We conclude that routine sampling of population- and species-level genetic diversity, combined with coalescent-based methods of species delineation, has great potential to become a standard procedure for the study of poorly known taxonomic groups.

Revisiting the Insect Mitochondrial Molecular Clock: The Mid-Aegean Trench Calibration

Phylogenetic trees in insects are frequently dated by applying a "standard" mitochondrial DNA (mtDNA) clock estimated at 2.3% My(-1), but despite its wide use reliable calibration points have been lacking. Here, we used a well-established biogeographic barrier, the mid-Aegean trench separating the western and eastern Aegean archipelago, to estimate substitution rates in tenebrionid beetles. Cytochrome oxidase I (cox1) for six codistributed genera across 28 islands (444 individuals) on both sides of the mid-Aegean trench revealed 60 independently coalescing entities delimited with a mixed Yule-coalescent model. One representative per entity was used for phylogenetic analysis of mitochondrial (cox1, 16S rRNA) and nuclear (Mp20, 28S rRNA) genes. Six nodes marked geographically congruent east-west splits whose separation was largely contemporaneous and likely to reflect the formation of the mid-Aegean trench at 9-12 Mya. Based on these "known" dates, a divergence rate of 3.54% My(-1) for the cox1 gene (2.69% when combined with the 16S rRNA gene) was obtained under the preferred partitioning scheme and substitution model selected using Bayes factors. An extensive survey suggests that discrepancies in mtDNA substitution rates in the entomological literature can be attributed to the use of different substitution models, the use of different mitochondrial gene regions, mixing of intraspecific with interspecific data, and not accounting for variance in coalescent times or postseparation gene flow. Different treatments of these factors in the literature confound estimates of mtDNA substitution rates in opposing directions and obscure lineage-specific differences in rates when comparing data from various sources.

DNA taxonomy in morphologically plastic taxa: Algorithmic species delimitation in the Boodlea complex (Chlorophyta: Cladophorales)

DNA-based taxonomy provides a convenient and reliable tool for species delimitation, especially in organisms in which morphological discrimination is difficult or impossible, such as many algal taxa. A group with a long history of confusing species circumscriptions is the morphologically plastic Boodlea complex, comprising the marine green algal genera Boodlea, Cladophoropsis, Phyllodictyon and Struveopsis. In this study, we elucidate species boundaries in the Boodlea complex by analysing nrDNA internal transcribed spacer sequences from 175 specimens collected from a wide geographical range. Algorithmic methods of sequence-based species delineation were applied, including statistical parsimony network analysis, and a maximum likelihood approach that uses a mixed Yule-coalescent model and detects species boundaries based on differences in branching rates at the level of species and populations. Sequence analyses resulted in the recognition of 13 phylogenetic species, although we failed to detect sharp species boundaries, possibly as a result of incomplete reproductive isolation. We found considerable conflict between traditional and phylogenetic species definitions. Identical morphological forms were distributed in different clades (cryptic diversity), and at the same time most of the phylogenetic species contained a mixture of different morphologies (indicating intraspecific morphological variation). Sampling outside the morphological range of the Boodlea complex revealed that the enigmatic, sponge-associated Cladophoropsis ( Spongocladia) vaucheriiformis, also falls within the Boodlea complex. Given the observed evolutionary complexity and nomenclatural problems associated with establishing a Linnaean taxonomy for this group, we propose to discard provisionally the misleading morphospecies and genus names, and refer to clade numbers within a single genus, Boodlea.

Accelerated Species Inventory on Madagascar Using Coalescent-Based Models of Species Delineation

High-throughput DNA sequencing has the potential to accelerate species discovery if it is able to recognize evolutionary entities from sequence data that are comparable to species. The general mixed Yule-coalescent (GMYC) model estimates the species boundary from DNA surveys by identifying independently evolving lineages as a transition from coalescent to speciation branching patterns on a phylogenetic tree. Applied here to 12 families from 4 orders of insects in Madagascar, we used the model to delineate 370 putative species from mitochondrial DNA sequence variation among 1614 individuals. These were compared with data from the nuclear genome and morphological identification and found to be highly congruent (98% and 94%). We developed a modified GMYC that allows for a variable transition from coalescent to speciation among lineages. This revised model increased the congruence with morphology (97%), suggesting that a variable threshold better reflects the clustering of sequence data into biological species. Local endemism was pronounced in all 5 insect groups. Most species (60-91%) and haplotypes (88-99%) were found at only 1 of the 5 study sites (40-1000 km apart). This pronounced endemism resulted in a 37% increase in species numbers using diagnostic nucleotides in a population aggregation analysis. Sample sizes between 7 and 10 individuals represented a threshold above which there was minimal increase in genetic diversity, broadly agreeing with coalescent theory and other empirical studies. Our results from > 1.4 Mb of empirical data suggest that the GMYC model captures species boundaries comparable to those from traditional methods without the need for prior hypotheses of population coherence. This provides a method of species discovery and biodiversity assessment using single-locus data from mixed or environmental samples while building a globally available taxonomic database for future identifications.

Comparative phylogeography of tenebrionid beetles in the Aegean archipelago: the effect of dispersal ability and habitat preference

Comparative phylogeographical studies in island archipelagos can reveal lineage-specific differential responses to the geological and climatic history. We analysed patterns of genetic diversity in six codistributed lineages of darkling beetles (Tenebrionidae) in the central Aegean archipelago which differ in wing development and habitat preferences. A total of 600 specimens from 30 islands and eight adjacent mainland regions were sequenced for mitochondrial cytochrome oxidase I and nuclear Muscular protein 20. Individual gene genealogies were assessed for the presence of groups that obey an independent coalescent process using a mixed Yule coalescent model. The six focal taxa differed greatly in the number of coalescent groups and depth of lineage subdivision, which was closely mirrored by the degree of geographical structuring. The most severe subdivision at both mitochondrial DNA and nuclear DNA level was found in flightless lineages associated with presumed stable compact-soil habitats (phrygana, maquis), in contrast to sand-obligate lineages inhabiting ephemeral coastal areas that displayed greater homogeneity across the archipelago. A winged lineage, although associated with stable habitats, showed no significant phylogenetic or geographical structuring. Patterns of nucleotide diversity and local genetic differentiation, as measured using PhiST and hierarchical AMOVA, were consistent with high levels of ongoing gene flow in the winged taxon; frequent local extinction and island recolonisation for flightless sand-obligate taxa; and very low gene flow and geographical structure largely defined by the palaeogeographical history of the region in flightless compact-soil taxa. These results show that differences in dispersal rate, mediated by habitat persistence, greatly influence the levels of phylogeographical subdivision in lineages that are otherwise subjected to the same geological events and palaeoclimatic changes.