Biodiversity of Philippine marine fishes: A DNA barcode reference library based on voucher specimens.

Accurate identification of fishes is essential for understanding their biology and to ensure food safety for consumers. DNA barcoding is an important tool because it can verify identifications of both whole and processed fishes that have had key morphological characters removed (e.g., filets, fish meal); however, DNA reference libraries are incomplete, and public repositories for sequence data contain incorrectly identified sequences. During a nine-year sampling program in the Philippines, a global biodiversity hotspot for marine fishes, we developed a verified reference library of cytochrome c oxidase subunit I (COI) sequences for 2,525 specimens representing 984 species. Specimens were primarily purchased from markets, with additional diversity collected using rotenone or fishing gear. Species identifications were verified based on taxonomic, phenotypic, and genotypic data, and sequences are associated with voucher specimens, live-color photographs, and genetic samples catalogued at Smithsonian Institution, National Museum of Natural History. The Biodiversity of Philippine Marine Fishes dataset is released herein to increase knowledge of species diversity and distributions and to facilitate accurate identification of market fishes.

Sampling multiple life stages significantly increases estimates of marine biodiversity.

Biodiversity assessments are critical for setting conservation priorities, understanding ecosystem function and establishing a baseline to monitor change. Surveys of marine biodiversity that rely almost entirely on sampling adult organisms underestimate diversity because they tend to be limited to habitat types and individuals that can be easily surveyed. Many marine animals have planktonic larvae that can be sampled from the water column at shallow depths. This life stage often is overlooked in surveys but can be used to relatively rapidly document diversity, especially for the many species that are rare or live cryptically as adults. Using DNA barcode data from samples of nemertean worms collected in three biogeographical regions-Northeastern Pacific, the Caribbean Sea and Eastern Tropical Pacific-we found that most species were collected as either benthic adults or planktonic larvae but seldom in both stages. Randomization tests show that this deficit of operational taxonomic units collected as both adults and larvae is extremely unlikely if larvae and adults were drawn from the same pool of species. This effect persists even in well-studied faunas. These results suggest that sampling planktonic larvae offers access to a different subset of species and thus significantly increases estimates of biodiversity compared to sampling adults alone. Spanish abstract is available in the electronic supplementary material.

Documenting decapod biodiversity in the Caribbean from DNA barcodes generated during field training in taxonomy.

DNA barcoding is a useful tool to identify the components of mixed or bulk samples, as well as to determine individuals that lack morphologically diagnostic features. However, the reference database of DNA barcode sequences is particularly sparsely populated for marine invertebrates and for tropical taxa. We used samples collected as part of two field courses, focused on graduate training in taxonomy and systematics, to generate DNA sequences of the barcode fragments of cytochrome c oxidase subunit I (COI) and mitochondrial ribosomal 16S genes for 447 individuals, representing at least 129 morphospecies of decapod crustaceans. COI sequences for 36% (51/140) of the species and 16S sequences for 26% (37/140) of the species were new to GenBank. Automatic Barcode Gap Discovery identified 140 operational taxonomic units (OTUs) which largely coincided with the morphospecies delimitations. Barcode identifications (i.e. matches to identified sequences) were especially useful for OTUs within Synalpheus, a group that is notoriously difficult to identify and rife with cryptic species, a number of which we could not identify to species, based on morphology. Non-concordance between morphospecies and barcode OTUs also occurred in a few cases of suspected cryptic species. As mitochondrial pseudogenes are particularly common in decapods, we investigate the potential for this dataset to include pseudogenes and discuss the utility of these sequences as species identifiers (i.e. barcodes). These results demonstrate that material collected and identified during training activities can provide useful incidental barcode reference samples for under-studied taxa.

Reference DNA barcodes and other mitochondrial markers for identifying Caribbean Octocorals.

DNA barcoding is a useful tool for documenting the diversity of metazoans. The most commonly used barcode markers, 16S and COI, are not considered suitable for species identification within some "basal" phyla of metazoans. Nevertheless metabarcoding studies of bulk mixed samples commonly use these markers and may obtain sequences for "basal" phyla. We sequenced mitochondrial DNA fragments of cytochrome oxidase c subunit I (COI), 16S ribosomal RNA (16S), NADH dehydrogenase subunits 2 (16S-ND2), 6 (ND6-ND3) and 4L (ND4L-MSH) for 27 species of Caribbean octocorals to create a reference barcode dataset and to compare the utility of COI and 16S to other markers more typically used for octocorals. The most common genera (Erythropodium, Ellisella, Briareum, Plexaurella, Muriceopsis and Pterogorgia) were effectively distinguished by small differences (5 or more substitutions or indels) in COI and 16S sequences. Gorgonia and Antillogorgia were effectively distinguished from each other by unique haplotypes, but the small genetic differences make distance approaches ineffective for these taxa. Plexaura, Pseudoplexaura and Eunicea were indistinguishable from each other but were generally effectively distinguished from other genera, further supporting the idea that these genera have undergone a rapid endemic radiation in the Caribbean.

Nuclear introns help unravel the diversification history of the Australo-Pacific Petroica robins.

Australo-Pacific Petroica robins are known for their striking variability in sexual plumage coloration. Molecular studies in recent years have revised the taxonomy of species and subspecies boundaries across the southwest Pacific and New Guinea. However, these studies have not been able to resolve phylogenetic relationships within Petroica owing to limited sampling of the nuclear genome. Here, we sequence five nuclear introns across all species for which fresh tissue was available. Nuclear loci offer support for major geographic lineages that were first inferred from mtDNA. We find almost no shared nuclear alleles between currently recognized species within the New Zealand and Australian lineages, whereas the Pacific robin radiation has many shared alleles. Multilocus coalescent species trees based on nuclear loci support a sister relationship between the Australian lineage and the Pacific robin radiation-a node that is poorly supported by mtDNA. We also find discordance in support for a sister relationship between the similarly plumaged Rose Robin (P. rosea) and Pink Robin (P. rodinogaster). Our nuclear data complement previous mtDNA studies in suggesting that the phenotypically cryptic eastern and western populations of Australia's Scarlet Robin (P. boodang) are genetically distinct lineages at the early stages of divergence and speciation.

DNA barcoding applied to ex situ tropical amphibian conservation programme reveals cryptic diversity in captive populations.

Amphibians constitute a diverse yet still incompletely characterized clade of vertebrates, in which new species are still being discovered and described at a high rate. Amphibians are also increasingly endangered, due in part to disease-driven threats of extinctions. As an emergency response, conservationists have begun ex situ assurance colonies for priority species. The abundance of cryptic amphibian diversity, however, may cause problems for ex situ conservation. In this study we used a DNA barcoding approach to survey mitochondrial DNA (mtDNA) variation in captive populations of 10 species of Neotropical amphibians maintained in an ex situ assurance programme at El Valle Amphibian Conservation Center (EVACC) in the Republic of Panama. We combined these mtDNA sequences with genetic data from presumably conspecific wild populations sampled from across Panama, and applied genetic distance-based and character-based analyses to identify cryptic lineages. We found that three of ten species harboured substantial cryptic genetic diversity within EVACC, and an additional three species harboured cryptic diversity among wild populations, but not in captivity. Ex situ conservation efforts focused on amphibians are therefore vulnerable to an incomplete taxonomy leading to misidentification among cryptic species. DNA barcoding may therefore provide a simple, standardized protocol to identify cryptic diversity readily applicable to any amphibian community.

Introduction to animal DNA barcoding protocols.

Procedures and protocols common to many DNA barcoding projects are summarized. Planning for any project should emphasize front-end procedures, especially the "genetic lockdown" of collected materials for downstream genetic procedures. Steps further into the DNA barcoding process chain, such as sequencing, data processing, and other back-end functions vary slightly, if at all, among projects and are presented elsewhere in the volume. Point-of-collection sample and tissue handling and data/metadata handling are stressed. Specific predictions of the future workflows and mechanics of DNA barcoding are difficult, so focus is on that which most or all future methods and technologies will surely share.

DNA barcoding fishes.

This chapter is an overview of the techniques for DNA barcoding of fishes from field collection to DNA sequence analysis. Recommendations for modifications of field protocols and best tissue sampling practices are made. A variety of DNA extraction protocols is provided, including high-throughput robot-assisted methods. A pair of well-tested forward and reverse primers for PCR amplification and sequencing are presented. These primers have been successfully used for DNA barcode on a wide array of marine fish taxa and also work well in most freshwater and cartilaginous fishes. Recipes and cycling protocols for both PCR amplification and sequencing and cleanup methods for the reaction products are provided. A method for the consistent production of high-quality DNA barcodes from DNA sequence data is given and stringent guidelines for judging the quality of raw sequence data are laid out.

Construction and analysis of phylogenetic trees using DNA barcode data.

The assembly of sequence data obtained from DNA barcodes into phylogenies or NJ trees has proven highly useful in estimating relatedness among species as well as providing a framework in which hypotheses regarding the evolution of traits or species distributions may be investigated. In this chapter, we outline the process by which DNA sequence data is assembled into a phylogenetically informative matrix, and then provide details on the methods to reconstruct NJ or phylogenetic trees that employ DNA barcode data, using only barcode data or combining barcodes with other data.

Multiple shifts between violet and ultraviolet vision in a family of passerine birds with associated changes in plumage coloration.

Colour vision in diurnal birds falls into two discrete classes, signified by the spectral sensitivity of the violet- (VS) or ultraviolet-sensitive (UVS) short wavelength-sensitive type 1 (SWS1) single cone. Shifts between sensitivity classes are rare; three or four are believed to have happened in the course of avian evolution, one forming UVS higher passerines. Such shifts probably affect the expression of shortwave-dominated plumage signals. We have used genomic DNA sequencing to determine VS or UVS affinity in fairy-wrens and allies, Maluridae, a large passerine family basal to the known UVS taxa. We have also spectrophotometrically analysed male plumage coloration as perceived by the VS and UVS vision systems. Contrary to any other investigated avian genus, Malurus (fairy-wrens) contains species with amino acid residues typical of either VS or UVS cone opsins. Three bowerbird species (Ptilonorhynchidae) sequenced for outgroup comparison carry VS opsin genes. Phylogenetic reconstructions render one UVS gain followed by one or more losses as the most plausible evolutionary scenario. The evolution of avian ultraviolet sensitivity is hence more complex, as a single shift no longer explains its distribution in Passeriformes. Character correlation analysis proposes that UVS vision is associated with shortwave-reflecting plumage, which is widespread in Maluridae.

A multigene phylogeny examining evolutionary and ecological relationships in the Australo-papuan wrens of the subfamily Malurinae (Aves).

Nucleotide sequences from four mitochondrial genes and three nuclear introns were used to examine phylogenetic relationships within the Australo-papuan fairy-wrens (Passeriformes: Maluridae: Malurinae). A well-resolved and well-supported phylogenetic hypothesis of all species in the subfamily was generated. The tree contained three clades corresponding to groups with similar plumages previously identified in earlier studies: the "bi-color," "blue," and "chestnut-shouldered" groups. The genus Malurus was not monophyletic -Malurusgrayi formed a clade with two New Guinean genera Sipodotus and Clytomyias. We recommend M. grayi be reclassified into the genus Chenorhamphus Oustalet 1898. One other taxonomic change is recommended based on the large genetic distance between the two subspecies of Chenorhamphus grayi - the elevation of C. g.campbelli to specific status (= C. campbelli). Although the family Maluridae appears to have had its origins in Australia, the DNA data supports a New Guinean origin for the Malurini (Sipodotus, Clytomyias, Chenorhamphus, Malurus).

Supertree bootstrapping methods for assessing phylogenetic variation among genes in genome-scale data sets.

Nonparamtric bootstrapping methods may be useful for assessing confidence in a supertree inference. We examined the performance of two supertree bootstrapping methods on four published data sets that each include sequence data from more than 100 genes. In "input tree bootstrapping," input gene trees are sampled with replacement and then combined in replicate supertree analyses; in "stratified bootstrapping," trees from each gene's separate (conventional) bootstrap tree set are sampled randomly with replacement and then combined. Generally, support values from both supertree bootstrap methods were similar or slightly lower than corresponding bootstrap values from a total evidence, or supermatrix, analysis. Yet, supertree bootstrap support also exceeded supermatrix bootstrap support for a number of clades. There was little overall difference in support scores between the input tree and stratified bootstrapping methods. Results from supertree bootstrapping methods, when compared to results from corresponding supermatrix bootstrapping, may provide insights into patterns of variation among genes in genome-scale data sets.

Prospects for building the tree of life from large sequence databases.

We assess the phylogenetic potential of approximately 300,000 protein sequences sampled from Swiss-Prot and GenBank. Although only a small subset of these data was potentially phylogenetically informative, this subset retained a substantial fraction of the original taxonomic diversity. Sampling biases in the databases necessitate building phylogenetic data sets that have large numbers of missing entries. However, an analysis of two "supermatrices" suggests that even data sets with as much as 92% missing data can provide insights into broad sections of the tree of life.

Phylogeny and evolution of the Australo-Papuan honeyeaters (Passeriformes, Meliphagidae).

We analyzed nucleotide variation at four loci for 75 species to produce a phylogenetic hypothesis for the Meliphagidae, and to examine the evolution and biogeographic history of the Meliphagidae. Both maximum parsimony and Bayesian methods of phylogenetic analysis were employed. The family was found to be monophyletic, though the genera Certhionyx, Anthochaera, and Phylidonyris were not. Four major clades were recovered and the spinebills (Acanthorhynchus) formed the sister clade to the remainder of the family in most analyses. The Australian endemic arid-adapted chats (Epthianura, Ashbyia) were found to be nested deeply within the family Meliphagidae. No evidence was found to support the hypothesis of separate New Guinean and Australian endemic radiations, nor of a close phylogenetic relationship between taxa from the New Guinea highlands and those from Australian northern rainforests.

The challenge of constructing large phylogenetic trees.

The amount of sequence data available to reconstruct the evolutionary history of genes and species has increased 20-fold in the past decade. Consequently the size of phylogenetic analyses has grown as well, and phylogenetic methods, algorithms and their implementations have struggled to keep pace. Computational and other challenges raised by this burgeoning database emerge at several stages of analysis, from the optimal assembly of large data matrices from sequence databases, to the efficient construction of trees from these large matrices and the piece-wise assembly of 'supertrees' from those trees in turn. A final challenge is posed by the difficulty of visualizing and making inferences from trees that might soon routinely contain thousands of species.

Obtaining maximal concatenated phylogenetic data sets from large sequence databases.

To improve the accuracy of tree reconstruction, phylogeneticists are extracting increasingly large multigene data sets from sequence databases. Determining whether a database contains at least k genes sampled from at least m species is an NP-complete problem. However, the skewed distribution of sequences in these databases permits all such data sets to be obtained in reasonable computing times even for large numbers of sequences. We developed an exact algorithm for obtaining the largest multigene data sets from a collection of sequences. The algorithm was then tested on a set of 100,000 protein sequences of green plants and used to identify the largest multigene ortholog data sets having at least 3 genes and 6 species. The distribution of sizes of these data sets forms a hollow curve, and the largest are surprisingly small, ranging from 62 genes by 6 species, to 3 genes by 65 species, with more symmetrical data sets of around 15 taxa by 15 genes. These upper bounds to sequence concatenation have important implications for building the tree of life from large sequence databases.