Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species

With millions of species and their life-stage transformations, the animal kingdom provides a challenging target for taxonomy. Recent work has suggested that a DNA-based identification system, founded on the mitochondrial gene, cytochrome c oxidase subunit 1 (COI), can aid the resolution of this diversity. While past work has validated the ability of COI sequences to diagnose species in certain taxonomic groups, the present study extends these analyses across the animal kingdom. The results indicate that sequence divergences at COI regularly enable the discrimination of closely allied species in all animal phyla except the Cnidaria. This success in species diagnosis reflects both the high rates of sequence change at COI in most animal groups and constraints on intraspecific mitochondrial DNA divergence arising, at least in part, through selective sweeps mediated via interactions with the nuclear genome.     

Calibrating the molecular clock beyond cytochrome b: assessing the evolutionary rate of COI in birds

Estimating the age of species or their component lineages based on sequence data is crucial for many studies in avian evolutionary biology. Although calibrations of the molecular clock in birds have been performed almost exclusively using cytochrome b (cyt b), they are commonly extrapolated to other mitochondrial genes. The existence of a large, standardized cytochrome c oxidase subunit I (COI) library generated as a result of the DNA barcoding initiative provides the opportunity to obtain a calibration for this mitochondrial gene in birds. In this study we compare the evolutionary rate of COI relative to cyt b across ten different avian orders. We obtained divergence estimates for both genes from nearly 300 phylogenetically independent pairs of species through the analysis of almost 5000 public sequences. For each pair of species we calculated the difference in divergence between COI and cyt b. Our results indicate that COI evolves on average 14% slower than cyt b, but also reveal considerable variation both among and within avian orders, precluding the use of this value as a standard adjustment for the COI molecular clock for birds. Our findings suggest that this variation is partially explained by a clear negative relationship between the difference in divergence in these genes and the age of species. Distances for cyt b are higher than those for COI for closely related species, but the values become similar as the divergence between the species increases. This appears to be the result of a stronger pattern of negative time‐dependency in the rate of cyt b than in that of COI, a difference that could be related to lower functional constraints on a small number of sites in cyt b that allow it to initially accumulate mutations more rapidly than COI.     

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.     

Intraspecific genetic variation in Paramecium revealed by mitochondrial cytochrome C oxidase I sequences

Studies of intraspecific genetic diversity of ciliates, such as population genetics and biogeography, are particularly hampered by the lack of suitable DNA markers. For example, sequences of the non-coding ribosomal internal transcribed spacer (ITS) regions are often too conserved for intraspecific analyses. We have therefore identified primers for the mitochondrial cytochrome c oxidase I (COI) gene and applied them for intraspecific investigations in Paramecium caudatum and Paramecium multimicronucleatum. Furthermore, we obtained sequences of the ITS regions from the same strains and carried out comparative sequence analyses of both data sets. The mitochondrial sequences revealed substantially higher variation in both Paramecium species, with intraspecific divergences up to 7% in P. caudatum and 9.5% in P. multimicronucleatum. Moreover, an initial survey of the population structure discovered different mitochondrial haplotypes of P. caudatum in one pond, thereby demonstrating the potential of this genetic marker for population genetic analyses. Our primers successfully amplified the COI gene of other Paramecium. This is the first report of intraspecific variation in free-living protozoans based on mitochondrial sequence data. Our results show that the high variation in mitochondrial DNA makes it a suitable marker for intraspecific and population genetic studies.     

Molecular population genetics of herbivore-induced protease inhibitor genes in European aspen (Populus tremula L., Salicaceae)

Plants defend themselves against the attack of natural enemies by using an array of both constitutively expressed and induced defenses. Long-lived woody perennials are overrepresented among plant species that show strong induced defense responses, whereas annual plants and crop species are underrepresented. However, most studies of plant defense genes have been performed on annual or short-lived perennial weeds or crop species. Here I use molecular population genetic methods to survey six wound-inducible protease inhibitors (PIs) in a long-lived woody, perennial plant species, the European aspen (Populus tremula), to evaluate the likelihood of either recurrent selective sweeps or balancing selection maintaining amino acid polymorphisms in these genes. The results show that none of the six PI genes have reduced diversities at synonymous sites, as would be expected in the presence of recurrent selective sweeps. However, several genes show some evidence of nonneutral evolution such as enhanced linkage disequilibrium and a large number of high-frequency-derived mutations. A group of at least four Kunitz trypsin inhibitor genes appear to have experienced elevated levels of nonsynonymous substitutions, indicating allelic turnover on an evolutionary timescale. One gene, TI1, has enhanced levels of intraspecific polymorphism at nonsynonymous sites and also has an unusual haplotype structure characterized by two divergent haplotypes occurring at roughly equal frequencies in the sample. One haplotype has very low levels of intraallelic nucleotide diversity, whereas the other haplotype has levels of diversity comparable to other genes in P. tremula. Patterns of sequence diversity at TI1 do not fit a simple model of either balancing selection or recurrent selective sweeps. This suggests that selection at TI1 is more complex, possibly involving allelic cycling.     

DNA barcodes provide new evidence of a recent radiation in the genus Sporophila (Aves: Passeriformes)

The capuchinos are a group of birds in the genus Sporophila that has apparently radiated recently, as evidenced by their lack of mitochondrial genetic diversity. We obtained cytochrome c oxidase I (COI) sequences (or DNA barcodes) for the 11 species of the group and various outgroups. We compared the patterns of COI variability of the capuchinos with those of the largest barcode data set from neotropical birds currently available (500 species representing 51% of avian richness in Argentina), and subjected COI sequences to neighbour-joining, maximum parsimony and Bayesian phylogenetic analyses as well as statistical parsimony network analysis. A clade within the capuchinos, the southern capuchinos, showed higher intraspecific and lower interspecific divergence than the remaining Argentine species. As most of the southern capuchinos shared COI haplotypes and pairwise distances within species were in many cases higher than distances between them, the phylogenetic affinities within the group remained unresolved. The observed genetic pattern is consistent with both incomplete lineage sorting and gene flow between species. The southern capuchinos constitute the only large group of species among the neotropical birds barcoded so far that are inseparable when using DNA barcodes, and one of few multispecies avian groups known to lack reciprocal monophyly. Extending the analysis to rapidly evolving nuclear and mitochondrial markers will be crucial to understanding this radiation. Apart from giving insights into the evolution of the capuchinos, this study shows how DNA barcoding can rapidly flag species or groups of species worthy of deeper study.     

Probing Evolutionary Patterns in Neotropical Birds through DNA Barcodes

Background

The Neotropical avifauna is more diverse than that of any other biogeographic region, but our understanding of patterns of regional divergence is limited. Critical examination of this issue is currently constrained by the limited genetic information available. This study begins to address this gap by assembling a library of mitochondrial COI sequences, or DNA barcodes, for Argentinian birds and comparing their patterns of genetic diversity to those of North American birds.

Methodology and Principal Findings

Five hundred Argentinian species were examined, making this the first major examination of DNA barcodes for South American birds. Our results indicate that most southern Neotropical bird species show deep sequence divergence from their nearest-neighbour, corroborating that the high diversity of this fauna is not based on an elevated incidence of young species radiations. Although species ages appear similar in temperate North and South American avifaunas, patterns of regional divergence are more complex in the Neotropics, suggesting that the high diversity of the Neotropical avifauna has been fueled by greater opportunities for regional divergence. Deep genetic splits were observed in at least 21 species, though distribution patterns of these lineages were variable. The lack of shared polymorphisms in species, even in species with less than 0.5M years of reproductive isolation, further suggests that selective sweeps could regularly excise ancestral mitochondrial polymorphisms.

Conclusions

These findings confirm the efficacy of species delimitation in birds via DNA barcodes, even when tested on a global scale. Further, they demonstrate how large libraries of a standardized gene region provide insight into evolutionary processes.     

COI barcoding of true bugs (Insecta, Heteroptera)

Several recent studies have proposed that partial DNA sequences of the cytochrome c oxidase I (COI) mitochondrial gene might serve as DNA barcodes for identifying and differentiating between animal species, such as birds, fish and insects. In this study, we tested the effectiveness of a COI barcode to identify true bugs from 139 species collected from Korea and adjacent regions (Japan, Northeastern China and Fareast Russia). All the species had a unique COI barcode sequence except for the genus Apolygus (Miridae), and the average interspecific genetic distance between closely related species was about 16 times higher than the average intraspecific genetic distance. DNA barcoding identified one probable new species of true bug and revealed identical or very recently divergent species that were clearly distinguished by morphological characteristics. Therefore, our results suggest that COI barcodes can reveal new cryptic true bug species and are able to contribute for the exact identification of the true bugs.     

Assessing the diversity of New Zealand freshwater harpacticoid copepods (Crustacea: Copepoda) using mitochondrial DNA (COI) barcodes

Taxonomic and ecological studies of freshwater harpacticoid copepods are limited globally by the ability to easily and accurately identify specimens. Here, we test the use of the mitochondrial cytochrome c oxidase subunit I (COI) gene locus as a tool for assessing the diversity of freshwater Harpacticoida. We obtained sequences from New Zealand harpacticoid copepods, representing two families, five genera and nine species, including the non-indigenous Elaphoidella sewelli. All species were delineated by the COI gene. However, high intraspecific diversity was evident among populations of Elaphoidella bidens (>12%), and between North and South Island populations of Bryocamptus pygmaeus (>18%), potentially indicating the presence of morphologically cryptic taxa. We suggest that mitochondrial DNA (COI) sequences can provide a useful tool for the routine identification of freshwater harpacticoid copepods. Applications of these data will include assessing species diversity and biogeography as well as assisting with the detection of non-indigenous species.     

Comprehensive DNA barcode coverage of North American birds

DNA barcoding seeks to assemble a standardized reference library for DNA-based identification of eukaryotic species. The utility and limitations of this approach need to be tested on well-characterized taxonomic assemblages. Here we provide a comprehensive DNA barcode analysis for North American birds including 643 species representing 93% of the breeding and pelagic avifauna of the USA and Canada. Most (94%) species possess distinct barcode clusters, with average neighbour-joining bootstrap support of 98%. In the remaining 6%, barcode clusters correspond to small sets of closely related species, most of which hybridize regularly. Fifteen (2%) currently recognized species are comprised of two distinct barcode clusters, many of which may represent cryptic species. Intraspecific variation is weakly related to census population size and species age. This study confirms that DNA barcoding can be effectively applied across the geographical and taxonomic expanse of North American birds. The consistent finding of constrained intraspecific mitochondrial variation in this large assemblage of species supports the emerging view that selective sweeps limit mitochondrial diversity.     

Improved COI barcoding primers for Southeast Asian perching birds (Aves: Passeriformes)

The All Birds Barcoding Initiative aims to assemble a DNA barcode database for all bird species, but the 648-bp 'barcoding' region of cytochrome c oxidase subunit I (COI) can be difficult to amplify in Southeast Asian perching birds (Aves: Passeriformes). Using COI sequences from complete mitochondrial genomes, we designed a primer pair that more reliably amplifies and sequences the COI barcoding region of Southeast Asian passerine birds. The 655-bp region amplified with these primers overlaps the COI region amplified with other barcoding primer pairs, enabling direct comparison of sequences with previously published DNA barcodes.     

Identification of Birds through DNA Barcodes

Short DNA sequences from a standardized region of the genome provide a DNA barcode for identifying species. Compiling a public library of DNA barcodes linked to named specimens could provide a new master key for identifying species, one whose power will rise with increased taxon coverage and with faster, cheaper sequencing. Recent work suggests that sequence diversity in a 648-bp region of the mitochondrial gene, cytochrome c oxidase I (COI), might serve as a DNA barcode for the identification of animal species. This study tested the effectiveness of a COI barcode in discriminating bird species, one of the largest and best-studied vertebrate groups. We determined COI barcodes for 260 species of North American birds and found that distinguishing species was generally straightforward. All species had a different COI barcode(s), and the differences between closely related species were, on average, 18 times higher than the differences within species. Our results identified four probable new species of North American birds, suggesting that a global survey will lead to the recognition of many additional bird species. The finding of large COI sequence differences between, as compared to small differences within, species confirms the effectiveness of COI barcodes for the identification of bird species. This result plus those from other groups of animals imply that a standard screening threshold of sequence difference (10× average intraspecific difference) could speed the discovery of new animal species. The growing evidence for the effectiveness of DNA barcodes as a basis for species identification supports an international exercise that has recently begun to assemble a comprehensive library of COI sequences linked to named specimens.     

Apparent Selection Intensity for the Cytochrome Oxidase Subunit I Gene Varies with Mode of Reproduction in Echinoderms

When most amino acid substitutions in protein-coding genes are slightly deleterious rather than selectively neutral, life history differences can potentially modify the effective population size or the selective regime, resulting in altered ratios of non-synonymous to synonymous substitutions among taxa. We studied substitution patterns for the mitochondrial cytochrome oxidase subunit I (COI) gene in a sea star genus (Leptasterias spp.) with an obligate brood-protecting mode of reproduction and small-scale population genetic subdivision, and compared the results to available COI sequences in nine other genera of echinoderms with pelagic larvae: three sea stars, five sea urchins and one brittle star. We predicted that this life history difference would be associated with differences in the ratio of non-synonymous (dN) to synonymous (dS) substitution rates. Leptasterias had a significantly greater dN/dS ratio (both between species and within species), a significantly smaller transition/transversion rate ratio, and a significantly lower average nucleotide diversity within species, than did the non-brooding genera. Other explanations for the results, such as altered mutation rates or selective sweeps, were not supported by the data analysis. These findings highlight the potential influence of reproductive traits and other life history factors on patterns of nucleotide substitution within and between species.     

Multiple cophylogenetic analyses reveal frequent cospeciation between pelecaniform birds and Pectinopygus lice

Lice in the genus Pectinopygus parasitize a single order of birds (Pelecaniformes). To examine the degree of congruence between the phylogenies of 17 Pectinopygus species and their pelecaniform hosts, sequences from mitochondrial 12S rRNA, 16S rRNA, COI, and nuclear wingless and EF1-alpha genes (2290 nucleotides) and from mitochondrial 12S rRNA, COI, and ATPases 8 and 6 genes (1755 nucleotides) were obtained for the lice and the birds, respectively. Louse data partitions were analyzed for evidence of incongruence and evidence of long-branch attraction prior to cophylogenetic analyses. Host-parasite coevolution was studied by different methods: TreeFitter, TreeMap, ParaFit, likelihood-ratio test, data-based parsimony method, and correlation of coalescence times. All methods agree that there has been extensive cospeciation in this host-parasite system, but the results are sensitive to the selection of different phylogenetic hypotheses and analytical methods for evaluating cospeciation. Perfect congruence between phylogenies is not found in this association, probably as a result of occasional host switching by the lice. Errors due to phylogenetic reconstruction methods, incorrect or incomplete taxon sampling, or to different loci undergoing different evolutionary histories cannot be rejected, thus emphasizing the need for improved cophylogenetic methodologies.     

Problematic barcoding in flatworms: A case-study on monogeneans and rhabdocoels (Platyhelminthes)

Some taxonomic groups are less amenable to mitochondrial DNA barcoding than others. Due to the paucity of molecular information of understudied groups and the huge molecular diversity within flatworms, primer design has been hampered. Indeed, all attempts to develop universal flatworm-specific COI markers have failed so far. We demonstrate how high molecular variability and contamination problems limit the possibilities for barcoding using standard COI-based protocols in flatworms. As a consequence, molecular identification methods often rely on other widely applicable markers. In the case of Monogenea, a very diverse group of platyhelminth parasites, and Rhabdocoela, representing one-fourth of all free-living flatworm taxa, this has led to a relatively high availability of nuclear ITS and 18S/28S rDNA sequences on GenBank. In a comparison of the effectiveness in species assignment we conclude that mitochondrial and nuclear ribosomal markers perform equally well. In case intraspecific information is needed, rDNA sequences can guide the selection of the appropriate (i.e. taxon-specific) COI primers if available.     

A cryptic, intergeneric cytochrome c oxidase I pseudogene in tyrant flycatchers (family: Tyrannidae)

Nuclear mitochondrial pseudogenes, or "numts", are nonfunctional copies of mitochondrial genes that have been translocated to the nuclear genome. Numts have been used to study differences in mutation rates between the nuclear and mitochondrial genomes, but have also been implicated as troublesome for phylogenetic studies and DNA-based species identification (i.e., DNA barcoding). In this study, a suspected numt discovered during a study of mitochondrial cytochrome c oxidase I (COI) diversity in North American birds was targeted and sequenced from tyrant flycatchers (family: Tyrannidae). In total, the numt was found in five taxa representing two genera. Substitution rates were compared between COI and numt sequences. None of the numt sequences harboured stop codons nor frameshift mutations, but phylogenetic analysis revealed they had accumulated more amino acid substitutions than the mitochondrial COI sequences. Mitochondrial COI appeared to be preferentially amplified in most cases, but methods for numt detection are discussed for cases like this where sequences lack obvious features for identification. Because of its persistence across a broad taxonomic lineage, this numt could form a valuable model system for studying evolution in numts. The full size of the numt and its location within the nuclear genome are yet to be determined.     

Genetic diversity and phylogenetic relationship among the western and the Asian honey bees based on two mitochondrial gene segments (COI and ND5)

The Asian honey bee species i.e., Apis cerana (the eastern honey bee), A. dorsata (the giant honey bee), and the western or European honey bee (A. mellifera) collected from Pakistan were studied using partial sequences from two mitochondrial genes (i) the Cytochrome c oxidase I (COI) and (ii) the mitochondrially encoded NADH dehydrogenase 5 (ND5) and then compared with other honey bees sequences (already submitted from different countries around the globe) obtained after the national center for biotechnology information (NCBI). DNA sequences were analyzed employing molecular evolutionary genetics analysis and Kimura 2-parameter model, neighbor-joining method was applied to investigate phylogenetic relationships, and DNA sequence polymorphism was applied to measure the genetic diversity within the genus Apis. The phylogenetic analyses yielded consistent results. Based on COI gene fragment in two Asian and European honey bee species from Pakistan and from other countries showed considerable genetic diversity levels and deviation among the species. While in contrast the phylogenetic analyses based on ND5 gene fragment in Asian and European honey bee species from Pakistan and other countries showed comparatively higher genetic diversity indices and variations than the COI gene. So, in the genus Apis, the mitochondrial ND5 region has shown the possibility to answer the interactions among species. A further detailed work (by linking the analysis of other genomic and mitochondrial genes) is required for good quality solution to establish the concise genetic diversity and interaction among the Apis species. The objective of this study was to explore the extent of genetic differences and phylogenetic links among the three kinds of honey bee species from Pakistan and comparing them with other bee species around the globe.     

Mitochondrial DNA sequence variation in the spruce budworm species complex (Choristoneura: Lepidoptera)

A combination of polymerase-chain-reaction amplification and automated DNA sequencing was used to survey variation in a species complex of pest insects, the spruce budworms (Choristoneura fumiferana species group), and an outgroup species, C. rosaceana. We sequenced an mtDNA region of 1,573 bp that extends from the middle of cytochrome oxidase subunit I (COI) through tRNA leucine (UUR) to the end of cytochrome oxidase subunit II. In addition, we examined levels of intraspecific variation within a 470-bp region of the COI gene. Choristoneura fumiferana clearly represented the oldest lineage within its species group, with 2.7%-2.9% sequence divergence from the other species. In contrast, the four remaining species (C. pinus, C. biennis, C. occidentalis, and C. orae) had closely related or identical mtDNA, with < 1% divergence among most of their haplotypes. Despite its older lineage and widespread geographic distribution, C. fumiferana showed significantly lower intraspecific genetic diversity than did C. occidentalis. Choristoneura orae shared haplotypes with C. occidentalis and C. biennis, and species-level separation of these three species was not supported. Two divergent, uncommon haplotypes were also found in C. occidentalis and C. biennis. The divergent haplotype in C. biennis had an unusually high number of inferred amino acid replacements, suggesting selective differences between mitochondrial DNA haplotypes. Transition:transversion ratios in Choristoneura paralleled those found in Drosophila; transition:transversion ratios were highest in closely related sequences but decreased with increasing sequence divergence. Nucleotide composition showed an A+T bias that was near the high end of the range known for insects. This work illustrates the potential utility of direct DNA sequencing in assessing population structures, species limits, and phylogenetic relationships among organisms that have not previously been subjected to DNA analysis.      

Mitochondrial DNA variation in the caramote prawn <Emphasis Type="Italic">Penaeus (Melicertus) kerathurus</Emphasis> across a transition zone in the Mediterranean Sea

In this study we analysed mitochondrial DNA variation in Penaeus kerathurus prawns collected from seven locations along a transect across the Siculo–Tunisian region in order to verify if any population structuring exists over a limited geographical scale and to delineate the putative transition zone with sufficient accuracy. Partial DNA sequences of COI and 16S genes were analysed. In contrast to the highly conservative 16S gene, the COI sequences exhibited sufficient diversity for population analysis. The COI gene revealed low levels of haplotype and nucleotide diversities. The size of the annual landings of this commercial species suggests large population sizes. Hence, the low genetic diversity detected in this study could indicate a possible reduction in effective population sizes in the past. We detected significant genetic differentiation between eastern and western populations likely due to restricted gene flow across the Siculo–Tunisian boundary. We discuss the different evolutionary forces that may have shaped the genetic variation and suggest that the genetic divide is probably maintained by present-day dispersal limitation. R. Zitari-Chatti and N. Chatti are contributed equally to the work.     

Applications of DNA barcoding to fish landings: authentication and diversity assessment

DNA barcoding methodologies are being increasingly applied not only for scientific purposes but also for diverse real-life uses. Fisheries assessment is a potential niche for DNA barcoding, which serves for species authentication and may also be used for estimating within-population genetic diversity of exploited fish. Analysis of single-sequence barcodes has been proposed as a shortcut for measuring diversity in addition to the original purpose of species identification. Here we explore the relative utility of different mitochondrial sequences (12S rDNA, COI, cyt b, and D-Loop) for application as barcodes in fisheries sciences, using as case studies two marine and two freshwater catches of contrasting diversity levels. Ambiguous catch identification from COI and cyt b was observed. In some cases this could be attributed to duplicated names in databases, but in others it could be due to mitochondrial introgression between closely related species that may obscure species assignation from mtDNA. This last problem could be solved using a combination of mitochondrial and nuclear genes. We suggest to simultaneously analyze one conserved and one more polymorphic gene to identify species and assess diversity in fish catches.