MSClust: A Multi-Seeds based Clustering algorithm for microbiome profiling using 16S rRNA sequence

Recent developments of next generation sequencing technologies have led to rapid accumulation of 16S rRNA sequences for microbiome profiling. One key step in data processing is to cluster short sequences into operational taxonomic units (OTUs). Although many methods have been proposed for OTU inferences, a major challenge is the balance between inference accuracy and computational efficiency, where inference accuracy is often sacrificed to accommodate the need to analyze large numbers of sequences. Inspired by the hierarchical clustering method and a modified greedy network clustering algorithm, we propose a novel multi-seeds based heuristic clustering method, named MSClust, for OTU inference. MSClust first adaptively selects multi-seeds instead of one seed for each candidate cluster, and the reads are then processed using a greedy clustering strategy. Through many numerical examples, we demonstrate that MSClust enjoys less memory usage, and better biological accuracy compared to existing heuristic clustering methods while preserving efficiency and scalability.     

Ecotypes of planktonic actinobacteria with identical 16S rRNA genes adapted to thermal niches in temperate, subtropical, and tropical freshwater habitats

Seven strains with identical 16S rRNA genes affiliated with the Luna2 cluster (Actinobacteria) were isolated from six freshwater habitats located in temperate (Austria and Australia), subtropical (People's Republic of China), and tropical (Uganda) climatic zones. The isolates had sequence differences at zero to five positions in a 2,310-nucleotide fragment of the ribosomal operon, including part of the intergenic spacer upstream of the 16S rRNA gene, the complete 16S rRNA gene, the complete 16S-23S internal transcribed spacer (ITS1), and a short part of the 23S rRNA gene. Most of the few sequence differences found were located in the internal transcribed spacer sequences. Two isolates obtained from habitats in Asia and Europe, as well as two isolates obtained from different habitats in the People's Republic of China, had identical sequences for the entire fragment sequenced. In spite of minimal sequence differences in the part of the ribosomal operon investigated, the strains exhibited significant differences in their temperature response curves (with one exception), as well as pronounced differences in their temperature optima (25.0 to 35.6 degrees C). The observed differences in temperature adaptation were generally in accordance with the thermal conditions in the habitats where the strains were isolated. Strains obtained from temperate zone habitats had the lowest temperature optima, strains from subtropical habitats had intermediate temperature optima, and a strain from a tropical habitat had the highest temperature optimum. Based on the observed temperature responses, we concluded that the strains investigated are well adapted to the thermal conditions in their home habitats. Consequently, these closely related strains represent different ecotypes adapted to different thermal niches.     

Towards a phylogeny of the clostridia based on 16S rRNA sequences

Abstract The genes hbl3, cpl1 and cpl7 coding for the pneumococcal phage lytic enzymes HBL3, CPL1 and CPL7, respectively, have been cloned into shuttle plasmids that can replicate in Streptococcus pneumoniae and Escherichia coli . All these genes were expressed in E. coli under the control of either the lyt^P promoter of the lytA gene, which codes for the major pneumococcal autolysin, or the promoter of the tetracycline-resistance gene (tet^P). In contrast, cpl1 and cpl7 genes that code for lysozymes were expressed in pneumococcus only under the control of tet^P, whereas the hbl3 gene that codes for an amidase can be expressed using either promoter. The phage lysozymes or amidase expressed in S. pneumoniae M31, a mutant deleted in the lyA gene coding for short chains, were placed under physiological control since these transformed bacteria grew as normal 'diplo' cells during the exponential phase and underwent autolysis only after long incubation at 37°C. The lysis genes appear to be expressed constitutively in the transformed pneumococci, since sharply defined lysis of these cultures could be induced prematurely during the exponential phase of growth by addition of sodium deoxycholate.     

Phylogenetic relationships of nonaxenic filamentous cyanobacterial strains based on 16S rRNA sequence analysis

In order to determine the nearly complete 16S rRNA gene sequences of cyanobacteria originating from nonaxenic cultures, a cyanobacterium-specific oligonucleotide probe was developed to distinguish polymerase chain reaction (PCR) products of the cyanobacterial rRNA operons from those resulting from amplification of contaminating bacteria. Using this screening method the 16S rRNA genes of four nonaxenic filamentous cyanobacterial strains belonging to the generaLeptolyngbya andOscillatoria were cloned and sequenced. For the genusLeptolyngbya, the 16S rRNA sequence of the axenic strain PCC 73110 was also determined. Phylogenetic trees were constructed based on complete and partial sequences. The results show that the strainsLeptolyngbya foveolarum Komárek 1964/112,Leptolyngbya sp. VRUC 135 Albertano 1985/1, andLeptolyngbya boryanum PCC 73110 belong to the same cluster. StrainOscillatoria cf.corallinae SAG 8.92, which contains the rare photosynthetic pigment CU-phycoerythrin, is not closely related to other CU-phycoerythrin-containing cyanobacteria.Oscillatoria agardhii CYA 18, which is a representative of planktonicOscillatoria species that form toxic blooms in Norwegian inland waters, has no close relatives in the tree.     

Cyanoseq: A database of cyanobacterial 16S rRNA gene sequences with curated taxonomy

Cyanobacteria are photosynthetic bacteria that occupy various habitats across the globe, playing critical roles in many of Earth's biogeochemical cycles both in both aquatic and terrestrial systems. Despite their well-known significance, their taxonomy remains problematic and is the subject of much research. Taxonomic issues of Cyanobacteria have consequently led to inaccurate curation within known reference databases, ultimately leading to problematic taxonomic assignment during diversity studies. Recent advances in sequencing technologies have increased our ability to characterize and understand microbial communities, leading to the generation of thousands of sequences that require taxonomic assignment. We herein propose CyanoSeq ( https://zenodo.org/record/7569105 ), a database of cyanobacterial 16S rRNA gene sequences with curated taxonomy. The taxonomy of CyanoSeq is based on the current state of cyanobacterial taxonomy, with ranks from the domain to genus level. Files are provided for use with common naive Bayes taxonomic classifiers, such as those included in DADA2 or the QIIME2 platform. Additionally, FASTA files are provided for creation of de novo phylogenetic trees with (near) full-length 16S rRNA gene sequences to determine the phylogenetic relationship of cyanobacterial strains and/or ASV/OTUs. The database currently consists of 5410 cyanobacterial 16S rRNA gene sequences along with 123 Chloroplast, Bacterial, and Vampirovibrionia (formally Melainabacteria) sequences.     

Comparing patterns of spatial autocorrelation of assemblages of benthic invertebrates in upland rivers in south-eastern Australia

Patterns of spatial autocorrelation of biota may reveal much about underlying ecological and biological forces responsible for generating the patterns. Operationally, ecological work and many applied problems (e.g., impact detection, ecosystem health assessment using reference sites) require statistical knowledge of autocorrelation patterns. Here, we report on assemblage-level autocorrelation in the benthic-invertebrate assemblages of riffles in two adjacent, relatively pristine rivers in south-eastern Victoria, Australia (40 km reaches of the Wellington and Wonnangatta rivers). The assemblages of the Wellington River were strongly autocorrelated, but those of the Wonnangatta River showed a distance-independent pattern. There was no effect of taxonomic resolution, rarity protocols or whole-assemblage surrogates on the inferred levels of autocorrelation. We conclude that there is little evidence that one can assume the pattern of spatial relationships among invertebrate faunas within a river, and this probably holds true for the usual set of taxonomic resolutions and subsets used to discern changes wrought by human impacts.     

Distribution patterns and phylogeny of marine stramenopiles in the north pacific ocean

Marine stramenopiles (MASTs) are a diverse suite of eukaryotic microbes found in marine environments. Several MAST lineages are thought to contain heterotrophic nanoflagellates. However, MASTs remain uncultured and data on distributions and trophic modes are limited. We investigated MASTs in provinces on the west and east sides of the North Pacific Subtropical Gyre, specifically the East China Sea (ECS) and the California Current system (CALC). For each province, DNA was sampled from three zones: coastal, mesotrophic transitional, and more oligotrophic euphotic waters. Along with diatoms, chrysophytes, and other stramenopiles, sequences were recovered from nine MAST lineages in the six ECS and four CALC 18S rRNA gene clone libraries. All but one of these libraries were from surface samples. MAST clusters 1, 3, 7, 8, and 11 were identified in both provinces, with MAST cluster 3 (MAST-3) being found the most frequently. Additionally, MAST-2 was detected in the ECS and MAST-4, -9, and -12 were detected in the CALC. Phylogenetic analysis indicated that some subclades within these lineages differ along latitudinal gradients. MAST-1A, -1B, and -1C and MAST-4 size and abundance estimates obtained using fluorescence in situ hybridization on 79 spring and summer ECS samples showed a negative correlation between size of MAST-1B and MAST-4 cells and temperature. MAST-1A was rarely detected, but MAST-1B and -1C and MAST-4 were abundant in summer and MAST-1C and MAST-4 were more so at the coast, with maximum abundances of 543 and 1,896 cells ml(-1), respectively. MAST-4 and Synechococcus abundances were correlated, and experimental work showed that MAST-4 ingests Synechococcus. Together with previous studies, this study helps refine hypotheses on distribution and trophic modes of MAST lineages.     

Communities of green sulfur bacteria in marine and saline habitats analyzed by gene sequences of 16S rRNA and Fenna-Matthews-Olson protein.

Communities of green sulfur bacteria were studied in selected marine and saline habitats on the basis of gene sequences of 16S rRNA and the Fenna- Matthews-Olson (FMO) protein. The availability of group-specific primers for both 16S rDNA and the fmoA gene, which is unique to green sulfur bacteria, has, for the first time, made it possible to analyze environmental communities of these bacteria by culture-independent methods using two independent genetic markers. Sequence results obtained with fmoA genes and with 16S rDNA were largely congruent to each other. All of the 16S rDNA and fmoA sequences from habitats of the Baltic Sea, the Mediterranean Sea, Sippewissett Salt Marsh (Massachusetts, USA), and Bad Water (Death Valley, California, USA) were found within salt-dependent phylogenetic lines of green sulfur bacteria established by pure culture studies. This strongly supports the existence of phylogenetic lineages of green sulfur bacteria specifically adapted to marine and saline environments and the exclusive occurrence of these bacteria in marine and saline habitats. The great majority of clone sequences belonged to different clusters of the Prosthecochloris genus and probably represent different species. Evidence for the occurrence of two new species of Prosthecochloris was also obtained. Different habitats were dominated by representatives from the Prosthecochloris group and different clusters or species of this genus were found either exclusively or as the clearly dominant green sulfur bacterium at different habitats.     

Identification of Vibrio anguillarum in fish by using partial 16S rRNA sequences and a specific 16S rRNA oligonucleotide probe

16S rRNA from seven different Vibrio anguillarum strains was partially sequenced and compared. From this sequence information we could design a 25-base-long oligonucleotide and use it as a specific probe for identification of V. anguillarum. This was determined by RNA-DNA colony hybridization and slot-blot hybridization. Strong, specific hybridization to the probe was observed for all V. anguillarum strains tested. Furthermore, no cross-hybridization could be seen against five other bacterial species. The detection limit was 5 x 10(3) bacteria per ml. It was even possible to detect V. anguillarum, by slot-blot hybridization, directly in a homogenized kidney from a fish that had died of vibriosis. The partial sequence information revealed small but significant differences between strains of the same species. These sequence differences are sufficiently significant to allow serotyping on the RNA level. Comparing strains of different serotypes revealed a 10-base and an 11-base difference in V. anguillarum serotypes O8 and O9, respectively, in a 122-base partial sequence.     

Identification of Vibrio anguillarum in fish by using partial 16S rRNA sequences and a specific 16S rRNA oligonucleotide probe.

16S rRNA from seven different Vibrio anguillarum strains was partially sequenced and compared. From this sequence information we could design a 25-base-long oligonucleotide and use it as a specific probe for identification of V. anguillarum. This was determined by RNA-DNA colony hybridization and slot-blot hybridization. Strong, specific hybridization to the probe was observed for all V. anguillarum strains tested. Furthermore, no cross-hybridization could be seen against five other bacterial species. The detection limit was 5 x 10(3) bacteria per ml. It was even possible to detect V. anguillarum, by slot-blot hybridization, directly in a homogenized kidney from a fish that had died of vibriosis. The partial sequence information revealed small but significant differences between strains of the same species. These sequence differences are sufficiently significant to allow serotyping on the RNA level. Comparing strains of different serotypes revealed a 10-base and an 11-base difference in V. anguillarum serotypes O8 and O9, respectively, in a 122-base partial sequence.     

16S rRNA gene phylogeny and tfdA gene analysis of 2,4-D-degrading bacteria isolated in China

Twenty-two 2,4-dichlorophenoxyacetic acid (2,4-D)-degrading bacterial isolates were collected from agricultural soils at three sites in China. Sequence analysis of the 16S rRNA genes indicated that the isolates were phylogenetically grouped into four categories: Ochrobactrum anthropi, in the Alpha- class of the phylum Proteobacteria (3 out of 22 isolates), Cupriavidus sp., of the Betaproteobacteria (3 out of 22), Pseudomonas sp. and Stenotrophomonas sp., which are Gammaproteobacteria (7 out of 22), and Bacillus sp., of the phylum Firmicutes (9 out of 22). Primers were designed to amplify the conserved domain of tfdA, which is known to be involved in the degradation of 2,4-D. Results showed that the tfdA genes of all 22 strains were most similar to that of Cupriavidus necator JMP134, which belongs to the 2,4-D/α-ketoglutarate dioxygenase TfdA protein family, indicating that the JMP134-type tfdA gene is likely to be almost universal among the 2,4-D-degrading bacteria isolated from China. Degradation abilities of these 22 strains were investigated in assays using 2,4-D as the sole source of carbon and energy. Thirteen strains degraded >60 % of the available 2,4-D (500 mg l^?1) over a 1-week incubation period, while a further nine Bacillus sp. strains degraded 50–81 % of the available 2,4-D. None of these nine strains degraded other selected herbicides, such as mecoprop, 2-methyl-4-chlorophenoxyacetic acid, quizalofop, and fluroxypyr. This is the first report of 2,4-D-degradation by Bacilli.     

Identification of Comamonas species using 16S rRNA gene sequence

A bacterial strain Bz02 was isolated from a water sample collected from river Gomti at the Indian city of Lucknow. We characterized the strain using 16S rRNA sequence. Phylogenetic analysis showed that the strain formed a monophyletic clade with members of the genus Comamonas. The closest phylogenetic relative was Comamonas testosteroni with 95% 16S rRNA gene sequence similarity. It is proposed that the identified strain Bz02 be assigned as the type strain of a species of the genus Comamonas (Comamonas sp Bz02) based on 16S rRNA gene sequence search in Ribosomal Database Project, small subunit rRNA and large subunit rRNA databases together with the phylogenetic tree analysis. The sequence is deposted in GenBank with the accession number {"type":"entrez-nucleotide","attrs":{"text":"FJ211417","term_id":"209164434"}}FJ211417.     

重庆市26个南亚果实蝇种群mtDNA 16S rRNA基因部分序列及其系统进化

对重庆市26个南亚果实蝇Bactrocera(Zeugodacus)tau(Walker)种群线粒体16S rRNA基因进行测序,获得长约350bp片段的序列。对获得的序列分析表明,A,T,C,G平均含量分别为35·0%,41·3%,7·2%,16·5%,其中保守位点数342个,变异位点数5个,简约信息位点2个,自裔位点2个,所有碱基转换总数为136,替换总数为50。利用MEGA2·1软件重建系统发生树,发现其中21个南亚果实蝇种群未出现分化,另外有5个南亚果实蝇种群出现了分化,但遗传分化程度小。     

The phylogeny of Mycoplasma bovis as determined by sequence analysis of the 16S rRNA gene

Abstract The nucleotide sequence of the 16S rRNA gene of Mycoplasma bovis has been determined. Comparisons with other 16S rRNA sequences of mycoplasmas showed that Mycoplasma agalactiae is phylogenetically the closet relative. In total, only eight nucleotides differed between the M. bovis and M. agalactiae 16S rRNA sequences. The phylogenetic position of M. bovis with respect to other mycoplasmas was determined by sequence comparisons and from features in the secondary structure of 16S rRNA.     

中国地鼠线粒体DNA 16S rRNA基因序列分析及分子系统发育研究

目的通过克隆分析中国地鼠16S基因的部分序列,对中国地鼠16S基因的结构和功能进行初步探索和揭示。方法从GenBank中已报道的啮齿动物16S基因保守区设计一对引物,进行PCR扩增,测序。用Blastn与GenBank中七种啮齿类动物的16S基因进行序列比较,分析其碱基组成及变异情况,并用邻接法(NJ)、非加权组平均法(UPGMA)构建分子系统树,在分子水平上探讨中国地鼠和其他啮齿类动物的进化关系,对中国地鼠的种属地位进行了进一步验证。结果获得了中国地鼠线粒体16S基因的部分序列,其碱基组成A、T、C、G分别为40.5%、24.5%、18.7%、16.3%,与其他七种啮齿类动物的碱基含量相比,各碱基含量基本相似。NJ进化树表明,中国地鼠、金黄地鼠与欧洲仓鼠先聚为一支,小鼠与大鼠先聚为一支,东方田鼠、台湾田鼠与东欧田鼠先聚为一支。结论中国地鼠和金黄地鼠的亲缘关系最近,与传统的分类地位基本吻合。     

Comparative phylogeny of the ammonia monooxygenase subunit A and 16S rRNA genes of ammonia-oxidizing bacteria

A fragment of the ammonia monooxygenase gene (amoA) from 31 strains of ammonia-oxidizing bacteria (AOB) was sequenced and analysed phylogenetically. The results were compared with the phylogeny of 16S rDNA from AOB. For most groups of AOB we found a high consistency between the phylogenetic trees based on the 16S rDNA and amoA sequences. Although it is not a phylogenetic marker, using the amoA as a probe when studying microbial diversity will probably reduce the amount of non-AOB detected, compared to using rDNA based probes. The data presented in this paper extend and improve the basis for application of amoA in studies of AOB in the environment.     

Benchmarking of 16S rRNA gene databases using known strain sequences

16S rRNA gene analysis is the most convenient and robust method for microbiome studies. Inaccurate taxonomic assignment of bacterial strains could have deleterious effects as all downstream analyses rely heavily on the accurate assessment of microbial taxonomy. The use of mock communities to check the reliability of the results has been suggested. However, often the mock communities used in most of the studies represent only a small fraction of taxa and are used mostly as validation of sequencing run to estimate sequencing artifacts. Moreover, a large number of databases and tools available for classification and taxonomic assignment of the 16S rRNA gene make it challenging to select the best-suited method for a particular dataset. In the present study, we used authentic and validly published 16S rRNA gene type strain sequences (full length, V3-V4 region) and analyzed them using a widely used QIIME pipeline along with different parameters of OTU clustering and QIIME compatible databases. Data Analysis Measures (DAM) revealed a high discrepancy in ratifying the taxonomy at different taxonomic hierarchies. Beta diversity analysis showed clear segregation of different DAMs. Limited differences were observed in reference data set analysis using partial (V3-V4) and full-length 16S rRNA gene sequences, which signify the reliability of partial 16S rRNA gene sequences in microbiome studies. Our analysis also highlights common discrepancies observed at various taxonomic levels using various methods and databases.     

Molecular phylogeny based on the 16S rRNA gene of elite rhizobial strains used in Brazilian commercial inoculants

Nitrogen is often a limiting nutrient, therefore the sustainability of food crops, forages and green manure legumes is mainly associated with their ability to establish symbiotic associations with stem and root-nodulating N2-fixing rhizobia. The selection, identification and maintenance of elite strains for each host are critical. Decades of research in Brazil resulted in a list of strains officially recommended for several legumes, but their genetic diversity is poorly known. This study aimed at gaining a better understanding of phylogenetic relationships of 68 rhizobial strains recommended for 64 legumes, based on the sequencing of the 16S rRNA genes. The strains were isolated from a wide range of legumes, including all three subfamilies and 17 tribes. Nine main phylogenetic branches were defined, seven of them related to the rhizobial species: Bradyrhizobium japonicum, B. elkanii, Rhizobium tropici, R. leguminosarum, Sinorhizobium meliloti/S. fredii, Mesorhizobium ciceri/M. loti, and Azorhizobium caulinodans. However, some strains differed by up to 35 nucleotides from the type strains, which suggests that they may represent new species. Two other clusters included bacteria showing similarity with the genera Methylobacterium and Burkholderia, and amplification with primers for nifH and/or nodC regions was achieved with these strains. Host specificity of several strains was very low, as they were capable of nodulating legumes of different tribes and subfamilies. Furthermore, host specificity was not related to 16S rRNA, therefore evolution of ribosomal and symbiotic genes may have been diverse. Finally, the great diversity observed in this study emphasizes that tropics are an important reservoir of N2-fixation genes.     

Direct detection of 16S rRNA in soil extracts by using oligonucleotide microarrays

We report on the development and validation of a simple microarray method for the direct detection of intact 16S rRNA from unpurified soil extracts. Total RNAs from Geobacter chapellei and Desulfovibrio desulfuricans were hybridized to an oligonucleotide array consisting of universal and species-specific 16S rRNA probes. PCR-amplified products from Geobacter and Desulfovibrio were easily and specifically detected under a range of hybridization times, temperatures, and buffers. However, reproducible, specific hybridization and detection of intact rRNA could be accomplished only by using a chaperone-detector probe strategy. With this knowledge, assay conditions were developed for rRNA detection using a 2-h hybridization time at room temperature. Hybridization specificity and signal intensity were enhanced using fragmented RNA. Formamide was required in the hybridization buffer in order to achieve species-specific detection of intact rRNA. With the chaperone detection strategy, we were able to specifically hybridize and detect G. chapellei 16S rRNA directly from a total-RNA soil extract, without further purification or removal of soluble soil constituents. The detection sensitivity for G. chapellei 16S rRNA in soil extracts was at least 0.5 microg of total RNA, representing approximately 7.5 x 10(6) Geobacter cell equivalents of RNA. These results suggest that it is now possible to apply microarray technology to the direct detection of microorganisms in environmental samples, without using PCR.