Utility of DNA taxonomy and barcoding for the inference of larval community structure in morphologically cryptic Chironomus (Diptera) species

The first record of Chironomusnuditarsis Keyl, 1961 from Sevan Lake (Armenia) confirmed by morphology, karyotype and COI gene sequence

Chironomusnuditarsis Keyl, 1961 is recorded from Sevan Lake for the first time. This species is widespread in Europe, the Caucasus, and Siberia. For species identification, we used a comprehensive approach that included morphological, cytogenetic and molecular genetic analyses. Morphological analysis showed a high similarity with the description. Nine chromosome banding sequences ndtA1, ndtA2, ndtB2, ndtC1, ndtD1, ndtE1, ndtF1, ndtG1, and ndtG2 were found. The banding sequences ndtA1, ndtA2, ndtG1, and ndtG2 are species-specific for C.nuditarsis and allow us to accurately distinguish it from the sibling species Ch.curabilis Belyanina, Sigareva et Loginova, 1990. Molecular-genetic analysis of the COI gene sequences has shown low genetic distances of 0.38-0.95% in the sibling species Ch.nuditarsis and Ch.curabilis complex and the insufficiency of using a single COI as a molecular marker for their separation.

Facing the infinity: tackling large samples of challenging Chironomidae (Diptera) with an integrative approach

Background

Integrative taxonomy is becoming ever more significant in biodiversity research as scientists are tackling increasingly taxonomically challenging groups. Implementing a combined approach not only guarantees more accurate species identification, but also helps overcome limitations that each method presents when applied on its own. In this study, we present one application of integrative taxonomy for the highly abundant and particularly diverse fly taxon Chironomidae (Diptera). Although non-biting midges are key organisms in merolimnic systems, they are often cast aside in ecological surveys because they are very challenging to identify and extremely abundant.

Methods

Here, we demonstrate one way of applying integrative methods to tackle this highly diverse taxon. We present a three-level subsampling method to drastically reduce the workload of bulk sample processing, then apply morphological and molecular identification methods in parallel to evaluate species diversity and to examine inconsistencies across methods.

Results

Our results suggest that using our subsampling approach, identifying less than 10% of a sample's contents can reliably detect >90% of its diversity. However, despite reducing the processing workload drastically, the performance of our taxonomist was affected by mistakes, caused by large amounts of material. We conducted misidentifications for 9% of vouchers, which may not have been recovered had we not applied a second identification method. On the other hand, we were able to provide species information in cases where molecular methods could not, which was the case for 14% of vouchers. Therefore, we conclude that when wanting to implement non-biting midges into ecological frameworks, it is imperative to use an integrative approach.

Phylogenetic Analysis of Wild Species and the Maternal Origin of Cultivars in the Genus Lilium Using 114 Plastid Genomes

Lilies are one of the most important ornamental flowers worldwide with approximately 100 wild species and numerous cultivars, but the phylogenetic relationships among wild species and their contributions to these cultivars are poorly resolved. We collected the major Lilium species and cultivars and assembled their plastome sequences. Our phylogenetic reconstruction using 114 plastid genomes, including 70 wild species representing all sections and 42 cultivars representing six hybrid divisions and two outgroups, uncovered well-supported genetic relationships within Lilium. The wild species were separated into two distinct groups (groups A and B) associated with geographical distribution, which further diversified into eight different clades that were phylogenetically well supported. Additional support was provided by the distributions of indels and single-nucleotide variants, which were consistent with the topology. The species of sections Archelirion, Sinomartagon III, and Leucolirion 6a and 6b were the maternal donors for Oriental hybrids, Asiatic hybrids, Trumpet hybrids, and Longiflorum hybrids, respectively. The maternal donors of the OT hybrids originated from the two sections Archelirion and Leucolirion 6a, and LA hybrids were derived from the two sections Leucolirion 6b and Sinomartagon. Our study provides an important basis for clarifying the infrageneric classification and the maternal origin of cultivars in Lilium.

Improved Chironomid Barcode Database Enhances Identification of Water Mite Dietary Content

Chironomids are one of the most biodiverse and abundant members of freshwater ecosystems. They are a food source for many organisms, including fish and water mites. The accurate identification of chironomids is essential for many applications in ecological research, including determining which chironomid species are present in the diets of diverse predators. Larval and adult chironomids from diverse habitats, including lakes, rivers, inland gardens, coastal vegetation, and nearshore habitats of the Great Lakes, were collected from 2012 to 2019. After morphological identification of chironomids, DNA was extracted and cytochrome oxidase I (COI) barcodes were PCR amplified and sequenced. Here we describe an analysis of biodiverse adult and larval chironomids in the Great Lakes region of North America based on new collections to improve chironomid identification by curating a chironomid DNA barcode database, thereby expanding the diversity and taxonomic specificity of DNA reference libraries for the Chironomidae family. In addition to reporting many novel chironomid DNA barcodes, we demonstrate here the use of this chironomid COI barcode database to improve the identification of DNA barcodes of prey in the liquefied diets of water mites. The species identifications of the COI barcodes of chironomids ingested by Lebertia davidcooki and L. quinquemaculosa are more diverse for L. davidcooki and include Parachironomus abortivus, Cryptochironomus ponderosus. Parachironomus tenuicaudatus, Glyptotendipes senilis, Dicrotendipes modestus, Chironomus riparius, Chironomus entis/plumosus, Chironomus maturus, Chironomus crassicaudatus, Endochironomus subtendens, Cricotopus sylvestris, Cricotopus festivellus, Orthocladius obumbratus, Tanypus punctipennis, Rheotanytarsus exiguus gr., and Paratanytarsus nr. bituberculatus.

The efficacy of DNA barcoding in the classification, genetic differentiation, and biodiversity assessment of benthic macroinvertebrates

Macroinvertebrates have been recognized as key ecological indicators of aquatic environment and are the most commonly used approaches for water quality assessment. However, species identification of macroinvertebrates (especially of aquatic insects) proves to be very difficult due to the lack of taxonomic expertise in some regions and can become time-consuming. In this study, we evaluated the feasibility of DNA barcoding for the classification of benthic macroinvertebrates and investigated the genetic differentiation in seven orders (Insecta: Ephemeroptera, Plecoptera, Trichoptera, Diptera, Hemiptera, Coleoptera, and Odonata) from four large transboundary rivers of northwest China and further explored its potential application to biodiversity assessment. A total of 1,144 COI sequences, belonging to 176 species, 112 genera, and 53 families were obtained and analyzed. The barcoding gap analysis showed that COI gene fragment yielded significant intra- and interspecific divergences and obvious barcoding gaps. NJ phylogenetic trees showed that all species group into monophyletic species clusters whether from the same population or not, except two species (Polypedilum. laetum and Polypedilum. bullum). The distance-based (ABGD) and tree-based (PTP and MPTP) methods were utilized for grouping specimens into Operational Taxonomic Units (OTUs) and delimiting species. The ABGD, PTP, and MPTP analysis were divided into 177 (p = .0599), 197, and 195 OTUs, respectively. The BIN analysis generated 186 different BINs. Overall, our study showed that DNA barcoding offers an effective framework for macroinvertebrate species identification and sheds new light on the biodiversity assessment of local macroinvertebrates. Also, the construction of DNA barcode reference library of benthic macroinvertebrates in Eurasian transboundary rivers provides a solid backup for bioassessment studies of freshwater habitats using modern high-throughput technologies in the near future.

Pond chironomid communities revealed by molecular species delimitation reflect eutrophication

Farm ponds, a valued habitat for freshwater organisms, are being negatively affected by the recent changes in the environment as well as anthropological activities. In these ponds, biodiversity researchers have tended to focus on species that prefer natural habitats and/or can be identified based on morphological characters. In contrast, this study focused on the insect family Chironomidae, which is widely distributed from clear to polluted waters of ponds, but is hard to identify morphologically as an aquatic larva. We adopted DNA barcoding and molecular species delimitation to identify every single specimen of quantitative collections. From bottom sediments of 17 ponds in summer in the Banshu Plain of Japan, a total of 62 species were delimited based on the DNA sequences of the mitochondrial COI region. Chironomid communities from these ponds were classified into four groups in a two-dimensional ordination of multivariate analysis (NMDS). One of the dimensions was well correlated with the gradient of eutrophication, while another dimension was not clearly assigned to any general feature of the environmental gradient, but rice cultivation could possibly be involved.

A Bioinformatics Tutorial for Comparative Development Genomics in Diverse Meiofauna

Miniaturization, which is a common feature in animals, is particularly manifest in meiofauna-animals sharing peculiar phenotypic features that evolved as adaptations to the highly specialized aquatic interstitial habitat. While revealing much about the extreme phyletic diversity of meiofauna, the genome structure of meiofaunal species could also characterize the phenotype of ancestral states as well as explain the origin and evolution of miniaturization. Here, we present a practical bioinformatics tutorial for genome assembly, genome comparison, and characterization of Hox clusters in meiofaunal species.

Molecular and morphometric divergence of four species of butterflies (Nymphalidae and Pieridae) from the Western Himalaya, India

Morphometric and molecular divergence among four butterfly species of the families Nymphalidae and Pieridae from the western Himalaya region were investigated using molecular tools, traditional morphometric measures and a truss network system. The considered species were Danaus chrysippus, Vanessa cardui, Pieris brassicae and Pieris canidia. Traditional taxonomy is sometimes unable to discriminate cryptic species or species that have close morphological features. Although taxonomists carefully examine external body features to differentiate the species; however, there is a risk for misidentification during a visual assessment of cryptic species. Therefore, we aimed to use the truss network system of 14 morphological landmarks interconnected to yield 90 variables about molecular taxonomy. Principal component analysis (PCA), discriminant function analysis (DFA) and cluster analysis (CA) were employed to determine morphometric variations. In the traditional analysis, 79, 68, 16 and 5 characters out of 90 were found significant (p < 0.05) for D. chrysippus, V. cardui, P. brassicae and P. canidia, respectively. One to seven principal components were extracted through PCA; they explained 87.5-100% of the total variance in samples. Notably, DFA correctly classified 100% of the original grouped cases and 100% of the cross-validated grouped cases. However, the variations were not the same for the two different methods (truss and traditional) employed for the analysis. We correctly identified all the species; the interspecies sequence divergence was between 0.1034 and 0.1398, and the intra-species sequence divergence range was 0.0001 to 0.0128 using the Cytochrome c oxidase subunit-I (COI) gene. The present study provides useful information about the application and complementary role of traditional with truss morphometric analysis for the precise identification and classification of the selected species.

DNA barcoding and molecular phylogeny of Dumasia (Fabaceae: Phaseoleae) reveals a cryptic lineage

Dumasia taxonomy and classification have long been problematic. Species within this genus have few morphological differences and plants without flowers or fruits are difficult to accurately identify. In this study, we evaluated the ability of six DNA barcoding sequences, one nuclear (ITS) and five chloroplast regions (trnH-psbA, matK, rbcL, trnL-trnF, psbB-psbF), to efficiently identify Dumasia species. Most single markers or their combinations identify obvious barcoding gaps between intraspecific and interspecific genetic variation. Most combined analyses including ITS showed good species resolution and identification efficiency. We therefore suggest that ITS alone or a combination of ITS with any cpDNA marker are most suitable for DNA barcoding of Dumasia. The phylogenetic analyses clearly indicated that Dumasia yunnanensis is not monophyletic and is separated as two independent branches, which may result from cryptic differentiation. Our results demonstrate that molecular data can deepen the comprehension of taxonomy of Dumasia and provide an efficient approach for identification of the species.

A DNA barcode library for 5,200 German flies and midges (Insecta: Diptera) and its implications for metabarcoding-based biomonitoring

This study summarizes results of a DNA barcoding campaign on German Diptera, involving analysis of 45,040 specimens. The resultant DNA barcode library includes records for 2,453 named species comprising a total of 5,200 barcode index numbers (BINs), including 2,700 COI haplotype clusters without species‐level assignment, so called “dark taxa.” Overall, 88 out of 117 families (75%) recorded from Germany were covered, representing more than 50% of the 9,544 known species of German Diptera. Until now, most of these families, especially the most diverse, have been taxonomically inaccessible. By contrast, within a few years this study provided an intermediate taxonomic system for half of the German Dipteran fauna, which will provide a useful foundation for subsequent detailed, integrative taxonomic studies. Using DNA extracts derived from bulk collections made by Malaise traps, we further demonstrate that species delineation using BINs and operational taxonomic units (OTUs) constitutes an effective method for biodiversity studies using DNA metabarcoding. As the reference libraries continue to grow, and gaps in the species catalogue are filled, BIN lists assembled by metabarcoding will provide greater taxonomic resolution. The present study has three main goals: (a) to provide a DNA barcode library for 5,200 BINs of Diptera; (b) to demonstrate, based on the example of bulk extractions from a Malaise trap experiment, that DNA barcode clusters, labelled with globally unique identifiers (such as OTUs and/or BINs), provide a pragmatic, accurate solution to the “taxonomic impediment”; and (c) to demonstrate that interim names based on BINs and OTUs obtained through metabarcoding provide an effective method for studies on species‐rich groups that are usually neglected in biodiversity research projects because of their unresolved taxonomy.

Hidden Diversity Hampers Conservation Efforts in a Highly Impacted Neotropical River System

Neotropical Rivers host a highly diverse ichthyofauna, but taxonomic uncertainty prevents appropriate conservation measures. The Doce River Basin (DRB), lying within two Brazilian threatened hotspots (Atlantic Forest and Brazilian Savanna) in south-east Brazil, faced the worst ever environmental accident reported for South American catchments, due to a dam collapse that spread toxic mining tailings along the course of its main river. Its ichthyofauna was known to comprise 71 native freshwater fish species, of which 13 endemic. Here, we build a DNA barcode library for the DRB ichthyofauna, using samples obtained before the 2015 mining disaster, in order to provide a more robust biodiversity record for this basin, as a baseline for future management actions. Throughout the whole DRB, we obtained a total of 306 barcodes, assigned to 69 putative species (with a mean of 4.54 barcodes per species), belonging to 45 genera, 18 families, and 5 orders. Average genetic distances within species, genus, and families were 2.59, 11.4, and 20.5%, respectively. The 69 species identified represent over 76% of the known DRB ichthyofauna, comprising 43 native (five endemic, of which three threatened by extinction), 13 already known introduced species, and 13 unknown species (such as Characidium sp., Neoplecostomus sp., and specimens identified only at the sub-family level Neoplecostominae, according to morphological identification provided by the museum collections). Over one fifth of all analyzed species (N = 16) had a mean intraspecific genetic divergence higher than 2%. An integrative approach, combining NND (nearest neighbor distance), BIN (barcode index number), ABGD (automatic barcode gap discovery), and bPTP (Bayesian Poisson Tree Processes model) analyses, suggested the occurrence of potential cryptic species, species complex, or historical errors in morphological identification. The evidence presented calls for a more robust, DNA-assisted cataloging of biodiversity-rich ecosystems, in order to enable effective monitoring and informed actions to preserve and restore these delicate habitats.

The genomic footprint of climate adaptation in Chironomus riparius

The gradual heterogeneity of climatic factors poses varying selection pressures across geographic distances that leave signatures of clinal variation in the genome. Separating signatures of clinal adaptation from signatures of other evolutionary forces, such as demographic processes, genetic drift and adaptation, to nonclinal conditions of the immediate local environment is a major challenge. Here, we examine climate adaptation in five natural populations of the harlequin fly Chironomus riparius sampled along a climatic gradient across Europe. Our study integrates experimental data, individual genome resequencing, Pool-Seq data and population genetic modelling. Common-garden experiments revealed significantly different population growth rates at test temperatures corresponding to the population origin along the climate gradient, suggesting thermal adaptation on the phenotypic level. Based on a population genomic analysis, we derived empirical estimates of historical demography and migration. We used an F_ST outlier approach to infer positive selection across the climate gradient, in combination with an environmental association analysis. In total, we identified 162 candidate genes as genomic basis of climate adaptation. Enriched functions among these candidate genes involved the apoptotic process and molecular response to heat, as well as functions identified in studies of climate adaptation in other insects. Our results show that local climate conditions impose strong selection pressures and lead to genomic adaptation despite strong gene flow. Moreover, these results imply that selection to different climatic conditions seems to converge on a functional level, at least between different insect species.

Population genetic structure and hybridization patterns in the cryptic sister species Chironomus riparius and Chironomus piger across differentially polluted freshwater systems

Chironomids are an integral and functionally important part of many freshwater ecosystems. Yet, to date, there is limited understanding of their microevolutionary processes under chemically polluted natural environments. In this study, we investigated the genetic variation within populations of the ecotoxicological model species Chironomus riparius and its cryptic sister species Chironomus piger at 18 metal-contaminated and reference sites in northwestern Portugal. Microsatellite analysis was conducted on 909 samples to answer if metal contamination affects genetic variation in natural chironomid populations as previously suggested from controlled laboratory experiments. Similarly high levels of genetic diversity and significant but weak genetic substructuring were found across all sites and temporal replicates, with no effects of metal contamination on the genetic variation or species' abundance, although C. piger tended to be less frequent at highly contaminated sites. Our results indicate that high levels of gene flow and population dynamic processes may overlay potential pollutant effects. At least for our study species, we conclude that the "genetic erosion hypothesis", which suggests that chemical pollution will reduce genome-wide genetic variability in affected populations, does not hold under natural conditions. Interestingly, our study provides evidence of successful hybridization between the two sister species under natural conditions.

Evolutionary consequences of historical metal contamination for natural populations of Chironomus riparius (Diptera: Chironomidae)

Populations inhabiting metal-impacted freshwater systems located nearby industrial and urban areas may be under intense selection. The present study aims to address two fundamental microevolutionary aspects of metal contamination in the midge Chironomus riparius (Meigen): Are populations inhabiting historically metal contaminated sites genetically adapted to metals? And, are populations from these sites genetically eroded? To answer these questions, C. riparius populations were sampled from three sites with well-known histories of metal contamination and three nearby-located references. Genetic adaptation to metals was investigated through acute and chronic exposures to cadmium (Cd), after rearing all populations for at least six generations under laboratory clean conditions. Genetic diversity was estimated based on the allelic variation of seven microsatellite markers. Results showed higher acute tolerance to Cd in populations originating from metal contaminated sites compared to their respective references and significant differences in two out of three pairwise comparisons. However, there was a mismatch between acute and chronic tolerance to Cd with results of the partial life-cycle tests suggesting fitness costs under control clean conditions in two metal-adapted populations. Despite no evidences of genetic erosion in populations sampled from metal contaminated sites, our results suggest genetically inherited tolerance to Cd in populations inhabiting historically contaminated sites. These findings lend support to the use of C. riparius as a model organism in evolutionary toxicology and highlight the importance of coupling measures of neutral genetic diversity with assessments of chemical tolerance of populations for a better understanding of contaminant-induced adaptation and evolutionary processes.

Improving the Conservation of Mediterranean Chondrichthyans: The ELASMOMED DNA Barcode Reference Library

Cartilaginous fish are particularly vulnerable to anthropogenic stressors and environmental change because of their K-selected reproductive strategy. Accurate data from scientific surveys and landings are essential to assess conservation status and to develop robust protection and management plans. Currently available data are often incomplete or incorrect as a result of inaccurate species identifications, due to a high level of morphological stasis, especially among closely related taxa. Moreover, several diagnostic characters clearly visible in adult specimens are less evident in juveniles. Here we present results generated by the ELASMOMED Consortium, a regional network aiming to sample and DNA-barcode the Mediterranean Chondrichthyans with the ultimate goal to provide a comprehensive DNA barcode reference library. This library will support and improve the molecular taxonomy of this group and the effectiveness of management and conservation measures. We successfully barcoded 882 individuals belonging to 42 species (17 sharks, 24 batoids and one chimaera), including four endemic and several threatened ones. Morphological misidentifications were found across most orders, further confirming the need for a comprehensive DNA barcoding library as a valuable tool for the reliable identification of specimens in support of taxonomist who are reviewing current identification keys. Despite low intraspecific variation among their barcode sequences and reduced samples size, five species showed preliminary evidence of phylogeographic structure. Overall, the ELASMOMED initiative further emphasizes the key role accurate DNA barcoding libraries play in establishing reliable diagnostic species specific features in otherwise taxonomically problematic groups for biodiversity management and conservation actions.

The role of genetic diversity and past-history selection pressures in the susceptibility of Chironomus riparius populations to environmental stress

Natural populations experiencing intense selection and genetic drift may exhibit limited potential to adapt to environmental change. The present study addresses the following aspects of the "genetic erosion" hypothesis in the midge Chironomus riparius: does long-term mercury (Hg) contamination affect the Hg tolerance of midge populations inhabiting such impacted areas? If so, is there any fitness cost under changing environmental conditions? And does genetic impoverishment influence the susceptibility of C. riparius to cope with environmental stressful conditions? For this end, we tested the acute and chronic tolerance to Hg and salinity in four C. riparius populations differing in their levels of genetic diversity (assessed through microsatellite markers) and past-histories of Hg exposure. Results showed that the midge population collected from a heavily Hg-contaminated site had higher Hg tolerance compared to the population collected from a closely-located reference site suggesting directional selection for Hg-tolerant traits in its native environment despite no genetic erosion in the field. No increased susceptibility under changing environmental conditions of salinity stress was observed. Moreover, results also showed that populations with higher genetic diversity performed better in the partial life-cycle assays providing evidence on the key role that genetic diversity plays as mediator of populations' susceptibility to environmental stress. Our findings are discussed in terms of the suitability of C. riparius as a model organism in evolutionary toxicology studies as well as the validity of ecotoxicological assessments using genetically eroded laboratory populations.

Support for the evolutionary speed hypothesis from intraspecific population genetic data in the non-biting midge Chironomus riparius

The evolutionary speed hypothesis (ESH) proposes a causal mechanism for the latitudinal diversity gradient. The central idea of the ESH is that warmer temperatures lead to shorter generation times and increased mutation rates. On an absolute time scale, both should lead to an acceleration of selection and drift. Based on the ESH, we developed predictions regarding the distribution of intraspecific genetic diversity: populations of ectothermic species with more generations per year owing to warmer ambient temperatures should be more differentiated from each other, accumulate more mutations and show evidence for increased mutation rates compared with populations in colder regions. We used the multivoltine insect species Chironomus riparius to test these predictions with cytochrome oxidase I (COI) sequence data and found that populations from warmer regions are indeed significantly more differentiated and have significantly more derived haplotypes than populations from colder regions. We also found a significant correlation of the annual mean temperature with the population mutation parameter θ that serves as a proxy for the per generation mutation rate under certain assumptions. This pattern could be corroborated with two nuclear loci. Overall, our results support the ESH and indicate that the thermal regime experienced may be crucially driving the evolution of ectotherms and may thus ultimately govern their speciation rate.

Assessing the performance of DNA barcoding using posterior predictive simulations

Accurate estimates of biodiversity are required for research in a broad array of biological subdisciplines including ecology, evolution, systematics, conservation and biodiversity science. The use of statistical models and genetic data, particularly DNA barcoding, has been suggested as an important tool for remedying the large gaps in our current understanding of biodiversity. However, the reliability of biodiversity estimates obtained using these approaches depends on how well the statistical models that are used describe the evolutionary process underlying the genetic data. In this study, we utilize data from the Barcode of Life Database and posterior predictive simulations to assess the performance of DNA barcoding under commonly used substitution models. We demonstrate that the success of DNA barcoding varies widely across DNA substitution models and that model choice has a substantial impact on the number of operational taxonomic units identified (changing results by ~4-31%). Additionally, we demonstrate that the widely followed practice of a priori assuming the Kimura 2-parameter model for DNA barcoding is statistically unjustified and should be avoided. Using both data-based and inference-based test statistics, we detect variation in model performance across taxonomic groups, clustering algorithms, genetic divergence thresholds and substitution models. Taken together, these results illustrate the importance of considering both model selection and model adequacy in studies quantifying biodiversity.

Integrated Taxonomy and DNA Barcoding of Alpine Midges (Diptera: Chironomidae)

Rapid and efficient DNA-based tools are recommended for the evaluation of the insect biodiversity of high-altitude streams. In the present study, focused principally on larvae of the genus Diamesa Meigen 1835 (Diptera: Chironomidae), the congruence between morphological/molecular delimitation of species as well as performances in taxonomic assignments were evaluated. A fragment of the mitochondrial cox1 gene was obtained from 112 larvae, pupae and adults (Diamesinae, Orthocladiinae and Tanypodinae) that were collected in different mountain regions of the Alps and Apennines. On the basis of morphological characters 102 specimens were attributed to 16 species, and the remaining ten specimens were identified to the genus level. Molecular species delimitation was performed using: i) distance-based Automatic Barcode Gap Discovery (ABGD), with no a priori assumptions on species identification; and ii) coalescent tree-based approaches as the Generalized Mixed Yule Coalescent model, its Bayesian implementation and Bayesian Poisson Tree Processes. The ABGD analysis, estimating an optimal intra/interspecific nucleotide distance threshold of 0.7%-1.4%, identified 23 putative species; the tree-based approaches, identified between 25-26 entities, provided nearly identical results. All species belonging to zernyi, steinboecki, latitarsis, bertrami, dampfi and incallida groups, as well as outgroup species, are recovered as separate entities, perfectly matching the identified morphospecies. In contrast, within the cinerella group, cases of discrepancy arose: i) the two morphologically separate species D. cinerella and D. tonsa are neither monophyletic nor diagnosable exhibiting low values of between-taxa nucleotide mean divergence (0.94%); ii) few cases of larvae morphological misidentification were observed. Head capsule color is confirmed to be a valid character able to discriminate larvae of D. zernyi, D. tonsa and D. cinerella, but it is here better defined as a color gradient between the setae submenti and genal setae. DNA barcodes performances were high: average accuracy was ~89% and precision of ~99%. On the basis of the present data, we can thus conclude that molecular identification represents a promising tool that could be effectively adopted in evaluating biodiversity of high-altitude streams.

Morphological identification and COI barcodes of adult flies help determine species identities of chironomid larvae (Diptera, Chironomidae)

Establishing reliable methods for the identification of benthic chironomid communities is important due to their significant contribution to biomass, ecology and the aquatic food web. Immature larval specimens are more difficult to identify to species level by traditional morphological methods than their fully developed adult counterparts, and few keys are available to identify the larval species. In order to develop molecular criteria to identify species of chironomid larvae, larval and adult chironomids from Western Lake Erie were subjected to both molecular and morphological taxonomic analysis. Mitochondrial cytochrome c oxidase I (COI) barcode sequences of 33 adults that were identified to species level by morphological methods were grouped with COI sequences of 189 larvae in a neighbor-joining taxon-ID tree. Most of these larvae could be identified only to genus level by morphological taxonomy (only 22 of the 189 sequenced larvae could be identified to species level). The taxon-ID tree of larval sequences had 45 operational taxonomic units (OTUs, defined as clusters with >97% identity or individual sequences differing from nearest neighbors by >3%; supported by analysis of all larval pairwise differences), of which seven could be identified to species or 'species group' level by larval morphology. Reference sequences from the GenBank and BOLD databases assigned six larval OTUs with presumptive species level identifications and confirmed one previously assigned species level identification. Sequences from morphologically identified adults in the present study grouped with and further classified the identity of 13 larval OTUs. The use of morphological identification and subsequent DNA barcoding of adult chironomids proved to be beneficial in revealing possible species level identifications of larval specimens. Sequence data from this study also contribute to currently inadequate public databases relevant to the Great Lakes region, while the neighbor-joining analysis reported here describes the application and confirmation of a useful tool that can accelerate identification and bioassessment of chironomid communities.

Similarity in the difference: changes in community functional features along natural and anthropogenic stress gradients

The effect of stressors on biodiversity can vary in relation to the degree to which biological communities have adapted over evolutionary time. We compared the responses of functional features of stream insect communities along chronic stress gradients with contrasting time persistence. Water salinity and land use intensification were used as examples of natural (long-term persistent) and anthropogenic (short-term persistent) stressors, respectively. A new trait-based approach was applied to quantify functional diversity components and functional redundancy within the same multidimensional space, using metrics at the taxon and community levels. We found similar functional responses along natural and anthropogenic stress gradients. In both cases, the mean taxon functional richness and functional similarity between taxa increased with stress, whereas community functional richness and functional redundancy decreased. Despite the differences in evolutionary persistence, both chronic stressors act as strong nonrandom environmental filters, producing convergent functional responses. These results can improve our ability to predict functional effects of novel stressors at ecoloiical and evolutionary scales.

Establishing a community-wide DNA barcode library as a new tool for arctic research

DNA sequences offer powerful tools for describing the members and interactions of natural communities. In this study, we establish the to-date most comprehensive library of DNA barcodes for a terrestrial site, including all known macroscopic animals and vascular plants of an intensively studied area of the High Arctic, the Zackenberg Valley in Northeast Greenland. To demonstrate its utility, we apply the library to identify nearly 20 000 arthropod individuals from two Malaise traps, each operated for two summers. Drawing on this material, we estimate the coverage of previous morphology-based species inventories, derive a snapshot of faunal turnover in space and time and describe the abundance and phenology of species in the rapidly changing arctic environment. Overall, 403 terrestrial animal and 160 vascular plant species were recorded by morphology-based techniques. DNA barcodes (CO1) offered high resolution in discriminating among the local animal taxa, with 92% of morphologically distinguishable taxa assigned to unique Barcode Index Numbers (BINs) and 93% to monophyletic clusters. For vascular plants, resolution was lower, with 54% of species forming monophyletic clusters based on barcode regions rbcLa and ITS2. Malaise catches revealed 122 BINs not detected by previous sampling and DNA barcoding. The insect community was dominated by a few highly abundant taxa. Even closely related taxa differed in phenology, emphasizing the need for species-level resolution when describing ongoing shifts in arctic communities and ecosystems. The DNA barcode library now established for Zackenberg offers new scope for such explorations, and for the detailed dissection of interspecific interactions throughout the community.

Two Mitochondrial Barcodes for one Biological Species: The Case of European Kuhl's Pipistrelles (Chiroptera)

The Kuhl's pipistrelle (Pipistrellus kuhlii) is a Western Palaearctic species of bat that exhibits several deeply divergent mitochondrial lineages across its range. These lineages could represent cryptic species or merely ancient polymorphism, but no nuclear markers have been studied so far to properly assess the taxonomic status of these lineages. We examined here two lineages occurring in Western Europe, and used both mitochondrial and nuclear markers to measure degrees of genetic isolation between bats carrying them. The sampling focused on an area of strict lineage sympatry in Switzerland but also included bats from further south, in North Africa. All individuals were barcoded for the COI gene to identify their mitochondrial lineages and five highly polymorphic microsatellite loci were used to cluster them according to their nuclear genotypes. Despite this low number of nuclear markers, all North African nuclear genotypes were grouped in a highly distinct subpopulation when compared with European samples sharing the same mitochondrial barcodes. The reverse situation prevailed in Switzerland where bats carrying distinct barcodes had similar nuclear genotypes. There was a weak east/west nuclear structure of populations, but this was independent of mitochondrial lineages as bats carrying either variant were completely admixed. Thus, the divergent mitochondrial barcodes present in Western Europe do not represent cryptic species, but are part of a single biological species. We argue that these distinct barcodes evolved in allopatry and came recently into secondary contact in an area of admixture north of the Alps. Historical records from this area and molecular dating support such a recent bipolar spatial expansion. These results also highlight the need for using appropriate markers before claiming the existence of cryptic species based on highly divergent barcodes.

How does molecular-assisted identification affect our estimation of α, β and γ biodiversity? An example from understory red seaweeds (Rhodophyta) of Laminaria kelp forests in Brittany, France

Using two distinct identification methods, one based on morphological characters only and the other combining morphological and molecular characters (integrative identification method), we investigated the differences in the biodiversity patterns of red seaweed communities associated with kelp forests at various spatial scales: the regional diversity of Brittany, France (γ-diversity), the local diversity at different Breton sites (α-diversity) and the differentiation in species diversity and abundances among those sites (β-diversity). To characterise α and β diversities, we conducted an initial survey in winter 2011 at 20 sites belonging to four different sub-regions, with specimens collected from six quadrats of 0.10 m(2) at each site, three in the tidal zone dominated by Laminaria digitata and three in the zone dominated by Laminaria hyperborea. To further characterise the regional diversity, we carried out another survey combining several sampling methods (quadrats and visual census) in different seasons (winter, spring and summer) and different years (2011 and 2012). In all, we collected 1990 specimens that were assigned to 76 taxa with the identification method based on morphological characters and 139 taxa using the integrative method. For γ and α diversity, the use of molecular characters revealed several cases of cryptic diversity and both increased the number of identified taxa and improved their taxonomic resolution. However, the addition of molecular characters for specimen identification only slightly affected estimates of β-diversity.

The Use of the Developmental Rate of the Aquatic Midge Chironomus riparius (Diptera, Chironomidae) in the Assessment of the Postsubmersion Interval

Nonbiting midges (Diptera, Chironomidae) are the most abundant members of the fauna associated with submerged carcasses, but their use in the medicolegal context is very restricted because of their complex ontogeny. In this case, the corpse of a woman was recovered in late spring from a river in Granada (Iberian Peninsula). It showed obvious signs of long permanence in the aquatic environment and, along with pulmonary and microscopical analyses, led to the conclusion that the cause of death was drowning. Several larvae-like specimens were sampled from the scalp and later identified by morphological external features as IV instar larvae of Chironomus riparius Meigen, 1804 (Diptera, Chironomidae). Sequencing of cytochrome oxidase subunit I was performed to confirm the identification. The knowledge of the biology of C. riparius at low temperatures was critical to assess a postsubmersion interval of 16-17 days.

A rapid method of species identification of wild chironomids (Diptera: Chironomidae) via electrophoresis of hemoglobin proteins in sodium dodecyl sulfate polyacrylamide gel (SDS-PAGE)

Studying aquatic benthic macroinvertebrates (BMIs) in the field requires accurate taxonomic identification, which can be difficult and time consuming. Conventionally, head capsule morphology has been used to identify wild larvae of Chironomidae. However, due to the number of species and possible damage and/or deformity of their head capsules, another supporting approach for identification is needed. Here, we provide hemoglobin (Hb) protein in hemolymph of chironomids as a new biomarker that may help resolve some of the ambiguities and difficulties encountered during taxonomic identification. Chironomids collected from two locations in Maine and New Jersey, USA were identified to the genus level and in some cases to the species-level using head capsule and body morphologies. The head capsule for a particular individual was then associated with a corresponding Hb protein profile generated from sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE). Distinct Hb profiles were observed from one group (Thienemannimyia) and four genera (Chironomus, Cricotopus, Dicrotendipes, and Glyptotendipes) of chironomids. Several species were polymorphic, having more than one Hb profile and/or having bands of the same size as those of other species. However, major bands and the combination of bands could distinguish individuals at the genus and sometimes species-level. Overall, this study showed that Hb profiles can be used in combination with head capsule morphology to identify wild chironomids.

'Direct PCR' optimization yields a rapid, cost-effective, nondestructive and efficient method for obtaining DNA barcodes without DNA extraction

Macroinvertebrates that are collected in large numbers pose major problems in basic and applied biodiversity research: identification to species via morphology is often difficult, slow and/or expensive. DNA barcodes are an attractive alternative or complementary source of information. Unfortunately, obtaining DNA barcodes from specimens requires many steps and thus time and money. Here, we promote a short cut to DNA barcoding, that is, a nondestructive PCR method that skips DNA extraction ('direct PCR') and that can be used for a broad range of invertebrate taxa. We demonstrate how direct PCR can be optimized for the larvae and adults of nonbiting midges (Diptera: Chironomidae), a typical invertebrate group that is abundant, contains important bioindicator species, but is difficult to identify based on morphological features. After optimization, direct PCR yields high PCR success rates (>90%), preserves delicate morphological features (e.g. details of genitalia, and larval head capsules) while allowing for the recovery of genomic DNA. We also document that direct PCR can be successfully optimized for a wide range of other invertebrate taxa that need routine barcoding (flies: Culicidae, Drosophilidae, Dolichopodidae, Sepsidae; sea stars: Oreasteridae). Key for obtaining high PCR success rates is optimizing (i) tissue quantity, (ii) body part, (iii) primer pair and (iv) type of Taq polymerase. Unfortunately, not all invertebrates appear suitable because direct PCR has low success rates for other taxa that were tested (e.g. Coleoptera: Dytiscidae, Copepoda, Hymenoptera: Formicidae and Odonata). It appears that the technique is less successful for heavily sclerotized insects and/or those with many exocrine glands.

Global haplotype analysis of the whitefly Bemisia tabaci cryptic species Asia I in Asia

The whitefly, Bemisia tabaci (Hemiptera: Aleyrodidiae), is a cryptic species complex comprising a minimum of 24 cryptic species. Some members of this complex are important agricultural pests, causing considerable damage to vegetable as well as ornamental and horticultural crops. Asia I, one of the cryptic species of B. tabaci, is widely distributed in Asia. One hundred and sixty mitochondrial cytochrome oxidase I (COI) sequences from eight countries have been analyzed to investigate the geographic origin and current genetic structure of this cryptic species. Sixty different haplotypes were identified, with levels of genetic distances ranging from 0.001 to 0.021. A sign of possible genetic differentiation emerges from the differential distribution of dominant haplotypes in Indonesia and India compared to China. A possible ancient separation between Asia I in India and Indonesia and secondary contact in China has been hypothesized.

Testing mitochondrial marker efficacy for DNA barcoding in spiders: a test case using the dwarf spider genus Oedothorax (Araneae : Linyphiidae : Erigoninae)

The present study focusses on comparatively assessing the efficacy for DNA barcoding of the two most commonly used mitochondrial markers (cox1 and 16S) in a genus of erigonine spiders. In total, 53 specimens representing five species, including four multi-sampled species, were sampled from several European localities. Initial evaluation of species monophyly was performed through parsimony and Bayesian phylogenetic analyses. Efficacy of mitochondrial markers was tested using operational (including distance-, tree-based measures and Barcode Gap) and evolutionary criteria (using the General Mixed Yule-coalescent Model) for species delimitation. We propose that the cox1 marker can potentially overestimate analyses of biodiversity and thus might not be the preferred marker for DNA species identification and delimitation methods in Oedothorax. Instead, our results suggest that the 16S marker appears to be a promising candidate for such endeavour. Evaluating the contribution and suitability of markers to the re-identification of species, measured by their recovery of well established morphological species, is critical for future studies and for reliable results in species identification in spiders.

Environmental monitoring using next generation sequencing: rapid identification of macroinvertebrate bioindicator species

Introduction

Invertebrate communities are central to many environmental monitoring programs. In freshwater ecosystems, aquatic macroinvertebrates are collected, identified and then used to infer ecosystem condition. Yet the key step of species identification is often not taken, as it requires a high level of taxonomic expertise, which is lacking in most organizations, or species cannot be identified as they are morphologically cryptic or represent little known groups. Identifying species using DNA sequences can overcome many of these issues; with the power of next generation sequencing (NGS), using DNA sequences for routine monitoring becomes feasible.

Results

In this study, we test if NGS can be used to identify species from field-collected samples in an important bioindicator group, the Chironomidae. We show that Cytochrome oxidase I (COI) and Cytochrome B (CytB) sequences provide accurate DNA barcodes for chironomid species. We then develop a NGS analysis pipeline to identifying species using megablast searches of high quality sequences generated using 454 pyrosequencing against comprehensive reference libraries of Sanger-sequenced voucher specimens. We find that 454 generated COI sequences successfully identified up to 96% of species in samples, but this increased up to 99% when combined with CytB sequences. Accurate identification depends on having at least five sequences for a species; below this level species not expected in samples were detected. Incorrect incorporation of some multiplex identifiers (MID’s) used to tag samples was a likely cause, and most errors could be detected when using MID tags on forward and reverse primers. We also found a strong quantitative relationship between the number of 454 sequences and individuals showing that it may be possible to estimate the abundance of species from 454 pyrosequencing data.

Conclusions

Next generation sequencing using two genes was successful for identifying chironomid species. However, when detecting species from 454 pyrosequencing data sets it was critical to include known individuals for quality control and to establish thresholds for detecting species. The NGS approach developed here can lead to routine species-level diagnostic monitoring of aquatic ecosystems.

Genomic basis of ecological niche divergence among cryptic sister species of non-biting midges

BACKGROUND

There is a lack of understanding the evolutionary forces driving niche segregation of closely related organisms. In addition, pinpointing the genes driving ecological divergence is a key goal in molecular ecology. Here, larval transcriptome sequences obtained by next-generation-sequencing are used to address these issues in a morphologically cryptic sister species pair of non-biting midges (Chironomus riparius and C. piger).

RESULTS

More than eight thousand orthologous open reading frames were screened for interspecific divergence and intraspecific polymorphisms. Despite a small mean sequence divergence of 1.53% between the sister species, 25.1% of 18,115 observed amino acid substitutions were inferred by α statistics to be driven by positive selection. Applying McDonald-Kreitman tests to 715 alignments of gene orthologues identified eleven (1.5%) genes driven by positive selection.

CONCLUSIONS

Three candidate genes were identified as potentially responsible for the observed niche segregation concerning nitrite concentration, habitat temperature and water conductivity. Additionally, signs of positive selection in the hydrogen sulfide detoxification pathway were detected, providing a new plausible hypothesis for the species' ecological differentiation. Finally, a divergently selected, nuclear encoded mitochondrial ribosomal protein may contribute to reproductive isolation due to cytonuclear coevolution.

trnL-F is a powerful marker for DNA identification of field vittarioid gametophytes (Pteridaceae)

BACKGROUND AND AIMS

The gametophyte phase of ferns plays an important role in habitat selection, dispersal, adaptation and evolution. However, ecological studies on fern gametophytes have been impeded due to the difficulty of species identification of free-living gametophytes. DNA barcoding provides an alternative approach to identifying fern gametophytes but is rarely applied to field studies. In this study, an example of field vittarioid gametophyte identification using DNA barcoding, which has not been done before, is given.

METHODS

A combination of distance-based and tree-based approaches was performed to evaluate the discriminating power of three candidate barcodes (matK, rbcL and trnL-F) on 16 vittarioid sporophytes. Sequences of the trnL-F region were generated from 15 fern gametophyte populations by tissue-direct PCR and were compared against the sporophyte dataset, using BLAST.

KEY RESULTS

trnL-F earns highest primer universality and discriminatory ability scores, whereas PCR success rates were very low for matK and rbcL regions (10·8 % and 41·3 %, respectively). BLAST analyses showed that all the sampled field gametophytes could be successfully identified to species level. Three gametophyte populations were also discovered to be living beyond the known occurrence of their sporophyte counterparts.

CONCLUSIONS

This study demonstrates that DNA barcoding (i.e. reference databasing, tissue-direct PCR and molecular analysis), especially the trnL-F region, is an efficient tool to identify field gametophytes, and has considerable potential in exploring the ecology of fern gametophytes.

Can DNA barcoding accurately discriminate megadiverse Neotropical freshwater fish fauna?

BACKGROUND

The megadiverse Neotropical freshwater ichthyofauna is the richest in the world with approximately 6,000 recognized species. Interestingly, they are distributed among only 17 orders, and almost 80% of them belong to only three orders: Characiformes, Siluriformes and Perciformes. Moreover, evidence based on molecular data has shown that most of the diversification of the Neotropical ichthyofauna occurred recently. These characteristics make the taxonomy and identification of this fauna a great challenge, even when using molecular approaches. In this context, the present study aimed to test the effectiveness of the barcoding methodology (COI gene) to identify the mega diverse freshwater fish fauna from the Neotropical region. For this purpose, 254 species of fishes were analyzed from the Upper Parana River basin, an area representative of the larger Neotropical region.

RESULTS

Of the 254 species analyzed, 252 were correctly identified by their barcode sequences (99.2%). The main K2P intra- and inter-specific genetic divergence values (0.3% and 6.8%, respectively) were relatively low compared with similar values reported in the literature, reflecting the higher number of closely related species belonging to a few higher taxa and their recent radiation. Moreover, for 84 pairs of species that showed low levels of genetic divergence (<2%), application of a complementary character-based nucleotide diagnostic approach proved useful in discriminating them. Additionally, 14 species displayed high intra-specific genetic divergence (>2%), pointing to at least 23 strong candidates for new species.

CONCLUSIONS

Our study is the first to examine a large number of freshwater fish species from the Neotropical area, including a large number of closely related species. The results confirmed the efficacy of the barcoding methodology to identify a recently radiated, megadiverse fauna, discriminating 99.2% of the analyzed species. The power of the barcode sequences to identify species, even with low interspecific divergence, gives us an idea of the distribution of inter-specific genetic divergence in these megadiverse fauna. The results also revealed hidden genetic divergences suggestive of reproductive isolation and putative cryptic speciation in some species (23 candidates for new species). Finally, our study constituted an important contribution to the international Barcoding of Life (iBOL.org) project, providing barcode sequences for use in identification of these species by experts and non-experts, and allowing them to be available for use in other applications.

Improving environmental and biodiversity monitoring in the Baltic Sea using DNA barcoding of Chironomidae (Diptera)

As for many other regions, environmental and biodiversity monitoring of the brackish Baltic Sea suffers from low species resolution for several taxa. One such case is the benthic larvae of midges Chironomidae (Diptera), which are estimated to constitute about 30% of the macrozoobenthos species of the Baltic Sea and are important indicators of environmental quality. We assessed the usefulness of COI (cytochrome oxidase I) gene barcoding to improve species resolution and its potential for implementation in monitoring programmes. Neighbour-Joining, Maximum parsimony and Bayesian-inference analyses all provided high congruency with morphological analyses of adult males for almost all 42 species studied. Barcoding was helpful to elucidate some cases of taxonomical difficulties, such as synonyms. In contrast to the high identification accuracy when using our local database, there were a number of cases where matching with GenBank and BOLD provided puzzling results. For reliable species identification at least 15-30 specimens from 5-10 well-distributed sites within the geographical range of the species might be needed in a database to adequately cover the intraspecific variability of chironomids. Implementation of DNA barcoding, as applied here, in monitoring would result in an increase from at present less than 10% to more than 90% successful chironomid species identification of Baltic Sea benthic samples, as it also would for many nearby lakes. Routine monitoring of benthic environmental samples based on Next-Generation sequencing techniques would provide a cost effective way to obtain a taxonomically much more complete assessment of environmental quality and biodiversity, as required by EU directives and national legislation.

Mitochondrial pseudogenes in insect DNA barcoding: differing points of view on the same issue

Molecular tools have been used in taxonomy for the purpose of identification and classification of living organisms. Among these, a short sequence of the mitochondrial DNA, popularly known as DNA barcoding, has become very popular. However, the usefulness and dependability of DNA barcodes have been recently questioned because mitochondrial pseudogenes, non-functional copies of the mitochondrial DNA incorporated into the nuclear genome, have been found in various taxa. When these paralogous sequences are amplified together with the mitochondrial DNA, they may go unnoticed and end up being analyzed as if they were orthologous sequences. In this contribution the different points of view regarding the implications of mitochondrial pseudogenes for entomology are reviewed and discussed. A discussion of the problem from a historical and conceptual perspective is presented as well as a discussion of strategies to keep these nuclear mtDNA copies out of sequence analyzes.

Integrative taxonomy and molecular phylogeny of genus Aplysina (Demospongiae: Verongida) from Mexican Pacific

Integrative taxonomy provides a major approximation to species delimitation based on integration of different perspectives (e.g. morphology, biochemistry and DNA sequences). The aim of this study was to assess the relationships and boundaries among Eastern Pacific Aplysina species using morphological, biochemical and molecular data. For this, a collection of sponges of the genus Aplysina from the Mexican Pacific was studied on the basis of their morphological, chemical (chitin composition), and molecular markers (mitochondrial COI and nuclear ribosomal rDNA: ITS1-5.8-ITS2). Three morphological species were identified, two of which are new to science. A. clathrata sp. nov. is a yellow to yellow-reddish or -brownish sponge, characterized by external clathrate-like morphology; A. revillagigedi sp. nov. is a lemon yellow to green, cushion-shaped sometimes lobate sponge, characterized by conspicuous oscules, which are slightly elevated and usually linearly distributed on rims; and A. gerardogreeni a known species distributed along the Mexican Pacific coast. Chitin was identified as the main structural component within skeletons of the three species using FTIR, confirming that it is shared among Verongida sponges. Morphological differences were confirmed by DNA sequences from nuclear ITS1-5.8-ITS2. Mitochondrial COI sequences showed extremely low but diagnostic variability for Aplysina revillagigedi sp. nov., thus our results corroborate that COI has limited power for DNA-barcoding of sponges and should be complemented with other markers (e.g. rDNA). Phylogenetic analyses of Aplysina sequences from the Eastern Pacific and Caribbean, resolved two allopatric and reciprocally monophyletic groups for each region. Eastern Pacific species were grouped in general accordance with the taxonomic hypothesis based on morphological characters. An identification key of Eastern Pacific Aplysina species is presented. Our results constitute one of the first approximations to integrative taxonomy, phylogeny and evolutionary biogeography of Eastern Pacific marine sponges; an approach that will significantly contribute to our better understanding of their diversity and evolutionary history.