Variable region analysis of 16S rRNA gene sequences is the most common tool in bacterial taxonomic studies. Although used for distinguishing bacterial species, its use remains limited due to the presence of variable copy numbers with sequence variation in the genomes. In this study, 16S rRNA gene sequences, obtained from completely assembled whole genome and Sanger electrophoresis sequencing of cloned PCR products from Serratia fonticola GS2, were compared. Sanger sequencing produced a combination of sequences from multiple copies of 16S rRNA genes. To determine whether the variant copies of 16S rRNA genes affected Sanger sequencing, two ratios (5:5 and 8:2) with different concentrations of cloned 16S rRNA genes were used; it was observed that the greater the number of copies with similar sequences the higher its chance of amplification. Effect of multiple copies for taxonomic classification of 16S rRNA gene sequences was investigated using the strain GS2 as a model. 16S rRNA copies with the maximum variation had 99.42% minimum pairwise similarity and this did not have an effect on species identification. Thus, PCR products from genomes containing variable 16S rRNA gene copies can provide sufficient information for species identification except from species which have high similarity of sequences in their 16S rRNA gene copies like the case of Bacillus thuringiensis and Bacillus cereus. In silico analysis of 1,616 bacterial genomes from long-read sequencing was also done. The average minimum pairwise similarity for each phylum was reported with their average genome size and average “unique copies” of 16S rRNA genes and we found that the phyla Proteobacteria and Firmicutes showed the highest amount of variation in their copies of their 16S rRNA genes. Overall, our results shed light on how the variations in the multiple copies of the 16S rRNA genes of bacteria can aid in appropriate species identification.

The intestinal tract is one of the most sensitive organs following irradiation. The protective effect of specific indigenous microbiota on irradiation-induced damage to intestinal epithelial cells has not been reported. Mice were irradiated with a single dose of 6 Gy of gamma rays. The intestinal damage was analyzed by histopathology. Intestinal stemness and differentiation were determined by intestinal organoid culture. Microbiota community was observed by high-throughput 16S rRNA gene sequencing and oligotyping analysis. We showed that distal small intestine was damaged by sublethal dose of gamma irradiation. Intestinal organoids derived from the irradiated mice showed defects in budding and mucin expression, suggesting the detrimental effect of irradiation on the intestinal stemness and differentiation. In addition, irradiation reduced intestinal immunoglobulin A level, concomitant with decreased microbiota diversity based on our high-throughput 16S rRNA gene sequencing data. Especially, the relative abundance of Lactobacillus was reduced at early time point post-irradiation; however, it was recovered at late time point. Oligotyping analysis within the Lactobacillus genus indicated that Lactobacillus-related oligotype 1 (OT1) including Lactobacillus acidophilus might drive recovery after irradiation as it was associated with increased long-term numbers post-exposure. We showed that treatment with heat-killed L. acidophilus rescued the budding-impaired organoids and induced sufficient differentiation in epithelial cells, and particularly mucin-producing cells, in intestinal organoids. This study provides the first evidence that the indigenous gut bacteria L. acidophilus enhance intestinal epithelial function with respect to irradiation-induced intestinal damage by improving intestinal stem cell function and cell differentiation.

Biliary atresia (BA) results in severe bile blockage and is caused by the absence of extrahepatic ducts. Even after successful hepatic portoenterostomy, a considerable number of patients are likely to show progressive deterioration in liver function. Recent studies show that mutations in protein-coding mitochondrial DNA (mtDNA) genes and/or mitochondrial genes in nuclear DNA (nDNA) are associated with hepatocellular dysfunction. This observation led us to investigate whether hepatic dysfunctions in BA is genetically associated with mtDNA mutations. We sequenced the mtDNA protein-coding genes in 14 liver specimens from 14 patients with BA and 5 liver specimens from 5 patients with choledochal cyst using next-generation sequencing. We found 34 common non-synonymous variations in mtDNA protein-coding genes in all patients examined. A systematic 3D structural analysis revealed the presence of several single nucleotide polymorphism-like mutations in critical regions of complexes I to V, that are involved in subunit assembly, proton-pumping activity, and/or supercomplex formation. The parameters of chronic hepatic injury and liver dysfunction in BA patients were also significantly correlated with the extent of hepatic failure, suggesting that the mtDNA mutations may aggravate hepatopathy. Therefore, mitochondrial mutations may underlie the pathological mechanisms associated with BA.

Quorum sensing (QS) enables bacteria to organize gene expression programs, thereby coordinating collective behaviors. It involves the production, release, and population-wide detection of extracellular signaling molecules. The cellular processes regulated by QS in bacteria are diverse and may be used in mutualistic coordination or in response to changing environmental conditions. Here, we focused on the influence of the QS-dependent genes of our model bacterial strain Serratia fonticola GS2 on potential plant growth promoting (PGP) activities including indole-3-acetic acid (IAA) production, 1-aminocyclopropane-1-carboxylate (ACC) deaminase activity, and biofilm formation. Based on genomic and phenotypic experimental data we identified and investigated the function of QS genes in the genome of the model strain. Our gene deletion study confirmed the biological functionality of the QS auto-inducer (gloI) and receptor (gloR) on potential PGP activities of GS2. A transcriptomic approach was also undertaken to understand the role of QS genes in regulation of genes primarily involved in PGP activities (IAA, ACC deaminase activity, and biofilm formation). Both transcriptomic and phenotypic data revealed that the QS-deletion mutants had considerably less PGP activities, as compared to the wild type. In addition, in vivo plant experiments showed that plants treated with GS2 had significantly higher growth rates than plants treated with the QS-deletion mutants. Overall, our results showed how QS-dependent genes regulate the potential PGP activities of GS2. This information may be helpful in understanding the relationship between QS-dependent genes and the PGP activity of bacteria, which aid in the production of practical bio-fertilizers for plant growth promotion.