SummaryTuberculosis (TB) is primarily associated with decline in immune health status. As gut microbiome (GM) is implicated in the regulation of host immunity and metabolism, here we investigate GM alteration in TB patients by 16S rRNA gene and whole-genome shotgun sequencing. The study group constituted of patients with pulmonary TB and their healthy household contacts as controls (HCs). Significant alteration of microbial taxonomic and functional capacity was observed in patients with active TB as compared to the HCs. We observed that Prevotella and Bifidobacterium abundance were associated with HCs, whereas butyrate and propionate-producing bacteria like Faecalibacterium, Roseburia, Eubacterium and Phascolarctobacterium were significantly enriched in TB patients. Functional analysis showed reduced biosynthesis of vitamins and amino acids in favour of enriched metabolism of butyrate and propionate in TB subjects. The TB subjects were also investigated during the course of treatment, to analyse the variation of GM. Although perturbation in microbial composition was still evident after a month's administration of anti-TB drugs, significant changes were observed in metagenome gene pool that pointed towards recovery in functional capacity. Therefore, the findings from this pilot study suggest that microbial dysbiosis may contribute to pathophysiology of TB by enhancing the anti-inflammatory milieu in the host.

No abstract

The outbreak of coronavirus disease 2019 (COVID-19) that started in Wuhan, China, in December 2019 has spread worldwide, emerging as a global pandemic. The severe respiratory pneumonia caused by novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has so far claimed more than 0.38 million lives and has impacted human lives worldwide. However, as the novel SARS-CoV-2 virus displays high transmission rates, the underlying genomic severity is required to be fully understood. We studied the complete genomes of 95 SARS-CoV-2 strains from different geographical regions worldwide to uncover the pattern of the spread of the virus. We show that there is no direct transmission pattern of the virus among neighboring countries, suggesting that its spread is a result of travel of infected humans to different countries. We revealed unique single nucleotide polymorphisms (SNPs) in nonstructural protein 13 (nsp13), nsp14, nsp15, and nsp16 (ORF1b polyproteins) and in the S-protein within 10 viral isolates from the United States. These viral proteins are involved in RNA replication and binding with the human receptors, indicating that the viral variants that are circulating in the population of the United States are different from those circulating in the populations of other countries. In addition, we found an amino acid addition in nsp16 (mRNA cap-1 methyltransferase) of a U.S. isolate (GenBank accession no. MT188341.1) leading to a shift in the amino acid frame from position 2540 onward. Through comparative structural analysis of the wild-type and mutant proteins, we showed that this addition of a phenylalanine residue renders the protein in the mutant less stable, which might affect mRNA cap-1 methyltransferase function. We further analyzed the SARS-CoV-2–human interactome, which revealed that the interferon signaling pathway is targeted by orf1ab during infection and that it also interacts with NF-κB-repressing factor (NKRF), which is a potential regulator of interleukin-8 (IL-8). We propose that targeting this interaction may subsequently improve the health condition of COVID-19 patients. Our analysis also emphasized that SARS-CoV-2 manipulates spliceosome machinery during infection; hence, targeting splicing might affect viral replication. In conclusion, the replicative machinery of SARS-CoV-2 is targeting interferon and the notch signaling pathway along with spliceosome machinery to evade host challenges. IMPORTANCE The COVID-19 pandemic continues to storm the world, with over 6.5 million cases worldwide. The severity of the disease varies with the territories and is mainly influenced by population density and age factor. In this study, we analyzed the transmission pattern of 95 SARS-CoV-2 genomes isolated from 11 different countries. Our study also revealed several nonsynonymous mutations in ORF1b and S-proteins and the impact on their structural stability. Our analysis showed the manipulation of host system by viral proteins through SARS-CoV-2–human protein interactome, which can be useful to understand the impact of virus on human health.

Bdellovibrio bacteriovorus are small Deltaproteobacteria that invade, kill and assimilate their prey. Metagenomic assembly analysis of the microbial mats of an arsenic rich, hot spring was performed to describe the genotypes of the predator Bdellovibrio and the ecogenetically adapted taxa Enterobacter. The microbial mats were enriched with Bdellovibrio (1.3%) and several Gram-negative bacteria including Bordetella (16%), Enterobacter (6.8%), Burkholderia (4.8%), Acinetobacter (2.3%) and Yersinia (1%). A high-quality (47 contigs, 25X coverage; 3.5 Mbp) draft genome of Bdellovibrio (strain ArHS; Arsenic Hot Spring) was reassembled, which lacked the marker gene Bd0108 associated with the usual method of prey interaction and invasion for this genus, while maintaining genes coding for the hydrolytic enzymes necessary for prey assimilation. By filtering microbial mat samples (< 0.45 μm) to enrich for small predatory cell sizes, we observed Bdellovibrio-like cells attached side-on to E. coli through electron microscopy. Furthermore, a draft pan-genome of the dominant potential host taxon, Enterobacter cloacae ArHS (4.8 Mb), along with three of its viral genotypes (n = 3; 42 kb, 49 kb and 50 kb), was assembled. These data were further used to analyse the population level evolutionary dynamics (taxonomical and functional) of reconstructed genotypes.

Aims/hypothesisEvaluation of the association of 31 common single nucleotide polymorphisms (SNPs) with fasting glucose, fasting insulin, HOMA-beta cell function (HOMA-β), HOMA-insulin resistance (HOMA-IR) and type 2 diabetes in the Indian population.MethodsWe genotyped 3,089 sib pairs recruited in the Indian Migration Study from four cities in India (Lucknow, Nagpur, Hyderabad and Bangalore) for 31 SNPs in 24 genes previously associated with type 2 diabetes in European populations. We conducted within-sib-pair analysis for type 2 diabetes and its related quantitative traits.ResultsThe risk-allele frequencies of all the SNPs were comparable with those reported in western populations. We demonstrated significant associations of CXCR4 (rs932206), CDKAL1 (rs7756992) and TCF7L2 (rs7903146, rs12255372) with fasting glucose, with β values of 0.007 (p = 0.05), 0.01 (p = 0.01), 0.007 (p = 0.05), 0.01 (p = 0.003) and 0.08 (p = 0.01), respectively. Variants in NOTCH2 (rs10923931), TCF-2 (also known as HNF1B) (rs757210), ADAM30 (rs2641348) and CDKN2A/B (rs10811661) significantly predicted fasting insulin, with β values of −0.06 (p = 0.04), 0.05 (p = 0.05), −0.08 (p = 0.01) and −0.08 (p = 0.02), respectively. For HOMA-IR, we detected associations with TCF-2, ADAM30 and CDKN2A/B, with β values of 0.05 (p = 0.04), −0.07 (p = 0.03) and −0.08 (p = 0.02), respectively. We also found significant associations of ADAM30 (β = −0.05; p = 0.01) and CDKN2A/B (β = −0.05; p = 0.03) with HOMA-β. THADA variant (rs7578597) was associated with type 2 diabetes (OR 1.5; 95% CI 1.04, 2.22; p = 0.03).Conclusions/interpretationWe validated the association of seven established loci with intermediate traits related to type 2 diabetes in an Indian population using a design resistant to population stratification.Electronic supplementary materialThe online version of this article (doi:10.1007/s00125-011-2355-6) contains peer-reviewed but unedited supplementary material, which is available to authorised users.

Few studies have investigated the association between genetic variation and obesity traits in Indian populations or the role of environmental factors as modifiers of these relationships. In the context of rapid urbanisation, resulting in significant lifestyle changes, understanding the aetiology of obesity is important. We investigated associations of FTO and MC4R variants with obesity traits in 3390 sibling pairs from four Indian cities, most of whom were discordant for current dwelling (rural or urban). The FTO variant rs9939609 predicted increased weight (0.09 Z-scores, 95% CI: 0.03, 0.15) and BMI (0.08 Z-scores, 95% CI: 0.02, 0.14). The MC4R variant rs17782313 was weakly associated with weight and hip circumference (P < .05). There was some indication that the association between FTO and weight was stronger in urban than that in rural dwellers (P for interaction = .03), but no evidence for effect modification by diet or physical activity. Further studies are needed to investigate ways in which urban environment may modify genetic risk of obesity.

SummaryThe aim of this study was to validate the single nucleotide polymorphisms (SNPs) of four candidate genes (TCF7L2, HHEX, KCNJ11, and ADIPOQ) related to type 2 diabetes (T2D) in an endogamous population of north India; the Aggarwal population, having 18-clans. This endogamous population model was heavily supported by recent land mark work and we also verified the homogeneity of this population by clan-based stratification analysis. Two SNPs (rs4506565; rs7903146) in TCF7L2 were found to be significant (p-value = 0.00191; p-value = 0.00179, respectively), and odds ratios of 2.1 (dominant-model) and 2.0 (recessive-model) respectively, were obtained for this population. The TTT haplotype in the TCF7L2 gene was significantly associated with T2D. Waist-Hip ratio (WHR), systolic blood pressure (SBP), and age were significant covariates for increasing risk of T2D. Single-SNP, combined-SNPs and haplotype analysis provides clear evidence that the causal mutation is near to or within the significant haplotype (TTT) of the TCF7L2 gene. In spite of a culturally-learned sedentary lifestyle and fat-enriched dietary habits, WHR rather than body-mass-index emerged as a robust predictor of risk for T2D in this population.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has resulted in 92 million cases in a span of 1 year. The study focuses on understanding population-specific variations attributing its high rate of infections in specific geographical regions particularly in the United States. Rigorous phylogenomic network analysis of complete SARS-CoV-2 genomes (245) inferred five central clades named a (ancestral), b, c, d, and e (subtypes e1 and e2). Clade d and subclade e2 were found exclusively comprised of U.S. strains. Clades were distinguished by 10 co-mutational combinations in Nsp3, ORF8, Nsp13, S, Nsp12, Nsp2, and Nsp6. Our analysis revealed that only 67.46% of single nucleotide polymorphism (SNP) mutations were at the amino acid level. T1103P mutation in Nsp3 was predicted to increase protein stability in 238 strains except for 6 strains which were marked as ancestral type, whereas co-mutation (P409L and Y446C) in Nsp13 were found in 64 genomes from the United States highlighting its 100% co-occurrence. Docking highlighted mutation (D614G) caused reduction in binding of spike proteins with angiotensin-converting enzyme 2 (ACE2), but it also showed better interaction with the TMPRSS2 receptor contributing to high transmissibility among U.S. strains. We also found host proteins, MYO5A, MYO5B, and MYO5C, that had maximum interaction with viral proteins (nucleocapsid [N], spike [S], and membrane [M] proteins). Thus, blocking the internalization pathway by inhibiting MYO5 proteins which could be an effective target for coronavirus disease 2019 (COVID-19) treatment. The functional annotations of the host-pathogen interaction (HPI) network were found to be closely associated with hypoxia and thrombotic conditions, confirming the vulnerability and severity of infection. We also screened CpG islands in Nsp1 and N conferring the ability of SARS-CoV-2 to enter and trigger zinc antiviral protein (ZAP) activity inside the host cell. IMPORTANCE In the current study, we presented a global view of mutational pattern observed in SARS-CoV-2 virus transmission. This provided a who-infect-whom geographical model since the early pandemic. This is hitherto the most comprehensive comparative genomics analysis of full-length genomes for co-mutations at different geographical regions especially in U.S. strains. Compositional structural biology results suggested that mutations have a balance of opposing forces affecting pathogenicity suggesting that only a few mutations are effective at the translation level. Novel HPI analysis and CpG predictions elucidate the proof of concept of hypoxia and thrombotic conditions in several patients. Thus, the current study focuses the understanding of population-specific variations attributing a high rate of SARS-CoV-2 infections in specific geographical regions which may eventually be vital for the most severely affected countries and regions for sharp development of custom-made vindication strategies.