Introducing W.A.T.E.R.S.: a Workflow for the Alignment,Taxonomy, and Ecology of Ribosomal Sequences

Background  

For more than two decades microbiologists have used a highly conserved microbial gene as a phylogenetic marker for bacteria and archaea. The small-subunit ribosomal RNA gene, also known as 16 S rRNA, is encoded by ribosomal DNA, 16 S rDNA, and has provided a powerful comparative tool to microbial ecologists. Over time, the microbial ecology field has matured from small-scale studies in a select number of environments to massive collections of sequence data that are paired with dozens of corresponding collection variables. As the complexity of data and tool sets have grown, the need for flexible automation and maintenance of the core processes of 16 S rDNA sequence analysis has increased correspondingly.     

Homogeneity of the 16S rDNA sequence among geographically disparate isolates of Taylorella equigenitalis

Background

At present, six accessible sequences of 16S rDNA from Taylorella equigenitalis (T. equigenitalis) are available, whose sequence differences occur at a few nucleotide positions. Thus it is important to determine these sequences from additional strains in other countries, if possible, in order to clarify any anomalies regarding 16S rDNA sequence heterogeneity. Here, we clone and sequence the approximate full-length 16S rDNA from additional strains of T. equigenitalis isolated in Japan, Australia and France and compare these sequences to the existing published sequences.

Results

Clarification of any anomalies regarding 16S rDNA sequence heterogeneity of T. equigenitalis was carried out. When cloning, sequencing and comparison of the approximate full-length 16S rDNA from 17 strains of T. equigenitalis isolated in Japan, Australia and France, nucleotide sequence differences were demonstrated at the six loci in the 1,469 nucleotide sequence. Moreover, 12 polymorphic sites occurred among 23 sequences of the 16S rDNA, including the six reference sequences.

Conclusion

High sequence similarity (99.5% or more) was observed throughout, except from nucleotide positions 138 to 501 where substitutions and deletions were noted.     

FastGroup: A program to dereplicate libraries of 16S rDNA sequences

Background

Ribosomal 16S DNA sequences are an essential tool for identifying and classifying microbes. High-throughput DNA sequencing now makes it economically possible to produce very large datasets of 16S rDNA sequences in short time periods, necessitating new computer tools for analyses. Here we describe FastGroup, a Java program designed to dereplicate libraries of 16S rDNA sequences. By dereplication we mean to: 1) compare all the sequences in a data set to each other, 2) group similar sequences together, and 3) output a representative sequence from each group. In this way, duplicate sequences are removed from a library.

Results

FastGroup was tested using a library of single-pass, bacterial 16S rDNA sequences cloned from coral-associated bacteria. We found that the optimal strategy for dereplicating these sequences was to: 1) trim ambiguous bases from the 5' end of the sequences and all sequence 3' of the conserved Bact517 site, 2) match the sequences from the 3' end, and 3) group sequences >=97% identical to each other.

Conclusions

The FastGroup program simplifies the dereplication of 16S rDNA sequence libraries and prepares the raw sequences for subsequent analyses.     

Phylogenetic analysis of the true water bugs (Insecta: Hemiptera: Heteroptera: Nepomorpha): evidence from mitochondrial genomes

Background  

The true water bugs are grouped in infraorder Nepomorpha (Insecta: Hemiptera: Heteroptera) and are of great economic importance. The phylogenetic relationships within Nepomorpha and the taxonomic hierarchies of Pleoidea and Aphelocheiroidea are uncertain. Most of the previous studies were based on morphological characters without algorithmic assessment. In the latest study, the molecular markers employed in phylogenetic analyses were partial sequences of 16S rDNA and 18S rDNA with a total length about 1 kb. Up to now, no mitochondrial genome of the true water bugs has been sequenced, which is one of the largest data sets that could be compared across animal taxa. In this study we analyzed the unresolved problems in Nepomorpha using evidence from mitochondrial genomes.     

Telonemia-specific environmental 18S rDNA PCR reveals unknown diversity and multiple marine-freshwater colonizations

Background  

Recent surveys of eukaryote 18S rDNA diversity in marine habitats have uncovered worldwide distribution of the heterotrophic eukaryote phylum Telonemia. Here we investigate the diversity and geographic distribution of Telonemia sequences by in-depth sequencing of several new 18S rDNA clone libraries from both marine and freshwater sites by using a Telonemia-specific PCR strategy.     

The biochemical differentiation of Enterobacter sakazakii genotypes

Background  

Enterobacter sakazakii is an emergent pathogen that has been associated with neonatal infections through contaminated powdered infant milk formula. The species was defined by Farmer et al. (1980) who described 15 biogroups according to the biochemical characterization of 57 strains. This present study compares genotypes (DNA cluster groups based on partial 16S rDNA sequence analysis) with the biochemical traits for 189 E. sakazakii strains.     

Analysis of 16S rRNA gene sequences and circulating cell-free DNA from plasma of chronic fatigue syndrome and non-fatigued subjects

Background  

The association of an infectious agent with chronic fatigue syndrome (CFS) has been difficult and is further complicated by the lack of a known lesion or diseased tissue. Cell-free plasma DNA could serve as a sentinel of infection and disease occurring throughout the body. This type of systemic sample coupled with broad-range amplification of bacterial sequences was used to determine whether a bacterial pathogen was associated with CFS. Plasma DNA from 34 CFS and 55 non-fatigued subjects was assessed to determine plasma DNA concentration and the presence of bacterial 16S ribosomal DNA (rDNA) sequences.     

Symbiotic functioning and bradyrhizobial biodiversity of cowpea (<Emphasis Type="Italic">Vigna unguiculata</Emphasis>L. Walp.) in Africa

Background  

Cowpea is the most important food grain legume in Sub-Saharan Africa. However, no study has so far assessed rhizobial biodiversity and/or nodule functioning in relation to strain IGS types at the continent level. In this study, 9 cowpea genotypes were planted in field experiments in Botswana, South Africa and Ghana with the aim of i) trapping indigenous cowpea root-nodule bacteria (cowpea "rhizobia") in the 3 countries for isolation, molecular characterisation using PCR-RFLP analysis, and sequencing of the 16S - 23S rDNA IGS gene, ii) quantifying N-fixed in the cowpea genotypes using the ¹⁵N natural abundance technique, and iii) relating the levels of nodule functioning (i.e. N-fixed) to the IGS types found inside nodules.     

The microbiology of Bandji,palm wine of Borassus akeassii from Burkina Faso: identification and genotypic diversity of yeasts,lactic acid and acetic acid bacteria

Aim

To investigate physicochemical characteristics and especially genotypic diversity of the main culturable micro‐organisms involved in fermentation of sap from Borassus akeassii, a newly identified palm tree from West Africa.

Methods and Results

Physicochemical characterization was performed using conventional methods. Identification of micro‐organisms included phenotyping and sequencing of: 26S rRNA gene for yeasts, 16S rRNA and gyrB genes for lactic acid bacteria (LAB) and acetic acid bacteria (AAB). Interspecies and intraspecies genotypic diversities of the micro‐organisms were screened respectively by amplification of the ITS1‐5.8S rDNA‐ITS2/16S‐23S rDNA ITS regions and repetitive sequence‐based PCR (rep‐PCR). The physicochemical characteristics of samples were: pH: 3·48–4·12, titratable acidity: 1·67–3·50 mg KOH g^?1, acetic acid: 0·16–0·37%, alcohol content: 0·30–2·73%, sugars (degrees Brix): 2·70–8·50. Yeast included mainly Saccharomyces cerevisiae and species of the genera Arthroascus, Issatchenkia, Candida, Trichosporon, Hanseniaspora, Kodamaea, Schizosaccharomyces, Trigonopsis and Galactomyces. Lactobacillus plantarum was the predominant LAB species. Three other species of Lactobacillus were also identified as well as isolates of Leuconostoc mesenteroides, Fructobacillus durionis and Streptococcus mitis. Acetic acid bacteria included nine species of the genus Acetobacter with Acetobacter indonesiensis as predominant species. In addition, isolates of Gluconobacter oxydans and Gluconacetobacter saccharivorans were also identified. Intraspecies diversity was observed for some species of micro‐organisms including four genotypes for Acet. indonesiensis, three for Candida tropicalis and Lactobacillus fermentum and two each for S. cerevisiae, Trichosporon asahii, Candida pararugosa and Acetobacter tropicalis.

Conclusion

fermentation of palm sap from B. akeassii involved multi‐yeast‐LAB‐AAB cultures at genus, species and intraspecies level.

Significance and Impact of the Study

First study describing microbiological and physicochemical characteristics of palm wine from B. akeassii. Genotypic diversity of palm wine LAB and AAB not reported before is demonstrated and this constitutes valuable information for better understanding of the fermentation which can be used to improve the product quality and develop added value by‐products.     

Characterization and flocculation properties of biopolymeric flocculant (glycosaminoglycan) produced by Cellulomonas sp. Okoh

Aims

Bioflocculant production potential of an actinobacteria isolated from a freshwater environment was evaluated and the bioflocculant characterized.

Methods and Results

16S rDNA nucleotide sequence and BLAST analysis was used to identify the actinobacteria and fermentation conditions, and nutritional requirements were evaluated for optimal bioflocculant production. Chemical analyses, FTIR, 1H NMR spectrometry and SEM imaging of the purified bioflocculant were carried out. The 16S rDNA nucleotide sequences showed 93% similarities to three Cellulomonas species (strain 794, Cellulomonas flavigena DSM 20109 and Cellulomonas flavigena NCIMB 8073), and the sequences was deposited in GenBank as Cellulomonas sp. Okoh (accession number HQ537132 ). Bioflocculant was optimally produced at an initial pH 7, incubation temperature 30°C, agitation speed of 160 rpm and an inoculum size of 2% (vol/vol) of cell density 1·5 × 10⁸ cfu ml^?1. Glucose (88·09% flocculating activity; yield: 4·04 ± 0·33 g l^?1), (NH₄)₂NO₃ (82·74% flocculating activity; yield: 4·47 ± 0·55 g l^?1) and MgCl₂ (90·40% flocculating activity; yield: 4·41 g l^?1) were the preferred nutritional source. Bioflocculant chemical analyses showed carbohydrate, protein and uronic acids in the proportion of 28·9, 19·3 and 18·7% in CPB and 31·4, 18·7 and 32·1% in PPB, respectively. FTIR and 1H NMR indicated the presence of carboxyl, hydroxyl and amino groups amongst others typical of glycosaminoglycan. SEM imaging revealed horizontal pleats of membranous sheets closely packed.

Conclusion

Cellulomonas sp. produces bioflocculant predominantly composed of glycosaminoglycan polysaccharides with high flocculation activity.

Significance and Impact of the Study

High flocculation activity suggests suitability for industrial applications; hence, it may serve to replace the hazardous flocculant used in water treatment.     

Repeated reunions and splits feature the highly dynamic evolution of 5S and 35S ribosomal RNA genes (rDNA) in the Asteraceae family

Background  

In flowering plants and animals the most common ribosomal RNA genes (rDNA) organisation is that in which 35S (encoding 18S-5.8S-26S rRNA) and 5S genes are physically separated occupying different chromosomal loci. However, recent observations established that both genes have been unified to a single 35S-5S unit in the genus Artemisia (Asteraceae), a genomic arrangement typical of primitive eukaryotes such as yeast, among others. Here we aim to reveal the origin, distribution and mechanisms leading to the linked organisation of rDNA in the Asteraceae by analysing unit structure (PCR, Southern blot, sequencing), gene copy number (quantitative PCR) and chromosomal position (FISH) of 5S and 35S rRNA genes in ~200 species representing the family diversity and other closely related groups.     

Comparing the usefulness of distance, monophyly and character-based DNA barcoding methods in species identification: a case study of neogastropoda

Background

DNA barcoding has recently been proposed as a promising tool for the rapid species identification in a wide range of animal taxa. Two broad methods (distance and monophyly-based methods) have been used. One method is based on degree of DNA sequence variation within and between species while another method requires the recovery of species as discrete clades (monophyly) on a phylogenetic tree. Nevertheless, some issues complicate the use of both methods. A recently applied new technique, the character-based DNA barcode method, however, characterizes species through a unique combination of diagnostic characters.

Methodology/Principal Findings

Here we analyzed 108 COI and 102 16S rDNA sequences of 40 species of Neogastropoda from a wide phylogenetic range to assess the performance of distance, monophyly and character-based methods of DNA barcoding. The distance-based method for both COI and 16S rDNA genes performed poorly in terms of species identification. Obvious overlap between intraspecific and interspecific divergences for both genes was found. The “10× rule” threshold resulted in lumping about half of distinct species for both genes. The neighbour-joining phylogenetic tree of COI could distinguish all species studied. However, the 16S rDNA tree could not distinguish some closely related species. In contrast, the character-based barcode method for both genes successfully identified 100% of the neogastropod species included, and performed well in discriminating neogastropod genera.

Conclusions/Significance

This present study demonstrates the effectiveness of the character-based barcoding method for species identification in different taxonomic levels, especially for discriminating the closely related species. While distance and monophyly-based methods commonly use COI as the ideal gene for barcoding, the character-based approach can perform well for species identification using relatively conserved gene markers (e.g., 16S rDNA in this study). Nevertheless, distance and monophyly-based methods, especially the monophyly-based method, can still be used to flag species.     

Comparison of genomes of eight species of sections <Emphasis Type="Italic">Linum</Emphasis> and <Emphasis Type="Italic">Adenolinum</Emphasis> from the genus<Emphasis Type="Italic"> Linum</Emphasis> based on chromosome banding,molecular markers and RAPD analysis

Karyotypes of species sects. Linum and Adenolinum have been studied using C/DAPI-banding, Ag-NOR staining, FISH with 5S and 26S rDNA and RAPD analysis. C/DAPI-banding patterns enabled identification of all homologous chromosome pairs in the studied karyotypes. The revealed high similarity between species L. grandiflorum (2n = 16) and L. decumbens by chromosome and molecular markers proved their close genome relationship and identified the chromosome number in L. decumbens as 2n = 16. The similarity found for C/DAPI-banding patterns between species with the same chromosome numbers corresponds with the results obtained by RAPD-analysis, showing clusterization of 16-, 18- and 30-chromosome species into three separate groups. 5S rDNA and 26S rDNA were co-localized in NOR-chromosome 1 in the genomes of all species investigated. In 30-chromosome species, there were three separate 5S rDNA sites in chromosomes 3, 8 and 13. In 16-chromosome species, a separate 5S rDNA site was also located in chromosome 3, whereas in 18-chromosome species it was found in the long arm of NOR-chromosome 1. Thus, the difference in localization of rDNA sites in species with 2n = 16, 2n = 30 and 2n = 18 confirms taxonomists opinion, who attributed these species to different sects. Linum and Adenolinum, respectively. The obtained results suggest that species with 2n = 16, 2n = 18 and 2n = 30 originated from a 16-chromosome ancestor.     

Inferring phylogenies with incomplete data sets: a 5-gene, 567-taxon analysis of angiosperms

Background  

Phylogenetic analyses of angiosperm relationships have used only a small percentage of available sequence data, but phylogenetic data matrices often can be augmented with existing data, especially if one allows missing characters. We explore the effects on phylogenetic analyses of adding 378 matK sequences and 240 26S rDNA sequences to the complete 3-gene, 567-taxon angiosperm phylogenetic matrix of Soltis et al.     

Phylogenetic placement of the enigmatic parasite, <Emphasis Type="Italic">Polypodium hydriforme</Emphasis>, within the Phylum Cnidaria

Background  

Polypodium hydriforme is a parasite with an unusual life cycle and peculiar morphology, both of which have made its systematic position uncertain. Polypodium has traditionally been considered a cnidarian because it possesses nematocysts, the stinging structures characteristic of this phylum. However, recent molecular phylogenetic studies using 18S rDNA sequence data have challenged this interpretation, and have shown that Polypodium is a close relative to myxozoans and together they share a closer affinity to bilaterians than cnidarians. Due to the variable rates of 18S rDNA sequences, these results have been suggested to be an artifact of long-branch attraction (LBA). A recent study, using multiple protein coding markers, shows that the myxozoan Buddenbrockia, is nested within cnidarians. Polypodium was not included in this study. To further investigate the phylogenetic placement of Polypodium, we have performed phylogenetic analyses of metazoans with 18S and partial 28S rDNA sequences in a large dataset that includes Polypodium and a comprehensive sampling of cnidarian taxa.     

On the value of nuclear and mitochondrial gene sequences for reconstructing the phylogeny of vanilloid orchids (Vanilloideae, Orchidaceae)

Background and Aims

Most molecular phylogenetic studies of Orchidaceae have relied heavily on DNA sequences from the plastid genome. Nuclear and mitochondrial loci have only been superficially examined for their systematic value. Since 40% of the genera within Vanilloideae are achlorophyllous mycoheterotrophs, this is an ideal group of orchids in which to evaluate non-plastid gene sequences.

Methods

Phylogenetic reconstructions for Vanilloideae were produced using independent and combined data from the nuclear 18S, 5·8S and 26S rDNA genes and the mitochondrial atpA gene and nad1b-c intron.

Key Results

These new data indicate placements for genera such as Lecanorchis and Galeola, for which plastid gene sequences have been mostly unavailable. Nuclear and mitochondrial parsimony jackknife trees are congruent with each other and previously published trees based solely on plastid data. Because of high rates of sequence divergence among vanilloid orchids, even the short 5·8S rDNA gene provides impressive levels of resolution and support.

Conclusions

Orchid systematists are encouraged to sequence nuclear and mitochondrial gene regions along with the growing number of plastid loci available.Key words: 26S rDNA, 18S rDNA, 5·8S rDNA, atpA, nad1, orchids, plastid, Vanilla, vanilloid orchids, Vanilloideae     

Genetic diversity of soybean-nodulating rhizobia in Nepal in relation to climate and soil properties

Backgroud and aims

This study was conducted to reveal the genetic diversity of soybean-nodulating rhizobia in Nepal in relation to climate and soil properties.

Method

A total of 102 bradyrhizobial strains were isolated from the root nodules of soybeans cultivated in 12 locations in Nepal varying in climate and soil properties, and their genetic diversity was examined based on 16S rDNA, ITS regions of 16S–23S rDNA, nodC and nifH. In vitro growth properties of some representative strains were examined to elucidate their characteristic distribution in Nepal.

Results

Four species of the genus Bradyrhizobium were isolated, and B. japonicum dominated at temperate locations, while in subtropical locations, B. elkanii, B. yuanmingense, and B. liaoningense dominated at acidic, moderately acidic, and slightly alkaline soils, respectively. The relative nodule occupancies could not be fully explained by their in vitro growth properties. Similar nodC and nifH genes among the strains suggested co-evolution of these genes also in Nepal, probably through horizontal gene transfer.

Conclusions

The influence of climate and soil pH on diversity at the sub-species level was revealed. It is concluded that the highly diverse climate and soils in Nepal might be conducive for the existence of diverse soybean rhizobial strains.     

Biocontrol potential of phylloplane bacterium Ochrobactrum anthropi BMO‐111 against blister blight disease of tea

Aims

The present study was carried out to screen the phylloplane bacteria from tea for antagonism against grey blight caused by Pestalotiopsis theae and blister bight caused by Exobasidium vexans and to further evaluate the efficient isolates for disease control potential under field condition.

Methods and Results

A total of 316 morphologically different phylloplane bacteria were isolated. Among the antagonists, the isolates designated as BMO‐075, BMO‐111 and BMO‐147 exhibited maximum inhibitory activity against both the pathogens under in vitro conditions and hence were selected for further evaluation under microplot field trial. Foliar application of 36‐h‐old culture of BMO‐111 (1 × 10⁸ colony‐forming units ml^?1) significantly reduced the blister blight disease incidence than the other isolates. The culture of BMO‐111 as well as its culture filtrate effectively inhibited the mycelial growth of various fungal plant pathogens. The isolate BMO‐111 was identified as Ochrobactrum anthropi based on the morphological and 16S rDNA sequence analyses.

Conclusions

It could be concluded that the biocontrol agent O. anthropi BMO‐111 was effective against blister blight disease of tea.

Significance and Impact of the Study

Further study is required to demonstrate the mechanism of its action and formulation for the biocontrol potential against blister blight disease of tea.     

IM-TORNADO: A Tool for Comparison of 16S Reads from Paired-End Libraries

Motivation

16S rDNA hypervariable tag sequencing has become the de facto method for accessing microbial diversity. Illumina paired-end sequencing, which produces two separate reads for each DNA fragment, has become the platform of choice for this application. However, when the two reads do not overlap, existing computational pipelines analyze data from read separately and underutilize the information contained in the paired-end reads.

Results

We created a workflow known as Illinois Mayo Taxon Organization from RNA Dataset Operations (IM-TORNADO) for processing non-overlapping reads while retaining maximal information content. Using synthetic mock datasets, we show that the use of both reads produced answers with greater correlation to those from full length 16S rDNA when looking at taxonomy, phylogeny, and beta-diversity.

Availability and Implementation

IM-TORNADO is freely available at http://sourceforge.net/projects/imtornado and produces BIOM format output for cross compatibility with other pipelines such as QIIME, mothur, and phyloseq.