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.

Exploring the human microbiome - A step forward for precision medicine in breast cancer.

BACKGROUND: The second most frequent cancer in the world and the most common malignancy in women is breast cancer. Breast cancer is a significant health concern in India with a high mortality-to-incidence ratio and presentation at a younger age. RECENT FINDINGS: Recent studies have identified gut microbiota as a significant factor that can have an influence on the development, treatment, and prognosis of breast cancer. This review article aims to describe the influence of microbial dysbiosis on breast cancer occurrence and the possible interactions between oncobiome and specific breast cancer molecular subtypes. The review further also discusses the role of epigenetics and diet/nutrition in the regulation of the gut and breast microbiome and its association with breast cancer prevention, therapy, and recurrence. Additionally, the recent technological advances in microbiome research, including next-generation sequencing (NGS) technologies, genome sequencing, single-cell sequencing, and microbial metabolomics along with recent advances in artificial intelligence (AI) have also been reviewed. This is an attempt to present a comprehensive status of the microbiome as a key cancer biomarker. CONCLUSION: We believe that correlating microbiome and carcinogenesis is important as it can provide insights into the mechanisms by which microbial dysbiosis can influence cancer development and progression, leading to the potential use of the microbiome as a tool for prognostication and personalized therapy.

Transcriptome analysis of human macrophages upon chikungunya virus (CHIKV) infection reveals regulation of distinct signaling and metabolic pathways during the early and late stages of infection.

Macrophages are efficient reservoirs for viruses that enable the viruses to survive over a longer period of infection. Alphaviruses such as chikungunya virus (CHIKV) are known to persist in macrophages even after the acute febrile phase. The viral particles replicate in macrophages at a very low level over extended period of time and are localized in tissues that are often less accessible by treatment. Comprehensive experimental studies are thus needed to characterize the CHIKV-induced modulation of host genes in these myeloid lineage cells and in one such pursuit, we obtained global transcriptomes of a human macrophage cell line infected with CHIKV, over its early and late timepoints of infection. We analyzed the pathways, especially immune related, perturbed over these timepoints and observed several host factors to be differentially expressed in infected macrophages in a time-dependent manner. We postulate that these pathways may play crucial roles in the persistence of CHIKV in macrophages.

Frequency spectrum of rare and clinically relevant markers in multiethnic Indian populations (ClinIndb): A resource for genomic medicine in India.

There have been concerted efforts toward cataloging rare and deleterious variants in different world populations using high-throughput genotyping and sequencing-based methods. The Indian population is underrepresented or its information with respect to clinically relevant variants is sparse in public data sets. The aim of this study was to estimate the burden of monogenic disease-causing variants in Indian populations. Toward this, we have assessed the frequency profile of monogenic phenotype-associated ClinVar variants. The study utilized a genotype data set (global screening array, Illumina) from 2795 individuals (multiple in-house genomics cohorts) representing diverse ethnic and geographically distinct Indian populations. Of the analyzed variants from Global Screening Array, ~9% were found to be informative and were either not known earlier or underrepresented in public databases in terms of their frequencies. These variants were linked to disorders, namely inborn errors of metabolism, monogenic diabetes, hereditary cancers, and various other hereditary conditions. We have also shown that our study cohort is genetically a better representative of the Indian population than its representation in the 1000 Genome Project (South Asians). We have created a database, ClinIndb, linked to the Leiden Open Variation Database, to help clinicians and researchers in diagnosis, counseling, and development of appropriate genetic screening tools relevant to the Indian populations and Indians living abroad.

Targeted exome sequencing in anti-factor H antibody negative HUS reveals multiple variations.

BACKGROUND: Genetic susceptibility to atypical hemolytic uremic syndrome (aHUS) may lie within genes regulating or activating the alternate complement and related pathways converging on endothelial cell activation. METHODS: We tested 32 Indian patients of aHUS negative for antibodies to complement factor H for genetic variations in a panel of 15 genes, i.e., CFH, CFHR1-5, CFI, CFB, C3, CD46, MASP2, DGKE, ADAMTS13, THBD and PLG using next-generation DNA sequencing and for copy number variation in CFHR1-3. RESULTS: Despite absence of a public database of exome variations in the Indian population and limited functional studies, we could establish a genetic diagnosis in 6 (18.8%) patients using a stringent scheme of prioritization. One patient carried a likely pathogenic variation. The number of patients carrying possibly pathogenic variation was as follows: 1 variation: 5 patients, 2 variations: 9 patients, 3 variations: 5 patients, 4 variations: 9 patients, 5 variations: 2 patients and 6 variations: 2 patients. Homozygous deletion of CFHR1-3 was present in five patients; none of these carried a diagnostic genetic variation. Patients with or without diagnostic variation did not differ significantly in terms of enrichment of genetic variations that were rare/novel or predicted deleterious, or for possible environmental triggers. CONCLUSION: We conclude that genetic testing for multiple genes in patients with aHUS negative for anti-FH antibodies reveals multiple candidate variations that require prioritization. Population data on variation frequency of the Indian population and supportive functional studies are likely to improve diagnostic yield.

Ayurgenomics for stratified medicine: TRISUTRA consortium initiative across ethnically and geographically diverse Indian populations.

BACKGROUND: Genetic differences in the target proteins, metabolizing enzymes and transporters that contribute to inter-individual differences in drug response are not integrated in contemporary drug development programs. Ayurveda, that has propelled many drug discovery programs albeit for the search of new chemical entities incorporates inter-individual variability "Prakriti" in development and administration of drug in an individualized manner. Prakriti of an individual largely determines responsiveness to external environment including drugs as well as susceptibility to diseases. Prakriti has also been shown to have molecular and genomic correlates. We highlight how integration of Prakriti concepts can augment the efficiency of drug discovery and development programs through a unique initiative of Ayurgenomics TRISUTRA consortium. METHODS: Five aspects that have been carried out are (1) analysis of variability in FDA approved pharmacogenomics genes/SNPs in exomes of 72 healthy individuals including predominant Prakriti types and matched controls from a North Indian Indo-European cohort (2) establishment of a consortium network and development of five genetically homogeneous cohorts from diverse ethnic and geo-climatic background (3) identification of parameters and development of uniform standard protocols for objective assessment of Prakriti types (4) development of protocols for Prakriti evaluation and its application in more than 7500 individuals in the five cohorts (5) Development of data and sample repository and integrative omics pipelines for identification of genomic correlates. RESULTS: Highlight of the study are (1) Exome sequencing revealed significant differences between Prakriti types in 28 SNPs of 11 FDA approved genes of pharmacogenomics relevance viz. CYP2C19, CYP2B6, ESR1, F2, PGR, HLA-B, HLA-DQA1, HLA-DRB1, LDLR, CFTR, CPS1. These variations are polymorphic in diverse Indian and world populations included in 1000 genomes project. (2) Based on the phenotypic attributes of Prakriti we identified anthropometry for anatomical features, biophysical parameters for skin types, HRV for autonomic function tests, spirometry for vital capacity and gustometry for taste thresholds as objective parameters. (3) Comparison of Prakriti phenotypes across different ethnic, age and gender groups led to identification of invariant features as well as some that require weighted considerations across the cohorts. CONCLUSION: Considering the molecular and genomics differences underlying Prakriti and relevance in disease pharmacogenomics studies, this novel integrative platform would help in identification of differently susceptible and drug responsive population. Additionally, integrated analysis of phenomic and genomic variations would not only allow identification of clinical and genomic markers of Prakriti for application in personalized medicine but also its integration in drug discovery and development programs.

Human brain harbors single nucleotide somatic variations in functionally relevant genes possibly mediated by oxidative stress.

Somatic variation in DNA can cause cells to deviate from the preordained genomic path in both disease and healthy conditions. Here, using exome sequencing of paired tissue samples, we show that the normal human brain harbors somatic single base variations measuring up to 0.48% of the total variations. Interestingly, about 64% of these somatic variations in the brain are expected to lead to non-synonymous changes, and as much as 87% of these represent G:C>T:A transversion events. Further, the transversion events in the brain were mostly found in the frontal cortex, whereas the corpus callosum from the same individuals harbors the reference genotype. We found a significantly higher amount of 8-OHdG (oxidative stress marker) in the frontal cortex compared to the corpus callosum of the same subjects (p<0.01), correlating with the higher G:C>T:A transversions in the cortex. We found significant enrichment for axon guidance and related pathways for genes harbouring somatic variations. This could represent either a directed selection of genetic variations in these pathways or increased susceptibility of some loci towards oxidative stress. This study highlights that oxidative stress possibly influence single nucleotide somatic variations in normal human brain.

Chromosome 22q11.2 microdeletion in monozygotic twins with discordant phenotype and deletion size.

We report on a pair of male monozygotic twins with 22q11.2 microdeletion, discordant phenotype and discordant deletion size. The second twin had findings suggestive of DiGeorge syndrome, while the first twin had milder anomalies without any cardiac malformation. The second twin had presented with intractable convulsion, cyanosis and cardiovascular failure in the fourth week of life and expired on the sixth week of life, whereas the first twin had some characteristic facial appearance with developmental delay but no other signs of the 22q11.2 microdeletion syndrome including cardiovascular malformation. The fluorescence in situ hybridization (FISH) analysis had shown a microdeletion on the chromosome 22q11.2 in both twins. The interphase FISH did not find any evidence for the mosaicism. The genomic DNA microarray analysis, using HumanCytoSNP-12 BeadChip (Illumina), was identical between the twins except different size of deletion of 22q11.2. The zygosity using HumanCytoSNP-12 BeadChip (Illumina) microarray analysis suggested monozygosity. This observation indicates that altered size of the deletion may be the underlying etiology for the discordance in phenotype in monozygotic twins. We think early post zygotic events (mitotic non-allelic homologous recombination) could have been played a role in the alteration of 22q11.2 deletion size and, thus phenotypic variability in the monozygotic twins.

Spectrum of large copy number variations in 26 diverse Indian populations: potential involvement in phenotypic diversity.

Copy number variations (CNVs) have provided a dynamic aspect to the apparently static human genome. We have analyzed CNVs larger than 100 kb in 477 healthy individuals from 26 diverse Indian populations of different linguistic, ethnic and geographic backgrounds. These CNVRs were identified using the Affymetrix 50K Xba 240 Array. We observed 1,425 and 1,337 CNVRs in the deletion and amplification sets, respectively, after pooling data from all the populations. More than 50% of the genes encompassed entirely in CNVs had both deletions and amplifications. There was wide variability across populations not only with respect to CNV extent (ranging from 0.04-1.14% of genome under deletion and 0.11-0.86% under amplification) but also in terms of functional enrichments of processes like keratinization, serine proteases and their inhibitors, cadherins, homeobox, olfactory receptors etc. These did not correlate with linguistic, ethnic, geographic backgrounds and size of populations. Certain processes were near exclusive to deletion (serine proteases, keratinization, olfactory receptors, GPCRs) or duplication (homeobox, serine protease inhibitors, embryonic limb morphogenesis) datasets. Populations having same enriched processes were observed to contain genes from different genomic loci. Comparison of polymorphic CNVRs (5% or more) with those cataloged in Database of Genomic Variants revealed that 78% (2473) of the genes in CNVRs in Indian populations are novel. Validation of CNVs using Sequenom MassARRAY revealed extensive heterogeneity in CNV boundaries. Exploration of CNV profiles in such diverse populations would provide a widely valuable resource for understanding diversity in phenotypes and disease.

Common variants of FTO and the risk of obesity and type 2 diabetes in Indians.

Common variants of fat mass and obesity-associated gene (FTO, fat mass- and obesity-associated gene) have been shown to be associated with obesity and type 2 diabetes in population of European and non-European ethnicity. However, studies in Indian population have provided inconsistent results. Here, we examined association of eight FTO variants (rs1421085, rs8050136, rs9939609, rs9930506, rs1861867, rs9926180, rs2540769 and rs708277) with obesity and type 2 diabetes in 5364 North Indians (2474 type 2 diabetes patients and 2890 non-diabetic controls) in two stages. None of the variants including previously reported intron 1 variants (rs1421085, rs8050136, rs9939609 and rs9930506) showed body mass index (BMI)-dependent/independent association with type 2 diabetes. However, rs1421085, rs8050136 and rs9939609 were associated with obesity status and measures of obesity (BMI, waist circumference and waist-to-hip ratio) in stage 2 and combined study population. Meta-analysis of the two study population results also revealed that rs1421085, rs8050136 and rs9939609 were significantly associated with BMI both under the random- and fixed-effect models (P (random/fixed)=0.02/0.0001, 0.004/0.0006 and 0.01/0.01, respectively). In conclusion, common variants of FTO were associated with obesity, but not with type 2 diabetes in North Indian population.

Genetic profile of patients with epilepsy on first-line antiepileptic drugs and potential directions for personalized treatment.

BACKGROUND: The first-line antiepileptic drugs, although affordable and effective in the control of seizures, are associated with adverse drug effects, and there is large interindividual variability in the appropriate dose at which patients respond favorably. This variability may partly be explained by functional consequences of genetic polymorphisms in the drug-metabolizing enzymes, such as the CYP450 family, microsomal epoxide hydrolase and UDP-glucuronosyltransferases, drug transporters, mainly ATP-binding cassette transporters, and drug targets, including sodium channels. The purpose of this study was to determine the allele and genotype frequencies of such genetic variants in patients with epilepsy from North India administered first-line antiepileptic drugs, such as phenobarbitone, phenytoin, carbamazepine and valproic acid, and compare them with worldwide epilepsy populations. MATERIALS & METHODS: SNP screening of 19 functional variants from 12 genes in 392 patients with epilepsy was carried out, and the patients were classified with respect to the metabolizing rate of their drug-metabolizing enzymes, efflux rate of drug transporters and sensitivity of drug targets. RESULTS: A total of 16 SNPs were found to be polymorphic, and the allelic frequencies for these SNPs were in conformance with Hardy-Weinberg equilibrium. Among all the polymorphisms studied, functional variants from genes encoding CYP2C19, EPHX1, ABCB1 and SCN1A were highly polymorphic in North Indian epilepsy patients, and might account for differential drug response to first-line antiepileptic drugs. CONCLUSION: Interethnic differences were elucidated for several polymorphisms that might be responsible for differential serum drug levels and optimal dose requirement for efficacious treatment.

Genetic determinants of immune-response to a polysaccharide vaccine for typhoid.

UNLABELLED: Differences in immunological response among vaccine recipients are determined both by their genetic differences and environmental factors. Knowledge of genetic determinants of immunological response to a vaccine can be used to design a vaccine that circumvents immunogenetic restrictions. The currently available vaccine for typhoid is a pure polysaccharide vaccine, immune response to which is T-cell independent. Little is known about whether genetic variation among vaccinees associates with variation in their antibody response to a polysaccharide vaccine. We conducted a study on 1,000 individuals resident in an area at high-risk for typhoid; vaccinated them with the typhoid vaccine, measured their antibody response to the vaccine, assayed >2,000 curated SNPs chosen from 283 genes that are known to participate in immune-response; and analyzed these data using a strategy to (a) minimize the statistical problems associated with testing of multiple hypotheses, and (b) internally cross-validate inferences, using a half-sample design, with little loss of statistical power. The first stage analysis, using the first half-sample, identified 54 SNPs in 43 genes to be significantly associated with immune response. In the second-stage, these inferences were cross-validated using the second half-sample. First-stage results of only 8 SNPs (out of 54) in 7 genes (out of 43) were cross-validated. We tested additional SNPs in these 7 genes, and found 8 more SNPs to be significantly associated. Haplotypes constructed with these SNPs in these 7 genes also showed significant association. These 7 genes are DEFB1, TLR1, IL1RL1, CTLA4, MAPK8, CD86 and IL17D. The overall picture that has emerged from this study is that (a) immune response to polysaccharide antigens is qualitatively different from that to protein antigens, and (b) polymorphisms in genes involved in polysaccharide recognition, signal transduction, inhibition of T-cell proliferation, pro-inflammatory signaling and eventual production of antimicrobial peptides are associated with antibody response to the polysaccharide vaccine for typhoid. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11568-010-9134-1) contains supplementary material, which is available to authorized users.