Longitudinal analysis of exposure to a low concentration of oxytetracycline on the zebrafish gut microbiome.

Oxytetracycline, a widely produced and administered antibiotic, is uncontrollably released in low concentrations in various types of environments. However, the impact of exposure to such low concentrations of antibiotics on the host remains poorly understood. In this study, we exposed zebrafish to a low concentration (5,000 ng/L) of oxytetracycline for 1 month, collected samples longitudinally (Baseline, and Days 3, 6, 9, 12, 24, and 30), and elucidated the impact of exposure on microbial composition, antibiotic resistance genes, mobile genetic elements, and phospholipid metabolism pathway through comparison of the sequenced data with respective sequence databases. We identified Pseudomonas aeruginosa, a well-known pathogen, to be significantly positively associated with the duration of oxytetracycline exposure (Adjusted P = 5.829e-03). Several tetracycline resistance genes (e.g., tetE) not only showed significantly higher abundance in the exposed samples but were also positively associated with the duration of exposure (Adjusted P = 1.114e-02). Furthermore, in the exposed group, the relative abundance of genes involved in phospholipid metabolism had also decreased. Lastly, we characterized the impact of exposure on zebrafish intestinal structure and found that the goblet cell counts were decreased (~82%) after exposure. Overall, our results show that a low concentration of oxytetracycline can increase the abundance of pathogenic bacteria and lower the abundance of key metabolic pathways in the zebrafish gut microbiome that can render them prone to bacterial infections and health-associated complications.

Genome-Resolved Characterization of Structure and Potential Functions of the Zebrafish Stool Microbiome.

Zebrafish have been used as a model organism for more than 50 years and are considered an excellent model for studying host-microbiome interactions. However, this largely depends on our understanding of the zebrafish gut microbiome itself. Despite advances in sequencing and data analysis methods, the zebrafish gut microbiome remains highly understudied. This study performed the de novo metagenome assembly and recovery of the metagenome-assembled genomes (MAGs) through genome binning (and refinement) of the contigs assembled from the zebrafish stool. The results indicate that majority of the MAGs had excellent quality i.e. high completeness (≥90%) and low contamination levels (≤5%). MAGs mainly belong to the taxa that are known to be members of the core zebrafish stool microbiome, including the phylum Proteobacteria, Fusobacteriota, and Actinobacteriota. However, most of the MAGs remained unclassified at the species level and reflected previously unexplored microbial taxa and their potential novelty. These MAGs also contained genes with predicted functions associated with diverse metabolic pathways that included carbohydrate, amino acid, and lipid metabolism pathways. Lastly, we performed a comparative analysis of Paucibacter MAGs and reference genomes that highlighted the presence of novel Paucibacter species and enriched metabolic potential in the recovered MAGs.

Patients with Primary and Secondary Bile Duct Stones Harbor Distinct Biliary Microbial Composition and Metabolic Potential.

Introduction: Cholelithiasis has a high incidence worldwide and limited treatment options due to its poorly understood pathogenesis. Furthermore, the role of biliary microbiota in cholelithiasis remains understudied. To address these questions, we performed microbial sequencing from biliary samples from primary bile duct stone (PBDS) and secondary bile duct stone (SBDS) patients. Results: We analyzed in total 45 biliary samples, including those from cholelithiasis patients with PBDS or SBDS and people with other digestive diseases. 16S rRNA sequencing showed the bacteria family Alcaligenaceae increased in relative abundance in the lithiasis group compared with the non-lithiasis group. In addition, the PBDS group showed significantly lower bacterial diversity than SBDS, with Propionibacteriaceae, Sphingomonadaceae, and Lactobacillaceae as the most significant bacteria families decreased in relative abundance. We further performed whole metagenomic shotgun sequencing (wMGS) and found increased ability of biofilm synthesis and the ability to sense external stimuli in PBDS based on functional annotation of mapped reads. From genome-resolved analysis of the samples, we identified 36 high-quality draft bacterial genome sequences with completion ≥70% and contamination ≤10%. Most of these genomes were classified into Proteobacteria, Firmicutes, or Actinobacteria. Conclusions: Our findings indicated that there is a subtle impact on biliary microbiome from cholelithiasis while the difference is more pronounced between the PBDS and SBDS. It was revealed that the diversity of biliary microbiota in PBDS is lower, while some metabolic pathways are up-regulated, including those linked to higher incidence of different types of cancer, providing new insights for the understanding of cholelithiasis with different origin.

Genome-resolved metagenomics using environmental and clinical samples.

Recent advances in high-throughput sequencing technologies and computational methods have added a new dimension to metagenomic data analysis i.e. genome-resolved metagenomics. In general terms, it refers to the recovery of draft or high-quality microbial genomes and their taxonomic classification and functional annotation. In recent years, several studies have utilized the genome-resolved metagenome analysis approach and identified previously unknown microbial species from human and environmental metagenomes. In this review, we describe genome-resolved metagenome analysis as a series of four necessary steps: (i) preprocessing of the sequencing reads, (ii) de novo metagenome assembly, (iii) genome binning and (iv) taxonomic and functional analysis of the recovered genomes. For each of these four steps, we discuss the most commonly used tools and the currently available pipelines to guide the scientific community in the recovery and subsequent analyses of genomes from any metagenome sample. Furthermore, we also discuss the tools required for validation of assembly quality as well as for improving quality of the recovered genomes. We also highlight the currently available pipelines that can be used to automate the whole analysis without having advanced bioinformatics knowledge. Finally, we will highlight the most widely adapted and actively maintained tools and pipelines that can be helpful to the scientific community in decision making before they commence the analysis.

Environmental concentrations of antibiotics alter the zebrafish gut microbiome structure and potential functions.

A paradoxical impact of high rates of production and consumption of antibiotics is their widespread release in the environment. Consequently, low concentrations of antibiotics and their byproducts have been routinely identified from various environmental settings especially from aquatic environments. However, the impact of such low concentrations of antibiotics on the exposed host especially in early life remains poorly understood. We exposed zebrafish to two different environmental concentrations of oxytetracycline and sulfamethoxazole, from larval stage to adulthood (∼120 days) and characterized their impact on the taxonomic diversity, antibiotic resistance genes, and metabolic pathways of the gut microbiome using metagenomic shotgun sequencing and analysis. Long term exposure of environmental concentrations of oxytetracycline and sulfamethoxazole significantly impacted the taxonomic composition and metabolic pathways of zebrafish gut microbiome. The antibiotic exposed samples exhibited significant enrichment of multiple flavobacterial species, including Flavobacterium sp. F52, Flavobacterium johnsoniae and Flavobacterium sp. Fl, which are well known pathogenic bacteria. The relative abundance of antibiotic resistance genes, especially several tetratcycline and sulfonamide resistance genes were significantly higher in the exposed samples and showed a linear correlation with the antibiotic concentrations. Furthermore, several metabolic pathways, including folate biosynthesis, oxidative phosphorylation, and biotin metabolism pathways, showed significant enrichment in the antibiotic exposed samples. Collectively, our results suggest that early life exposure of the environmental concentrations of antibiotics can increase the abundance of unfavorable bacteria, antibiotic resistance genes and associated pathways in the gut microbiome of zebrafish.

Prediction and analysis of metagenomic operons via MetaRon: a pipeline for prediction of Metagenome and whole-genome opeRons.

BACKGROUND: Efficient regulation of bacterial genes in response to the environmental stimulus results in unique gene clusters known as operons. Lack of complete operonic reference and functional information makes the prediction of metagenomic operons a challenging task; thus, opening new perspectives on the interpretation of the host-microbe interactions. RESULTS: In this work, we identified whole-genome and metagenomic operons via MetaRon (Metagenome and whole-genome opeRon prediction pipeline). MetaRon identifies operons without any experimental or functional information. MetaRon was implemented on datasets with different levels of complexity and information. Starting from its application on whole-genome to simulated mixture of three whole-genomes (E. coli MG1655, Mycobacterium tuberculosis H37Rv and Bacillus subtilis str. 16), E. coli c20 draft genome extracted from chicken gut and finally on 145 whole-metagenome data samples from human gut. MetaRon consistently achieved high operon prediction sensitivity, specificity and accuracy across E. coli whole-genome (97.8, 94.1 and 92.4%), simulated genome (93.7, 75.5 and 88.1%) and E. coli c20 (87, 91 and 88%,), respectively. Finally, we identified 1,232,407 unique operons from 145 paired-end human gut metagenome samples. We also report strong association of type 2 diabetes with Maltose phosphorylase (K00691), 3-deoxy-D-glycero-D-galacto-nononate 9-phosphate synthase (K21279) and an uncharacterized protein (K07101). CONCLUSION: With MetaRon, we were able to remove two notable limitations of existing whole-genome operon prediction methods: (1) generalizability (ability to predict operons in unrelated bacterial genomes), and (2) whole-genome and metagenomic data management. We also demonstrate the use of operons as a subset to represent the trends of secondary metabolites in whole-metagenome data and the role of secondary metabolites in the occurrence of disease condition. Using operonic data from metagenome to study secondary metabolic trends will significantly reduce the data volume to more precise data. Furthermore, the identification of metabolic pathways associated with the occurrence of type 2 diabetes (T2D) also presents another dimension of analyzing the human gut metagenome. Presumably, this study is the first organized effort to predict metagenomic operons and perform a detailed analysis in association with a disease, in this case type 2 diabetes. The application of MetaRon to metagenomic data at diverse scale will be beneficial to understand the gene regulation and therapeutic metagenomics.

Characterization of the Gastrointestinal and Reproductive Tract Microbiota in Fertile and Infertile Pakistani Couples.

The human microbiota is recognized as a vital "virtual" organ of the human body that influences human health, metabolism, and physiology. While the microbiomes of the gut, oral cavity, and skin have been extensively studied in the literature, relatively little work has been done on characterizing the microbiota of the human reproductive tract organs, and specifically on investigating its association to fertility. Here, we implemented a 16S ribosomal RNA (rRNA) amplicon sequencing approach to sequence and characterize the gut and genital tract microbiomes from several married Pakistani couples. The recruited individuals included 31 fertile and 35 infertile individuals, with ages ranging from 19-45 years. We identified several fluctuations in the diversity and composition of the gut and genital microbiota among fertile and infertile samples. For example, measures of α-diversity varied significantly between the genital samples donated by fertile and infertile men and there was overall greater between-sample variability in genital samples regardless of gender. In terms of taxonomic composition, Actinobacteria, Bacteroidetes, and Firmicutes fluctuated significantly between the gut microbiomes of fertile and infertile samples. Finally, biomarker analyses identified features (genera and molecular functions and pathways) that differed significantly between the fertile and infertile samples and in the past have been associated with bacterial vaginosis. However, we emphasize that 16S amplicon data alone has no bearing on individual health and is merely representative of microbial taxonomic differences that could also arise due to multiple other factors. Our findings, however, represent the first effort to characterize the microbiome associated with fertile and infertile couples in Pakistan and will hopefully pave the way for more comprehensive and broad-scale investigations in the future.

DNA phosphorothioate modification facilitates the dissemination of mcr-1 and blaNDM-1 in drinking water supply systems.

The mechanism driving the dissemination of antibiotic resistance genes (ARGs) in drinking water supply systems (DWSSs) with multiple barriers remains poorly understood despite several recent efforts. Phosphorothioate (PT) modifications, governed by dndABCDE genes, occur naturally in various bacteria and involve the incorporation of sulfur into the DNA backbone. PT is regarded as a mild antioxidant in vivo and is known to provide protection against bacterial genomes. We combined quantitative polymerase chain reaction, metagenomic, and network analyses for the water treatment process and laboratory-scale experiments for chlorine treatment using model strains to determine if DNA PT modification occurred in DWSS and facilitated the dissemination of mobilized colistin resistance-1 (mcr-1) and New Delhi metallo-ß-lactamase-1 (blaNDM-1) in DWSS. Our results indicated that the relative abundance of dndB increased in the effluent, compared with the influent, in the water treatment plants. Presence of dndB copies had a positive correlation with the concentration of chloramine disinfectant. Network analysis revealed Bdellovibrio as a potential host for MCR genes, NDM genes, and dndB in the DWSS. E. coli DH10B (Wild-type with the dndABCDE gene cluster and ΔdndB) model strains were used to investigate resistance to chlorine treatment at the concentration range of 0.5-3 mg/L. The resistance of the wild-type strain increased with increasing concentration of chlorine. DNA PT modification protected MCR- and NDM-carrying bacteria from chloramine disinfection during the water treatment process. The higher relative abundance of ARGs in the effluent of the water treatment plants may be due to the resistance of DNA PT modification to chloramine disinfection, thereby causing the enrichment of genera carrying MCR, NDM, and dndB. This study provides a new understanding on the mechanism of ARG dissemination in DWSS, which will help to improve the performance of drinking water treatment to control the risk associated with antibiotic-resistant bacteria.

Dynamics of the Salivary Microbiome During Different Phases of Crohn's Disease.

Crohn's disease is a chronic disorder that typically affects the gastrointestinal tract. The increased incidence in the recent years, especially in Asian countries, prompts for performing studies and gain newer insights into the etiology and pathogenesis of the disease. Among other causative factors, gut microbiome and its cross-talk with the salivary microbiome is a known factor that has a plausible role in the pathogenesis of Crohn's disease. The gut microbiome has been extensively studied, however, the salivary microbiome and its dynamics during different phases of this disease remain understudied. In this study, we obtained saliva samples from the patients during active and remission phases of the disease and compared them with control samples and highlighted the differences in taxonomic as well as predicted functional pathways among them. Our results indicated that the α and ß diversities were significantly lower during the active phase in contrast with remission phase and healthy samples. In general, Firmicutes were most abundant among the three sample groups, followed by Bacteroidetes and Proteobacteria. Genus level distribution highlighted Streptococcus, Neisseria, Prevotella, Haemophilus, and Veillonella as the five most abundant taxa. Differential abundance analysis of the three sample groups identified significant enrichment of 30 bacterial taxa in the active phase that included g_Prevotella, f_Prevotellaceae, and p_Bacteroidetes. Furthermore, remission phase and control also exhibited significant enrichment of 24 and 22 bacterial taxa, respectively. Eleven differentially abundant pathways were also identified, four were significantly enriched in healthy controls whereas other seven were significantly enriched in active phase of the disease. Several important pathways, such as ribosome biogenesis and Energy metabolism were depleted in the active phase. Our study has highlighted several taxa and functional categories that could be implicated with the onset of Crohn's disease and thus have the potential to serve as biomarkers of the active disease. However, these findings require further validation through functional studies in the future.

Metagenomic analysis of basal ice from an Alaskan glacier.

BACKGROUND: Glaciers cover ~ 10% of land but are among the least explored environments on Earth. The basal portion of glaciers often harbors unique aquatic microbial ecosystems in the absence of sunlight, and knowledge on the microbial community structures and their metabolic potential is very limited. Here, we provide insights into the microbial lifestyle present at the base of the Matanuska Glacier, Alaska. RESULTS: DNA and RNA were extracted from samples of the Matanuska Glacier basal ice. Using Illumina MiSeq and HiSeq sequencing, we investigated the microbial diversity with the metagenomic shotgun reads and 16S ribosomal RNA data. We further assembled 9 partial and draft bacterial genomes from the metagenomic assembly, and identified key metabolic pathways such as sulfur oxidation and nitrification. Collectively, our analyses suggest a prevalence of lithotrophic and heterotrophic metabolisms in the subglacial microbiome. CONCLUSION: Our results present the first metagenomic assembly and bacterial draft genomes for a subglacial environment. These results extend our understanding of the chemical and biological processes in subglacial environments critically influenced by global climate change.

Comprehensive analysis of circRNA expression profiles in humans by RAISE.

Circular RNAs (circRNAs) are pervasively expressed circles of non­coding RNAs. Even though many circRNAs have been reported in humans, their expression patterns and functions remain poorly understood. In this study, we employed a pipeline named RAISE to detect circRNAs in RNA­seq data. RAISE can fully characterize circRNA structure and abundance. We evaluated inter-individual variations in circRNA expression in humans by applying this pipeline to numerous non­poly(A)-selected RNA­seq data. We identified 59,128 circRNA candidates in 61 human liver samples, with almost no overlap in the circRNA of the recruited samples. Approximately 89% of the circRNAs were detected in one or two samples. In comparison, 10% of the linear mRNAs and non­coding RNAs were detected in each sample. We estimated the variation in other tissues, especially the circRNA high-abundance tissues, in advance. Only 0.5% of the 50,631 brain circRNA candidates were shared among the 30 recruited brain samples, which is similar to the proportion in liver. Moreover, we found inter- and intra-individual diversity in circRNAs expression in the granulocyte RNA­seq data from seven individuals sampled 3 times at one-month intervals. Our findings suggest that careful consideration of inter-individual diversity is required when extensively identifying human circRNAs or proposing their use as potential biomarkers and therapeutic targets in disease.

The essential mycobacterial amidotransferase GatCAB is a modulator of specific translational fidelity.

Although regulation of translation fidelity is an essential process1-7, diverse organisms and organelles have differing requirements of translational accuracy8-15, and errors in gene translation serve an adaptive function under certain conditions16-20. Therefore, optimal levels of fidelity may vary according to context. Most bacteria utilize a two-step pathway for the specific synthesis of aminoacylated glutamine and/or asparagine tRNAs, involving the glutamine amidotransferase GatCAB21-25, but it had not been appreciated that GatCAB may play a role in modulating mistranslation rates. Here, by using a forward genetic screen, we show that the mycobacterial GatCAB enzyme complex mediates the translational fidelity of glutamine and asparagine codons. We identify mutations in gatA that cause partial loss of function in the holoenzyme, with a consequent increase in rates of mistranslation. By monitoring single-cell transcription dynamics, we demonstrate that reduced gatCAB expression leads to increased mistranslation rates, which result in enhanced rifampicin-specific phenotypic resistance. Consistent with this, strains with mutations in gatA from clinical isolates of Mycobacterium tuberculosis show increased mistranslation, with associated antibiotic tolerance, suggesting a role for mistranslation as an adaptive strategy in tuberculosis. Together, our findings demonstrate a potential role for the indirect tRNA aminoacylation pathway in regulating translational fidelity and adaptive mistranslation.

Transmitted Extended-Spectrum Extensively Drug-Resistant Tuberculosis in Beijing, China, with Discordant Whole-Genome Sequencing Analysis Results.

Drug resistance to tuberculosis remains a major public health threat. Here, we report two cases of extended-spectrum extensively drug-resistant (XXDR) tuberculosis showing resistance to most first- and second-line agents. The results of a correlation of whole-genome sequencing (WGS) and phenotypic testing were discordant, suggesting that overreliance on WGS may miss clinically relevant resistance in extensively drug-resistant disease.

Role of miRNAs in breast cancer.

miRNAs belong to an important class of endogenous molecules which are present in a wide range of organisms including animals, plants and viruses. They are involved in regulating expression of several genes inside a cell due to presence of complementary region against specific mRNA molecules. Altered expression patterns cause progression of multiple diseases inside an organism. They have also been confirmed to be involved in different cancers including breast cancer. In this review, we discuss role of miRNAs with respect to uncontrolled division of cells, promotion, progression and metastasis in breast cancer.