A cross-sectional study on the nasopharyngeal microbiota of individuals with SARS-CoV-2 infection across three COVID-19 waves in India.

Background: Multiple variants of the SARS-CoV-2 virus have plagued the world through successive waves of infection over the past three years. Independent research groups across geographies have shown that the microbiome composition in COVID-19 positive patients (CP) differs from that of COVID-19 negative individuals (CN). However, these observations were based on limited-sized sample-sets collected primarily from the early days of the pandemic. Here, we study the nasopharyngeal microbiota in COVID-19 patients, wherein the samples have been collected across the three COVID-19 waves witnessed in India, which were driven by different variants of concern. Methods: The nasopharyngeal swabs were collected from 589 subjects providing samples for diagnostics purposes at the Centre for Cellular and Molecular Biology (CSIR-CCMB), Hyderabad, India and subjected to 16s rRNA gene amplicon - based sequencing. Findings: We found variations in the microbiota of symptomatic vs. asymptomatic COVID-19 patients. CP showed a marked shift in the microbial diversity and composition compared to CN, in a wave-dependent manner. Rickettsiaceae was the only family that was noted to be consistently depleted in CP samples across the waves. The genera Staphylococcus, Anhydrobacter, Thermus, and Aerococcus were observed to be highly abundant in the symptomatic CP patients when compared to the asymptomatic group. In general, we observed a decrease in the burden of opportunistic pathogens in the host microbiota during the later waves of infection. Interpretation: To our knowledge, this is the first analytical cross-sectional study of this scale, which was designed to understand the relation between the evolving nature of the virus and the changes in the human nasopharyngeal microbiota. Although no clear signatures were observed, this study shall pave the way for a better understanding of the disease pathophysiology and help gather preliminary evidence on whether interventions to the host microbiota can help in better protection or faster recovery.

Sensing Host Health: Insights from Sensory Protein Signature of the Metagenome.

The human microbiota, which comprises an ensemble of taxonomically and functionally diverse but often mutually cooperating microorganisms, benefits its host by shaping the host immunity, energy harvesting, and digestion of complex carbohydrates as well as production of essential nutrients. Dysbiosis in the human microbiota, especially the gut microbiota, has been reported to be linked to several diseases and metabolic disorders. Recent studies have further indicated that tracking these dysbiotic variations could potentially be exploited as biomarkers of disease states. However, the human microbiota is not geography agnostic, and hence a taxonomy-based (microbiome) biomarker for disease diagnostics has certain limitations. In comparison, (microbiome) function-based biomarkers are expected to have a wider applicability. Given that (i) the host physiology undergoes certain changes in the course of a disease and (ii) host-associated microbial communities need to adapt to this changing microenvironment of their host, we hypothesized that signatures emanating from the abundance of bacterial proteins associated with the signal transduction system (herein referred to as sensory proteins [SPs]) might be able to distinguish between healthy and diseased states. To test this hypothesis, publicly available metagenomic data sets corresponding to three diverse health conditions, namely, colorectal cancer, type 2 diabetes mellitus, and schizophrenia, were analyzed. Results demonstrated that SP signatures (derived from host-associated metagenomic samples) indeed differentiated among healthy individual and patients suffering from diseases of various severities. Our finding was suggestive of the prospect of using SP signatures as early biomarkers for diagnosing the onset and progression of multiple diseases and metabolic disorders. IMPORTANCE The composition of the human microbiota, a collection of host-associated microbes, has been shown to differ among healthy and diseased individuals. Recent studies have investigated whether tracking these variations could be exploited for disease diagnostics. It has been noted that compared to microbial taxonomies, the ensemble of functional proteins encoded by microbial genes are less likely to be affected by changes in ethnicity and dietary preferences. These functions are expected to help the microbe adapt to changing environmental conditions. Thus, healthy individuals might harbor a different set of genes than diseased individuals. To test this hypothesis, we analyzed metagenomes from healthy and diseased individuals for signatures of a particular group of proteins called sensory proteins (SP), which enable the bacteria to sense and react to changes in their microenvironment. Results demonstrated that SP signatures indeed differentiate among healthy individuals and those suffering from diseases of various severities.

Can machines learn the mutation signatures of SARS-CoV-2 and enable viral-genotype guided predictive prognosis?

MOTIVATION: Continuous emergence of new variants through appearance/accumulation/disappearance of mutations is a hallmark of many viral diseases. SARS-CoV-2 variants have particularly exerted tremendous pressure on global healthcare system owing to their life threatening and debilitating implications. The sheer plurality of variants and huge scale of genomic data have added to the challenges of tracing the mutations/variants and their relationship to infection severity (if any). RESULTS: We explored the suitability of virus-genotype guided machine-learning in infection prognosis and identification of features/mutations-of-interest. Total 199,519 outcome-traced genomes, representing 45,625 nucleotide-mutations, were employed. Among these, post data-cleaning, Low and High severity genomes were classified using an integrated model (employing virus genotype, epitopic-influence and patient-age) with consistently high ROC-AUC (Asia:0.97 ± 0.01, Europe:0.94 ± 0.01, N.America:0.92 ± 0.02, Africa:0.94 ± 0.07, S.America:0.93 ± 03). Although virus-genotype alone could enable high predictivity (0.97 ± 0.01, 0.89 ± 0.02, 0.86 ± 0.04, 0.95 ± 0.06, 0.9 ± 0.04), the performance was not found to be consistent and the models for a few geographies displayed significant improvement in predictivity when the influence of age and/or epitope was incorporated with virus-genotype (Wilcoxon p_BH < 0.05). Neither age or epitopic-influence or clade information could out-perform the integrated features. A sparse model (6 features), developed using patient-age and epitopic-influence of the mutations, performed reasonably well (>0.87 ± 0.03, 0.91 ± 0.01, 0.87 ± 0.03, 0.84 ± 0.08, 0.89 ± 0.05). High-performance models were employed for inferring the important mutations-of-interest using Shapley Additive exPlanations (SHAP). The changes in HLA interactions of the mutated epitopes of reference SARS-CoV-2 were then subsequently probed. Notably, we also describe the significance of a 'temporal-modeling approach' to benchmark the models linked with continuously evolving pathogens. We conclude that while machine learning can play a vital role in identifying relevant mutations and factors driving the severity, caution should be exercised in using the genotypic signatures for predictive prognosis.

Does immune recognition of SARS-CoV2 epitopes vary between different ethnic groups?

The SARS-CoV2 mediated Covid-19 pandemic has impacted humankind at an unprecedented scale. While substantial research efforts have focused towards understanding the mechanisms of viral infection and developing vaccines/ therapeutics, factors affecting the susceptibility to SARS-CoV2 infection and manifestation of Covid-19 remain less explored. Given that the Human Leukocyte Antigen (HLA) system is known to vary among ethnic populations, it is likely to affect the recognition of the virus, and in turn, the susceptibility to Covid-19. To understand this, we used bioinformatic tools to probe all SARS-CoV2 peptides which could elicit T-cell response in humans. We also tried to answer the intriguing question of whether these potential epitopes were equally immunogenic across ethnicities, by studying the distribution of HLA alleles among different populations and their share of cognate epitopes. Results indicate that the immune recognition potential of SARS-CoV2 epitopes tend to vary between different ethnic groups. While the South Asians are likely to recognize higher number of CD8-specific epitopes, Europeans are likely to identify higher number of CD4-specific epitopes. We also hypothesize and provide clues that the newer mutations in SARS-CoV2 are unlikely to alter the T-cell mediated immunogenic responses among the studied ethnic populations. The work presented herein is expected to bolster our understanding of the pandemic, by providing insights into differential immunological response of ethnic populations to the virus as well as by gaging the possible effects of mutations in SARS-CoV2 on efficacy of potential epitope-based vaccines through evaluating ∼40,000 viral genomes.

Trans-ethnic gut microbiota signatures of type 2 diabetes in Denmark and India.

BACKGROUND: Type 2 diabetes (T2D), a multifactorial disease influenced by host genetics and environmental factors, is the most common endocrine disease. Several studies have shown that the gut microbiota as a close-up environmental mediator influences host physiology including metabolism. The aim of the present study is to examine the compositional and functional potential of the gut microbiota across individuals from Denmark and South India with a focus on T2D. Many earlier studies have investigated the microbiome aspects of T2D, and it has also been anticipated that such microbial associations would be dependent on diet and ethnic origin. However, there has been no large scale trans-ethnic microbiome study earlier in this direction aimed at evaluating any "universal" microbiome signature of T2D. METHODS: 16S ribosomal RNA gene amplicon sequencing was performed on stool samples from 279 Danish and 294 Indian study participants. Any differences between the gut microbiota of both populations were explored using diversity measures and negative binomial Wald tests. Study samples were stratified to discover global and country-specific microbial signatures for T2D and treatment with the anti-hyperglycemic drug, metformin. To identify taxonomical and functional signatures of the gut microbiota for T2D and metformin treatment, we used alpha and beta diversity measures and differential abundances analysis, comparing metformin-naive T2D patients, metformin-treated T2D patients, and normoglycemic individuals. RESULTS: Overall, the gut microbial communities of Danes and Indians are compositionally very different. By analyzing the combined study materials, we identify microbial taxonomic and functional signatures for T2D and metformin treatment. T2D patients have an increased relative abundance of two operational taxonomic units (OTUs) from the Lachnospiraceae family, and a decreased abundance of Subdoligranulum and Butyricicoccus. Studying each population per se, we identified T2D-related microbial changes at the taxonomic level within the Danish population only. Alpha diversity indices show that there is no significant difference between normoglycemic individuals and metformin-naive T2D patients, whereas microbial richness is significantly decreased in metformin-treated T2D patients compared to metformin-naive T2D patients and normoglycemic individuals. Enrichment of two OTUs from Bacteroides and depletion of Faecalibacterium constitute a trans-ethnic signature of metformin treatment. CONCLUSIONS: We demonstrate major compositional differences of the gut microbiota between Danish and South Indian individuals, some of which may relate to differences in ethnicity, lifestyle, and demography. By comparing metformin-naive T2D patients and normoglycemic individuals, we identify T2D-related microbiota changes in the Danish and Indian study samples. In the present trans-ethnic study, we confirm that metformin changes the taxonomic profile and functional potential of the gut microbiota.

Trans-ethnic gut microbial signatures of prediabetic subjects from India and Denmark.

BACKGROUND: Recent studies have indicated an association of gut microbiota and microbial metabolites with type 2 diabetes mellitus (T2D). However, large-scale investigation of the gut microbiota of "prediabetic" (PD) subjects has not been reported. Identifying robust gut microbiome signatures of prediabetes and characterizing early prediabetic stages is important for the understanding of disease development and could be crucial in early diagnosis and prevention. METHODS: The current study performed amplification and sequencing on the variable regions (V1-V5) of the 16S rRNA genes to profile and compare gut microbiota of prediabetic individuals (N = 262) with normoglycemic individuals (N = 275) from two cohorts in India and Denmark. Similarly, fasting serum inflammatory biomarkers were profiled from the study participants. RESULTS: After correcting for strong country-specific cohort effect, 16 operational taxonomic units (OTUs) including members from the genera Prevotella9, Phascolarctobacterium, Barnesiella, Flavonifractor, Tyzzerella_4, Bacteroides, Faecalibacterium, and Agathobacter were identified as enriched in normoglycaemic subjects with respect to the subjects with prediabetes using a negative binomial Wald test. We also identified 144 OTUs enriched in the prediabetic subjects, which included members from the genera Megasphaera, Streptococcus, Prevotella9, Alistipes, Mitsuokella, Escherichia/Shigella, Prevotella2, Vibrio, Lactobacillus, Alloprevotella, Rhodococcus, and Klebsiella. Comparative analyses of relative abundance of bacterial taxa revealed that the Streptococcus, Escherichia/Shigella, Prevotella2, Vibrio, and Alloprevotella OTUs exhibited more than fourfold enrichment in the gut microbiota of prediabetic subjects. When considering subjects from the two geographies separately, we were able to identify additional gut microbiome signatures of prediabetes. The study reports a probable association of Megasphaera OTU(s) with impaired glucose tolerance, which is significantly pronounced in Indian subjects. While the overall results confirm a state of proinflammation as early as in prediabetes, the Indian cohort exhibited a characteristic pattern of abundance of inflammatory markers indicating low-grade intestinal inflammation at an overall population level, irrespective of glycemic status. CONCLUSIONS: The results present trans-ethnic gut microbiome and inflammation signatures associated with prediabetes, in Indian and Danish populations. The identified associations may be explored further as potential early indicators for individuals at risk of dysglycemia.

Can Targeting Non-Contiguous V-Regions With Paired-End Sequencing Improve 16S rRNA-Based Taxonomic Resolution of Microbiomes?: An In Silico Evaluation.

Background: Next-generation sequencing (NGS) technologies have enabled probing of microbial diversity in different environmental niches with unprecedented sequencing depth. However, due to read-length limitations of popular NGS technologies, 16S amplicon sequencing-based microbiome studies rely on targeting short stretches of the 16S rRNA gene encompassing a selection of variable (V) regions. In most cases, such a short stretch constitutes a single V-region or a couple of V-regions placed adjacent to each other on the 16S rRNA gene. Given that different V-regions have different resolving ability with respect to various taxonomic groups, selecting the optimal V-region (or a combination thereof) remains a challenge. Methods: The accuracy of taxonomic profiles generated from sequences encompassing 1) individual V-regions, 2) adjacent V-regions, and 3) pairs of non-contiguous V-regions were assessed and compared. Subsequently, the discriminating capability of different V-regions with respect to different taxonomic lineages was assessed. The possibility of using paired-end sequencing protocols to target combinations of non-adjacent V-regions was finally evaluated with respect to the utility of such an experimental design in providing improved taxonomic resolution. Results: Extensive validation with simulated microbiome datasets mimicking different environmental and host-associated microbiome samples suggest that targeting certain combinations of non-contiguously placed V-regions might yield better taxonomic classification accuracy compared to conventional 16S amplicon sequencing targets. This work also puts forward a novel in silico combinatorial strategy that enables creation of consensus taxonomic profiles from experiments targeting multiple pair-wise combinations of V-regions to improve accuracy in taxonomic classification. Conclusion: The study suggests that targeting non-contiguous V-regions with paired-end sequencing can improve 16S rRNA-based taxonomic resolution of microbiomes. Furthermore, employing the novel in silico combinatorial strategy can improve taxonomic classification without any significant additional experimental costs and/or efforts. The empirical observations obtained can potentially serve as a guideline for future 16S microbiome studies, and facilitate researchers in choosing the optimal combination of V-regions for a specific experiment/sampled environment.

Correction: Alterations in the gut bacterial microbiome in fungal Keratitis patients.

[This corrects the article DOI: 10.1371/journal.pone.0199640.].

Alterations in the gut bacterial microbiome in fungal Keratitis patients.

Dysbiosis in the gut microbiome has been implicated in several diseases including auto-immune diseases, inflammatory diseases, cancers and mental disorders. Keratitis is an inflammatory disease of the eye significantly contributing to corneal blindness in the developing world. It would be worthwhile to investigate the possibility of dysbiosis in the gut microbiome being associated with Keratitis. Here, we have analyzed fungal and bacterial populations in stool samples through high-throughput sequencing of the ITS2 region for fungi and V3-V4 region of 16S rRNA gene for bacteria in healthy controls (HC, n = 31) and patients with fungal keratitis (FK, n = 32). Candida albicans (2 OTUs), Aspergillus (1 OTU) and 3 other denovo-OTUs were enriched in FK samples and an unclassified denovo-OTU was enriched in HC samples. However, the overall abundances of these 'discriminatory' OTUs were very low (< 0.001%) and not indicative of significant dysbiosis in the fungal community inhabiting the gut of FK patients. In contrast, the gut bacterial richness and diversity in FK patients was significantly decreased when compared to HC. 52 OTUs were significantly enriched in HC samples whereas only 5 OTUs in FK. The OTUs prominently enriched in HC were identified as Faecalibacterium prausnitzii, Bifidobacterium adolescentis, Lachnospira, Mitsuokella multacida, Bacteroides plebeius, Megasphaera and Lachnospiraceae. In FK samples, 5 OTUs affiliated to Bacteroides fragilis, Dorea, Treponema, Fusobacteriaceae, and Acidimicrobiales were significantly higher in abundance. The functional implications are that Faecalibacterium prausnitzii, an anti-inflammatory bacterium and Megasphaera, Mitsuokella multacida and Lachnospira are butyrate producers, which were enriched in HC patients, whereas Treponema and Bacteroides fragilis, which are pathogenic were abundant in FK patients, playing a potential pro-inflammatory role. Heatmap, PCoA plots and functional profiles further confirm the distinct patterns of gut bacterial composition in FK and HC samples. Our study demonstrates dysbiosis in the gut bacterial microbiomes of FK patients compared to HC. Further, based on inferred functions, it appears that dysbiosis in the gut of FK subjects is strongly associated with the disease phenotype with decrease in abundance of beneficial bacteria and increase in abundance of pro-inflammatory and pathogenic bacteria.