Taxonomic Classification of Bacterial 16S rRNA Genes Using Short Sequencing Reads: Evaluation of Effective Study Designs

Massively parallel high throughput sequencing technologies allow us to interrogate the microbial composition of biological samples at unprecedented resolution. The typical approach is to perform high-throughout sequencing of 16S rRNA genes, which are then taxonomically classified based on similarity to known sequences in existing databases. Current technologies cause a predicament though, because although they enable deep coverage of samples, they are limited in the length of sequence they can produce. As a result, high-throughout studies of microbial communities often do not sequence the entire 16S rRNA gene. The challenge is to obtain reliable representation of bacterial communities through taxonomic classification of short 16S rRNA gene sequences. In this study we explored properties of different study designs and developed specific recommendations for effective use of short-read sequencing technologies for the purpose of interrogating bacterial communities, with a focus on classification using naïve Bayesian classifiers. To assess precision and coverage of each design, we used a collection of ∼8,500 manually curated 16S rRNA gene sequences from cultured bacteria and a set of over one million bacterial 16S rRNA gene sequences retrieved from environmental samples, respectively. We also tested different configurations of taxonomic classification approaches using short read sequencing data, and provide recommendations for optimal choice of the relevant parameters. We conclude that with a judicious selection of the sequenced region and the corresponding choice of a suitable training set for taxonomic classification, it is possible to explore bacterial communities at great depth using current technologies, with only a minimal loss of taxonomic resolution.     

16S Classifier: A Tool for Fast and Accurate Taxonomic Classification of 16S rRNA Hypervariable Regions in Metagenomic Datasets

The diversity of microbial species in a metagenomic study is commonly assessed using 16S rRNA gene sequencing. With the rapid developments in genome sequencing technologies, the focus has shifted towards the sequencing of hypervariable regions of 16S rRNA gene instead of full length gene sequencing. Therefore, 16S Classifier is developed using a machine learning method, Random Forest, for faster and accurate taxonomic classification of short hypervariable regions of 16S rRNA sequence. It displayed precision values of up to 0.91 on training datasets and the precision values of up to 0.98 on the test dataset. On real metagenomic datasets, it showed up to 99.7% accuracy at the phylum level and up to 99.0% accuracy at the genus level. 16S Classifier is available freely at http://metagenomics.iiserb.ac.in/16Sclassifier and http://metabiosys.iiserb.ac.in/16Sclassifier.     

Prokaryotic Genome Size and SSU rDNA Copy Number: Estimation of Microbial Relative Abundance from a Mixed Population

Abstract Determination of the relative abundance of a specific prokaryote in an environmental sample is of major interest in applied and environmental microbiology. Relative abundance can be calculated using knowledge of SSU rDNA copy number, amount of SSU rDNA in the sample, and a weighted average estimate of the genome sizes for organisms in the original sample. By surveying the literature, we provide estimates of genome size and SSU rDNA copy number for 303 and 101 prokaryotes, respectively. This compilation can be used to make reasonable estimates for a wide range of organisms in the calculation of relative abundance. A statistical analysis suggests that no correlation exists between genome size and SSU rDNA copy number. A phylogenetic analysis is used to offer insights into the evolution of both genome size and SSU rDNA copy number. Received: 29 January 1999; Accepted: 29 April 1999     

Phylogeny and Classification of the Trapdoor Spider Genus Myrmekiaphila: An Integrative Approach to Evaluating Taxonomic Hypotheses

Background

Revised by Bond and Platnick in 2007, the trapdoor spider genus Myrmekiaphila comprises 11 species. Species delimitation and placement within one of three species groups was based on modifications of the male copulatory device. Because a phylogeny of the group was not available these species groups might not represent monophyletic lineages; species definitions likewise were untested hypotheses. The purpose of this study is to reconstruct the phylogeny of Myrmekiaphila species using molecular data to formally test the delimitation of species and species-groups. We seek to refine a set of established systematic hypotheses by integrating across molecular and morphological data sets.

Methods and Findings

Phylogenetic analyses comprising Bayesian searches were conducted for a mtDNA matrix composed of contiguous 12S rRNA, tRNA-val, and 16S rRNA genes and a nuclear DNA matrix comprising the glutamyl and prolyl tRNA synthetase gene each consisting of 1348 and 481 bp, respectively. Separate analyses of the mitochondrial and nuclear genome data and a concatenated data set yield M. torreya and M. millerae paraphyletic with respect to M. coreyi and M. howelli and polyphyletic fluviatilis and foliata species groups.

Conclusions

Despite the perception that molecular data present a solution to a crisis in taxonomy, studies like this demonstrate the efficacy of an approach that considers data from multiple sources. A DNA barcoding approach during the species discovery process would fail to recognize at least two species (M. coreyi and M. howelli) whereas a combined approach more accurately assesses species diversity and illuminates speciation pattern and process. Concomitantly these data also demonstrate that morphological characters likewise fail in their ability to recover monophyletic species groups and result in an unnatural classification. Optimizations of these characters demonstrate a pattern of “Dollo evolution” wherein a complex character evolves only once but is lost multiple times throughout the group''s history.     

Taxonomic Studies of the Diplazium mettenianum Complex (Woodsiaceae; Pteridophyta) in Japan: Morphology, Cytology and Taxonomy of Plants with Normal-Shaped Spores

Diplazium mettenianum complex, including two polymorphic species D. mettenianum and D. griffithii, with normal-shaped spores collected from 62 sites. An analysis of 20 qualitative morphological characters showed that the complex can be divided into five forms. A statistical analysis of 16 quantitative morphological characters supported the distinctness of five forms which are regarded as independent species: D. mettenianum, D. fauriei, D. deciduum nom. nov. (=D. mettenianum var. tenuifolium), D. griffithii, and D. hayatamae sp. nov. A cytological study on 150 plants showed that D. deciduum and plants in a single population of D. fauriei are sexual hexaploids (2n=6x=246=123ll), and all others are sexual tetraploids (2n=4x=164=82ll). Received 18 May 1998/ Accepted in revised form 3 December 1998     

Resolving Genetic Functions within Microbial Populations: In Situ Analyses Using rRNA and mRNA Stable Isotope Probing Coupled with Single-Cell Raman-Fluorescence In Situ Hybridization

Prokaryotes represent one-half of the living biomass on Earth, with the vast majority remaining elusive to culture and study within the laboratory. As a result, we lack a basic understanding of the functions that many species perform in the natural world. To address this issue, we developed complementary population and single-cell stable isotope (¹³C)-linked analyses to determine microbial identity and function in situ. We demonstrated that the use of rRNA/mRNA stable isotope probing (SIP) recovered the key phylogenetic and functional RNAs. This was followed by single-cell physiological analyses of these populations to determine and quantify in situ functions within an aerobic naphthalene-degrading groundwater microbial community. Using these culture-independent approaches, we identified three prokaryote species capable of naphthalene biodegradation within the groundwater system: two taxa were isolated in the laboratory (Pseudomonas fluorescens and Pseudomonas putida), whereas the third eluded culture (an Acidovorax sp.). Using parallel population and single-cell stable isotope technologies, we were able to identify an unculturable Acidovorax sp. which played the key role in naphthalene biodegradation in situ, rather than the culturable naphthalene-biodegrading Pseudomonas sp. isolated from the same groundwater. The Pseudomonas isolates actively degraded naphthalene only at naphthalene concentrations higher than 30 μM. This study demonstrated that unculturable microorganisms could play important roles in biodegradation in the ecosystem. It also showed that the combined RNA SIP-Raman-fluorescence in situ hybridization approach may be a significant tool in resolving ecology, functionality, and niche specialization within the unculturable fraction of organisms residing in the natural environment.     

From Genus to Phylum: Large-Subunit and Internal Transcribed Spacer rRNA Operon Regions Show Similar Classification Accuracies Influenced by Database Composition

We compared the classification accuracy of two sections of the fungal internal transcribed spacer (ITS) region, individually and combined, and the 5′ section (about 600 bp) of the large-subunit rRNA (LSU), using a naive Bayesian classifier and BLASTN. A hand-curated ITS-LSU training set of 1,091 sequences and a larger training set of 8,967 ITS region sequences were used. Of the factors evaluated, database composition and quality had the largest effect on classification accuracy, followed by fragment size and use of a bootstrap cutoff to improve classification confidence. The naive Bayesian classifier and BLASTN gave similar results at higher taxonomic levels, but the classifier was faster and more accurate at the genus level when a bootstrap cutoff was used. All of the ITS and LSU sections performed well (>97.7% accuracy) at higher taxonomic ranks from kingdom to family, and differences between them were small at the genus level (within 0.66 to 1.23%). When full-length sequence sections were used, the LSU outperformed the ITS1 and ITS2 fragments at the genus level, but the ITS1 and ITS2 showed higher accuracy when smaller fragment sizes of the same length and a 50% bootstrap cutoff were used. In a comparison using the larger ITS training set, ITS1 and ITS2 had very similar accuracy classification for fragments between 100 and 200 bp. Collectively, the results show that any of the ITS or LSU sections we tested provided comparable classification accuracy to the genus level and underscore the need for larger and more diverse classification training sets.     

Evolution and Taxonomic Classification of Alphapapillomavirus 7 Complete Genomes: HPV18, HPV39, HPV45, HPV59, HPV68 and HPV70

Background

The species Alphapapillomavirus 7 (alpha-7) contains human papillomavirus genotypes that account for 15% of invasive cervical cancers and are disproportionately associated with adenocarcinoma of the cervix. Complete genome analyses enable identification and nomenclature of variant lineages and sublineages.

Methods

The URR/E6 region was sequenced to screen for novel variants of HPV18, 39, 45, 59, 68, 70, 85 and 97 from 1147 cervical samples obtained from multiple geographic regions that had previously been shown to contain an alpha-7 HPV isolate. To study viral heterogeneity, the complete 8 kb genome of 128 isolates, including 109 sequenced for this analysis, were annotated and analyzed. Viral evolution was characterized by constructing phylogenic trees using maximum-likelihood and Bayesian algorithms. Global and pairwise alignments were used to calculate total and ORF/region nucleotide differences; lineages and sublineages were assigned using an alphanumeric system. The prototype genome was assigned to the A lineage or A1 sublineage.

Results

The genomic diversity of alpha-7 HPV types ranged from 1.1% to 6.7% nucleotide sequence differences; the extent of genome-genome pairwise intratype heterogeneity was 1.1% for HPV39, 1.3% for HPV59, 1.5% for HPV45, 1.6% for HPV70, 2.1% for HPV18, and 6.7% for HPV68. ME180 (previously a subtype of HPV68) was designated as the representative genome for HPV68 sublineage C1. Each ORF/region differed in sequence diversity, from most variable to least variable: noncoding region 1 (NCR1) / noncoding region 2 (NCR2) > upstream regulatory region (URR) > E6 / E7 > E2 / L2 > E1 / L1.

Conclusions

These data provide estimates of the maximum viral genomic heterogeneity of alpha-7 HPV type variants. The proposed taxonomic system facilitates the comparison of variants across epidemiological and molecular studies. Sequence diversity, geographic distribution and phylogenetic topology of this clinically important group of HPVs suggest an independent evolutionary history for each type.     

One Size Doesn't Fit All - RefEditor: Building Personalized Diploid Reference Genome to Improve Read Mapping and Genotype Calling in Next Generation Sequencing Studies

With rapid decline of the sequencing cost, researchers today rush to embrace whole genome sequencing (WGS), or whole exome sequencing (WES) approach as the next powerful tool for relating genetic variants to human diseases and phenotypes. A fundamental step in analyzing WGS and WES data is mapping short sequencing reads back to the reference genome. This is an important issue because incorrectly mapped reads affect the downstream variant discovery, genotype calling and association analysis. Although many read mapping algorithms have been developed, the majority of them uses the universal reference genome and do not take sequence variants into consideration. Given that genetic variants are ubiquitous, it is highly desirable if they can be factored into the read mapping procedure. In this work, we developed a novel strategy that utilizes genotypes obtained a priori to customize the universal haploid reference genome into a personalized diploid reference genome. The new strategy is implemented in a program named RefEditor. When applying RefEditor to real data, we achieved encouraging improvements in read mapping, variant discovery and genotype calling. Compared to standard approaches, RefEditor can significantly increase genotype calling consistency (from 43% to 61% at 4X coverage; from 82% to 92% at 20X coverage) and reduce Mendelian inconsistency across various sequencing depths. Because many WGS and WES studies are conducted on cohorts that have been genotyped using array-based genotyping platforms previously or concurrently, we believe the proposed strategy will be of high value in practice, which can also be applied to the scenario where multiple NGS experiments are conducted on the same cohort. The RefEditor sources are available at https://github.com/superyuan/refeditor.

This is a PLOS Computational Biology Software Article.

     

Complete mitochondrial genome of the bullhead torrent catfish, Liobagrus obesus (Siluriformes, Amblycipididae): Genome description and phylogenetic considerations inferred from the Cyt b and 16S rRNA genes

Mitochondrial DNA (mtDNA) from the bullhead torrent catfish, Liobagrus obesus, was isolated by long-polymerase chain reaction (Long-PCR) with universal primers and was fully sequenced by primer working using flanking sequences. The complete mtDNA from L. obesus was 16,531 bp in length and contained 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and a control region, demonstrating a structure very similar to that of other bony fish. An analysis of the protein-coding genes revealed a statistically substantiated bias in (T+C): (A+G) content, supporting earlier findings regarding this peculiarity. As indicated by a chi-square test, the observed scores for pyrimidine and purine content were different from those expected assuming a 50:50 ratio: chi(2)=41.63, d.f.=5, p<0.000001 for three categories, including the 1st, 2nd, and 3rd codon positions. Further, there was a difference in nucleotide content between ND6 and the other 12 protein-coding genes in L. obesus. The values of p-distances, as summarized for different scales of evolutionary history at the Cyt b gene, revealed a clear pattern of increased nucleotide diversity at four levels: (1) intraspecies, (2) intragenus, (3) intrafamily, and (4) intraorder. Scores of average p-distances of the four categories in catfish were (1) 1.59+/-0.54%, (2) 5.28+/-1.72% (3) 16.37+/-1.26%, and (4) 19.81+/-0.14%, respectively. These data support the hypothesis that speciation in the order Siluriformes, in most cases, follows a geographic mode through the accumulation of a numerous small genetic changes over a long time period. A phylogenetic tree for the bullhead torrent catfish and several other fish species belonging to the order Siluriformes was developed on the basis of respective Cyt b sequences (1138 bp); the analysis revealed a monophyletic origin for the five examined families. A species-specific clustering of sequences from single species was obtained, supporting additionally basic phylogenetic information for the catfish and the barcoding suitability of Cyt b sequence data. Lastly, one of the well-supported properties of our phylogenetic tree (99% repetition level in our analysis) was the monophyletic placement of all catfish (order Siluriformes) among other ray-finned fish of the class Actinopterigii. Also discussed herein are the aspects of phylogeny based on the 16S rRNA gene.     

Taxonomic redescriptions of two ciliates, Protogastrostyla pulchra n. g., n. comb. and Hemigastrostyla enigmatica (ciliophora: spirotrichea, stichotrichia), with phylogenetic analyses based on 18S and 28S rRNA gene sequences

The morphology and infraciliature of two stichotrichid ciliates, Gastrostyla pulchra(Perejaslawzewa 1886) Kahl, 1932 and Hemigastrostyla enigmatica(Dragesco and Dragesco-Kernéis 1986) Song & Wilbert, 1997, collected from marine and brackish sediments, were investigated by using living observations and protargol impregnations. Both 18S and 28S rRNA genes of these two species were sequenced. The 18S rDNA show high similarities (98.4%-99.7%) among populations of each species. There is about 94% similarity in 18S rDNA genes between G. pulchra and Gastrostyla steinii, the type species of the genus, which has been confirmed to be an oxytrichid by previous studies. In the phylogenetic trees of 18S, 28S, and combined 18S and 28S rDNA, both G. pulchra and H. enigmatica are consistently placed outside the well-established oxytrichid clade. Based on our analyses and previous ontogenetic data, we conclude that these two species may represent some lower groups in the subclass Stichotrichia, and that G. pulchra should represent a new genus, Protogastrostyla n. g. This new genus, which is morphologically similar to Gastrostyla, differs in its morphogenesis: the apical part of the old AZM is retained combining with the newly built membranelles that develop from the proter's oral primordium; the primary primordia of the dorsal kinety; and marginal primordia commence de novo without a definite contribution from the old structure.     

Taxonomic rearrangement of the genus Schizochytrium sensu lato based on morphology, chemotaxonomic characteristics, and 18S rRNA gene phylogeny (Thraustochytriaceae, Labyrinthulomycetes): emendation for Schizochytrium and erection of Aurantiochytrium and Oblongichytrium gen. nov.

The genus Schizochytrium sensu lato has been characterized by successive binary division of its vegetative cells. However, the molecular phylogeny strongly suggests that this genus is not a natural taxon, because the original and recorded strains that have been identified as Schizochytrium spp. separately form three well-supported monophyletic groups in the 18S rRNA gene tree. These three groups are clearly distinguishable by their combined morphological characteristics and the profiles of the polyunsaturated fatty acids and carotenoid pigments they contain, although these are hard to distinguish using only a single feature. Therefore, three different genera are proposed to accommodate these three groups, i.e., Schizochytrium sensu stricto, Aurantiochytrium, and Oblongichytrium gen. nov.     

The Shine and Dalgarno hypothesis for termination: the 3' terminus of the 16S rRNA of the Escherichia coli ribosome can be modified or base-paired with a complementary oligonucleotide without affecting termination in vitro

The occurrence of the nucleotides "...CCUUAOH" at the 3' terminus of the 16S rRNA of the small subunit of the Escherichia coli ribosome led to the suggestion that they may have a direct base pairing with the termination codon in the termination event of protein biosynthesis (Shine and Dalgarno 1974). We have examined this concept with two approaches, firstly using a 30S subunit whose 16S rRNA has been modified with a fluorescein moiety on the terminal adenosine together with the antibody against the moiety, and secondly with an oligonucleotide, UAAGG, complementary to the terminal pentanucleotide sequence of the rRNA. Collectively the data suggest that the nucleotides at the 3' terminus of 16S rRNA are not critically involved in base pairing during termination codon recognition.     

Taxonomic rearrangement of the genus Ulkenia sensu lato based on morphology, chemotaxonomical characteristics, and 18S rRNA gene phylogeny (Thraustochytriaceae, Labyrinthulomycetes): emendation for Ulkenia and erection of Botryochytrium, Parietichytrium, and Sicyoidochytrium gen. nov.

The genus Ulkenia is characterized by the naked protoplast stage within its life cycle. However, the 18S rRNA gene tree clearly shows that this genus is not a natural taxon, because our own isolates and reported strains separately form four well-supported monophyletic groups. These four groups are clearly distinguishable by their profiles of polyunsaturated fatty acids and carotenoid pigments and cell and colony morphology, e.g., persistence of sporangial wall, manner of the cell cleavage at the zoospore formation, and development of the ectoplasmic nets. Therefore, the four groups are assigned to four genera including three new genera, i.e., Ulkenia sensu stricto, Botryochytrium, Parietichytrium, and Sicyoidochytrium gen. nov.     

Treatment Extension of Pegylated Interferon Alpha and Ribavirin Does Not Improve SVR in Patients with Genotypes 2/3 without Rapid Virological Response (OPTEX Trial): A Prospective,Randomized, Two-Arm,Multicentre Phase IV Clinical Trial

Although sofosbuvir has been approved for patients with genotypes 2/3 (G2/3), many parts of the world still consider pegylated Interferon alpha (P) and ribavirin (R) as standard of care for G2/3. Patients with rapid virological response (RVR) show response rates >80%. However, SVR (sustained virological response) in non-RVR patients is not satisfactory. Longer treatment duration may be required but evidence from prospective trials are lacking. A total of 1006 chronic HCV genotype 2/3 patients treated with P/R were recruited into a German HepNet multicenter screening registry. Of those, only 226 patients were still HCV RNA positive at week 4 (non-RVR). Non-RVR patients with ongoing response after 24 weeks P-2b/R qualified for OPTEX, a randomized trial investigating treatment extension of additional 24 weeks (total 48 weeks, Group A) or additional 12 weeks (total 36 weeks, group B) of 1.5 μg/kg P-2b and 800-1400 mg R. Due to the low number of patients without RVR, the number of 150 anticipated study patients was not met and only 99 non-RVR patients (n=50 Group A, n=49 Group B) could be enrolled into the OPTEX trial. Baseline factors did not differ between groups. Sixteen patients had G2 and 83 patients G3. Based on the ITT (intention-to-treat) analysis, 68% [55%; 81%] in Group A and 57% [43%; 71%] in Group B achieved SVR (p= 0.31). The primary endpoint of better SVR rates in Group A compared to a historical control group (SVR 70%) was not met. In conclusion, approximately 23% of G2/3 patients did not achieve RVR in a real world setting. However, subsequent recruitment in a treatment-extension study was difficult. Prolonged therapy beyond 24 weeks did not result in higher SVR compared to a historical control group.

Trial Registration

ClinicalTrials.gov NCT00803309