Design of phylum-specific hybrid primers for DNA barcoding: addressing the need for efficient COI amplification in the Echinodermata

Recent research has shown the usefulness of the Folmer region of the cytochrome oxidase I (COI) as a genetic barcode to assist in species delimitation of echinoderms. However, amplification of COI is often challenging in echinoderms (low success or pseudogenes). We present a method that allows the design of phylum-specific hybrid primers, and use this to develop COI primers for the Echinodermata. We aligned COI sequences from 310 echinoderm species and designed all possible primers along the consensus sequence with two methods (standard degenerate and hybrid). We found much lower degeneracy for hybrid primers (4-fold degeneracy) than for standard degenerate primers (≥48-fold degeneracy). We then designed the most conserved hybrid primers to amplify a >500-bp region within COI. These primers successfully amplified this gene region in all tested taxa (123 species across all echinoderm classes). Sequencing of 30 species among these confirmed both the quality of the sequences (>500 bp, no pseudogenes) and their utility as a DNA barcode. This method should be useful for developing primers for other mitochondrial genes and other phyla. The method will also be of interest for the development of future projects involving both community-based genetic assessments on macroorganisms and biodiversity assessment of environmental samples using high-throughput sequencing.     

Development of internal COI primers to improve and extend barcoding of fruit flies (Diptera: Tephritidae: Dacini)

Accurate species-level identifications underpin many aspects of basic and applied biology;however,identifications can be hampered by a lack of discriminating morphological characters,taxonomic expertise or time.Molecular approaches,such as DNA"barcoding"of the cytochrome c oxidase(COI)gene,are argued to overcome these issues.However,nuclear encoding of mitochondrial genes(numts)and poor amplification success of suboptimally preserved specimens can lead to erroneous identifications.One insect group for which these molecular and morphological problems are significant are the dacine fruit flies(Diptera:Tephritidae:Dacini).We addressed these issues associated with COI barcoding in the dacines by first assessing several"universal"COI primers against public mitochondrial genome and numt sequences for dacine taxa.We then modified a set of four primers that more closely matched true dacine COI sequence and amplified two overlapping portions of the COI barcode region.Our new primers were tested alongside universal primers on a selection of dacine species,including both fresh preserved and decades-old dry specimens.Additionally,Bactrocera tiyoni mitochondrial and nuclear genomes were compared to identify putative numts.Four numt clades were identified,three of which were amplified using existing universal primers.In contrast,our new primers preferentially amplified the"true"mitochondrial COI barcode in all dacine species tested.The new primers also successfully amplified partial barcodes from dry specimens for which full length barcodes were unobtainable.Thus we recommend these new primers be incorporated into the suites of primers used by diagnosticians and quarantine labs for the accurate identification of dacine species.     

Environmental metabarcodes for insects: in silico PCR reveals potential for taxonomic bias

Studies of insect assemblages are suited to the simultaneous DNA‐based identification of multiple taxa known as metabarcoding. To obtain accurate estimates of diversity, metabarcoding markers ideally possess appropriate taxonomic coverage to avoid PCR‐amplification bias, as well as sufficient sequence divergence to resolve species. We used in silico PCR to compare the taxonomic coverage and resolution of newly designed insect metabarcodes (targeting 16S) with that of existing markers [16S and cytochrome oxidase c subunit I (COI)] and then compared their efficiency in vitro. Existing metabarcoding primers amplified in silico <75% of insect species with complete mitochondrial genomes available, whereas new primers targeting 16S provided >90% coverage. Furthermore, metabarcodes targeting COI appeared to introduce taxonomic PCR‐amplification bias, typically amplifying a greater percentage of Lepidoptera and Diptera species, while failing to amplify certain orders in silico. To test whether bias predicted in silico was observed in vitro, we created an artificial DNA blend containing equal amounts of DNA from 14 species, representing 11 insect orders and one arachnid. We PCR‐amplified the blend using five primer sets, targeting either COI or 16S, with high‐throughput amplicon sequencing yielding more than 6 million reads. In vitro results typically corresponded to in silico PCR predictions, with newly designed 16S primers detecting 11 insect taxa present, thus providing equivalent or better taxonomic coverage than COI metabarcodes. Our results demonstrate that in silico PCR is a useful tool for predicting taxonomic bias in mixed template PCR and that researchers should be wary of potential bias when selecting metabarcoding markers.     

Performance comparison of genetic markers for high‐throughput sequencing‐based biodiversity assessment in complex communities

Metabarcode surveys of DNA extracted from environmental samples are increasingly popular for biodiversity assessment in natural communities. Such surveys rely heavily on robust genetic markers. Therefore, analysis of PCR efficiency and subsequent biodiversity estimation for different types of genetic markers and their corresponding primers is important. Here, we test the PCR efficiency and biodiversity recovery potential of three commonly used genetic markers – nuclear small subunit ribosomal DNA (18S), mitochondrial cytochrome c oxidase subunit I (COI) and 16S ribosomal RNA (mt16S) – using 454 pyrosequencing of a zooplankton community collected from Hamilton Harbour, Ontario. We found that biodiversity detection power and PCR efficiency varied widely among these markers. All tested primers for COI failed to provide high‐quality PCR products for pyrosequencing, but newly designed primers for 18S and 16S passed all tests. Furthermore, multiple analyses based on large‐scale pyrosequencing (i.e. 1/2 PicoTiter plate for each marker) showed that primers for 18S recover more (38 orders) groups than 16S (10 orders) across all taxa, and four vs. two orders and nine vs. six families for Crustacea. Our results showed that 18S, using newly designed primers, is an efficient and powerful tool for profiling biodiversity in largely unexplored communities, especially when amplification difficulties exist for mitochondrial markers such as COI. Universal primers for higher resolution markers such as COI are still needed to address the possible low resolution of 18S for species‐level identification.     

Effect of marker choice and thermal cycling protocol on zooplankton DNA metabarcoding studies

DNA metabarcoding is a promising approach for rapidly surveying biodiversity and is likely to become an important tool for measuring ecosystem responses to environmental change. Metabarcoding markers need sufficient taxonomic coverage to detect groups of interest, sufficient sequence divergence to resolve species, and will ideally indicate relative abundance of taxa present. We characterized zooplankton assemblages with three different metabarcoding markers (nuclear 18S rDNA, mitochondrial COI, and mitochondrial 16S rDNA) to compare their performance in terms of taxonomic coverage, taxonomic resolution, and correspondence between morphology‐ and DNA‐based identification. COI amplicons sequenced on separate runs showed that operational taxonomic units representing >0.1% of reads per sample were highly reproducible, although slightly more taxa were detected using a lower annealing temperature. Mitochondrial COI and nuclear 18S showed similar taxonomic coverage across zooplankton phyla. However, mitochondrial COI resolved up to threefold more taxa to species compared to 18S. All markers revealed similar patterns of beta‐diversity, although different taxa were identified as the greatest contributors to these patterns for 18S. For calanoid copepod families, all markers displayed a positive relationship between biomass and sequence reads, although the relationship was typically strongest for 18S. The use of COI for metabarcoding has been questioned due to lack of conserved primer‐binding sites. However, our results show the taxonomic coverage and resolution provided by degenerate COI primers, combined with a comparatively well‐developed reference sequence database, make them valuable metabarcoding markers for biodiversity assessment.     

Comparing diversity levels in environmental samples: DNA sequence capture and metabarcoding approaches using 18S and COI genes

Environmental DNA studies targeting multiple taxa using metabarcoding provide remarkable insights into levels of species diversity in any habitat. The main drawbacks are the presence of primer bias and difficulty in identifying rare species. We tested a DNA sequence‐capture method in parallel with the metabarcoding approach to reveal possible advantages of one method over the other. Both approaches were performed using the same eDNA samples and the same 18S and COI regions, followed by high throughput sequencing. Metabarcoded eDNA libraries were PCR amplified with one primer pair from 18S and COI genes. DNA sequence‐capture libraries were enriched with 3,639 baits targeting the same gene regions. We tested amplicon sequence variants (ASVs) and operational taxonomic units (OTUs) in silico approaches for both markers and methods, using for this purpose the metabarcoding data set. ASVs methods uncovered more species for the COI gene, whereas the opposite occurred for the 18S gene, suggesting that clustering reads into OTUs could bias diversity richness especially using 18S with relaxed thresholds. Additionally, metabarcoding and DNA sequence‐capture recovered 80%–90% of the control sample species. DNA sequence‐capture was 8x more expensive, nonetheless it identified 1.5x more species for COI and 13x more genera for 18S than metabarcoding. Both approaches offer reliable results, sharing ca. 40% species and 72% families and retrieve more taxa when nuclear and mitochondrial markers are combined. eDNA metabarcoding is quite well established and low‐cost, whereas DNA‐sequence capture for biodiversity assessment is still in its infancy, is more time‐consuming but provides more taxonomic assignments.     

Meeting the challenge of DNA barcoding Neotropical amphibians: polymerase chain reaction optimization and new COI primers

Amphibians are one of the most threatened vertebrate classes, yet at the same time new species are being described every year, demonstrating that the number of existing species is grossly underestimated. In groups such as amphibians, with high extinction rates and poorly known species boundaries, DNA barcoding is a tool that can rapidly assess genetic diversity and estimate species richness for prioritizing conservation decisions. However, reliable recovery of the 5′ region of the cytochrome c oxidase subunit 1 (COI) gene is critical for the ongoing effort to gather DNA barcodes for all amphibian species. Here, we provide new PCR conditions and tested new primers that increase the efficiency of barcode recovery in amphibians. We found that a low extension temperature for PCR cycles significantly improves the efficiency of amplification for all combinations of primers. Combining low PCR extension temperature and primers AnF1 + AnR1, we were able to recover COI sequences for 100% of the species analysed (N = 161), encompassing ~15% of the species known from Brazil (representing 77 genera and 23 families), which is an important improvement over previous studies. The preliminary assessment of species diversity suggested that number of species might be underestimated by about 25%. We conclude that DNA barcoding is an efficient, simple, and standardized protocol for identifying cryptic diversity in amphibians and advocate for its use in biodiversity inventories and across widespread populations within known species.     

High mitochondrial diversity in geographically widespread butterflies of Madagascar: A test of the DNA barcoding approach

The standardized use of mitochondrial cytochrome c oxidase subunit I (COI) gene sequences as DNA barcodes has been widely promoted as a high-throughput method for species identification and discovery. Species delimitation has been based on the following criteria: (1) monophyletic association and less frequently (2) a minimum 10× greater divergence between than within species. Divergence estimates, however, can be inflated if sister species pairs are not included and the geographic extent of variation within any given taxon is not sampled comprehensively. This paper addresses both potential biases in DNA divergence estimation by sampling range-wide variation in several morphologically distinct, endemic butterfly species in the genus Heteropsis, some of which are sister taxa. We also explored the extent to which mitochondrial DNA from the barcode region can be used to assess the effects of historical rainforest fragmentation by comparing genetic variation across Heteropsis populations with an unrelated forest-associated taxon Saribia tepahi. Unexpectedly, generalized primers led to the inadvertent amplification of the endosymbiont Wolbachia, undermining the use of universal primers and necessitating the design of genus-specific COI primers alongside a Wolbachia-specific PCR assay. Regardless of the high intra-specific genetic variation observed, most species satisfy DNA barcoding criteria and can be differentiated in the nuclear phylogeny. Nevertheless, two morphologically distinguishable candidate species fail to satisfy the barcoding 10× genetic distance criterion, underlining the difficulties of applying a standard distance threshold to species delimitation. Phylogeographic analysis of COI data suggests that forest fragmentation may have played an important role in the recent evolutionary diversification of these butterflies. Further work on other Malagasy taxa using both mitochondrial and nuclear data will provide better insight into the role of historical habitat fragmentation in species diversification and may potentially contribute to the identification of priority areas for conservation.     

Escaping introns in COI through cDNA barcoding of mushrooms: Pleurotus as a test case

DNA barcoding involves the use of one or more short, standardized DNA fragments for the rapid identification of species. A 648‐bp segment near the 5′ terminus of the mitochondrial cytochrome c oxidase subunit I (COI) gene has been adopted as the universal DNA barcode for members of the animal kingdom, but its utility in mushrooms is complicated by the frequent occurrence of large introns. As a consequence, ITS has been adopted as the standard DNA barcode marker for mushrooms despite several shortcomings. This study employed newly designed primers coupled with cDNA analysis to examine COI sequence diversity in six species of Pleurotus and compared these results with those for ITS. The ability of the COI gene to discriminate six species of Pleurotus, the commonly cultivated oyster mushroom, was examined by analysis of cDNA. The amplification success, sequence variation within and among species, and the ability to design effective primers was tested. We compared ITS sequences to their COI cDNA counterparts for all isolates. ITS discriminated between all six species, but some sequence results were uninterpretable, because of length variation among ITS copies. By comparison, a complete COI sequences were recovered from all but three individuals of Pleurotus giganteus where only the 5′ region was obtained. The COI sequences permitted the resolution of all species when partial data was excluded for P. giganteus. Our results suggest that COI can be a useful barcode marker for mushrooms when cDNA analysis is adopted, permitting identifications in cases where ITS cannot be recovered or where it offers higher resolution when fresh tissue is. The suitability of this approach remains to be confirmed for other mushrooms.     

DNA metabarcoding and the cytochrome c oxidase subunit I marker: not a perfect match

DNA metabarcoding enables efficient characterization of species composition in environmental DNA or bulk biodiversity samples, and this approach is making significant and unique contributions in the field of ecology. In metabarcoding of animals, the cytochrome c oxidase subunit I (COI) gene is frequently used as the marker of choice because no other genetic region can be found in taxonomically verified databases with sequences covering so many taxa. However, the accuracy of metabarcoding datasets is dependent on recovery of the targeted taxa using conserved amplification primers. We argue that COI does not contain suitably conserved regions for most amplicon-based metabarcoding applications. Marker selection deserves increased scrutiny and available marker choices should be broadened in order to maximize potential in this exciting field of research.     

DNA barcoding reveals 24 distinct lineages as cryptic bird species candidates in and around the Japanese Archipelago

DNA barcoding using a partial region (648 bp) of the cytochrome c oxidase I (COI) gene is a powerful tool for species identification and has revealed many cryptic species in various animal taxa. In birds, cryptic species are likely to occur in insular regions like the Japanese Archipelago due to the prevention of gene flow by sea barriers. Using COI sequences of 234 of the 251 Japanese‐breeding bird species, we established a DNA barcoding library for species identification and estimated the number of cryptic species candidates. A total of 226 species (96.6%) had unique COI sequences with large genetic divergence among the closest species based on neighbour‐joining clusters, genetic distance criterion and diagnostic substitutions. Eleven cryptic species candidates were detected, with distinct intraspecific deep genetic divergences, nine lineages of which were geographically separated by islands and straits within the Japanese Archipelago. To identify Japan‐specific cryptic species from trans‐Paleartic birds, we investigated the genetic structure of 142 shared species over an extended region covering Japan and Eurasia; 19 of these species formed two or more clades with high bootstrap values. Excluding six duplicated species from the total of 11 species within the Japanese Archipelago and 19 trans‐Paleartic species, we identified 24 species that were cryptic species candidates within and surrounding the Japanese Archipelago. Repeated sea level changes during the glacial and interglacial periods may be responsible for the deep genetic divergences of Japanese birds in this insular region, which has led to inconsistencies in traditional taxonomies based on morphology.     

Performance of amplicon and shotgun sequencing for accurate biomass estimation in invertebrate community samples

New applications of DNA and RNA sequencing are expanding the field of biodiversity discovery and ecological monitoring, yet questions remain regarding precision and efficiency. Due to primer bias, the ability of metabarcoding to accurately depict biomass of different taxa from bulk communities remains unclear, while PCR‐free whole mitochondrial genome (mitogenome) sequencing may provide a more reliable alternative. Here, we used a set of documented mock communities comprising 13 species of freshwater macroinvertebrates of estimated individual biomass, to compare the detection efficiency of COI metabarcoding (three different amplicons) and shotgun mitogenome sequencing. Additionally, we used individual COI barcoding and de novo mitochondrial genome sequencing, to provide reference sequences for OTU assignment and metagenome mapping (mitogenome skimming), respectively. We found that, even though both methods occasionally failed to recover very low abundance species, metabarcoding was less consistent, by failing to recover some species with higher abundances, probably due to primer bias. Shotgun sequencing results provided highly significant correlations between read number and biomass in all but one species. Conversely, the read–biomass relationships obtained from metabarcoding varied across amplicons. Specifically, we found significant relationships for eight of 13 (amplicons B1FR‐450 bp, FF130R‐130 bp) or four of 13 (amplicon FFFR, 658 bp) species. Combining the results of all three COI amplicons (multiamplicon approach) improved the read–biomass correlations for some of the species. Overall, mitogenomic sequencing yielded more informative predictions of biomass content from bulk macroinvertebrate communities than metabarcoding. However, for large‐scale ecological studies, metabarcoding currently remains the most commonly used approach for diversity assessment.     

“Forms” of water mites (Acari: Hydrachnidia): intraspecific variation or valid species?

In many groups of organisms, especially in the older literature, it has been common practice to recognize sympatrically occurring phenotypic variants of a species as “forms”. However, what these forms really represent often remains unclear, especially in poorly studied groups. With new algorithms for DNA‐based species delimitation, the status of forms can be explicitly tested with molecular data. In this study, we test a number of what is now recognized as valid species of water mites (Hydrachnidia), but have in the past been treated as forms sympatrically occurring with their nominate species. We also test a form without prior taxonomical status, using DNA and morphometrics. The barcoding fragment of COI, nuclear 28S and quantitative analyses of morphological data were used to test whether these taxa merit species status, as suggested by several taxonomists. Our results confirm valid species. Genetic distances between the form and nominate species (Piona dispersa and Piona variabilis, COI 11%), as well as likelihood ratio tests under the general mixed‐Yule coalescent model, supported that these are separately evolving lineages as defined by the unified species concept. In addition, they can be diagnosed with morphological characters. The study also reveals that some taxa genetically represent more than one species. We propose that P. dispersa are recognized as valid taxa at the species level. Unionicola minor (which may consist of several species), Piona stjordalensis, P. imminuta s. lat., and P. rotundoides are confirmed as species using this model. The results also imply that future studies of other water mite species complexes are likely to reveal many more genetically and morphologically distinct species.     

Molecular Phylogeny and Systematics of Anoplocephaline Cestodes in Rodents and Lagomorphs

A molecular phylogenetic hypothesis is presented for the anoplocephaline cestodes of placental mammals based on sequence data from the mitochondrial cytochrome c oxidase I (COI) gene, the nuclear-encoded 28S rRNA gene and the internal transcribed spacer region I of rRNA (ITS1). The material consists of 35 species representing nine genera of cestodes, with emphasis on taxa parasitising rodents and lagomorphs in the Holarctic region. The resulting phylogenies show considerable disagreement with earlier systematic and phylogenetic hypotheses derived from morphology. Specifically, the results contradict the view of uterine morphology being the primary determinant of deeper phylogenetic splits within Anoplocephalinae. Also, the role of genital duplication as a means of generic divergence was not found to follow consistently the pattern suggested by earlier hypotheses. Colonisation of novel host lineages has evidently been the predominant mode of diversification in anoplocephaline cestodes of placental mammals; evidence for phyletic co-evolution was obscure. The phylogenies consistently distinguished a large monophyletic group including all species from arvicoline rodents (voles and lemmings), primarily representing the genera Anoplocephaloides Baer, 1923 and Paranoplocephala Lühe, 1910. Phylogenetic relationships within the “arvicoline clade” of cestodes were generally poorly resolved. Consistent support for nodes above and below the unresolved polytomy indicates a rapid radiation involving a nearly simultaneous diversification of many lineages, a scenario also proposed for the arvicoline hosts.     

A review of cross‐backed grasshoppers of the genus Dociostaurus Fieber (Orthoptera: Acrididae) from the western Mediterranean: insights from phylogenetic analyses and DNA‐based species delimitation

Phylogenetic analyses and species delimitation methods are powerful tools for understanding patterns of species diversity. Given the current biodiversity crisis, such approaches are invaluable for urgent assessment and delimitation of truthful species, particularly of endangered and morphologically cryptic taxa from vulnerable areas submitted to strong climate change and progressive human intervention such as the M editerranean region. In this study, we applied two DNA ‐based species delimitation methods and performed a B ayesian phylogenetic reconstruction using three mitochondrial gene fragments (12S , 16S and COI) to solve several taxonomic uncertainties among species of cross‐backed grasshoppers (genus Dociostaurus F ieber) from the western M editerranean. P hylogenetic analyses demonstrate the polyphyletic character of subgenera Dociostaurus , Kazakia B ey‐B ienko and Stauronotulus T arbinsky and, thus, the necessity of revising the currently accepted taxonomic subgenera within the genus Dociostaurus . We propose the split of closely related taxa with allopatric distributions such as D. (S.) kraussi and D. (S.) crassiusculus , considering the later a distinct species limited to the I berian P eninsula and excluding the name crassiusculus from other forms of D. (S.) kraussi from E ast E urope and A sia. Estimates of divergence times indicate that diversification of Dociostaurus probably happened during the M iocene–P liocene (3–7 Ma), and the split of the studied pairs of sister taxa took place during the middle and late P leistocene (1–2 Ma). This study highlights the need for more molecular studies on the genus and their different species for a better understanding of their evolution, genetic variation and population dynamics in order to prioritize strategies for their adequate conservation and management.     

Evaluation of species boundaries in sympatric and parapatric populations of Mesoamerican toads

Approaches that integrate multiple independent, yet complimentary, lines of evidence have been effectively utilized to identify and evaluate species diversity. Integrative approaches are especially useful in taxa that exhibit cryptic diversity and are highly morphologically conserved, as well as organisms whose distributions may be sympatric or parapatric. The Incilius coccifer complex in Honduras is comprised of three putative taxa: I. coccifer, I. ibarrai and I. porteri. The taxonomy of the I. coccifer complex has been a source of debate among specialists, with some recognizing three species, while others choose to recognize one widespread taxon. To assess species boundaries and evaluate the taxonomic structure for the I. coccifer complex, we utilized a combination of comprehensive field sampling, molecular phylogenetics and macroecological modelling. Using 58 samples representing all three putative taxa, we generated sequence data from the mitochondrial loci 16S and COI in order to assess genetic diversity and phylogenetic relationships, and tested putative species boundaries using General Mixed Yule‐Coalescent models. To evaluate macroecological differences in the distribution of putative taxa, we utilized maximum entropy modelling and identified areas of suitable and non‐suitable habitat, as well as identifying potential areas of overlap between species habitats. We recovered three clades that broadly correspond to the three named taxa that, while being monophyletic, are separated by relatively small genetic distances. Species distribution models revealed that I. coccifer is macroecologically different than the other two taxa, but that I. ibarrai and I. porteri are highly similar. We uncovered cases of sympatry between pairs of species in at least three localities in Honduras, suggesting the potential for hybridization in these closely related lineages.     

Phylogenomic species delimitation and host‐symbiont coevolution in the fungus‐farming ant genus Sericomyrmex Mayr (Hymenoptera: Formicidae): ultraconserved elements (UCEs) resolve a recent radiation

Ants in the Neotropical genus Sericomyrmex Mayr cultivate fungi for food. Both ants and fungi are obligate, coevolved symbionts. The taxonomy of Sericomyrmex is problematic because the morphology of the worker caste is generally homogeneous across all of the species within the genus, species limits are vague, and the relationships between them are unknown. We used ultraconserved elements (UCEs) as genome‐scale markers to reconstruct evolutionary history and to infer species boundaries in Sericomyrmex. We recovered an average of ～990 UCE loci for 88 Sericomyrmex samples from across the geographical range of the genus as well as for five outgroup taxa. Using maximum likelihood and species‐tree approaches, we recovered nearly identical topologies across datasets with 50–95% matrix completeness. We identify nine species‐level lineages in Sericomyrmex, including two new species. This is less than the previously described 19 species, even accounting for two species for which we had no UCE samples, which brings the total number of Sericomyrmex species to 11. Divergence‐dating analyses recovered 4.3 Ma as the crown‐group age estimates for Sericomyrmex, indicating a recent, rapid radiation. We also sequenced mitochondrial cytochrome oxidase subunit I (COI) for 125 specimens. Resolution and support for clades in our COI phylogeny are weak, indicating that COI is not an appropriate species‐delimitation tool. However, taxa within species consistently cluster together, suggesting that COI is useful as a species identification (‘DNA barcoding’) tool. We also sequenced internal transcribed spacer (ITS) and large subunit (LSU) for 32 Sericomyrmex fungal cultivars. The fungal phylogeny confirms that Sericomyrmex fungi are generalized higher‐attine cultivars, interspersed with Trachymyrmex‐associated fungal species, indicating cultivar sharing and horizontal transfer between these two genera. Our results indicate that UCEs offer immense potential for delimiting and resolving relationships of problematic, recently diverged species.     

Comparing COI and ITS as DNA barcode markers for mushrooms and allies (Agaricomycotina)

DNA barcoding is an approach to rapidly identify species using short, standard genetic markers. The mitochondrial cytochrome oxidase I gene (COI) has been proposed as the universal barcode locus, but its utility for barcoding in mushrooms (ca. 20,000 species) has not been established. We succeeded in generating 167 partial COI sequences (~450 bp) representing ~100 morphospecies from ~650 collections of Agaricomycotina using several sets of new primers. Large introns (~1500 bp) at variable locations were detected in ~5% of the sequences we obtained. We suspect that widespread presence of large introns is responsible for our low PCR success (~30%) with this locus. We also sequenced the nuclear internal transcribed spacer rDNA regions (ITS) to compare with COI. Among the small proportion of taxa for which COI could be sequenced, COI and ITS perform similarly as a barcode. However, in a densely sampled set of closely related taxa, COI was less divergent than ITS and failed to distinguish all terminal clades. Given our results and the wealth of ITS data already available in public databases, we recommend that COI be abandoned in favor of ITS as the primary DNA barcode locus in mushrooms.     

Comparative analysis of the mitochondrial cytochrome c oxidase subunit I (COI) gene in ciliates (Alveolata,Ciliophora) and evaluation of its suitability as a biodiversity marker

The mitochondrial cytochrome c oxidase subunit 1 (COI) gene of ciliates was first successfully sequenced in species of the genera Tetrahymena and Paramecium (Class Oligohymenophorea). The sequence of the COI gene is extremely divergent from other eukaryotes and includes an insert, which is over 300 nucleotides long. In this study, we designed a primer pair that successfully amplified the COI gene of ciliates from five different classes: Heterotrichea, Spirotrichea, Oligohymenophorea, Nassophorea and Colpodea. These classes represent the diversity of the phylum Ciliophora very well, since they are widely distributed on the ciliate small subunit rRNA tree. The amplified region is approximately 850 nucleotides long and corresponds to the general barcoding region; it also includes the insert region. In this study, 58 new COI sequences from over 38 species in 13 orders are analysed and compared, and distance trees are constructed. While the COI gene shows high divergence within ciliates, the insert region, which is present in all classes, is even more divergent. Genetic distances calculated with and without the insert region remain in the same range at the intraspecific level, but they differ considerably at or above genus level. This suggests that the entire barcoding region is under similar selective constraints and that the evolutionary rate of the ciliate COI is extremely high and shows unequal rate variation. Although many problems still remain regarding standardization of barcoding methods in ciliates, the development of a universal or almost universal primer combination for the Phylum Ciliophora represents important progress. As shown in four examples, the resolution of COI at the intraspecific level is much greater than that of any nuclear genes and shows great potential to (1) identify species based on molecular data if a reliable database exists, and (2) resolve the relationships of closely related ciliate taxa and uncover cryptic species.