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

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

Primers for identification and polymorphism assessment of Vespertilionid bats in the Pacific Northwest

Species‐specific primers for portions of the mitochondrial DNA 16S ribosomal subunit gene were designed to identify bats of the family Vespertilionidae (Mammalia). Two fragments, each of 190 base pairs, were found to contain sufficient genetic variability to consistently resolve 10 of the 14 Pacific Northwest species included in this study. The remaining four species could be resolved into two sets of paired species that were unique from the other 10 species. These primers are demonstrably useful for purposes of surveys designed for species identification of bats, including using DNA extracted from guano collected from roost sites when bats are absent.     

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

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

Modular tagging of amplicons using a single PCR for high‐throughput sequencing

High‐throughput sequencing (HTS) of PCR amplicons is becoming the method of choice to sequence one or several targeted loci for phylogenetic and DNA barcoding studies. Although the development of HTS has allowed rapid generation of massive amounts of DNA sequence data, preparing amplicons for HTS remains a rate‐limiting step. For example, HTS platforms require platform‐specific adapter sequences to be present at the 5′ and 3′ end of the DNA fragment to be sequenced. In addition, short multiplex identifier (MID) tags are typically added to allow multiple samples to be pooled in a single HTS run. Existing methods to incorporate HTS adapters and MID tags into PCR amplicons are either inefficient, requiring multiple enzymatic reactions and clean‐up steps, or costly when applied to multiple samples or loci (fusion primers). We describe a method to amplify a target locus and add HTS adapters and MID tags via a linker sequence using a single PCR. We demonstrate our approach by generating reference sequence data for two mitochondrial loci (COI and 16S) for a diverse suite of insect taxa. Our approach provides a flexible, cost‐effective and efficient method to prepare amplicons for HTS.     

Ribosomal protein-dependent orientation of the 16 S rRNA environment of S15

Ribosomal protein S15 binds specifically to the central domain of 16 S ribosomal RNA (16 S rRNA) and directs the assembly of four additional proteins to this domain. The central domain of 16 S rRNA along with these five proteins form the platform of the 30 S subunit. Previously, directed hydroxyl radical probing from Fe(II)-S15 in small ribonucleoprotein complexes was used to study assembly of the central domain of 16 S rRNA. Here, this same approach was used to understand the 16 S rRNA environment of Fe(II)-S15 in 30 S subunits and to determine the ribosomal proteins that are involved in forming the mature S15-16 S rRNA environment. We have identified additional sites of Fe(II)-S15-directed cleavage in 30S subunits compared to the binary complex of Fe(II)-S15/16 S rRNA. Along with novel targets in the central domain, sites within the 5' and 3' minor domains are also cleaved. This suggests that during the course of 30S subunit assembly these elements are positioned in the vicinity of S15. Besides the previously determined role for S8, roles for S5, S6+S18, and S16 in altering the 16 S rRNA environment of S15 were established. These studies reveal that ribosomal proteins can alter the assembly of regions of the 30 S subunit from a considerable distance and influence the overall conformation of this ribonucleoprotein particle.     

Assembly of the central domain of the 30S ribosomal subunit: roles for the primary binding ribosomal proteins S15 and S8

Assembly of the 30S ribosomal subunit occurs in a highly ordered and sequential manner. The ordered addition of ribosomal proteins to the growing ribonucleoprotein particle is initiated by the association of primary binding proteins. These proteins bind specifically and independently to 16S ribosomal RNA (rRNA). Two primary binding proteins, S8 and S15, interact exclusively with the central domain of 16S rRNA. Binding of S15 to the central domain results in a conformational change in the RNA and is followed by the ordered assembly of the S6/S18 dimer, S11 and finally S21 to form the platform of the 30S subunit. In contrast, S8 is not part of this major platform assembly branch. Of the remaining central domain binding proteins, only S21 association is slightly dependent on S8. Thus, although S8 is a primary binding protein that extensively contacts the central domain, its role in assembly of this domain remains unclear. Here, we used directed hydroxyl radical probing from four unique positions on S15 to assess organization of the central domain of 16S rRNA as a consequence of S8 association. Hydroxyl radical probing of Fe(II)-S15/16S rRNA and Fe(II)-S15/S8/16S rRNA ribonucleoprotein particles reveal changes in the 16S rRNA environment of S15 upon addition of S8. These changes occur predominantly in helices 24 and 26 near previously identified S8 binding sites. These S8-dependent conformational changes are consistent with 16S rRNA folding in complete 30S subunits. Thus, while S8 binding is not absolutely required for assembly of the platform, it appears to affect significantly the 16S rRNA environment of S15 by influencing central domain organization.     

metaxa2: improved identification and taxonomic classification of small and large subunit rRNA in metagenomic data

The ribosomal rRNA genes are widely used as genetic markers for taxonomic identification of microbes. Particularly the small subunit (SSU; 16S/18S) rRNA gene is frequently used for species‐ or genus‐level identification, but also the large subunit (LSU; 23S/28S) rRNA gene is employed in taxonomic assignment. The metaxa software tool is a popular utility for extracting partial rRNA sequences from large sequencing data sets and assigning them to an archaeal, bacterial, nuclear eukaryote, mitochondrial or chloroplast origin. This study describes a comprehensive update to metaxa – metaxa 2 – that extends the capabilities of the tool, introducing support for the LSU rRNA gene, a greatly improved classifier allowing classification down to genus or species level, as well as enhanced support for short‐read (100 bp) and paired‐end sequences, among other changes. The performance of metaxa 2 was compared to other commonly used taxonomic classifiers, showing that metaxa 2 often outperforms previous methods in terms of making correct predictions while maintaining a low misclassification rate. metaxa 2 is freely available from http://microbiology.se/software/metaxa2/ .     

Identifying the true oysters (Bivalvia: Ostreidae) with mitochondrial phylogeny and distance-based DNA barcoding

Oysters (family Ostreidae), with high levels of phenotypic plasticity and wide geographic distribution, are a challenging group for taxonomists and phylogenetics. As a useful tool for molecular species identification, DNA barcoding might offer significant potential for oyster identification and taxonomy. This study used two mitochondrial fragments, cytochrome c oxidase I (COI) and the large ribosomal subunit (16S rDNA), to assess whether oyster species could be identified by phylogeny and distance-based DNA barcoding techniques. Relationships among species were estimated by the phylogenetic analyses of both genes, and then pairwise inter- and intraspecific genetic divergences were assessed. Species forming well-differentiated clades in the molecular phylogenies were identical for both genes even when the closely related species were included. Intraspecific variability of 16S rDNA overlapped with interspecific divergence. However, average intra- and interspecific genetic divergences for COI were 0-1.4% (maximum 2.2%) and 2.6-32.2% (minimum 2.2%), respectively, indicating the existence of a barcoding gap. These results confirm the efficacy of species identification in oysters via DNA barcodes and phylogenetic analysis.     

Phylogenetic study of Geitlerinema and Microcystis (Cyanobacteria) using PC‐IGS and 16S–23S ITS as markers: investigation of horizontal gene transfer

Selection of genes that have not been horizontally transferred for prokaryote phylogenetic inferences is regarded as a challenging task. The markers internal transcribed spacer of ribosomal genes (16S–23S ITS) and phycocyanin intergenic spacer (PC‐IGS), based on the operons of ribosomal and phycocyanin genes respectively, are among the most used markers in cyanobacteria. The region of the ribosomal genes has been considered stable, whereas the phycocyanin operon may have undergone horizontal transfer. To investigate the occurrence of horizontal transfer of PC‐IGS, phylogenetic trees of Geitlerinema and Microcystis strains were generated using PC‐IGS and 16S–23S ITS and compared. Phylogenetic trees based on the two markers were mostly congruent for Geitlerinema and Microcystis, indicating a common evolutionary history among ribosomal and phycocyanin genes with no evidence for horizontal transfer of PC‐IGS. Thus, PC‐IGS is a suitable marker, along with 16S–23S ITS for phylogenetic studies of cyanobacteria.     

A new plant protein interacts with eIF3 and 60S to enhance virus‐activated translation re‐initiation

The plant viral re‐initiation factor transactivator viroplasmin (TAV) activates translation of polycistronic mRNA by a re‐initiation mechanism involving translation initiation factor 3 (eIF3) and the 60S ribosomal subunit (60S). QJ;Here, we report a new plant factor—re‐initiation supporting protein (RISP)—that enhances TAV function in re‐initiation. RISP interacts physically with TAV in vitro and in vivo. Mutants defective in interaction are less active, or inactive, in transactivation and viral amplification. RISP alone can serve as a scaffold protein, which is able to interact with eIF3 subunits a/c and 60S, apparently through the C‐terminus of ribosomal protein L24. RISP pre‐bound to eIF3 binds 40S, suggesting that RISP enters the translational machinery at the 43S formation step. RISP, TAV and 60S co‐localize in epidermal cells of infected plants, and eIF3–TAV–RISP–L24 complex formation can be shown in vitro. These results suggest that RISP and TAV bridge interactions between eIF3‐bound 40S and L24 of 60S after translation termination to ensure 60S recruitment during repetitive initiation events on polycistronic mRNA; RISP can thus be considered as a new component of the cell translation machinery.     

Illumina-based analysis of microbial community diversity

Microbes commonly exist in milieus of varying complexity and diversity. Although cultivation-based techniques have been unable to accurately capture the true diversity within microbial communities, these deficiencies have been overcome by applying molecular approaches that target the universally conserved 16S ribosomal RNA gene. The recent application of 454 pyrosequencing to simultaneously sequence thousands of 16S rDNA sequences (pyrotags) has revolutionized the characterization of complex microbial communities. To date, studies based on 454 pyrotags have dominated the field, but sequencing platforms that generate many more sequence reads at much lower costs have been developed. Here, we use the Illumina sequencing platform to design a strategy for 16S amplicon analysis (iTags), and assess its generality, practicality and potential complications. We fabricated and sequenced paired-end libraries of amplified hyper-variable 16S rDNA fragments from sets of samples that varied in their contents, ranging from a single bacterium to highly complex communities. We adopted an approach that allowed us to evaluate several potential sources of errors, including sequencing artifacts, amplification biases, non-corresponding paired-end reads and mistakes in taxonomic classification. By considering each source of error, we delineate ways to make biologically relevant and robust conclusions from the millions of sequencing reads that can be readily generated by this technology.     

Phylogeny and adaptive radiation in the Onychopoda (Crustacea,Cladocera): evidence from multiple gene sequences

Members of the order Cladocera show remarkable morphological and ecological diversity. One of the most spectacular adaptive radiations in this group has involved species of the suborder Onychopoda, which have adopted a novel feeding strategy, predation, and have colonized habitats with a broad range of salinities. In order to evaluate the origins and systematics of this group, we derived a molecular phylogeny for its three component families including nine of 10 recognized genera based on three mitochondrial (mt) gene sequences: cytochrome c oxidase subunit I (COI), the ribosomal small and large subunits (12S and 16S) and one nuclear gene sequence: the small ribosomal subunit (18S). Maximum‐parsimony, maximum‐likelihood and neighbour‐joining phylogenetic analyses were largely congruent, supporting the monophyly of the suborder and each of its families. Comparative analyses of data based on total evidence and the conditional combination of the ribosomal genes produced relatively congruent patterns of phylogenetic affinity. By contrast, analyses of single gene results were inconsistent in recovering the monophyletic groups identified by the multigene analyses. Based on the reconstructed phylogeny, we discriminate among the existing hypotheses regarding the evolutionary history of the onychopods. We identify a prolonged episode of speciation from the Miocene to the Pleistocene with two pulses of diversification. We discuss our results with reference to the geological history of the Ponto‐Caspian basin, the region which fostered the onychopod radiation.     

Identification and Molecular Characterization of ‘Candidatus Phytoplasma mali’ Isolates in North‐western Italy

Apple proliferation (AP) is an important disease and is prevalent in several European countries. The causal agent of AP is ‘Candidatus Phytoplasma mali’ (‘Ca. Phytoplasma mali’). In this work, isolates of ‘Ca. Phytoplasma mali’ were detected and characterized through polymerase chain reaction (PCR) and restriction fragment length polymorphism (RFLP) analyses of 16S rRNA gene and non‐ribosomal DNA fragment. The presence of three AP subtypes (AT‐1, AT‐2 and AP‐15) was identified in 31 symptomatic apple trees and two samples each constituted by a pool of five insects, collected in north‐western Italy, where AT‐1 is a dominant subtype. Subsequent nucleotide sequence analysis of the PCR‐amplified 1.8 kb (P1/P7) fragment, containing the 16S rDNA, the 16S–23S intergenic ribosomal region and the 5′‐end of the 23S rDNA, revealed the presence of at least two phytoplasmal genetic lineages within the AT‐1 subtype, designed AT‐1a and AT‐1b. Moreover, in silico single nucleotide polymorphism (SNP) analysis based on 16S rDNA sequence can differentiate AT‐1 subtype from AT‐2 and AP‐15 subtypes. Our data showed a high degree of genetic diversity among ‘Ca. Phytoplasma mali’ population in north‐western Italy and underlined the possible use of the 16S rDNA analysis for the identification and the geographical origin assignation of isolates of AP phytoplasma. Molecular markers on 16S rDNA, here identified, could be useful for studying the epidemiology of AP disease.     

Analysis of conformational changes in 16 S rRNA during the course of 30 S subunit assembly

Ribosome biogenesis involves an integrated series of binding events coupled with conformational changes that ultimately result in the formation of a functional macromolecular complex. In vitro, Escherichia coli 30 S subunit assembly occurs in a cooperative manner with the ordered addition of 20 ribosomal proteins (r-proteins) with 16 S rRNA. The assembly pathway for 30 S subunits has been dissected in vitro into three steps, where specific r-proteins associate with 16 S rRNA early in 30 S subunit assembly, followed by a mid-assembly conformational rearrangement of the complex that then enables the remaining r-proteins to associate in the final step. Although the three steps of 30 S subunit assembly have been known for some time, few details have been elucidated about changes that occur as a result of these three specific stages. Here, we present a detailed analysis of the concerted early and late stages of small ribosomal subunit assembly. Conformational changes, roles for base-pairing and r-proteins at specific stages of assembly, and a polar nature to the assembly process have been revealed. This work has allowed a more comprehensive and global view of E.coli 30 S ribosomal subunit assembly to be obtained.     

Universal COI primers for DNA barcoding amphibians

DNA barcoding is a proven tool for the rapid and unambiguous identification of species, which is essential for many activities including the vouchering tissue samples in the genome 10K initiative, genealogical reconstructions, forensics and biodiversity surveys, among many other applications. A large‐scale effort is underway to barcode all amphibian species using the universally sequenced DNA region, a partial fragment of mitochondrial cytochrome oxidase subunit I COI. This fragment is desirable because it appears to be superior to 16S for barcoding, at least for some groups of salamanders. The barcoding of amphibians is essential in part because many species are now endangered. Unfortunately, existing primers for COI often fail to achieve this goal. Herein, we report two new pairs of primers (?, ?) that in combination serve to universally amplify and sequence all three orders of Chinese amphibians as represented by 36 genera. This taxonomic diversity, which includes caecilians, salamanders and frogs, suggests that the new primer pairs will universally amplify COI for the vast majority species of amphibians.     

Molecular identification and phylogenetic relationships of Coleoidea (Mollusca: Cephalopoda) from the Persian Gulf and Oman Sea reveals a case of cryptic diversity

The Persian Gulf and Oman Sea constitute one of the most important marine ecosystems and have many economically important aquatic species, including several coleoid cephalopods. Some coleoids are difficult to identify using traditional morphological characteristics. In this study, two mitochondrial fragments, cytochrome oxidase I (COI) and the large ribosomal subunit (16S rRNA), were used for identification of coleoid species in four regions in the northern Persian Gulf and Oman Sea. The study led to the identification of potential cryptic species of Sepia, Amphioctopus and Uroteuthis. Furthermore, Euprymna hyllebergi was reported for the first time from the Persian Gulf. A high diversity of Coeloidea was found in the study area. Mean intraspecific and interspecific nucleotide distances for COI were 0%–2% and 2%–7%, respectively, while these values for 16S rRNA sequences were 0%–1% and 1%–4%. Given the uncertainty about species identity and the high levels of intraspecific genetic diversity reported for some species in GenBank, a comprehensive global study will be needed to resolve the taxonomic status of several coleoid species.     

Taxon-specific PCR primers to detect two inconspicuous arbuscular mycorrhizal fungi from temperate agricultural grassland

Taxon-specific polymerase chain reaction (PCR) primers enable detection of arbuscular mycorrhizal fungi (AMF, Glomeromycota) in plant roots where the fungi lack discriminative morphological and biochemical characters. We designed and validated pairs of new PCR primers targeted to the flanking regions of the variable domain 1 of the nuclear ribosomal large subunit RNA gene to specifically detect Acaulospora paulinae and an undescribed member of the Diversisporaceae. These two fungal taxa, sporulating late in soil-trap cultures and showing small, faintly coloured spores and weakly staining intraradical structures, were frequently found in roots of Trifolium repens from a high-input agricultural grassland. The newly developed PCR primers may thus enable studies on two inconspicuous AMF taxa that appear to have been overlooked in previous molecular AMF community analyses and for which no specific PCR primers have been published.     

SELECTIVE RECOVERY OF MICROALGAE FROM DIVERSE HABITATS USING “PHYTO‐SPECIFIC” 16S rDNA PRIMERS1

Environmental PCR is a common tool for surveying aquatic microalgae; however, universal primers generally employed are not specific to phytoplankton and typically recover nonphotosynthetic bacteria at high frequencies. Using a 16S rDNA “phyto‐specific” primer, we were able to selectively amplify sequences of photosynthetic species from several mixed aquatic samples, even when large numbers of nonphotosynthetic microorganisms were present. We identified 21 microalgal sequences from three different habitats: salt marshes in Virginia, river basins in North Carolina, and sea ice in Alaska. In contrast, universal 16S primers recovered a majority of nonphotosynthetic organisms from some of the same samples. Our results indicate that phytoplankton‐specific primers are efficient in selectively amplifying a broad diversity of microalgae in mixed environmental samples and, therefore, can reduce the noise from extraneous species that often dominates molecular surveys of aquatic samples.     

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

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