GMYC Model Research Articles

Leaf gall diversity in the neotropical plant Coccoloba barbadensis (Polygonaceae) and its associated insect species richness

Galls are an abnormal growth of plant tissue in response to the presence generally of an inducing insect, which ensures food and protection during specific periods of its life. Besides gall formers, a vast community of arthropods are associated with galls, including inquilines and parasitoids. Few studies have assessed the gall diversity and its associated insect community in Neotropical vascular plants. Here, we characterised the leaf gall diversity of Coccoloba barbadensis Jacq. (Polygonaceae) in a Mexican tropical dry forest, as well as their associated entomofauna based on morphology and DNA barcoding. Five different gall morphotypes were observed during both dry (April-June) and rainy (November) seasons. A total of 34 and 38 species of Diptera, Hymenoptera, Coleoptera, and Lepidoptera were delimited with the 2% divergence criterion and the GMYC model, respectively. Based on our rearing observations and literature, Cecidomyiidae (Diptera) species might induce all leaf gall morphotypes, whereas hymenopterans are represented by parasitoid and probably inquiline species of the families Braconidae, Eulophidae, Eupelmidae, Platygastridae and Torymidae. Our results highlight the importance of performing integrative species delineation studies of arthropods present in galls to have an accurate knowledge of their diversity and trophic interactions.

Molecular phylogeny and species delimitation of the genus Tonkinacris (Orthoptera, Acrididae, Melanoplinae) from China.

Tonkinacris is a small group in Acrididae. While a few species were occasionally sampled in some previous molecular studies, there is no revisionary research devoted to the genus. In this study, we explored the phylogeny of and the relationships among Chinese species of the genus Tonkinacris using the mitochondrial COI barcode and the complete sequences of ITS1 and ITS2 of the nuclear ribosomal DNA. The phylogeny was reconstructed in maximum likelihood and Bayesian inference frameworks, respectively. The overlap range between intraspecific variation and interspecific divergence was assessed via K2P distances. Species boundaries were delimitated using phylogenetic species concept, NJ tree, K2P distance, the statistical parsimony network as well as the GMYC model. The results demonstrate that the Chinese Tonkinacris species is a monophyletic group and the phylogenetic relationship among them is (T. sinensis, (T. meridionalis, (T. decoratus, T. damingshanus))). While T. sinensis, T. meridionalis and T. decoratus were confirmed being good independent species strongly supported by both morphological and molecular evidences, the validity of T. damingshanus was not perfectly supported by molecular evidence in this study.

Two ways to be endemic. Alps and Apennines are different functional refugia during climatic cycles.

Endemics co-occur because they evolved in situ and persist regionally or because they evolved ex situ and later dispersed to shared habitats, generating evolutionary or ecological endemicity centres, respectively. We investigate whether different endemicity centres can intertwine in the region ranging from Alps to Sicily, by studying their butterfly fauna. We gathered an extensive occurrence data set for butterflies of the study area (27,123 records, 269 species, in cells of 0.5×0.5 degrees of latitude-longitude). We applied molecular-based delimitation methods (GMYC model) to 26,557 cytochrome c oxidase subunit 1 (COI) sequences of Western Palearctic butterflies. We identified entities based on molecular delimitations and/or the checklist of European butterflies and objectively attributed occurrences to their most probable entity. We obtained a zoogeographic regionalisation based on the 69 endemics of the area. Using phylogenetic ANOVA we tested if endemics from different centres differ from each other and from nonendemics for key ecological traits and divergence time. Endemicity showed high incidence in the Alps and Southern Italy. The regionalisation separated the Alps from the Italian Peninsula and Sicily. The endemics of different centres showed a high turnover and differed in phylogenetic distances, phenology and distribution traits. Endemics are on average younger than nonendemics and the Peninsula-Sicily endemics also have lower variance in divergence than those from the Alps. The observed variation identifies Alpine endemics as paleoendemics, now occupying an ecological centre, and the Peninsula-Sicily ones as neoendemics, that diverged in the region since the Pleistocene. The results challenge the common view of the Alpine-Apennine area as a single "Italian refugium".

Natural enemies of armored scales (Hemiptera: Diaspididae) and soft scales (Hemiptera: Coccidae) in Chile: Molecular and morphological identification.

Scale insects (Hemiptera: Sternorrhyncha: Coccomorpha) are key pests of agricultural crops and ornamental plants worldwide. Their populations are difficult to control, even with insecticides, due to their cryptic habits. Moreover, there is growing concern over the use of synthetic pesticides for their control, due to deleterious environmental effects and the emergence of resistant populations of target pests. In this context, biological control may be an effective and sustainable approach. Hymenoptera Chalcidoidea includes natural enemies of scale insects that have been successfully used in many biological control programs. However, the correct identification of pest scale species and their natural enemies is particularly challenging because these insects are very small and highly specialized. Integrative taxonomy, coupling DNA barcoding and morphological analysis, has been successfully used to characterize pests and natural enemy species. In this study, we performed a survey of parasitoids and predators of armored and soft scales in Chile, based on 28S and COI barcodes. Fifty-three populations of Diaspididae and 79 populations of Coccidae were sampled over the entire length of the country, from Arica (18°S) to Frutillar (41°S), between January 2015 and February 2016. The phylogenetic relationships obtained by Bayesian inference from multilocus haplotypes revealed 41 putative species of Chalcidoidea, five Coccinellidae and three Neuroptera. Species delimitation was confirmed using ABGD, GMYC and PTP model. In Chalcidoidea, 23 species were identified morphologically, resulting in new COI barcodes for 12 species and new 28S barcodes for 14 species. Two predator species (Rhyzobius lophantae and Coccidophilus transandinus) were identified morphologically, and two parasitoid species, Chartocerus niger and Signiphora bifasciata, were recorded for the first time in Chile.

A Multi-Gene Analysis and Potential Spatial Distribution of Species of the Strodei Subgroup of the Genus Nyssorhynchus (Diptera: Culicidae).

Nyssorhynchus strodei (Root) is an understudied species of potential epidemiological importance, having been found naturally infected in Brazil with Plasmodium falciparum Welch, Plasmodium vivax Grassi & Feletti, and Plasmodium malariae Feletti & Grassi. It belongs to the strodei subgroup that includes Nyssorhynchus albertoi (Unti), Nyssorhynchus arthuri (Unti), Nyssorhynchus rondoni (Neiva & Pinto), Nyssorhynchus striatus (Sant'Ana & Sallum), and three unnamed species, Nyssorhynchus arthuri B, Nyssorhynchus arthuri C, and Nyssorhynchus arthuri D. As the accurate identification of vector species is of fundamental importance for public health entomology, the aim of the study was to provide additional information for the presence of seven species that had been previously misidentified as Ny. strodei. Bayesian phylogenetic analyses using DNA sequences of the genes COI mtDNA, white, CAD and CAT nuclear genes confirmed Ny. albertoi, Ny. arthuri, Ny. strodei, and Ny. striatus as distinct groups within the strodei subgroup and corroborated the presence of three undescribed species under the name Ny. arthuri. Results of the GMYC model analysis corroborated Ny. arthuri B, Ny. arthuri C, and Ny. arthuri D; however, they did not distinguish between Ny. strodei and Ny. albertoi. Predicted distribution of seven species based on maximum entropy in MaxEnt showed that each species has its specific ecological niche suitability.

Genetic and Morphological Analyses Demonstrate That Schizolecis guntheri (Siluriformes: Loricariidae) Is Likely to Be a Species Complex.

Schizolecis is a monotypic genus of Siluriformes widely distributed throughout isolated coastal drainages of southeastern Brazil. Previous studies have shown that fish groups found in isolated river basins tend to differentiate over time in the absence of gene flow, resulting in allopatric speciation. In this study, we used partial sequences of the mitochondrial gene COI with the analysis of the General Mixed Yule Coalescent model (GMYC) and the Automatic Barcode Gap Discovery (ABGD) for single locus species delimitation, and a Principal Component Analysis (PCA) of external morphology to test the hypothesis that Schizolecis guntheri is a complex of species. We analyzed 94 samples of S. guntheri for GMYC and ABGD, and 82 samples for PCA from 22 coastal rivers draining to the Atlantic in southeastern Brazil from the Paraná State to the north of the Rio de Janeiro State. As a result, the GMYC model and the ABGD delimited five operational taxonomy units (OTUs – a nomenclature referred to in the present study of the possible new species delimited for the genetic analysis), a much higher number compared to the traditional alfa taxonomy that only recognizes S. guntheri across the isolated coastal rivers of Brazil. Furthermore, the PCA analysis suggests that S. guntheri is highly variable in aspects of external body proportions, including dorsal-fin spine length, pectoral-fin spine length, pelvic-fin spine length, lower caudal-fin spine length, caudal peduncle depth, anal width and mandibular ramus length. However, no exclusive character was found among the isolated populations that could be used to describe a new species of Schizolecis. Therefore, we can conclude, based on our results of PCA contrasting with the results of GMYC and ABGD, that S. guntheri represents a complex of species.

Temporal patterns and processes of genetic differentiation of the Brachionus calyciflorus (Rotifera) complex in a subtropical shallow lake

Sympatric distribution and temporal overlap of cryptic zooplankton species pose a challenge to the framework of the niche differentiation theory and the mechanisms underlying the coexistence of potentially strong competitors. In this study, we collected Brachionus calyciflorus Pallas monthly from Lake Tingtang over the course of 1 year and clonally cultured them in the laboratory. Finally, 204 clones were established, and their nuDNA ITS1 and mtCOI sequences were analyzed. The phylogenetic relationships and GMYC model analyses revealed that B. calyciflorus was a species complex composed of three cryptic species. Cryptic species BcI-W, and BcII-S and Bc-SW underwent a clear seasonal succession. However, BcII-S and Bc-SW overlapped in June and July, and all the three cryptic species coexisted in November. Different preferences to water temperature might contribute to the seasonal succession. More investment in sexual reproduction accelerated the disappearance of cryptic species from the water column during the seasonal succession process. Differential responses in rate of population growth of the two overlapping cryptic species to increasing food level might reduce the intensity of competition between them and promote their coexistence over time. Differences in population growth rate partially resulted in the different relative abundances of the two overlapping cryptic species.

Neoplecostomus canastra, a new catfish (Teleostei: Siluriformes) species from upper Rio Paraná basin

In this study, we describe a new species of Neoplecostomus from córrego Tamborete, municipality of Capitólio, Minas Gerais state, upper rio Paraná basin. The new species is distinguished from its congeners by the presence of conspicuous dark spots all over the body. Moreover, the new species can be diagnosed by the absence of moderate keels along each lateral series of plates, presence of developed adipose fin, absence of enlarged, fleshy folds between dentaries, dorsal-fin spinelet larger than dorsal-fin spine, lacking enlarged odontodes along lateral margins of snout and along ridges before eyes, having reduced number of dentary teeth, having two irregular and conspicuous rows of large and transversally-flattened papillae, just posterior to dentary teeth, smaller mandibullary width/head length, smaller Interorbital width/mandibullary width, smaller number of lateral-line plates, higher and lower number of dentary teeth depending of the species compared, absence of a conspicuous sexual dimorphism, absence of a conspicuous horseshoe-shaped light blotch posterior to the supraoccipital, with a central dark area which rarely contacts the edge of the light area, smaller mandibulary width/head length, smaller Interorbital width/Mandibullary width and a greater Caudal peduncle length/Caudal peduncle depth. Furthermore, we used DNA barcoding techniques of 2% threshold criteria and GMYC model to distinguish the new species from its congeners. These molecular techniques were unable to distinguish the new species from N. langeanii and N. selenae from N. botucatu, that were differentiated using external morphology.

Species delimitation of the Hyphydrus ovatus complex in western Palaearctic with an update of species distributions (Coleoptera, Dytiscidae).

The species status of Hyphydrus anatolicus Guignot, 1957 and H. sanctus Sharp, 1882, previously often confused with the widespread H. ovatus (Linnaeus, 1760), are tested with molecular and morphological characters. Cytochrome c oxidase subunit 1 (CO1) was sequenced for 32 specimens of all three species. Gene-trees were inferred with parsimony, time-free bayesian and strict clock bayesian analyses. The GMYC model was used to estimate species limits. All three species were reciprocally monophyletic with CO1 and highly supported. The GMYC species delimitation analysis unequivocally delimited the three species with no other than the three species solution included in the confidence interval. A likelihood ratio test rejected the one-species null model. Important morphological characters distinguishing the species are provided and illustrated. New distributional data are given for the following species: Hyphydrus anatolicus from Slovakia and Ukraine, and H. aubei Ganglbauer, 1891, and H. sanctus from Turkey.

Genetic and morphological analyses indicate that the Australian endemic scorpion Urodacus yaschenkoi (Scorpiones: Urodacidae) is a species complex.

BackgroundAustralian scorpions have received far less attention from researchers than their overseas counterparts. Here we provide the first insight into the molecular variation and evolutionary history of the endemic Australian scorpion Urodacus yaschenkoi. Also known as the inland robust scorpion, it is widely distributed throughout arid zones of the continent and is emerging as a model organism in biomedical research due to the chemical nature of its venom.MethodsWe employed Bayesian Inference (BI) methods for the phylogenetic reconstructions and divergence dating among lineages, using unique haplotype sequences from two mitochondrial loci (COXI, 16S) and one nuclear locus (28S). We also implemented two DNA taxonomy approaches (GMYC and PTP/dPTP) to evaluate the presence of cryptic species. Linear Discriminant Analysis was used to test whether the linear combination of 21 variables (ratios of morphological measurements) can predict individual’s membership to a putative species.ResultsGenetic and morphological data suggest that U. yaschenkoi is a species complex. High statistical support for the monophyly of several divergent lineages was found both at the mitochondrial loci and at a nuclear locus. The extent of mitochondrial divergence between these lineages exceeds estimates of interspecific divergence reported for other scorpion groups. The GMYC model and the PTP/bPTP approach identified major lineages and several sub-lineages as putative species. Ratios of several traits that approximate body shape had a strong predictive power (83–100%) in discriminating two major molecular lineages. A time-calibrated phylogeny dates the early divergence at the onset of continental-wide aridification in late Miocene and Pliocene, with finer-scale phylogeographic patterns emerging during the Pleistocene. This structuring dynamics is congruent with the diversification history of other fauna of the Australian arid zones.DiscussionOur results indicate that the taxonomic status of U. yaschenkoi requires revision, and we provide recommendations for such future efforts. A complex evolutionary history and extensive diversity highlights the importance of conserving U. yaschenkoi populations from different Australian arid zones in order to preserve patterns of endemism and evolutionary potential.

Multilocus Phylogeography and Species Delimitation in the Cumberland Plateau Salamander, Plethodon kentucki: Incongruence among Data Sets and Methods.

Species are a fundamental unit of biodiversity, yet can be challenging to delimit objectively. This is particularly true of species complexes characterized by high levels of population genetic structure, hybridization between genetic groups, isolation by distance, and limited phenotypic variation. Previous work on the Cumberland Plateau Salamander, Plethodon kentucki, suggested that it might constitute a species complex despite occupying a relatively small geographic range. To examine this hypothesis, we sampled 135 individuals from 43 populations, and used four mitochondrial loci and five nuclear loci (5693 base pairs) to quantify phylogeographic structure and probe for cryptic species diversity. Rates of evolution for each locus were inferred using the multidistribute package, and time calibrated gene trees and species trees were inferred using BEAST 2 and *BEAST 2, respectively. Because the parameter space relevant for species delimitation is large and complex, and all methods make simplifying assumptions that may lead them to fail, we conducted an array of analyses. Our assumption was that strongly supported species would be congruent across methods. Putative species were first delimited using a Bayesian implementation of the GMYC model (bGMYC), Geneland, and Brownie. We then validated these species using the genealogical sorting index and BPP. We found substantial phylogeographic diversity using mtDNA, including four divergent clades and an inferred common ancestor at 14.9 myr (95% HPD: 10.8–19.7 myr). By contrast, this diversity was not corroborated by nuclear sequence data, which exhibited low levels of variation and weak phylogeographic structure. Species trees estimated a far younger root than did the mtDNA data, closer to 1.0 myr old. Mutually exclusive putative species were identified by the different approaches. Possible causes of data set discordance, and the problem of species delimitation in complexes with high levels of population structure and introgressive hybridization, are discussed.

Extraordinary micro-endemism in Australian desert spring amphipods

Increasing pressure for water in the Australian arid zone is placing enormous stress on the diverse endemic communities inhabiting desert springs. Detailed information about the evolutionary processes occurring within and between individual endemic species will help to develop effective and biologically relevant management strategies this fragile ecosystem. To help determine conservation priorities, we documented the genetic structure of the endemic freshwater amphipod populations in springs fed by the Great Artesian Basin in central Australia. Phylogenetic and phylogeographic history and genetic diversity measures were examined using nuclear and mitochondrial DNA from approximately 500 chiltoniid amphipods across an entire group of springs. Pronounced genetic diversity was identified, demonstrating that levels of endemism have been grossly underestimated in these amphipods. Using the GMYC model, 13 genetically divergent lineages were recognized as Evolutionarily Significant Units (ESUs), all of which could be considered as separate species. The results show that due to the highly fragmented ecosystem, these taxa have highly restricted distributions. Many of the identified ESUs are endemic to a very small number of already degraded springs, with the rarest existing in single springs. Despite their extraordinarily small ranges, most ESUs showed relative demographic stability and high levels of genetic diversity, and genetic diversity was not directly linked to habitat extent. The relatively robust genetic health of ESUs does not preclude them from endangerment, as their limited distributions ensure they will be highly vulnerable to future water extraction.

Phylogenetic estimation error can decrease the accuracy of species delimitation: a Bayesian implementation of the general mixed Yule-coalescent model

BackgroundSpecies are considered the fundamental unit in many ecological and evolutionary analyses, yet accurate, complete, accessible taxonomic frameworks with which to identify them are often unavailable to researchers. In such cases DNA sequence-based species delimitation has been proposed as a means of estimating species boundaries for further analysis. Several methods have been proposed to accomplish this. Here we present a Bayesian implementation of an evolutionary model-based method, the general mixed Yule-coalescent model (GMYC). Our implementation integrates over the parameters of the model and uncertainty in phylogenetic relationships using the output of widely available phylogenetic models and Markov-Chain Monte Carlo (MCMC) simulation in order to produce marginal probabilities of species identities.ResultsWe conducted simulations testing the effects of species evolutionary history, levels of intraspecific sampling and number of nucleotides sequenced. We also re-analyze the dataset used to introduce the original GMYC model. We found that the model results are improved with addition of DNA sequence and increased sampling, although these improvements have limits. The most important factor in the success of the model is the underlying phylogenetic history of the species under consideration. Recent and rapid divergences result in higher amounts of uncertainty in the model and eventually cause the model to fail to accurately assess uncertainty in species limits.ConclusionOur results suggest that the GMYC model can be useful under a wide variety of circumstances, particularly in cases where divergences are deeper, or taxon sampling is incomplete, as in many studies of ecological communities, but that, in accordance with expectations from coalescent theory, rapid, recent radiations may yield inaccurate results. Our implementation differs from existing ones in two ways: it allows for the accounting for important sources of uncertainty in the model (phylogenetic and in parameters specific to the model) and in the specification of informative prior distributions that can increase the precision of the model. We have incorporated this model into a user-friendly R package available on the authors’ websites.

Geographical and seasonal evidence of cryptic diversity in the Baetis rhodani complex (Ephemeroptera, Baetidae) revealed by means of DNA taxonomy

Previous phylogenetic investigations on the mayfly Baetis rhodani Pictet from several European countries, excluding Italy, strongly suggested the presence of cryptic species. Our paper reports a DNA-taxonomy phylogenetic analysis of B. rhodani with additional populations coming from Italian and UK sites, and aims to identify potential cryptic species with a coalescent-based method (GMYC model) and to understand the mechanisms of local coexistence of cryptic species. Twenty-five haplotypes of Italian samples and five haplotypes of UK samples were identified and added to a large European dataset. A total of 11 potential cryptic species have been recognised, and three of them co-occured in one Italian area. Such cryptic species seem to be phylogenetically over-dispersed on the tree and temporally segregated, and the seasonal substitution pattern of cryptic species could explain the apparently widespread distribution of the B. rhodani complex and its ability to adapt to different temperatures and food resources, justifying some of the differences observed in the relationship between water temperature, growth rates and phenology documented from field studies.

Endemism and evolutionary history in conflict over Madagascar’s freshwater conservation priorities

Regional-scale biodiversity indicators provide important criteria for the selection of protected areas in conservation, but their application is often hindered by a lack of taxonomic knowledge. Moreover, different indicators include different types of information, sometimes leading to divergent conservation priorities. Madagascar tops the world list of biodiversity hotspots and much conservation effort has been directed toward its threatened plants and vertebrates. In contrast, its highly diverse freshwater invertebrate fauna has received comparatively little conservation attention. We conducted an inventory of Malagasy adephagan water beetles (Coleoptera, Dytiscidae, Noteridae, Gyrinidae, Haliplidae) using a combined morphological and molecular approach. In total, 2043 beetles from 153 sites were sequenced for cytochrome oxidase subunit I ( cox1), and species delimitation was carried out using the coalescent-based GMYC model. Phylogenetic relationships of the resulting entities were established using cox1 combined with partial 16S rRNA and 28S rRNA sequences. Ten national parks were assessed for their species richness, phylogenetic diversity (PD) and endemism. We were particularly interested in the contribution of endemic species to PD. Congruence between molecular and taxonomic identifications was high (91%), with 69% of sampled species endemic to Madagascar. Interestingly, we found that PD at a site was negatively correlated to the proportion of endemic species, most likely because endemics are the result of recent radiations with relatively little branch-length contribution to the measure of PD. This suggests that ranking sites for conservation priority based solely on PD potentially disfavor endemic species by underrating areas where the evolutionary process is most active.

Patterns and processes in the genetic differentiation of the Brachionus calyciflorus complex, a passively dispersing freshwater zooplankton

Elucidating the evolutionary patterns and processes of extant species is an important objective of any research program that seeks to understand population divergence and, ultimately, speciation. The island-like nature and temporal fluctuation of limnetic habitats create opportunities for genetic differentiation in rotifers through space and time. To gain further understanding of spatio-temporal patterns of genetic differentiation in rotifers other than the well-studied Brachionus plicatilis complex in brackish water, a total of 318 nrDNA ITS sequences from the B. calyciflorus complex in freshwater were analysed using phylogenetic and phylogeographic methods. DNA taxonomy conducted by both the sequence divergence and the GMYC model suggested the occurrence of six potential cryptic species, supported also by reproductive isolation among the tested lineages. The significant genetic differentiation and non-significant correlation between geographic and genetic distances existed in the most abundant cryptic species, BcI-W and Bc-SW. The large proportion of genetic variability for cryptic species Bc-SW was due to differences between sampling localities within seasons, rather than between different seasons. Nested Clade Analysis suggested allopatric or past fragmentation, contiguous range expansion and long-distance colonization possibly coupled with subsequent fragmentation as the probable main forces shaping the present-day phylogeographic structure of the B. calyciflorus species complex.