Pieces of a puzzle

Why DNA barcoding works as well as it does is an unsolved scientific puzzle. It is long observed that mitochondrial DNA differences within animal species are generally much smaller than those among species and, in the landscape of phylogenetic trees, mitochondrial DNA sequences of most species form single clusters distinct from those of other species. As a result “mtDNA data and traditional taxonomic assignments tend to converge on what may be “real” biotic units in nature”  (Avise and Walker 1999 Proc Natl Acad Sci USA 96:992). Although Avise and Walker’s original observation was largely based on terrestrial, temperate zone vertebrates, growing barcode libraries demonstrate similar patterning in diverse invertebrates, vertebrates, and protists in marine and terrestrial environments, and in tropical and temperate zones, and in at least some fungi and plants (see last week’s post on COI barcodes in red algae

Smith et al 2004 Science 305:371What underlies the usual patterning of small differences within and large differences among most animal species? The unsolved puzzle is how to reconcile these two findings. Large differences among closely-related species indicates mitochondrial DNA undergoes rapid sequence evolution, and there are reasonable mechanistic explanations for why this might be so. On the other hand, rapid sequence evolution should also lead to accumulation of sequence diversity within species over time and in those with large populations. Instead the data shows a relative absence of variation within most species, including those thought to be ancient and those with enormous population sizes. I will set aside two of the usual suspects: population bottlenecks and small effective population size. Population bottlenecks are implausible given the diversity of species showing this pattern. Postulating a small effective population size is a restatement of the finding of absence of variation, not an explanation. 

This table-napkin analysis leads me to selective sweeps as pruning mitochondrial diversity within species (eg Bazin et al 2006 Science 312:570, see also editorial and reader commentary). If selective sweeps restrict mitochondrial diversity, then the question becomes what is being selected for? Environmental adaptation seems unlikely, as restricted variation is seen in species that are as best one can tell morphologically and ecologically unchanged (eg see earlier posts on horseshoe crabs, salamanders). It might be there is little tolerance for genetic variation due to interactions of mitochondrial proteins with other cellular components, but if so there should be species with genetic stasis in mitochondrial DNA, just as there are many species with apparent morphologic stasis. However, in simple distance trees most species show roughly similar genetic distances. 

I am intrigued by a time series of influenza A hemagglutinin gene evolution which reflects competition between virus and host and wonder if there might be some kind of competition that helps drive mitochondrial sequence evolution forward and at the same time suppresses variation. It is exciting there will be an EMBO workshop “Molecular Biodiversity and DNA Barcodes” May 2007 in Rome which may help answer scientific questions posed by DNA barcode data.

COI outperforms in red algae, oldest known multicellular eukaryotes

Red seaweeds, kingdom Rhodophyta, are “weird, wonderful, and extremely ancient” organisms distantly related to plants (Tudge 2000 The Variety of Life). Multicellular red algae arose at least 1.2 billion years ago, predating the earliest multicellular animals by 600 million years. Visual identification is challenging, as “morphology can be highly variable within and between species, and conspicuous features with which they can be readily identified are often lacking. In addition, highly convergent morphology is commonly encountered. …Identification is further compounded by the complexities of red algal life histories, many of which have a heteromorphic alternation of generations. Different life history stages of species have frequently been described as separate species and have only been linked through observations of life histories in culture and use of molecular techniques” (Robba et al 2006 Am J Botany 93:1101). Red seaweeds call out for DNA-based identification.

In a preliminary appraisal, COI barcodes “resolved accurately and unequivocally species identities” in three morphologically confusing species complexes of red algae (Saunders 2005 Phil Trans Royal Soc B 360:1879). In September 2006 Am J Botany, researchers from Natural History Museum, London, and University of Bristol, analyzed COI sequences of 79 specimens from 37 red algae, chosen to span 6 orders and represent species “extremely difficult to distinguish morphologically”. All species could be discriminated and intraspecific variation was generally much lower (0.0-0.7%) than interspecific differences (5.2-27.3%). Large intraspecific distances suggested the presence of cryptic species or incipient speciation in two cases. COI outperformed a commonly used genetic marker, the plastid Rubisco spacer, in revealing cryptic diversity. As is generally observed, “a good taxonomic foundation coupled with extensive sampling of taxa is essential for the development of an effective identification system.”  A holiday wish for red algae identifiers: a well-curated DNA barcode library.

Labelling specimens and species with standardized DNA sequences

Labelling specimens helps make natural history museum collections valuable. Once specimens have been carefully examined by an expert, a species label can be applied. Labelled specimens can then be re-found and re-examined. Unsorted specimens lacking species designations are more raw material than scientific resource.  

As in all sciences, taxonomic knowledge undergoes continuous revision. Changes in names and understanding of species boundaries may mean that specimen labels need to be updated in light of current knowledge. Partly to compensate, there is enormous effort to link current taxonomic understanding to historical species descriptions. Nonetheless, the accuracy of specimen labels may decay over time. In one study, a revision of predatory flies (Euscelidia) (Meier and Dikow Conservation Biol 2004 18:478), 83% of 1383 specimens from 19 collections were found to be incorrectly identified.

DNA sequences offer a simple approach to help keep specimen labels up to date. DNA sequences are an intrinsic, unvarying characteristic of a specimen. A DNA sequence from a standardized locus (ie a DNA barcode) can serve as a permanent “name” for a specimen. Results so far with more than 20,000 invertebrate and vertebrate species show it is generally straightforward to use a COI barcode to assign specimens to known species. Future taxonomic revisions may change species names or boundaries, but that will not change DNA barcodes of specimens or the clustering patterns of barcode sequences. Thus it should be simple to use a specimen’s barcode sequence “name” to search a regularly revised public database for the current species name it corresponds to. A public database of sequences, specimens, and associated data as is BOLD can undergo continuous refinement, whereas revising labels in every collection around the world is impractical. Given that assigning specimens to species involves sorting among millions of species names, this approach does not have to give 100% resolution to species level to be valuable. It will be useful both for those species with unique barcode clusters or characters and for those sets of closely-related species with overlapping or indistinguishable barcodes.

Finally, on a related note I observe that diagnostic keys are unwieldy and not easily amenable to computerization. For example, Dragonflies of North America (Gainesville: Scientific Publishers, 2000) by Needham, Westfall, and May covers 350 species in 939 pages, or about 2.7 pages per species, which by my count is typical for diagnostic keys. At this rate, a world key for the 5500 known Odonates would run to 15,000 pages. The recondite language that is required to describe morphologic detail in this and other keys makes the tools inacessible except to highly-trained persons.  On the other hand, anyone can interpret a DNA sequence. It is exciting that taxonomists and others are increasingly taking up the challenge of translating taxonomic knowledge into a much more widely-accessible format, namely DNA barcode libraries.

CoML Showcase Public report

The Census of Marine Life released its annual highlights, now called the “Showcase Public Report” (5 MB). Jesse Ausubel and Paul Waggoner prepared the report with Darlene Crist and Darrell McIntyre of the University of Rhode Island, drawing on the remarkable work of the thousands of contributors to the CoML. The CoML highlights 2006 Site also enables visitors to:

  • Download the Press Release
  • View images & download hi-res versions
  • View additional images & hi-res versions
  • Download PDF version of Report
  • Download PDF version of Report COVER
  • Download Showcase Public Report. This report is the same text and images as the Highlights Report laid out to print on 8.5″ x 11″ paper with additional hyperlinks to individual projects, more information, images, or news releases.
  • View videos
  • Download PDFs of the translated text for the 2006 CoML Highlights Report in German, Portuguese, Chinese, and French

Single locus mtDNA resolves freshwater snail species

Freshwater snails are intermediate hosts for schistosomiasis and flukes, trematode parasites that infect approximately 10% of world’s human population. Freshwater snails are also indicator species for water quality. Snail identification is essential for reducing disease burden and monitoring water quality.

Researchers at the University of Frankfurt (November 2006 BMC Evol Biol 6:100) compared efficacy of morphologic and DNA-based taxonomy in freshwater snails in the genus Radix. Regarding Radix species in northwestern Europe, “species determination by shell morphology is difficult [and] unreliable…intraspecific variability of the putatively distinctive anatomical measurements largely overlaps among species” and identifications are “further complicated by recent nomenclatorial revisions”. 

In their report, Pfenninger, Cordellier and Streit analyze morphology, mitochondrial COI and nuclear ITS-1 sequences, and describe breeding experiments with Radix snails collected at 60 sites throughout Europe. Using mtCOI sequences they found five MOTU (molecular operational taxonomic units), defined as “terminal clades with bootstrap support of 90% or more”. Populations of these MOTU overlapped broadly in geographic range and none corresponded to described species. Nuclear ITS sequences analyzed in a subset of specimens produced MOTU congruent with those generated by mtCOI.

ALL crosses between individuals from the same MOTU population were viable, whereas NONE of crosses between individuals from different MOTU produced eggs. In morphometric analysis, Radix MOTU overlapped as shown at left, and in rearing experiments, shell shape changed in 4 of 5 populations, demonstrating phenotypic plasticity of putative morphologic characters. In northwestern European Radix snails, DNA trumps morphology.

This work follows what might be a “best practices” pathway for single-locus mtDNA species discovery, aka DNA barcoding applied to species discovery:

1. COI sequence clusters (MOTU), found in analyzing  multiple specimens from geographically widespread locations, are proposed as putative species.

2. COI clustering is congruent with nuclear sequence data.

3. COI clusters show corresponding biological differences, such as morphologic characters, behavioral differences, or breeding incompatibility.

In some cases a virus or bacteria is recognized to be the causative agent even though not all of Koch’s postulates have been fulfilled. In a similar way in some cases it might be desirable to recognize mtDNA clusters as representing species without fulfilling all of the above criteria.

Criminal Feminism

Feminism involves the entry of women into occupations from which they were excluded.  A report on the imprisonment of women in the US provides data on convicted women that presumably reflect a rise in women pursuing criminal careers.  Plotted by our colleague Cesare Marchetti, the data fit a logistic curve beautifully and thus suggest the coherence of the feminist movement as a whole and its likeness to thousands of other biological and social growth processes.  The midpoint of the growth, which quadrupled women convicts, was 1992.  Have women hit a ceiling in their penetration of criminal occupations or will another pulse follow?  For more information on logistic growth processes, see Loglet Lab.

Visualizing large data sets

Growing barcode libraries challenge understanding. There are already about 200,000 mtCOI barcodes from about 25,000 species in BOLD, the Barcode of Life Data Systems Database. The burgeoning data sets hint at insights into biological diversity, revealed by looking at many species at once. A map of counties of the United States shows both large and small scale patterns, shaped by history, geography, and politics.  Viewed through the lens of mitochondrial variation, what would a map of species show? Are differences among and within species similar in birds and butterflies? Do species boundaries differ in marine vs terrestrial species, or in tropical vs. temperate zones? 

 

   

Here I offer one possible way of visualizing differences in barcode data sets using as an example the BOLD  “Hesperiidae of the ACG 1” Project containing 2,185 COI sequences from 355 species of skipper butterfly in ACG conservation area in Costa Rica (Hajibabei et al Proc Natl Acad Sci USA 2006 103:968). The BOLD-generated neighbor-joining tree of ACG Hesperiid COI sequences shown at left offers a traditional way of comparing sequences and is an essential step in looking at individual species and their close genetic neighbors. However the NJ tree contains only about 100 sequences from 20 species per page, and so runs to 22 pages. In the future it will likely be desirable to compare much larger data sets from, say, all 3700 known species of world skipper butterflies. 

For DNA barcoding, the essential information is differences among and within species. The higher-level groupings of species which are inevitably generated by a tree are of less interest. (In the following analysis distances are used simply to examine patterns of variation, NOT to determine whether they are sufficient for diagnosing species.) 

One useful approach is to generate histograms of differences within and between species. BOLD has a “Nearest Neighbor” analytic function which generates a table of mean and maximum variation within each species, “nearest neighbor” distance to the next closest species, and histogram summaries of the results.  

The histograms quickly show distances within most species are small and minimum distances between species are generally larger. Histograms are summaries with unlimited capacity. However, one might want to know more about individual species. For example, do species with higher intraspecific distances also show greater interspecific distances? One also wonders about the variation below 1% in both panels. In Beautiful Evidence, Edward Tufte points out histograms display relatively small amounts of data, usually 1 value per column. How to generate something with more information, more like the US Counties map, but not 22 pages long?

The graph at left uses the same 2 essential parameters: distance within each species and distance to nearest genetic neighbor. Because the usual distances within and between species are very different, plotting on a logarithmic scale allows one to inspect the variation in each set simultaneously. The results with 278 of the ACG skipper species (all those for which more than one individual was sampled, thereby generating a mean intraspecific distance) are shown. For each species, there is a black dot showing intraspecific distance and a red dot directly above or below showing distance to nearest neighbor. Sorting by intra- and interspecific distance allows the relative distances for each species to be seen. This graph highlights the relatively few species with nearest neighbor distances less than the mean intraspecific distance for that species. A line drawn at 1% appears to separate most of the intraspecific from interspecific values.

This graph is remarkably compressible, as shown by the small inset in the US county map above and in the figure at right. Here this is used to compare variation in Costa Rican skippers (278 species in 1 Family) to that in Australian fish (172 species in 1 Class) (Ward et al 2006 Phil Trans Royal Soc B 360:1471). The distribution of intraspecific variation seems quite similar while the nearest neighbor distances in fish are higher, presumably reflecting less dense sampling of a larger taxon. In the Fish paper, the red dots at bottom were thought to be ID errors, so perhaps some of the those in the skipper data set showing zero distance between species are taxonomic errors as well. This graphic approach could be useful in comparing patterning of intra- and inter-specific variation in marine vs terrestrial, tropical vs temperate, and allopatric vs sympatric species.

 

The spider and the fly: learning about applying COI to species identification

Two recent articles suggest how and how not to learn about applying mtCOI sequences to identifying species.  In Zoologica Scripta 2006 35:441 researchers from Koenig Zoological Research Museum, Bonn, analyze 113 specimens of 61 morphologically-defined species of pholcid (daddy long-legs) spiders. Important for this analysis and for future study, collection locations are given and voucher numbers are provided for each specimen and DNA extract.

(Some pholcid spiders vibrate in their webs when disturbed, moving so rapidly they become invisible; here is a wonderful video)

16s and COI sequences were successfully amplified using a single primer pair for each gene from 79% and 80% of specimens, respectively. It is striking that strong clustering within species was observed despite using short segments of mtDNA (COI, 312 bp; 16s 287 bp), which are less than half as long as the standard 648 bp COI barcode. In NJ trees with either mtDNA sequence, all morphologic conspecifics grouped together and were reciprocally monophyletic (ie no overlaps between species). Likely splits based on large intraspecific distances and differing geographic distributions were observed in 6 (25%) of 24 multiply-sampled morphospecies.

The authors go on to propose graphic and statistical metrics to calibrate how well simple distances can define species limits. They find that mtDNA distances will often diagnose species: “tree-based taxon clustering and statistical taxon analysis indicate that molecular evidence does coincide with morphological hypotheses” and “we disagree with [Meyer and Paulay’s] point that independently of the group of organisms studied, a “barcoding gap” between interspecific and intraspecific distance values would likely disappear in studies featuring both dense within-species sampling and closely related species”, ie distance-based clustering often corresponds to species limits. 

This study uses vouchered specimens from known locations and accurate modern sequencing technology, focuses on a relatively small clade (959 known species pholcid spiders), and analyzes in a positive way how distance measures might be used to define species, helping us learn about DNA barcoding as a tool for species identification.   

In another recent study Syst Biol 55:715 2006 researchers from National University of Singapore examine COI sequences deposited in GenBank from Diptera (flies, mosquitos, and gnats). They found 449 of the 150,000 known species of dipterans represented, with multiple sequences from 127 species, and analyze these to “test two key claims of molecular taxonomy”. The scientists found that there were often large differences in COI within species and also frequent overlaps between species such that some sequences were more closely related to those of another species than to conspecifics. The litany of failures is quite long, including “even when two COI sequences are identical, there is a 6% chance they belong to different species”.

I do not understand why the authors put so much effort into analyzing such a heterogeneous set of data, except that they are worried about molecular taxonomy in general and DNA barcoding in particular. To my reading this study suggests that many GenBank records contain errors, either because current morphologic taxonomy is incorrect (for example, study cited above suggests probable splits in 25% of pholcid spiders), specimens used for GenBank records are incorrectly identified, or because DNA sequences in GenBank contain errors due to human factors or older sequencing technology. There must be some limitations to COI barcode identification of dipteran species, mostly presumably closely-related young species, but this study has not shown where such problems might lie. 

I hope that future studies will use more of the “best practices” demonstrated in Astrin et al’s study of pholcid spiders and so help us learn more about how to apply COI sequences to species identification.    

Alexander (Sandy) Mather

We note with great sadness the passing of Alexander (Sandy) Mather on Tuesday, November 14, 2006, aged 63 years, Professor, Department of Geography and Environment, University of Aberdeen. Sandy passed away one day after the publication in PNAS of our joint paper on “Returning forests analyzed with the forest identity.” Our article in PNAS is a memorial to Sandy’s work. We first learned of Sandy while writing our papers about encroachment of farming on forests. His concept of the Forest Transition inspired us. Whenever anyone writes forest transition, they remember Sandy, pictured here leading a group in the Costa Rican forest. His passing coincides with the high-profile validation of his work over many years.