Performance comparison of genetic markers for high‐throughput sequencing‐based biodiversity assessment in complex communities

Metabarcode surveys of DNA extracted from environmental samples are increasingly popular for biodiversity assessment in natural communities. Such surveys rely heavily on robust genetic markers. Therefore, analysis of PCR efficiency and subsequent biodiversity estimation for different types of genetic markers and their corresponding primers is important. Here, we test the PCR efficiency and biodiversity recovery potential of three commonly used genetic markers – nuclear small subunit ribosomal DNA (18S), mitochondrial cytochrome c oxidase subunit I (COI) and 16S ribosomal RNA (mt16S) – using 454 pyrosequencing of a zooplankton community collected from Hamilton Harbour, Ontario. We found that biodiversity detection power and PCR efficiency varied widely among these markers. All tested primers for COI failed to provide high‐quality PCR products for pyrosequencing, but newly designed primers for 18S and 16S passed all tests. Furthermore, multiple analyses based on large‐scale pyrosequencing (i.e. 1/2 PicoTiter plate for each marker) showed that primers for 18S recover more (38 orders) groups than 16S (10 orders) across all taxa, and four vs. two orders and nine vs. six families for Crustacea. Our results showed that 18S, using newly designed primers, is an efficient and powerful tool for profiling biodiversity in largely unexplored communities, especially when amplification difficulties exist for mitochondrial markers such as COI. Universal primers for higher resolution markers such as COI are still needed to address the possible low resolution of 18S for species‐level identification.     

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.     

Universal Primers for Amplification of Mitochondrial Small Subunit Ribosomal RNA-Encoding Gene in Scleractinian Corals

We describe the construction of polymerase chain reaction primers designed to amplify a portion of the mitochondrial (mt) small subunit ribosomal (SSU) RNA-encoding genes in scleractinian corals. Combinations of cloning and sequencing show that the amplified fragments are between 694 and 896 bp in length. Alignment of the amplified DNA sequences to the published mt SSU rRNA genes of Metridium senile and Sarcophyton glaucum indicates several conserved regions among actiniarian, corallimorpharian, octocorallian, and scleractinians, suggesting this primer set can successfully amplify over 80% of the mt SSU rDNA region of scleractinian corals. Surveys of sequence variation and estimation of the rate of evolution show an extremely slow divergence of the SSU rRNA gene in the family Acroporidae. Received June 11, 1999; accepted October 4, 1999.     

Mitogenome polymorphism in a single branch sample revealed by SOLiD deep sequencing of the Lophelia pertusa coral genome

We present an initial genomic analysis of the non-symbiotic scleractinian coral Lophelia pertusa, the dominant cold-water reef-building coral species in the North Atlantic Ocean. A significant fraction of the deep sequencing reads was of mitochondrial and microbial origins. SOLiD deep sequencing reads from fragment library experiments of total DNA and PCR amplified mitogenome generated about 21,000 times and 136,000 times coverage, respectively, of the 16,150bp mitogenome. Five polymorphic sites that include two non-synonymous sites in the NADH dehydrogenase subunit 5 genes were detected in both experiments. This observation is surprising since anthozoans in general exhibit very low mtDNA sequence variation at intraspecific level compared to nuclear sequences. More than fifty bacterial species associated with the coral isolate were also sequence detected, representing at least ten complete genomes. Most reads, however, were predicted to originate from the Lophelia nuclear genome.     

Polymorphism of the Gene for Subunit 6 of the NADH Dehydrogenase Complex in Ethnic Russians of Russia

A sample of ethnic Russians of Russia was tested for polymorphism of the NADH dehydrogenase subunit 6 (ND6) gene mapping to the mtDNA region 14,170–14,569. Genetic diversity of ND6 haplotypes was estimated at 0.406, and probability of haplotype random match, at 0.598. Combined with typing the mtDNA control region, analysis of the ND6 gene polymorphism was assumed to improve the reliability of forensic identification. Several point substitutions in the ND6 gene region proved to be associated with particular transitions in the mtDNA control region; the association was characterized with the coefficient.     

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.     

16S rRNA序列分析法在医学微生物鉴定中的应用

１６Ｓ ｒＲＮＡ序列分析作为微生物系统分类的主要依据已得到了广泛认同,随着微生物核糖体数据库的日益完善,该技术成为细菌分类和鉴定的一个有力工具。本文概述了 １６５ ｒＲＮＡ序列分析法的技术步骤以及该技术在医学微生物研究中的应用,总结了目前文献报导的各种致病微生物种属特异性 １６５ ｒＲＮＡ引物和探针序列,同时分析了该技术在应用中存在的一些问题。     

Microbial community dynamics based on 16S rRNA gene profiles in a Pacific Northwest estuary and its tributaries

We analyzed bacterioplankton community structure in Tillamook Bay, Oregon and its tributaries to evaluate phylogenetic variability and its relation to changes in environmental conditions along an estuarine gradient. Using eubacterial primers, we amplified 16S rRNA genes from environmental DNA and analyzed the PCR products by length heterogeneity polymerase chain reaction (LH-PCR), which discriminates products based on naturally occurring length differences. Analysis of LH-PCR profiles by multivariate ordination methods revealed differences in community composition along the estuarine gradient that were correlated with changes in environmental variables. Microbial community differences were also detected among different rivers. Using partial 16S rRNA sequences, we identified members of dominant or unique gene fragment size classes distributed along the estuarine gradient. Gammaproteobacteria and Betaproteobacteria and members of the Bacteroidetes dominated in freshwater samples, while Alphaproteobacteria, Cyanobacteria and chloroplast genes dominated in marine samples. Changes in the microbial communities correlated most strongly with salinity and dissolved silicon, but were also strongly correlated with precipitation. We also identified specific gene fragments that were correlated with inorganic nutrients. Our data suggest that there is a significant and predictable change in microbial species composition along an estuarine gradient, shifting from a more complex community structure in freshwater habitats to a community more typical of open ocean samples in the marine-influenced sites. We also demonstrate the resolution and power of LH-PCR and multivariate analyses to provide a rapid assessment of major community shifts, and show how these shifts correlate with environmental variables.     

Molecular phylogeny of parasitic Platyhelminthes based on sequences of partial 28S rDNA D1 and mitochondrial cytochrome c oxidase subunit I

The phylogenic relationships existing among 14 parasitic Platyhelminthes in the Republic of Korea were investigated via the use of the partial 28S ribosomal DNA (rDNA) D1 region and the partial mitochondrial cytochrome c oxidase subunit 1 (mCOI) DNA sequences. The nucleotide sequences were analyzed by length, G + C %, nucleotide differences and gaps in order to determine the analyzed phylogenic relationships. The phylogenic patterns of the 28S rDNA D1 and mCOI regions were closely related within the same class and order as analyzed by the PAUP 4.0 program, with the exception of a few species. These findings indicate that the 28S rDNA gene sequence is more highly conserved than are the mCOI gene sequences. The 28S rDNA gene may prove useful in studies of the systematics and population genetic structures of parasitic Platyhelminthes.     

Identification of mycorrhizal fungi from single pelotons of Dactylorhiza majalis (Orchidaceae) using single-strand conformation polymorphism and mitochondrial ribosomal large subunit DNA sequences

The mitochondrial ribosomal large subunit (Ls) DNA was used to identify the orchid mycorrhizal fungi found in roots of Dactylorhiza majalis. The gene was amplified using DNA extracted from single pelotons obtained from fresh and silica gel dried roots. Furthermore, sequencing a variety of well-characterized orchid isolates expanded the fungal database of the mitochondrial ribosomal LsDNA. Polymerase chain reaction product length variants present in D. majalis were sequenced and identified using the expanded database. These analyses revealed two different peloton-forming fungi in samples from D. majalis, which sometimes occurred together as a single two-taxa peloton within the same cortex cell. The first taxon belonged to the genus Tulasnella and the second taxon was distantly related to Laccaria.     

Optimizing the use of shed feathers for genetic analysis

Shed feathers obtained by noninvasive genetic sampling (NGS) are a valuable source of DNA for genetic studies of birds. They can be collected across a large geographical range and facilitate research on species that would otherwise be extremely difficult to study. A limitation of this approach is uncertainty concerning the quality of the extracted DNA. Here we investigate the relationship between feather type, feather condition and DNA quality (amplification success) in order to provide a simple, cost-effective method for screening samples prior to genetic analysis. We obtained 637 shed feathers of the powerful owl (Ninox strenua) from across its range in southeastern Australia. The extracted DNA was amplified using polymerase chain reaction for a range of markers including mitochondrial DNA, ND3 and nuclear DNA, a simple sequence repeat (Nst02) and a portion of the CHD-1 gene (P2/P8). We found that feather condition significantly influenced the amplification success of all three loci, with feathers characterized as 'good' having greater success. Feather type was found to be of lower importance, with good quality feathers of all types consistently producing high success for all three loci. We also found that the successful amplification of multilocus genotypes was dependant on the condition of the starting material and was highly correlated with successful amplification of the sex-linked CHD-1 locus. Samples with low DNA quality have a higher probability of amplification failure and are more likely to produce incorrect genotypes; therefore, identifying samples with high DNA quality can save substantial time and cost associated with the genetic analysis of NGS. As a result, we propose a method for screening shed feathers in order to provide a subset of samples which will have a greater probability of containing high quality DNA suitable for the amplification of multilocus genotypes.     

Development of a selective and sensitive polymerase chain reaction assay for the detection of Mycoplasma pirum

Abstract A new assay using the polymerase chain reaction to amplify a 173-nucleotide DNA fragment within the 16S ribosomal RNA gene of Mycoplasma pirum has been developed. The assay selectively amplified DNA from all strains of M. pirum tested with a high level of sensitivity, even in a context of human DNA. DNA from other mollicute species, including those closely related to M. pirum , from bacteria phylogenetically close to mollicutes ( Clostridium innocuum, C. ramosum and Bacillus subtilis ), from Escherichia coli and from human peripheral blood mononuclear cells, did not produce the amplified DNA product specific for M. pirum .     

Vector soup: high‐throughput identification of Neotropical phlebotomine sand flies using metabarcoding

Phlebotomine sand flies are haematophagous dipterans of primary medical importance. They represent the only proven vectors of leishmaniasis worldwide and are involved in the transmission of various other pathogens. Studying the ecology of sand flies is crucial to understand the epidemiology of leishmaniasis and further control this disease. A major limitation in this regard is that traditional morphological‐based methods for sand fly species identifications are time‐consuming and require taxonomic expertise. DNA metabarcoding holds great promise in overcoming this issue by allowing the identification of multiple species from a single bulk sample. Here, we assessed the reliability of a short insect metabarcode located in the mitochondrial 16S rRNA for the identification of Neotropical sand flies, and constructed a reference database for 40 species found in French Guiana. Then, we conducted a metabarcoding experiment on sand flies mixtures of known content and showed that the method allows an accurate identification of specimens in pools. Finally, we applied metabarcoding to field samples caught in a 1‐ha forest plot in French Guiana. Besides providing reliable molecular data for species‐level assignations of phlebotomine sand flies, our study proves the efficiency of metabarcoding based on the mitochondrial 16S rRNA for studying sand fly diversity from bulk samples. The application of this high‐throughput identification procedure to field samples can provide great opportunities for vector monitoring and eco‐epidemiological studies.     

The influence of diet on faecal DNA amplification and sex identification in brown bears (Ursus arctos)

To evaluate the influence of diet on faecal DNA amplification, 11 captive brown bears (Ursus arctos) were placed on six restricted diets: grass (Trifolium spp., Haplopappus hirtus and Poa pratensis), alfalfa (Lupinus spp.), carrots (Daucus spp.), white-tailed deer (Odocoileus virginianus), blueberries (Vaccinium spp.) and salmon (Salmo spp.). DNA was extracted from 50 faecal samples of each restricted diet, and amplification of brown bear DNA was attempted for a mitochondrial DNA (mtDNA) locus and nuclear DNA (nDNA) locus. For mtDNA, no significant differences were observed in amplification success rates across diets. For nDNA, amplification success rates for salmon diet extracts were significantly lower than all other diet extracts (P < 0.001). To evaluate the accuracy of faecal DNA sex identification when female carnivores consume male mammalian prey, female bears were fed male white-tailed deer. Four of 10 extracts amplified, and all extracts were incorrectly scored as male due to amplification of X and Y-chromosome fragments. The potential biases highlighted in this study have broad implications for researchers using faecal DNA for individual and sex identification, and should be evaluated in other species.     

Rapid identification of marine bioluminescent bacteria by amplified 16S ribosomal RNA gene restriction analysis

To rapidly identify natural isolates of marine bioluminescent bacteria, we developed amplified ribosomal DNA restriction analysis (ARDRA) methods. ARDRA, which is based on the restriction patterns of 16S rRNA gene digested with five enzymes (EcoRI, DdeI, HhaI, HinfI, RsaI), clearly distinguished the 14 species of marine bioluminescent bacteria currently known, which belong to the genera Vibrio, Photobacterium, and Shewanella. When we applied ARDRA to 129 natural isolates from two cruises in Sagami Bay, Japan, 127 were grouped into six ARDRA types with distinctive restriction patterns; these isolates represented the bioluminescent species, P. angustum, P. leiognathi, P. phosphoreum, S. woodyi, V. fischeri, and V. harveyi. The other two isolates showing unexpected ARDRA patterns turned out to have 16S rRNA gene sequences similar to P. leiognathi and P. phosphoreum. Nevertheless, ARDRA provides a simple and fairly robust means for rapid identification of the natural isolates of marine bioluminescent bacteria, and is therefore useful in studying their diversity.     

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

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

Characterization of the predominant bacterial population of different mangrove rhizosphere soils using 16S rRNA gene-based single-strand conformation polymorphism (SSCP)

Variations in chemical parameters and bacterial populations in mangrove rhizosphere samples were noted for different sites. The C, N, P and K contents as well as pH, EC and salinity showed variation between sites. Significant differences in soil properties were also found in sampling sites. Two types of soil were noted among sites. Guesthouse had significantly higher organic matter and nutrient content (N) than other three sites suggesting that human discharges, litter deposition and surface runoff were major nutrient inputs. This contaminated site was located at the landward edges. Positive correlations between organic matter, N, P and K contents were found suggesting that these nutrients were from similar input sources. Effects of sampling sites on microbial diversity were also analyzed via SSCP. Porteresia coarctata and Rhizophora mucronata did not show any variation in the banding patterns among replicates sampled in short distance within site. But Sonneratia apetala showed variation among replicates sampled in distance within site. A significant variation was noted in the SSCP profile among replicates between sites. The majority of dominant SSCP band sequences were related to bacterial genera of root and root-free soil environments, namely Bacillus, Planococcus, Planomicrobium, low G+C Gram-positive bacterium, glacial ice bacterium and unidentified bacteria. In the analysis of 16S rRNA sequences, members belonging to the phylum Firmicutes dominated the sequence collection. The phylogenetic analysis of 16S rRNA gene sequences showed close relationships to a wide range of clones or bacterial species of phylum Firmicutes and unidentified bacteria.