Revisiting bacterial phylogeny: Natural and experimental evidence for horizontal gene transfer of 16S rRNA

Harnessing the intragenomic variability of rRNA operons to improve differentiation of Vibrio species

Although the 16S rRNA gene is frequently used as a phylogenetic marker in analysis of environmental DNA, this marker often fails to distinguish closely related species, including those in the genus Vibrio. Here, we investigate whether inclusion and analysis of 23S rRNA sequence can help overcome the intrinsic weaknesses of 16S rRNA analyses for the differentiation of Vibrio species. We construct a maximum likelihood 16S rRNA gene tree to assess the use of this gene to identify clades of Vibrio species. Within the 16S rRNA tree, we identify the putative informative bases responsible for polyphyly, and demonstrate the association of these positions with tree topology. We demonstrate that concatenation of 16S and 23S rRNA genes increases the number of informative nucleotide positions, thereby overcoming ambiguities in 16S rRNA-based phylogenetic reconstructions. Finally, we experimentally demonstrate that this approach considerably improves the differentiation and identification of Vibrio species in environmental samples.

Isolation and Genomics of Futiania mangrovii gen. nov., sp. nov., a Rare and Metabolically Versatile Member in the Class Alphaproteobacteria

Mangrove microorganisms are a major part of the coastal ecosystem and are directly associated with nutrient cycling. Despite their ecological significance, the collection of culturable mangrove microbes is limited due to difficulties in isolation and cultivation. Here, we report the isolation and genome sequence of strain FT118^T, the first cultured representative of a previously uncultivated order UBA8317 within Alphaproteobacteria, based on the combined results of 16S rRNA gene similarity, phylogenomic, and average amino acid identity analyses. We propose Futianiales ord. nov. and Futianiaceae fam. nov. with Futiania as the type genus, and FT118^T represents the type species with the name Futiania mangrovii gen. nov, sp. nov. The 16S rRNA gene sequence comparison reveals that this novel order is a rare member but has a ubiquitous distribution across various habitats worldwide, which is corroborated by the experimental confirmation that this isolate can physiologically adapt to a wide range of oxygen levels, temperatures, pH and salinity levels. Biochemical characterization, genomic annotation, and metatranscriptomic analysis of FT118^T demonstrate that it is metabolically versatile and active in situ. Genomic analysis reveals adaptive features of Futianiales to fluctuating mangrove environments, including the presence of high- and low-affinity terminal oxidases, N-type ATPase, and the genomic capability of producing various compatible solutes and polyhydroxybutyrate, which possibly allow for the persistence of this novel order across various habitats. Collectively, these results expand the current culture collection of mangrove microorganisms, providing genomic insights of how this novel taxon adapts to fluctuating environments and the culture reference to unravel possible microbe-environment interactions. IMPORTANCE The rare biosphere constitutes an essential part of the microbial community and may drive nutrient cycling and other geochemical processes. However, the difficulty in microbial isolation and cultivation has hampered our understanding of the physiology and ecology of uncultured rare lineages. In this study, we successfully isolated a novel alphaproteobacterium, designated as FT118^T, and performed a combination of phenotypic, phylogenetic, and phylogenomic analyses, confirming that this isolate represents the first cultured member of a previously uncultivated order UBA8317 within Alphaproteobacteria. It is a rare species with a ubiquitous distribution across different habitats. Genomic and metatranscriptomic analyses demonstrate that it is metabolically versatile and active in situ, suggesting its potential role in nutrient cycling despite being scarce. This work not only expands the current phylogeny of isolated Alphaproteobacteria but also provides genomic and culture reference to unravel microbial adaptation strategies in mangrove sediments and possible microbe-environment interactions.

Best practice for wildlife gut microbiome research: A comprehensive review of methodology for 16S rRNA gene investigations

Extensive research in well-studied animal models underscores the importance of commensal gastrointestinal (gut) microbes to animal physiology. Gut microbes have been shown to impact dietary digestion, mediate infection, and even modify behavior and cognition. Given the large physiological and pathophysiological contribution microbes provide their host, it is reasonable to assume that the vertebrate gut microbiome may also impact the fitness, health and ecology of wildlife. In accordance with this expectation, an increasing number of investigations have considered the role of the gut microbiome in wildlife ecology, health, and conservation. To help promote the development of this nascent field, we need to dissolve the technical barriers prohibitive to performing wildlife microbiome research. The present review discusses the 16S rRNA gene microbiome research landscape, clarifying best practices in microbiome data generation and analysis, with particular emphasis on unique situations that arise during wildlife investigations. Special consideration is given to topics relevant for microbiome wildlife research from sample collection to molecular techniques for data generation, to data analysis strategies. Our hope is that this article not only calls for greater integration of microbiome analyses into wildlife ecology and health studies but provides researchers with the technical framework needed to successfully conduct such investigations.

Intact ribosomal DNA arrays of Potentilla origin detected in Erythronium nucleus suggest recent eudicot-to-monocot horizontal transfer

During our initial phylogenetic study of the monocot genus Erythronium (Liliaceae), we observed peculiar eudicot-type internal transcribed spacer (ITS) sequences in a dataset derived from genomic DNA of Erythronium dens-canis. This raised the possibility of horizontal transfer of a eudicot alien ribosomal DNA (rDNA) into the Erythronium genome. In this work we aimed to support this hypothesis by carrying out genomic, molecular, and cytogenetic analyses. Genome skimming coupled by PacBio HiFi sequencing of a bacterial artificial chromosome clone derived from flow-sorted nuclei was used to characterise the alien 45S rDNA. Integration of alien rDNA in the recipient genome was further proved by Southern blotting and fluorescence in situ hybridization using specific probes. Alien rDNA, nested among Potentilla species in phylogenetic analysis, likely entered the Erythronium lineage in the common ancestor of E. dens-canis and E. caucasicum. Transferred eudicot-type rDNA preserved its tandemly arrayed feature on a single chromosome and was found to be transcribed in the monocot host, albeit much less efficiently than the native counterpart. This study adds a new example to the rarely documented nuclear-to-nuclear jumps of DNA between eudicots and monocots while holding the scientific community continually in suspense about the mode of DNA transfer.

Microbial genomic trait evolution is dominated by frequent and rare pulsed evolution

On the macroevolutionary time scale, does trait evolution proceed gradually or by rapid bursts (pulses) separated by prolonged periods of stasis or slow evolution? Although studies have shown that pulsed evolution is prevalent in animals, our knowledge about the tempo and mode of evolution across the tree of life is very limited. This long-standing debate calls for a test in bacteria and archaea, the most ancient and diverse forms of life with unique population genetic properties. Using a likelihood-based framework, we show that pulsed evolution is not only present but also prevalent and predominant in microbial genomic trait evolution. We detected two distinct types of pulsed evolution (small frequent and large rare jumps) that are predicted by the punctuated equilibrium and quantum evolution theories. Our findings suggest that major bacterial lineages could have originated in quick bursts and that pulsed evolution is a common theme across the tree of life.

Phylogenies of the 16S rRNA gene and its hypervariable regions lack concordance with core genome phylogenies

BACKGROUND

The 16S rRNA gene is used extensively in bacterial phylogenetics, in species delineation, and now widely in microbiome studies. However, the gene suffers from intragenomic heterogeneity, and reports of recombination and an unreliable phylogenetic signal are accumulating. Here, we compare core gene phylogenies to phylogenies constructed using core gene concatenations to estimate the strength of signal for the 16S rRNA gene, its hypervariable regions, and all core genes at the intra- and inter-genus levels. Specifically, we perform four intra-genus analyses (Clostridium, n = 65; Legionella, n = 47; Staphylococcus, n = 36; and Campylobacter, n = 17) and one inter-genus analysis [41 core genera of the human gut microbiome (31 families, 17 orders, and 12 classes), n = 82].

RESULTS

At both taxonomic levels, the 16S rRNA gene was recombinant and subject to horizontal gene transfer. At the intra-genus level, the gene showed one of the lowest levels of concordance with the core genome phylogeny (50.7% average). Concordance for hypervariable regions was lower still, with entropy masking providing little to no benefit. A major factor influencing concordance was SNP count, which showed a positive logarithmic association. Using this relationship, we determined that 690 ± 110 SNPs were required for 80% concordance (average 16S rRNA gene SNP count was 254). We also found a wide range in 16S-23S-5S rRNA operon copy number among genomes (1-27). At the inter-genus level, concordance for the whole 16S rRNA gene was markedly higher (73.8% - 10th out of 49 loci); however, the most concordant hypervariable regions (V4, V3-V4, and V1-V2) ranked in the third quartile (62.5 to 60.0%).

CONCLUSIONS

Ramifications of a poor phylogenetic performance for the 16S rRNA gene are far reaching. For example, in addition to incorrect species/strain delineation and phylogenetic inference, it has the potential to confound community diversity metrics if phylogenetic information is incorporated - for example, with popular approaches such as Faith's phylogenetic diversity and UniFrac. Our results highlight the problematic nature of these approaches and their use (along with entropy masking) is discouraged. Lastly, the wide range in 16S rRNA gene copy number among genomes also has a strong potential to confound diversity metrics. Video Abstract.

The Antibacterial, Antitumor Activities, and Bioactive Constituents’ Identification of Alectra sessiliflora Bacterial Endophytes

Due to increased antimicrobial resistance against current drugs, new alternatives are sought. Endophytic bacteria associated with medicinal plants are recognized as valuable sources of novel secondary metabolites possessing antimicrobial, antitumor, insecticidal, and antiviral activities. In this study, five bacterial endophytes were isolated and identified from the medicinal plant, Alectra sessiliflora, and their antibacterial and antitumor activities were investigated. In addition, the crude extracts of the endophytes were analyzed using gas chromatography (GC) coupled with time-of-flight mass spectrometry (TOF-MS). The identified bacterial endophytes belong to three genera viz Lysinibacillus, Peribacillus, and Bacillus, with the latter as the dominant genus with three species. Ethyl acetate extracts from the endophytes were used for antimicrobial activity against eleven pathogenic strains through minimum inhibitory concentration (MIC). The antitumor activity against the Hela cervical, Hek 293 kidney, and A549 lung carcinoma cells was determined by the MTS [3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxy-phenyl)-2-(4-sulfophenyl)-2H-tetrazolium] assay. Lysinibacillus sp. strain AS_1 exhibited broad antibacterial activity against the pathogenic strains with MIC values ranging from 4 to 8 mg/ml, while Bacillus sp. strain AS_3 displayed MIC of 0.25 mg/ml. Crude extracts of Lysinibacillus sp. strain AS_1, Peribacillus sp. strain AS_2, and Bacillus sp. strain AS_3 showed growth inhibition of more than 90% against all the cancer cell lines at a concentration of 1,000 μg/ml. Untargeted secondary metabolite profiling of the crude extracts revealed the presence of compounds with reported biological activity, such as antimicrobial, antioxidant, anti-inflammatory, antitumor, and antidiabetic properties. This study reported for the first time, bacterial endophytes associated with A. sessiliflora with antibacterial and antitumor activities.

Relation between two evolutionary clocks reveal new insights in bacterial evolution

New insights in evolution are available thanks to next-generation sequencing technologies in recent years. However, due to the network of complex relations between species, caused by the intensive horizontal gene transfer (HGT) between different bacterial species, it is difficult to discover bacterial evolution. This difficulty leads to new research in the field of phylogeny, including the gene-based phylogeny, in contrast to sequence-based phylogeny. In previous articles, we presented evolutionary insights of Synteny Index (SI) study on a large biological dataset. We showed that the SI approach naturally clusters 1133 species into 39 cliques of closely related species. In addition, we presented a model that enables calculation of the number of translocation events between genomes based on their SI distance. Here, these two studies are combined together and lead to new insights. A principal result is the relation between two evolutionary clocks: the well-known sequence-based clock influenced by point mutations, and SI distance clock influenced by translocation events. A surprising linear relation between these two evolutionary clocks rising for closely related species across all genus. In other words, these two different clocks are ticking at the same rate inside the genus level. Conversely, a phase-transition manner discovered between these two clocks across non-closely related species. This may suggest a new genus definition based on an analytic approach, since the phase-transition occurs where each gene, on average, undergoes one translocation event. In addition, rare cases of HGT among highly conserved genes can be detected as outliers from the phase-transition pattern.

Screening and Molecular Identification of Bacteria from the Midgut of Amphimallon solstitiale Larvae Exhibiting Antagonistic Activity against Bacterial Symbionts of Entomopathogenic Nematodes

Entomopathogenic nematodes (Rhabditida: Steinernematidae and Heterorhabditidae) are a group of organisms capable of infecting larvae of insects living in soil, including representatives of the family Scarabaeidae. Their insecticidal activity is related to the presence of symbiotic bacteria Xenorhabdus spp. or Photorhabdus spp. in the alimentary tract, which are released into the insect body, leading to its death caused by bacterial toxins and septicemia. Although the antibacterial activities of symbionts of entomopathogenic nematodes have been well described, there is insufficient knowledge of the interactions between these bacteria and microorganisms that naturally inhabit the alimentary tract of insects infested by nematodes. In this study, 900 bacterial strains isolated from midgut samples of Amphimallon solstitiale larvae were tested for their antagonistic activity against the selected five Xenorhabdus and Photorhabdus species. Cross-streak tests showed significant antibacterial activity of 20 isolates. These bacteria were identified as Bacillus [Brevibacterium] frigoritolerans, Bacillus toyonensis, Bacillus wiedmannii, Chryseobacterium lathyri, Chryseobacterium sp., Citrobacter murliniae, Enterococcus malodoratus, Paenibacillus sp., Serratia marcescens and Serratia sp. Since some representatives of the intestinal microbiota of A. solstitiale are able to inhibit the growth of Xenorhabdus and Photorhrhabdus bacteria in vitro, it can be assumed that this type of bacterial interaction may occur at certain stages of insect infection by Steinernema or Heterorhabditis nematodes.

Agnostic Framework for the Classification/Identification of Organisms Based on RNA Post-Transcriptional Modifications

We propose a novel approach for building a classification/identification framework based on the full complement of RNA post-transcriptional modifications (rPTMs) expressed by an organism at basal conditions. The approach relies on advanced mass spectrometry techniques to characterize the products of exonuclease digestion of total RNA extracts. Sample profiles comprising identities and relative abundances of all detected rPTM were used to train and test the capabilities of different machine learning (ML) algorithms. Each algorithm proved capable of identifying rigorous decision rules for differentiating closely related classes and correctly assigning unlabeled samples. The ML classifiers resolved different members of the Enterobacteriaceae family, alternative Escherichia coli serotypes, a series of Saccharomyces cerevisiae knockout mutants, and primary cells of the Homo sapiens central nervous system, which shared very similar genetic backgrounds. The excellent levels of accuracy and resolving power achieved by training on a limited number of classes were successfully replicated when the number of classes was significantly increased to escalate complexity. A dendrogram generated from ML-curated data exhibited a hierarchical organization that closely resembled those afforded by established taxonomic systems. Finer clustering patterns revealed the extensive effects induced by the deletion of a single pivotal gene. This information provided a putative roadmap for exploring the roles of rPTMs in their respective regulatory networks, which will be essential to decipher the epitranscriptomics code. The ubiquitous presence of RNA in virtually all living organisms promises to enable the broadest possible range of applications, with significant implications in the diagnosis of RNA-related diseases.

Current Classification of the Bacillus pumilus Group Species, the Rubber-Pathogenic Bacteria Causing Trunk Bulges Disease in Malaysia as Assessed by MLSA and Multi rep-PCR Approaches

Bacillus pumilus is the causal agent of trunk bulges disease affecting rubber and rubberwood quality and yield production. In this study, B. pumilus and other closely related species were included in B. pumilus group, as they shared over 99.5% similarity from 16S rRNA analysis. Multilocus sequence analysis (MLSA) of five housekeeping genes and repetitive elements-based polymerase chain reaction (rep-PCR) using REP, ERIC, and BOX primers conducted to analyze the diversity and systematic relationships of 20 isolates of B. pumilus group from four rubber tree plantations in Peninsular Malaysia (Serdang, Tanah Merah, Baling, and Rawang). Multi rep-PCR results revealed the genetic profiling among the B. pumilus group isolates, while MLSA results showed 98-100% similarity across the 20 isolates of B. pumilus group species. These 20 isolates, formerly established as B. pumilus, were found not to be grouped with B. pumilus. However, being distributed within distinctive groups of the B. pumilus group comprising of two clusters, A and B. Cluster A contained of 17 isolates close to B. altitudinis, whereas Cluster B consisted of three isolates attributed to B. safensis. This is the first MLSA and rep-PCR study on B. pumilus group, which provides an in-depth understanding of the diversity of these rubber-pathogenic isolates in Malaysia.

The Antibacterial Activity of Crude Extracts of Secondary Metabolites from Bacterial Endophytes Associated with Dicoma anomala

Endophytic bacteria isolated from medicinal plants are recognized valuable sources of novel bioactive compounds with various activities such as antimicrobial, anticancer, and antiviral. In this study, eleven bacterial endophytes were isolated from surface sterilized roots and leave tissues, of medicinal plant Dicoma anomala. The bacterial endophytes were identified by sequencing the 16S rRNA gene, and belong to five genera viz Bacillus, Staphylococcus, Stenotrophomonas, Enterobacter, and Pantoea. The dominant genera were Bacillus with five strains, Staphylococcus with two strains, and Stenotrophomonas with two strains. The crude extracts of seven selected bacterial endophytes indicated antimicrobial activity against five pathogenic strains Escherichia coli (ATCC 25922), Bacillus cereus (ATCC 10876), Staphylococcus aureus (NCTC 6571), Pseudomonas aeruginosa (ATCC 27853), and Klebsiella oxytoca (ATCC 13182), with significant inhibition concentration ranging from 0.312 mg/ml to 0.625 mg/ml. Finally, based on the data analysis of the crude extracts of the endophytes, we identified bioactive secondary metabolites with reported biological activities such as antimicrobial, anti-inflammatory, and antioxidant properties with biotechnological applications in medicine, agriculture, and other industries. This study reported for the first time bacterial endophytes associated with D. anomala, with antimicrobial activity against bacterial pathogens.

16S rRNA Gene Copy Number Normalization Does Not Provide More Reliable Conclusions in Metataxonomic Surveys

Sequencing 16S rRNA gene amplicons is the gold standard to uncover the composition of prokaryotic communities. The presence of multiple copies of this gene makes the community abundance data distorted and gene copy normalization (GCN) necessary for correction. Even though GCN of 16S data provided a picture closer to the metagenome before, it should also be compared with communities of known composition due to the fact that library preparation is prone to methodological biases. Here, we process 16S rRNA gene amplicon data from eleven simple mock communities with DADA2 and estimate the impact of GCN. In all cases, the mock community composition derived from the 16S sequencing differs from those expected, and GCN fails to improve the classification for most of the analysed communities. Our approach provides empirical evidence that GCN does not improve the 16S target sequencing analyses in real scenarios. We therefore question the use of GCN for metataxonomic surveys until a more comprehensive catalogue of copy numbers becomes available.

Orthogonal translation enables heterologous ribosome engineering in E. coli

The ribosome represents a promising avenue for synthetic biology, but its complexity and essentiality have hindered significant engineering efforts. Heterologous ribosomes, comprising rRNAs and r-proteins derived from different microorganisms, may offer opportunities for novel translational functions. Such heterologous ribosomes have previously been evaluated in E. coli via complementation of a genomic ribosome deficiency, but this method fails to guide the engineering of refractory ribosomes. Here, we implement orthogonal ribosome binding site (RBS):antiRBS pairs, in which engineered ribosomes are directed to researcher-defined transcripts, to inform requirements for heterologous ribosome functionality. We discover that optimized rRNA processing and supplementation with cognate r-proteins enhances heterologous ribosome function for rRNAs derived from organisms with ≥76.1% 16S rRNA identity to E. coli. Additionally, some heterologous ribosomes undergo reduced subunit exchange with E. coli-derived subunits. Cumulatively, this work provides a general framework for heterologous ribosome engineering in living cells.

Reconstructing the Phylogeny of Corynebacteriales while Accounting for Horizontal Gene Transfer

Horizontal gene transfer is a common mechanism in Bacteria that has contributed to the genomic content of existing organisms. Traditional methods for estimating bacterial phylogeny, however, assume only vertical inheritance in the evolution of homologous genes, which may result in errors in the estimated phylogenies. We present a new method for estimating bacterial phylogeny that accounts for the presence of genes acquired by horizontal gene transfer between genomes. The method identifies and corrects putative transferred genes in gene families, before applying a gene tree-based summary method to estimate bacterial species trees. The method was applied to estimate the phylogeny of the order Corynebacteriales, which is the largest clade in the phylum Actinobacteria. We report a collection of 14 phylogenetic trees on 360 Corynebacteriales genomes. All estimated trees display each genus as a monophyletic clade. The trees also display several relationships proposed by past studies, as well as new relevant relationships between and within the main genera of Corynebacteriales: Corynebacterium, Mycobacterium, Nocardia, Rhodococcus, and Gordonia. An implementation of the method in Python is available on GitHub at https://github.com/UdeS-CoBIUS/EXECT (last accessed April 2, 2020).

Diversity, Phylogeny and Plant Growth Promotion Traits of Nodule Associated Bacteria Isolated from Lotus parviflorus

Lotus spp. are widely used as a forage to improve pastures, and inoculation with elite rhizobial strains is a common practice in many countries. However, only a few Lotus species have been studied in the context of plant-rhizobia interactions. In this study, forty highly diverse bacterial strains were isolated from root nodules of wild Lotus parviflorus plants growing in two field locations in Portugal. However, only 10% of these isolates could nodulate one or more legume hosts tested, whereas 90% were thought to be opportunistic nodule associated bacteria. Phylogenetic studies place the nodulating isolates within the Bradyrhizobium genus, which is closely related to B. canariense and other Bradyrhizobium sp. strains isolated from genistoid legumes and Ornithopus spp. Symbiotic nodC and nifH gene phylogenies were fully consistent with the taxonomic assignment and host range. The non-nodulating bacteria isolated were alpha- (Rhizobium/Agrobacterium), beta- (Massilia) and gamma-proteobacteria (Pseudomonas, Lysobacter, Luteibacter, Stenotrophomonas and Rahnella), as well as some bacteroidetes from genera Sphingobacterium and Mucilaginibacter. Some of these nodule-associated bacteria expressed plant growth promotion (PGP) traits, such as production of lytic enzymes, antagonistic activity against phytopathogens, phosphate solubilization, or siderophore production. This argues for a potential beneficial role of these L. parviflorus nodule-associated bacteria.

First-tier detection of intragenomic 16S rRNA gene variation in culturable endophytic bacteria from cacao seeds

Background

Intragenomic variability in 16S rDNA is a limiting factor for taxonomic and diversity characterization of Bacteria, and studies on its occurrence in natural/environmental populations are scarce. In this work, direct DNA amplicon sequencing coupled with frequent-cutter restriction analysis allowed detection of intragenomic 16S rDNA variation in culturable endophytic bacteria from cacao seeds in a fast and attractive manner.

Methods

Total genomic DNA from 65 bacterial strains was extracted and the 16S rDNA hyper variable V5-V9 regions were amplified for enzyme digestion and direct Sanger-type sequencing. The resulting electropherograms were visually inspected and compared to the corresponding AluI-restriction profiles, as well as to complete genome sequences in databases. Restriction analysis were employed to substitute the need of amplicon cloning and re-sequencing. A specifically improved polyacrylamide-gradient electrophoresis allowed to resolve 5-bp differences in restriction fragment sizes. Chi-square analysis on 2 × 2 contingency table tested for the independence between the 'number of AluI bands' and 'type of eletropherogram'.

Results

Two types of electropherograms were obtained: unique template, with single peaks per base (clean chromatograms), and heterogeneous template, with various levels of multiple peaks per base (mixed chromatograms). Statistics revealed significant interaction between number of restriction fragments and type of electropherogram for the same amplicons: clean or mixed ones associated to ≤5 or ≥6 bands, respectively. The mixed-template pattern combined with the AluI-restriction profiles indicated a high proportion of 49% of the culturable endophytes from a tropical environment showing evidence of intragenomic 16S rDNA heterogeneity.

Conclusion

The approach presented here was useful for a rapid, first-tier detection of intragenomic variation in culturable isolates, which can be applied in studies of other natural populations; a preliminary view of intragenomic heterogeneity levels can complement culture-dependent and -independent methods. Consequences of these findings in taxonomic and diversity studies in complex bacterial communities are discussed.

Identification, Classification and Screening for γ-Amino-butyric Acid Production in Lactic Acid Bacteria from Cambodian Fermented Foods

Screening for various types of lactic acid bacteria (LAB) that form the biological agent γ-amino-butyric acid (GABA) is important to produce different kinds of GABA-containing fermented foods. So far, no GABA-producing LAB have been reported from Cambodian fermented foods. Most small-scale fermentations and even some industrial processes in this country still rely on indigenous LAB. The application of GABA-producing autochthonous starters would allow the production of Cambodian fermented foods with an additional nutritional value that meet the population's dietary habits and that are also more attractive for the international food market. Matrix-assisted laser desorption/ionizing time-of-flight mass spectrometry (MALDI-TOF MS) and partial 16S rDNA sequencing were used to identify 68 LAB isolates from Cambodian fermented foods. These isolates were classified and grouped with (GTG)5 rep-PCR, resulting in 50 strains. Subsequently, all strains were investigated for their ability to produce GABA by thin layer chromatography. GABA-positive strains were further analyzed by the GABase assay. Of the six GABA-positive LAB strains-one Lactobacillus futsaii, two Lactobacillus namurensis, and three Lactobacillus plantarum strains-two Lactobacillus plantarum strains produced high amounts of GABA (20.34 mM, 16.47 mM). These strains should be further investigated for their potential application as GABA-producing starter cultures in the food applications.

Occurrence of randomly recombined functional 16S rRNA genes in Thermus thermophilus suggests genetic interoperability and promiscuity of bacterial 16S rRNAs

Based on the structural complexity of ribosomes, 16S rRNA genes are considered species-specific and hence used for bacterial phylogenetic analysis. However, a growing number of reports suggest the occurrence of horizontal gene transfer, raising genealogical questions. Here we show the genetic interoperability and promiscuity of 16S rRNA in the ribosomes of an extremely thermophilic bacterium, Thermus thermophilus. The gene in this thermophile was systematically replaced with a diverse array of heterologous genes, resulting in the discovery of various genes that supported growth, some of which were from different phyla. Moreover, numerous functional chimeras were spontaneously generated. Remarkably, cold-adapted mutants were obtained carrying chimeric or full-length heterologous genes, indicating that horizontal gene transfer promoted adaptive evolution. The ribosome may well be understood as a patchworked supramolecule comprising patchworked components. We here propose the “random patch model” for ribosomal evolution.

A 16S rDNA PCR-based theoretical to actual delta approach on culturable mock communities revealed severe losses of diversity information

BACKGROUND

Subunits of ribosomal RNA genes (rDNAs) characterized by PCR-based protocols have been the proxy for studies in microbial taxonomy, phylogenetics, evolution and ecology. However, relevant factors have shown to interfere in the experimental outputs in a variety of systems. In this work, a 'theoretical' to 'actual' delta approach was applied to data on culturable mock bacterial communities (MBCs) to study the levels of losses in operational taxonomic units (OTUs) detectability. Computational and lab-bench strategies based on 16S rDNA amplification by 799F and U1492R primers were employed, using a fingerprinting method with highly improved detectability of fragments as a case-study tool. MBCs were of two major types: in silico MBCs, assembled with database-retrieved sequences, and in vitro MBCs, with AluI digestions of PCR data generated from culturable endophytes isolated from cacao trees.

RESULTS

Interfering factors for the 16 s rDNA amplifications, such as the type of template, direct and nested PCR, proportion of chloroplast DNA from a tropical plant source (Virola officinalis), and biased-amplification by the primers resulted in altered bacterial 16S rDNA amplification, both on MBCs and V. officinalis leaf-extracted DNA. For the theoretical data, the maximum number of fragments for in silico and in vitro cuts were not significantly different from each other. Primers' preferences for certain sequences were detected, depending on the MBCs' composition prior to PCR. The results indicated overall losses from 2.3 up to 8.2 times in the number of OTUs detected from actual AluI digestions of MBCs when compared to in silico and in vitro theoretical data.

CONCLUSIONS

Due to all those effects, the final amplification profile of the bacterial community assembled was remarkably simplified when compared to the expected number of detectable fragments known to be present in the MBC. From these findings, the scope of hypotheses generation and conclusions from experiments based on PCR amplifications of bacterial communities was discussed.

Taxonomy annotation and guide tree errors in 16S rRNA databases

Sequencing of the 16S ribosomal RNA (rRNA) gene is widely used to survey microbial communities. Specialized 16S rRNA databases have been developed to support this approach including Greengenes, RDP and SILVA. Most taxonomy annotations in these databases are predictions from sequence rather than authoritative assignments based on studies of type strains or isolates. In this work, I investigated the taxonomy annotations and guide trees provided by these databases. Using a blinded test, I estimated that the annotation error rate of the RDP database is ∼10%. The branching orders of the Greengenes and SILVA guide trees were found to disagree at comparable rates with each other and with taxonomy annotations according to the training set (authoritative reference) provided by RDP, indicating that the trees have comparable quality. Pervasive conflicts between tree branching order and type strain taxonomies strongly suggest that the guide trees are unreliable guides to phylogeny. I found 249,490 identical sequences with conflicting annotations in SILVA v128 and Greengenes v13.5 at ranks up to phylum (7,804 conflicts), indicating that the annotation error rate in these databases is ∼17%.

Multiple Ribosomal RNA Operons in Bacteria; Their Concerted Evolution and Potential Consequences on the Rate of Evolution of Their 16S rRNA

Bacterial species differ greatly in the number and location of the rRNA operons which may be present in the bacterial chromosomes and plasmids. Most bacterial species contain more than one ribosomal RNA operon copy in their genomes, with some species containing up to 15 such copies. We review the number and location of the rRNA operons and discuss evolution of 16S rRNA (rrs) genes -which are considered as ultimate chronometers for phylogenetic classification- in bacteria with multiple copies of these genes. In these bacterial species, the rrs genes must evolve in concert and sequence changes generated by mutation or horizontal gene transfer must be either erased or spread to every gene copy to avoid divergence, as it occurs when they are present in different species. Analysis of polymorphic sites in intra-genomic rrs copies identifies putative conversion events and demonstrates that sequence conversion is patchy and occurs in small conversion tracts. Sequence conversion probably arises by a non-reciprocal transfer between two or more copies where one copy contributes only a small contiguous segment of DNA, whereas the other copy contributes the rest of the genome in a fairly well understood molecular process. Because concerted evolution implies that a mutation in any of the rrs copies is either eliminated or transferred to every rrs gene in the genome, this process should slow their evolution rate relative to that of single copy genes. However, available data on the rrs genes in bacterial genomes do not show a clear relationship between their evolution rates and the number of their copies in the genome.

Functional metagenomic approach to identify overlooked antibiotic resistance mutations in bacterial rRNA

Our knowledge as to how bacteria acquire antibiotic resistance is still fragmented, especially for the ribosome-targeting drugs. In this study, with the aim of finding novel mechanisms that render bacteria resistant to the ribosome-targeting antibiotics, we developed a general method to systematically screen for antibiotic resistant 16 S ribosomal RNAs (rRNAs), which are the major target for multiple antibiotics (e.g. spectinomycin, tetracycline, and aminoglycosides), and identify point mutations therein. We used Escherichia coli ∆7, a null mutant of the rrn (ribosomal RNA) operons, as a surrogate host organism to construct a metagenomic library of 16 S rRNA genes from the natural (non-clinical) environment. The library was screened for spectinomycin resistance to obtain four resistant 16 S rRNA genes from non-E. coli bacterial species. Bioinformatic analysis and site-directed mutagenesis identified three novel mutations - U1183C (the first mutation discovered in a region other than helix 34), and C1063U and U1189C in helix 34 - as well as three well-described mutations (C1066U, C1192G, and G1193A). These results strongly suggest that uncharacterized antibiotic resistance mutations still exist, even for traditional antibiotics.

Phylogenetic Network Analysis Revealed the Occurrence of Horizontal Gene Transfer of 16S rRNA in the Genus Enterobacter

Horizontal gene transfer (HGT) is a ubiquitous genetic event in bacterial evolution, but it seldom occurs for genes involved in highly complex supramolecules (or biosystems), which consist of many gene products. The ribosome is one such supramolecule, but several bacteria harbor dissimilar and/or chimeric 16S rRNAs in their genomes, suggesting the occurrence of HGT of this gene. However, we know little about whether the genes actually experience HGT and, if so, the frequency of such a transfer. This is primarily because the methods currently employed for phylogenetic analysis (e.g., neighbor-joining, maximum likelihood, and maximum parsimony) of 16S rRNA genes assume point mutation-driven tree-shape evolution as an evolutionary model, which is intrinsically inappropriate to decipher the evolutionary history for genes driven by recombination. To address this issue, we applied a phylogenetic network analysis, which has been used previously for detection of genetic recombination in homologous alleles, to the 16S rRNA gene. We focused on the genus Enterobacter, whose phylogenetic relationships inferred by multi-locus sequence alignment analysis and 16S rRNA sequences are incompatible. All 10 complete genomic sequences were retrieved from the NCBI database, in which 71 16S rRNA genes were included. Neighbor-joining analysis demonstrated that the genes residing in the same genomes clustered, indicating the occurrence of intragenomic recombination. However, as suggested by the low bootstrap values, evolutionary relationships between the clusters were uncertain. We then applied phylogenetic network analysis to representative sequences from each cluster. We found three ancestral 16S rRNA groups; the others were likely created through recursive recombination between the ancestors and chimeric descendants. Despite the large sequence changes caused by the recombination events, the RNA secondary structures were conserved. Successive intergenomic and intragenomic recombination thus shaped the evolution of 16S rRNA genes in the genus Enterobacter.

A Robust Framework for Microbial Archaeology

Microbial archaeology is flourishing in the era of high-throughput sequencing, revealing the agents behind devastating historical plagues, identifying the cryptic movements of pathogens in prehistory, and reconstructing the ancestral microbiota of humans. Here, we introduce the fundamental concepts and theoretical framework of the discipline, then discuss applied methodologies for pathogen identification and microbiome characterization from archaeological samples. We give special attention to the process of identifying, validating, and authenticating ancient microbes using high-throughput DNA sequencing data. Finally, we outline standards and precautions to guide future research in the field.

Functional Constraints on Replacing an Essential Gene with Its Ancient and Modern Homologs

Genes encoding proteins that carry out essential informational tasks in the cell, in particular where multiple interaction partners are involved, are less likely to be transferable to a foreign organism. Here, we investigated the constraints on transfer of a gene encoding a highly conserved informational protein, translation elongation factor Tu (EF-Tu), by systematically replacing the endogenous tufA gene in the Escherichia coli genome with its extant and ancestral homologs. The extant homologs represented tuf variants from both near and distant homologous organisms. The ancestral homologs represented phylogenetically resurrected tuf sequences dating from 0.7 to 3.6 billion years ago (bya). Our results demonstrate that all of the foreign tuf genes are transferable to the E. coli genome, provided that an additional copy of the EF-Tu gene, tufB, remains present in the E. coli genome. However, when the tufB gene was removed, only the variants obtained from the gammaproteobacterial family (extant and ancestral) supported growth which demonstrates the limited functional interchangeability of E. coli tuf with its homologs. Relative bacterial fitness correlated with the evolutionary distance of the extant tuf homologs inserted into the E. coli genome. This reduced fitness was associated with reduced levels of EF-Tu and reduced rates of protein synthesis. Increasing the expression of tuf partially ameliorated these fitness costs. In summary, our analysis suggests that the functional conservation of protein activity, the amount of protein expressed, and its network connectivity act to constrain the successful transfer of this essential gene into foreign bacteria.IMPORTANCE Horizontal gene transfer (HGT) is a fundamental driving force in bacterial evolution. However, whether essential genes can be acquired by HGT and whether they can be acquired from distant organisms are very poorly understood. By systematically replacing tuf with ancestral homologs and homologs from distantly related organisms, we investigated the constraints on HGT of a highly conserved gene with multiple interaction partners. The ancestral homologs represented phylogenetically resurrected tuf sequences dating from 0.7 to 3.6 bya. Only variants obtained from the gammaproteobacterial family (extant and ancestral) supported growth, demonstrating the limited functional interchangeability of E. coli tuf with its homologs. Our analysis suggests that the functional conservation of protein activity, the amount of protein expressed, and its network connectivity act to constrain the successful transfer of this essential gene into foreign bacteria.

Toxic genes present a unique phylogenetic signature

Horizontal gene transfer (HGT) is a major part of the evolution of Archaea and Bacteria, to the extent that the validity of the Tree of Life concept for prokaryotes has been seriously questioned. The patterns and routes of HGT remain a subject of intense study and debate. It was discovered that while several genes exhibit rampant HGT across the whole prokaryotic tree of life, others are lethal to certain organisms and therefore cannot be successfully transferred to them. We distinguish between these two classes of genes and show analytically that genes found to be toxic to a specific species (E. coli) also resist HGT in general. Several tools we employ show evidence to support that claim. One of those tools is the quartet plurality distribution (QPD), a mathematical tool that measures tendency to HGT over a large set of genes and species. When aggregated over a collection of genes, it can reveal important properties of this collection. We conclude that evidence of toxicity of certain genes to a wide variety of prokaryotes are revealed using the new tool of quartet plurality distribution.

A Multilocus Sequence Analysis Scheme for Phylogeny of Thioclava Bacteria and Proposal of Two Novel Species

A multilocus sequence analysis (MLSA) was established and performed on the genus Thioclava, including 23 strains isolated from diverse marine environments, with the aim of better differentiation of strains and species within this genus. The study was based on sequences of 16S rRNA gene and five protein-coding housekeeping genes, gyrB, rpoD, dnaK, trpB, and recA. In contrast to 16S rRNA gene-based tree that was unable to separate some species within this genus, each tree based on a single housekeeping gene and MLSA had consistently defined seven clades, corresponding to the five established ones and two novel ones. The digital DNA-DNA hybridization and average nucleotide identity analyses based on genome sequences of the representative strains reconfirmed the validity of the MLSA analysis, and recommended a 97.3% MLSA similarity as the soft species threshold and nine species representing the five known and four putative novel species. Two of the four new species were identified as Thioclava sediminum sp. nov. (type strain TAW-CT134^T = MCCC 1A10143^T = LMG 29615^T) and Thioclava marinus sp. nov. (type strain 11.10-0-13^T = MCCC 1A03502^T = LMG 29618^T) by using a polyphasic taxonomic approach. Taken together, the newly established MLSA in this study first described the variability and phylogeny of the genus Thioclava which contributes to better understanding its ecology and evolution.

Marker genes that are less conserved in their sequences are useful for predicting genome-wide similarity levels between closely related prokaryotic strains

BACKGROUND

The 16s rRNA gene is so far the most widely used marker for taxonomical classification and separation of prokaryotes. Since it is universally conserved among prokaryotes, it is possible to use this gene to classify a broad range of prokaryotic organisms. At the same time, it has often been noted that the 16s rRNA gene is too conserved to separate between prokaryotes at finer taxonomic levels.

RESULTS

In this paper, we examine how well levels of similarity of 16s rRNA and 73 additional universal or nearly universal marker genes correlate with genome-wide levels of gene sequence similarity. We demonstrate that the percent identity of 16s rRNA predicts genome-wide levels of similarity very well for distantly related prokaryotes, but not for closely related ones. In closely related prokaryotes, we find that there are many other marker genes for which levels of similarity are much more predictive of genome-wide levels of gene sequence similarity. Finally, we show that the identities of the markers that are most useful for predicting genome-wide levels of similarity within closely related prokaryotic lineages vary greatly between lineages. However, the most useful markers are always those that are least conserved in their sequences within each lineage.

CONCLUSIONS

Our results show that by choosing markers that are less conserved in their sequences within a lineage of interest, it is possible to better predict genome-wide gene sequence similarity between closely related prokaryotes than is possible using the 16s rRNA gene. We point readers towards a database we have created (POGO-DB) that can be used to easily establish which markers show lowest levels of sequence conservation within different prokaryotic lineages.

Bacterial clade with the ribosomal RNA operon on a small plasmid rather than the chromosome

rRNA is essential for life because of its functional importance in protein synthesis. The rRNA (rrn) operon encoding 16S, 23S, and 5S rRNAs is located on the "main" chromosome in all bacteria documented to date and is frequently used as a marker of chromosomes. Here, our genome analysis of a plant-associated alphaproteobacterium, Aureimonas sp. AU20, indicates that this strain has its sole rrn operon on a small (9.4 kb), high-copy-number replicon. We designated this unusual replicon carrying the rrn operon on the background of an rrn-lacking chromosome (RLC) as the rrn-plasmid. Four of 12 strains close to AU20 also had this RLC/rrn-plasmid organization. Phylogenetic analysis showed that those strains having the RLC/rrn-plasmid organization represented one clade within the genus Aureimonas. Our finding introduces a previously unaddressed viewpoint into studies of genetics, genomics, and evolution in microbiology and biology in general.

Rare Events of Intragenus and Intraspecies Horizontal Transfer of the 16S rRNA Gene

Horizontal gene transfer (HGT) of operational genes has been widely reported in prokaryotic organisms. However, informational genes such as those involved in transcription and translation processes are very difficult to be horizontally transferred, as described by Woese’s complexity hypothesis. Here, we analyzed all of the completed prokaryotic genome sequences (2,143 genomes) in the NCBI (National Center for Biotechnology Information) database, scanned for genomes with high intragenomic heterogeneity of 16S rRNA gene copies, and explored potential HGT events of ribosomal RNA genes based on the phylogeny, genomic organization, and secondary structures of the ribosomal RNA genes. Our results revealed 28 genomes with relatively high intragenomic heterogeneity of multiple 16S rRNA gene copies (lowest pairwise identity <98.0%), and further analysis revealed HGT events and potential donors of the heterogeneous copies (such as HGT from Chlamydia suis to Chlamydia trachomatis) and mutation events of some heterogeneous copies (such as Streptococcus suis JS14). Interestingly, HGT of the 16S rRNA gene only occurred at intragenus or intraspecies levels, which is quite different from the HGT of operational genes. Our results improve our understanding regarding the exchange of informational genes.

Comparing K-mer based methods for improved classification of 16S sequences

Background

The need for precise and stable taxonomic classification is highly relevant in modern microbiology. Parallel to the explosion in the amount of sequence data accessible, there has also been a shift in focus for classification methods. Previously, alignment-based methods were the most applicable tools. Now, methods based on counting K-mers by sliding windows are the most interesting classification approach with respect to both speed and accuracy. Here, we present a systematic comparison on five different K-mer based classification methods for the 16S rRNA gene. The methods differ from each other both in data usage and modelling strategies. We have based our study on the commonly known and well-used naïve Bayes classifier from the RDP project, and four other methods were implemented and tested on two different data sets, on full-length sequences as well as fragments of typical read-length.

Results

The difference in classification error obtained by the methods seemed to be small, but they were stable and for both data sets tested. The Preprocessed nearest-neighbour (PLSNN) method performed best for full-length 16S rRNA sequences, significantly better than the naïve Bayes RDP method. On fragmented sequences the naïve Bayes Multinomial method performed best, significantly better than all other methods. For both data sets explored, and on both full-length and fragmented sequences, all the five methods reached an error-plateau.

Conclusions

We conclude that no K-mer based method is universally best for classifying both full-length sequences and fragments (reads). All methods approach an error plateau indicating improved training data is needed to improve classification from here. Classification errors occur most frequent for genera with few sequences present. For improving the taxonomy and testing new classification methods, the need for a better and more universal and robust training data set is crucial.

Evolution of the tRNALeu (UAA) Intron and Congruence of Genetic Markers in Lichen-Symbiotic Nostoc

The group I intron interrupting the tRNALeu UAA gene (trnL) is present in most cyanobacterial genomes as well as in the plastids of many eukaryotic algae and all green plants. In lichen symbiotic Nostoc, the P6b stem-loop of trnL intron always involves one of two different repeat motifs, either Class I or Class II, both with unresolved evolutionary histories. Here we attempt to resolve the complex evolution of the two different trnL P6b region types. Our analysis indicates that the Class II repeat motif most likely appeared first and that independent and unidirectional shifts to the Class I motif have since taken place repeatedly. In addition, we compare our results with those obtained with other genetic markers and find strong evidence of recombination in the 16S rRNA gene, a marker widely used in phylogenetic studies on Bacteria. The congruence of the different genetic markers is successfully evaluated with the recently published software Saguaro, which has not previously been utilized in comparable studies.

Pre_GI: a global map of ontological links between horizontally transferred genomic islands in bacterial and archaeal genomes

The Predicted Genomic Islands database (Pre_GI) is a comprehensive repository of prokaryotic genomic islands (islands, GIs) freely accessible at http://pregi.bi.up.ac.za/index.php. Pre_GI, Version 2015, catalogues 26 744 islands identified in 2407 bacterial/archaeal chromosomes and plasmids. It provides an easy-to-use interface which allows users the ability to query against the database with a variety of fields, parameters and associations. Pre_GI is constructed to be a web-resource for the analysis of ontological roads between islands and cartographic analysis of the global fluxes of mobile genetic elements through bacterial and archaeal taxonomic borders. Comparison of newly identified islands against Pre_GI presents an alternative avenue to identify their ontology, origin and relative time of acquisition. Pre_GI aims to aid research on horizontal transfer events and materials through providing data and tools for holistic investigation of migration of genes through ecological niches and taxonomic boundaries.

Extending the conserved phylogenetic core of archaea disentangles the evolution of the third domain of life

Initial studies of the archaeal phylogeny relied mainly on the analysis of the RNA component of the small subunit of the ribosome (SSU rRNA). The resulting phylogenies have provided interesting but partial information on the evolutionary history of the third domain of life because SSU rRNA sequences do not contain enough phylogenetic signal to resolve all nodes of the archaeal tree. Thus, many relationships, and especially the most ancient ones, remained elusive. Moreover, SSU rRNA phylogenies can be heavily biased by tree reconstruction artifacts. The sequencing of complete genomes allows using a variety of protein markers as an alternative to SSU rRNA. Taking advantage of the recent burst of archaeal complete genome sequences, we have carried out an in-depth phylogenomic analysis of this domain. We have identified 200 new protein families that, in addition to the ribosomal proteins and the subunits of the RNA polymerase, form a conserved phylogenetic core of archaeal genes. The accurate analysis of these markers combined with desaturation approaches shed new light on the evolutionary history of Archaea and reveals that several relationships recovered in recent analyses are likely the consequence of tree reconstruction artifacts. Among others, we resolve a number of important relationships, such as those among methanogens Class I, and we propose the definition of two new superclasses within the Euryarchaeota: Methanomada and Diaforarchaea.

Identification of horizontally transferred genes in the genus Colletotrichum reveals a steady tempo of bacterial to fungal gene transfer

BACKGROUND

Horizontal gene transfer (HGT) is the stable transmission of genetic material between organisms by means other than vertical inheritance. HGT has an important role in the evolution of prokaryotes but is relatively rare in eukaryotes. HGT has been shown to contribute to virulence in eukaryotic pathogens. We studied the importance of HGT in plant pathogenic fungi by identifying horizontally transferred genes in the genomes of three members of the genus Colletotrichum.

RESULTS

We identified eleven HGT events from bacteria into members of the genus Colletotrichum or their ancestors. The HGT events include genes involved in amino acid, lipid and sugar metabolism as well as lytic enzymes. Additionally, the putative minimal dates of transference were calculated using a time calibrated phylogenetic tree. This analysis reveals a constant flux of genes from bacteria to fungi throughout the evolution of subphylum Pezizomycotina.

CONCLUSIONS

Genes that are typically transferred by HGT are those that are constantly subject to gene duplication and gene loss. The functions of some of these genes suggest roles in niche adaptation and virulence. We found no evidence of a burst of HGT events coinciding with major geological events. In contrast, HGT appears to be a constant, albeit rare phenomenon in the Pezizomycotina, occurring at a steady rate during their evolution.