Assessment of loci for DNA barcoding in the genus Thrips (Thysanoptera:Thripidae)

Identification of the juveniles of economically important thrips species on imports by morphology alone can be challenging and culturing is usually required. In the case of EU quarantine species such as Thrips palmi, rapid and accurate identification is essential. DNA barcoding using the Cytochrome oxidase I (COI) gene has become a popular technique for species identification; however, in some invertebrate genera COI has been shown to provide insufficient variability for species discrimination. This study presents a comparison of five different loci to investigate their ability to discriminate a small number of Thrips species. All five loci discriminated the species by neighbour-joining tree and varying degrees of discrimination were determined upon further investigation of the intraspecific and interspecific distances. Two distinct COI clades were observed for T. Palmi and judged to be COI haplotypes when data from the other four additional loci and geographical collection data were taken into consideration. COI was shown to provide sufficient variation to be used in future DNA barcoding efforts within the genus Thrips.     

DNA barcoding and genetic diversity of phyllostomid bats from the Yucatan Peninsula with comparisons to Central America

The mitochondrial cytochrome c oxidase subunit I gene is the standard DNA barcoding region used for species identification and discovery. We examined the variation of COI (454 bp) to discriminate 20 species of bats in the family Phyllostomidae that are found in the Yucatan Peninsula of southeastern Mexico and northern Guatemala and compared them genetically to other samples from Central America. The majority of these species had low intraspecific variation (mean = 0.75%), but some taxa had intraspecific variation ranging to 8.8%, suggesting the possibility of cryptic species (i.e. Desmodus rotundus and Artibeus jamaicensis). There was a recurring biogeographic pattern in eight species with a separation of northern and southern Middle American localities. The Yucatan Peninsula was a discrete area identified in four species, whereas Panama was recovered in five species of phyllostomid bats. Our study establishes a foundation for further molecular work incorporating broader taxonomic and geographic coverage to better understand the phylogeography and genetic diversity that have resulted from the ecological constraints in this region and the remarkable differentiation of bats in the Neotropics.     

Patterns of DNA Barcode Variation in Canadian Marine Molluscs

Background

Molluscs are the most diverse marine phylum and this high diversity has resulted in considerable taxonomic problems. Because the number of species in Canadian oceans remains uncertain, there is a need to incorporate molecular methods into species identifications. A 648 base pair segment of the cytochrome c oxidase subunit I gene has proven useful for the identification and discovery of species in many animal lineages. While the utility of DNA barcoding in molluscs has been demonstrated in other studies, this is the first effort to construct a DNA barcode registry for marine molluscs across such a large geographic area.

Methodology/Principal Findings

This study examines patterns of DNA barcode variation in 227 species of Canadian marine molluscs. Intraspecific sequence divergences ranged from 0–26.4% and a barcode gap existed for most taxa. Eleven cases of relatively deep (>2%) intraspecific divergence were detected, suggesting the possible presence of overlooked species. Structural variation was detected in COI with indels found in 37 species, mostly bivalves. Some indels were present in divergent lineages, primarily in the region of the first external loop, suggesting certain areas are hotspots for change. Lastly, mean GC content varied substantially among orders (24.5%–46.5%), and showed a significant positive correlation with nearest neighbour distances.

Conclusions/Significance

DNA barcoding is an effective tool for the identification of Canadian marine molluscs and for revealing possible cases of overlooked species. Some species with deep intraspecific divergence showed a biogeographic partition between lineages on the Atlantic, Arctic and Pacific coasts, suggesting the role of Pleistocene glaciations in the subdivision of their populations. Indels were prevalent in the barcode region of the COI gene in bivalves and gastropods. This study highlights the efficacy of DNA barcoding for providing insights into sequence variation across a broad taxonomic group on a large geographic scale.     

Cryptic diversity in flathead fishes (Scorpaeniformes: Platycephalidae) across the Indo‐West Pacific uncovered by DNA barcoding

Identification of taxonomical units underpins most biological endeavours ranging from accurate biodiversity estimates to the effective management of sustainably harvested, protected or endangered species. Successful species identification is now frequently based on a combination of approaches including morphometrics and DNA markers. Sequencing of the mitochondrial COI gene is an established methodology with an international campaign directed at barcoding all fishes. We employed COI sequencing alongside traditional taxonomic identification methods and uncovered instances of deep intraspecific genetic divergences among flathead species. Sixty‐five operational taxonomic units (OTUs) were observed across the Indo‐West Pacific from just 48 currently recognized species. The most comprehensively sampled taxon, Platycephalus indicus, exhibited the highest levels of genetic diversity with eight lineages separated by up to 16.37% genetic distance. Our results clearly indicate a thorough reappraisal of the current taxonomy of P. indicus (and its three junior synonyms) is warranted in conjunction with detailed taxonomic work on the other additional Platycephalidae OTUs detected by DNA barcoding.     

Deep sympatric mtDNA divergence in the autumnal moth (Epirrita autumnata)

Deep sympatric intraspecific divergence in mtDNA may reflect cryptic species or formerly distinct lineages in the process of remerging. Preliminary results from DNA barcoding of Scandinavian butterflies and moths showed high intraspecific sequence variation in the autumnal moth, Epirrita autumnata. In this study, specimens from different localities in Norway and some samples from Finland and Scotland, with two congeneric species as outgroups, were sequenced with mitochondrial and nuclear markers to resolve the discrepancy found between mtDNA divergence and present species‐level taxonomy. We found five COI sub‐clades within the E. autumnata complex, most of which were sympatric and with little geographic structure. Nuclear markers (ITS2 and Wingless) showed little variation and gave no indications that E. autumnata comprises more than one species. The samples were screened with primers for Wolbachia outer surface gene (wsp) and 12% of the samples tested positive. Two Wolbachia strains were associated with different mtDNA sub‐clades within E. autumnata, which may indicate indirect selection/selective sweeps on haplotypes. Our results demonstrate that deep mtDNA divergences are not synonymous with cryptic speciation and this has important implications for the use of mtDNA in species delimitation, like in DNA barcoding.     

Using DNA barcodes for assessing diversity in the?family Hybotidae (Diptera,Empidoidea)

Empidoidea is one of the largest extant lineages of flies, but phylogenetic relationships among species of this group are poorly investigated and global diversity remains scarcely assessed. In this context, one of the most enigmatic empidoid families is Hybotidae. Within the framework of a pilot study, we barcoded 339 specimens of Old World hybotids belonging to 164 species and 22 genera (plus two Empis as outgroups) and attempted to evaluate whether patterns of intra- and interspecific divergences match the current taxonomy. We used a large sampling of diverse Hybotidae. The material came from the Palaearctic (Belgium, France, Portugal and Russian Caucasus), the Afrotropic (Democratic Republic of the Congo) and the Oriental realms (Singapore and Thailand). Thereby, we optimized lab protocols for barcoding hybotids. Although DNA barcodes generally well distinguished recognized taxa, the study also revealed a number of unexpected phenomena: e.g., undescribed taxa found within morphologically very similar or identical specimens, especially when geographic distance was large; some morphologically distinct species showed no genetic divergence; or different pattern of intraspecific divergence between populations or closely related species. Using COI sequences and simple Neighbour-Joining tree reconstructions, the monophyly of many species- and genus-level taxa was well supported, but more inclusive taxonomical levels did not receive significant bootstrap support. We conclude that in hybotids DNA barcoding might be well used to identify species, when two main constraints are considered. First, incomplete barcoding libraries hinder efficient (correct) identification. Therefore, extra efforts are needed to increase the representation of hybotids in these databases. Second, the spatial scale of sampling has to be taken into account, and especially for widespread species or species complexes with unclear taxonomy, an integrative approach has to be used to clarify species boundaries and identities.     

Cryptic diversity in Indo-Australian rainbowfishes revealed by DNA barcoding: implications for conservation in a biodiversity hotspot candidate

The rainbowfishes of the family Melanotaeniidae represent one of the largest radiations of freshwater fishes from the Indo-Australian archipelago. A total of 75 nominal species have been described, among which several have become very popular among tropical fish hobbyists because of their tendency to form large schools of colourful individuals. Facing habitat loss and competition or predation by introduced species, this group has become a priority in the conservation of ornamental fishes in Indonesia. In this context, several expeditions have been conducted between 2007 and 2010 in Indonesian Papua with the aim to initiate a large-scale survey of the genetic resources in this group. We assessed the diversity of the Papua rainbowfishes with DNA barcoding. We sequenced the mitochondrial COI gene for 350 specimens belonging to 53 nominal species throughout the Indo-Australian archipelago. Unexpected levels of cryptic diversity and endemism were detected since additional cryptic lineages were detected in several watersheds from the Vogelkop and the Lengguru massif. DNA barcoding supports the presence of nearly 30 evolutionary lineages among the 15 nominal species sampled in the Vogelkop and all these lineages are endemic to a single lake or watershed. This result highlights that the diversity of the family has been largely underestimated and urges for the identification of conservation priorities in Papua.     

DNA barcoding of African fruit bats (Mammalia, Pteropodidae). The mitochondrial genome does not provide a reliable discrimination between Epomophorus gambianus and Micropteropus pusillus

Sequences of the mitochondrial cytochrome c oxidase subunit I (COI) gene have been shown to be useful for species identification in various groups of animals. However, the DNA barcoding approach has never been tested on African fruit bats of the family Pteropodidae (Mammalia, Chiroptera). In this study, the COI gene was sequenced from 120 bats collected in the Central African Republic and belonging to either Epomophorus?gambianus or Micropteropus?pusillus, two species easily diagnosed on the basis of morphological characters, such as body size, skull shape and palatal ridges. Two additional molecular markers were used for comparisons: the complete mitochondrial cytochrome b gene and the intron 7 of the nuclear β-fibrinogen (FGB) gene. Our results reveal an unexpected discordance between mitochondrial and nuclear genes. The nuclear FGB signal agrees with our morphological identifications, as the three alleles detected for E.?gambianus are divergent from the fourteen alleles found for M.?pusillus. By contrast, this taxonomic distinction is not recovered with the analyses of mitochondrial genes, which support rather a polyphyletic pattern for both species. The conflict between molecular markers is explained by multiple mtDNA introgression events from M.?pusillus into E.?gambianus or, alternatively, by incomplete lineage sorting of mtDNA haplotypes associated with positive selection on FGB alleles of M.?pusillus. Our work shows the failure of DNA barcoding to discriminate between two morphologically distinct fruit bat species and highlights the importance of using both mitochondrial and nuclear markers for taxonomic identification.     

DNA barcoding of Neotropical bats: species identification and discovery within Guyana

Sequence diversity in the cytochrome c oxidase subunit 1 gene has been shown to be an effective tool for species identification and discovery in various groups of animals, but has not been extensively tested in mammals. We address this gap by examining the performance of DNA barcodes in the discrimination of 87 species of bats from Guyana. Eighty‐one of these species showed both low intraspecific variation (mean = 0.60%), and clear sequence divergence from their congeners (mean = 7.80%), while the other six showed deeply divergent intraspecific lineages suggesting that they represent species complexes. Although further work is needed to examine patterns of sequence diversity at a broader geographical scale, the present study validates the effectiveness of barcoding for the identification of regional bat assemblages, even highly diverse tropical faunas.     

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.     

DNA Barcoding for the Identification of Sand Fly Species (Diptera,Psychodidae, Phlebotominae) in Colombia

Sand flies include a group of insects that are of medical importance and that vary in geographic distribution, ecology, and pathogen transmission. Approximately 163 species of sand flies have been reported in Colombia. Surveillance of the presence of sand fly species and the actualization of species distribution are important for predicting risks for and monitoring the expansion of diseases which sand flies can transmit. Currently, the identification of phlebotomine sand flies is based on morphological characters. However, morphological identification requires considerable skills and taxonomic expertise. In addition, significant morphological similarity between some species, especially among females, may cause difficulties during the identification process. DNA-based approaches have become increasingly useful and promising tools for estimating sand fly diversity and for ensuring the rapid and accurate identification of species. A partial sequence of the mitochondrial cytochrome oxidase gene subunit I (COI) is currently being used to differentiate species in different animal taxa, including insects, and it is referred as a barcoding sequence. The present study explored the utility of the DNA barcode approach for the identification of phlebotomine sand flies in Colombia. We sequenced 700 bp of the COI gene from 36 species collected from different geographic localities. The COI barcode sequence divergence within a single species was <2% in most cases, whereas this divergence ranged from 9% to 26.6% among different species. These results indicated that the barcoding gene correctly discriminated among the previously morphologically identified species with an efficacy of nearly 100%. Analyses of the generated sequences indicated that the observed species groupings were consistent with the morphological identifications. In conclusion, the barcoding gene was useful for species discrimination in sand flies from Colombia.     

Identification of Channa species using the partial cytochrome c oxidase subunit I (COI) gene as a DNA barcoding marker

The snakehead fish of the genus Channa are an important food fish in China. However, the molecular identification and phylogeny of this genus is poorly understood. Here, we present the utility of partial sequences of the COI gene for use in DNA barcoding for the identification of Channa individuals, which includes four species: Channa argus, Channa maculata, Channa asiatica, and Channa striata. A total of 19 haplotypes were identified in this study. The interspecific K2P distances were higher than intraspecific distances. The lowest interspecific distance (0.091) was between C. argus and C. maculata while the highest interspecific distance (0.219) was between C. argus and C. striata. No intraspecific–interspecific distance overlaps were observed, and a distinct barcoding gap was found between intraspecific and interspecific distances in each species. Our results showed that the partial COI gene is an effective DNA barcoding marker for identifying Channa species.     

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.     

Deciphering amphibian diversity through DNA barcoding: chances and challenges

Amphibians globally are in decline, yet there is still a tremendous amount of unrecognized diversity, calling for an acceleration of taxonomic exploration. This process will be greatly facilitated by a DNA barcoding system; however, the mitochondrial population structure of many amphibian species presents numerous challenges to such a standardized, single locus, approach. Here we analyse intra- and interspecific patterns of mitochondrial variation in two distantly related groups of amphibians, mantellid frogs and salamanders, to determine the promise of DNA barcoding with cytochrome oxidase subunit I (cox1) sequences in this taxon. High intraspecific cox1 divergences of 7-14% were observed (18% in one case) within the whole set of amphibian sequences analysed. These high values are not caused by particularly high substitution rates of this gene but by generally deep mitochondrial divergences within and among amphibian species. Despite these high divergences, cox1 sequences were able to correctly identify species including disparate geographic variants. The main problems with cox1 barcoding of amphibians are (i) the high variability of priming sites that hinder the application of universal primers to all species and (ii) the observed distinct overlap of intraspecific and interspecific divergence values, which implies difficulties in the definition of threshold values to identify candidate species. Common discordances between geographical signatures of mitochondrial and nuclear markers in amphibians indicate that a single-locus approach can be problematic when high accuracy of DNA barcoding is required. We suggest that a number of mitochondrial and nuclear genes may be used as DNA barcoding markers to complement cox1.     

DNA barcoding reveals taxonomic uncertainty in Salminus (Characiformes)

Salminus is a genus composed of four species of migratory fishes and top predators. Although this group has great economic and ecological importance, the species level diversity of Salminus is not yet completely clarified. Our goal was to detect if this taxonomic problem is the consequence of lineage divergence within species, and, if so, whether these divergences are sufficient to flag potentially undescribed taxa. We employed the standard DNA barcoding analyses and a generalized mixed Yule-coalescent model (GMYC) using one mitochondrial (COI) marker and Bayesian Inference (BI) reconstruction for one nuclear (RAG2) marker for all currently recognized species of Salminus, sampled across different hydrographic basins. Eight MOTUs (Molecular Operational Taxonomic Units) were determined by distance and model-based analyses, and recovered with BI analyses for COI. Only Salminus affinis and Salminus franciscanus formed monophyletic haplogroups. Salminus brasiliensis and Salminus hilarii had two and four distinct mitochondrial lineages, respectively, and higher intraspecific K2P distances than the adopted optimum threshold. The RAG2 gene tree supported two lineages of S. hilarii (S. hilarii Amazon and S. hilarii Araguaia), while the other mitochondrial lineages of S. hilarii and S. brasiliensis were not supported. All lineages of both species, corresponded to morphological variation described in previous studies. We suggest, based on the DNA barcoding analysis, a new taxonomic scenario and conservation polices for Salminus in the Brazilian territory.     

中国蛛缘蝽科物种DNA条形码研究(英文)

【目的】本研究旨在探讨DNA条形码对中国蛛缘蝽科(半翅目：缘蝽总科)物种界定的适用性。【方法】对中国蛛缘蝽科13属23种207个样本的线粒体COI基因DNA条形码序列进行扩增,并扩增稻缘蝽属Leptocorisa 3个物种的31条内转录间隔区1(ITS-1)序列作为辅助标记。使用MEGA 11软件计算种间和种内遗传距离(Kimura 2-parameter, K2P);采用邻接法(neighbor-joining, NJ)进行物种聚类分析;利用中介邻接网络算法构建单倍型网络图。【结果】基于线粒体COI DNA条形码序列得出测试的中国蛛缘蝽科所有23个种的种内平均K2P距离在2%以下,种间K2P距离在0.98%~23.98%之间(平均17.50%)。多数物种彼此能够被较好地分开,且支持率较高。其中,中稻缘蝽Leptocorisa chinensis和大稻缘蝽L. oratoria共享部分COI单倍型,造成COI条形码无法区分二者,可通过ITS-1序列在单倍型网络分析中将二者区分。【结论】本研究得出的中国蛛缘蝽科中绝大部分物种的DNA条形码数据分析结果与基于形态特征的分类单元一致。然而,对于其中亲缘关系极近的物种,单靠线粒体数据尤其是COI条形码序列无法进行准确界定,需引入其他DNA序列或其他类型数据进行区分。     

DNA barcodes for globally threatened marine turtles: a registry approach to documenting biodiversity

DNA barcoding is a global initiative that provides a standardized and efficient tool to catalogue and inventory biodiversity, with significant conservation applications. Despite progress across taxonomic realms, globally threatened marine turtles remain underrepresented in this effort. To obtain DNA barcodes of marine turtles, we sequenced a segment of the cytochrome c oxidase subunit I (COI) gene from all seven species in the Atlantic and Pacific Ocean basins (815 bp; n = 80). To further investigate intraspecific variation, we sequenced green turtles (Chelonia mydas) from nine additional Atlantic/Mediterranean nesting areas (n = 164) and from the Eastern Pacific (n = 5). We established character-based DNA barcodes for each species using unique combinations of character states at 76 nucleotide positions. We found that no haplotypes were shared among species and the mean of interspecific variation ranged from 1.68% to 13.0%, and the mean of intraspecific variability was relatively low (0–0.90%). The Eastern Pacific green turtle sequence was identical to an Australian haplotype, suggesting that this marker is not appropriate for identifying these phenotypically distinguishable populations. Analysis of COI revealed a north–south gradient in green turtles of Western Atlantic/Mediterranean nesting areas, supporting a hypothesis of recent dispersal from near equatorial glacial refugia. DNA barcoding of marine turtles is a powerful tool for species identification and wildlife forensics, which also provides complementary data for conservation genetic research.     

Resolving genetic diversity in Australasian Culex mosquitoes: Incongruence between the mitochondrial cytochrome c oxidase I and nuclear acetylcholine esterase 2

Insects that vector pathogens are under constant surveillance in Australasia although the repertoire of genetic markers to distinguish what are often cryptic mosquito species remains limited. We present a comparative assessment of the second exon–intron region of the acetylcholine esterase 2 gene (ace-2) and the mitochondrial DNA cytochrome c oxidase I (COI) using two closely related Australasia mosquitoes Culex annulirostris and Culex palpalis. The COI revealed eight divergent lineages of which four were confirmed with the ace-2. We dissect out the nuclear chromosomal haplotypes of the ace-2 as well as the exon–intron regions by assessing the protein’s tertiary structure to reveal a hypervariable 5’-exon that forms part of an external protein loop and displays a higher polymorphic rate than the intron. We retrace the evolutionary history of these mosquitoes by phylogenetic inference and by testing different evolutionary hypotheses. We conclude that DNA barcoding using COI may overestimate the diversity of Culex mosquitoes in Australasia and should be applied cautiously with support from the nuclear DNA such as the ace-2. Together the COI and ace-2 provide robust evidence for distinct cryptic Culex lineages—one of which correlates exactly with the southern limit of Japanese encephalitis virus activity in Australasia.     

Population genetic structure and intraspecific genetic distance of Periplaneta americana (Blattodea: Blattidae) based on mitochondrial and nuclear DNA markers

The American cockroach (Periplaneta americana) is a globally invasive pest that can cause significant economic loss and threaten human health. Although it is abundant and lives in close proximity to humans, few studies have investigated the genetic diversity of P. americana. Our study analyzed 1,053 P. americana and other Periplaneta species' samples from different locations in China and the United States. A traditional tree‐based method using 17 unique mitochondrial COI haplotypes of P. americana and 20 haplotypes of the other Periplaneta species accurately identified P. americana with a barcoding threshold of 5.1%. To identify the population genetic structure of P. americana, we investigated wingless gene and pooled them with obtained mtDNA data for a combined analysis. Although the genetic diversity of the USA group was relatively higher than the China group, the number of haplotypes and alleles of both groups was small. The analysis of molecular variance (AMOVA), intraspecific phylogeny, and haplotype networks indicated that P. americana had very little global genetic differentiation. The weak geographic genetic structure might reflect the human‐mediated dispersal of P. americana. Despite no apparent phylogeographic assignment of mtDNA and nuclear lineages was observed in both BI trees, the integrated COI sequence data identified four distinct P. americana haplotype groups, showing four ancient maternal lineages of P. americana in China and the United States.     

Opening a can of worms: unprecedented sympatric cryptic diversity within British lumbricid earthworms

Earthworms play a major role in many aspects of soil fertility, food web ecology and ecosystem functioning, and hence are frequently the subjects of, for example, ecological and toxicological research. Our aim was to examine the genetic structure of common earthworm species, to identify cryptic lineages or species that may be distinct ecotypes or biotypes (and hence confound current research based upon morphotypes) and to try to explain the massive cryptic diversity that eventually emerged. We demonstrated that species such as Allolobophora chlorotica, Aporrectodea longa, Aporrectodea rosea and Lumbricus rubellus all comprise highly divergent lineages with species-level divergence at the mitochondrial cytochrome oxidase I (COI) gene. In Allo. chlorotica alone, we found 55 haplotypes for COI, with 35 of these being found in pink and 20 in green morph worms. There were no cases of the two colour morphs sharing COI haplotypes. Phylogenetic analyses of mitochondrial COI and 16S genes showed the presence of five highly divergent lineages, suggesting the presence of multiple cryptic species within Allo. chlorotica. There was no clear geographical pattern to lineage distribution and many populations were polymorphic for both mitochondrial DNA lineage and colour morph. Amplified fragment length polymorphism results, based on two primer combinations, were broadly congruent with mitochondrial DNA results with one significant exception. Despite showing over 14% divergence at COI, amplified fragment length polymorphism markers showed that the two green morph lineages may be interbreeding and therefore represent a single taxon. The cryptic diversity revealed by these results has profound consequences for all areas of earthworm research.