Extreme constitution "Prakriti" types of Ayurveda exhibit systemic physiological attributes. Our earlier genetic study has revealed differences in EGLN1, key modulator of hypoxia axis between Prakriti types. This was associated with differences in high altitude adaptation and susceptibility to high altitude pulmonary edema (HAPE). In this study we investigate other molecular differences that contribute to systemic attributes of Prakriti that would be relevant in predictive marker discovery.

Non-coding RNAs from transposable elements of human genome are gaining prominence in modulating transcriptome dynamics. Alu elements, as exonized, edited and antisense components within same transcripts could create novel regulatory switches in response to different transcriptional cues. We provide the first evidence for co-occurrences of these events at transcriptome-wide scale through integrative analysis of data sets across diverse experimental platforms and tissues. This involved the following: (i) positional anchoring of Alu exonization events in the UTRs and CDS of 4663 transcript isoforms from RefSeq mRNAs and (ii) mapping on to them A→I editing events inferred from ∼7 million ESTs from dbEST and antisense transcripts identified from virtual serial analysis of gene expression tags represented in Cancer Genome Anatomy Project next-generation sequencing data sets across 20 tissues. We observed significant enrichment of these events in the 3'UTR as well as positional preference within the embedded Alus. More than 300 genes had co-occurrence of all these events at the exon level and were significantly enriched in apoptosis and lysosomal processes. Further, we demonstrate functional evidence of such dynamic interactions between Alu-mediated events in a time series data from Integrated Personal Omics Profiling during recovery from a viral infection. Such 'single transcript-multiple fate' opportunity facilitated by Alu elements may modulate transcriptional response, especially during stress.

The Indian Genome Variation Consortium (IGVC) project, an initiative of the Council for Scientific and Industrial Research, has been the first large-scale comprehensive study of the Indian population. One of the major aims of the project is to study and catalog the variations in nearly thousand candidate genes related to diseases and drug response for predictive marker discovery, founder identification and also to address questions related to ethnic diversity, migrations, extent and relatedness with other world population. The Phase I of the project aimed at providing a set of reference populations that would represent the entire genetic spectrum of India in terms of language, ethnicity and geography and Phase II in providing variation data on candidate genes and genome wide neutral markers on these reference set of populations. We report here development of the IGVBrowser that provides allele and genotype frequency data generated in the IGVC project. The database harbors 4229 SNPs from more than 900 candidate genes in contrasting Indian populations. Analysis shows that most of the markers are from genic regions. Further, a large fraction of genes are implicated in cardiovascular, metabolic, cancer and immune system-related diseases. Thus, the IGVC data provide a basal level variation data in Indian population to study genetic diseases and pharmacology. Additionally, it also houses data on ∼50,000 (Affy 50 K array) genome wide neutral markers in these reference populations. In IGVBrowser one can analyze and compare genomic variations in Indian population with those reported in HapMap along with annotation information from various primary data sources. Database URL: http://igvbrowser.igib.res.in.

Human mitochondrial DNA (mtDNA) variations have been implicated in a broad spectrum of diseases. With over 3000 mtDNA variations reported across databases, establishing pathogenicity of variations in mtDNA is a major challenge. We have designed and developed a comprehensive weighted scoring system (MtSNPscore) for identification of mtDNA variations that can impact pathogenicity and would likely be associated with disease. The criteria for pathogenicity include information available in the literature, predictions made by various in silico tools and frequency of variation in normal and patient datasets. The scoring scheme also assigns scores to patients and normal individuals to estimate the cumulative impact of variations. The method has been implemented in an automated pipeline and has been tested on Indian ataxia dataset (92 individuals), sequenced in this study, and other publicly available mtSNP dataset comprising of 576 mitochondrial genomes of Japanese individuals from six different groups, namely, patients with Parkinson's disease, patients with Alzheimer's disease, young obese males, young non-obese males, and type-2 diabetes patients with or without severe vascular involvement. MtSNPscore, for analysis can extract information from variation data or from mitochondrial DNA sequences. It has a web-interface http://bioinformatics.ccmb.res.in/cgi-bin/snpscore/Mtsnpscore.pl that provides flexibility to update/modify the parameters for estimating pathogenicity.

The recent release of SARS-CoV-2 genomic data from several countries has provided clues into the potential antigenic drift of the coronavirus population. In particular, the genomic instability observed in the spike protein necessitates immediate action and further exploration in the context of viral-host interactions. By temporally tracking 527,988 SARS-CoV-2 genomes, we identified invariant and hypervariable regions within the spike protein. We evaluated combination of mutations from SARS-CoV-2 lineages and found that maximum number of lineage-defining mutations were present in the N-terminal domain (NTD). Based on mutant 3D-structural models of known Variants of Concern (VOCs), we found that structural properties such as accessibility, secondary structural type, and intra-protein interactions at local mutation sites are greatly altered. Further, we observed significant differences between intra-protein networks of wild-type and Delta mutant, with the latter showing dense intra-protein contacts. Extensive molecular dynamics simulations of D614G mutant spike structure with hACE2 further revealed dynamic features with 47.7% of mutations mapping on flexible regions of spike protein. Thus, we propose that significant changes within spike protein structure have occurred that may impact SARS-CoV-2 pathogenesis, and repositioning of vaccine candidates is required to contain the spread of COVID-19 pathogen.

Semi-autonomous functioning of mitochondria in eukaryotic cell necessitates coordination with nucleus. Several RNA species fine-tune mitochondrial processes by synchronizing with the nuclear program, however the involved components remain enigmatic. In this study, we identify a widely conserved dually localized protein Myg1, and establish its role as a 3'-5' RNA exonuclease. We employ mouse melanoma cells, and knockout of the Myg1 ortholog in Saccharomyces cerevisiae with complementation using human Myg1 to decipher the conserved role of Myg1 in selective RNA processing. Localization of Myg1 to nucleolus and mitochondrial matrix was studied through imaging and confirmed by sub-cellular fractionation studies. We developed Silexoseqencing, a methodology to map the RNAse trail at single-nucleotide resolution, and identified in situ cleavage by Myg1 on specific transcripts in the two organelles. In nucleolus, Myg1 processes pre-ribosomal RNA involved in ribosome assembly and alters cytoplasmic translation. In mitochondrial matrix, Myg1 processes 3'-termini of the mito-ribosomal and messenger RNAs and controls translation of mitochondrial proteins. We provide a molecular link to the possible involvement of Myg1 in chronic depigmenting disorder vitiligo. Our study identifies a key component involved in regulating spatially segregated organellar RNA processing and establishes the evolutionarily conserved ribonuclease as a coordinator of nucleo-mitochondrial crosstalk.

Admixture mapping has been enormously resourceful in identifying genetic variations linked to phenotypes, adaptation, and diseases. In this study through analysis of copy number variable regions (CNVRs), we report extensive restructuring in the genomes of the recently admixed African-Indian population (OG-W-IP) that inhabits a highly saline environment in Western India. The study included subjects from OG-W-IP (OG), five different Indian and three HapMap populations that were genotyped using Affymetrix version 6.0 arrays. Copy number variations (CNVs) detected using Birdsuite were used to define CNVRs. Population structure with respect to CNVRs was delineated using random forest approach. OG genomes have a surprising excess of CNVs in comparison to other studied populations. Individual ancestry proportions computed using STRUCTURE also reveals a unique genetic component in OGs. Population structure analysis with CNV genotypes indicates OG to be distant from both the African and Indian ancestral populations. Interestingly, it shows genetic proximity with respect to CNVs to only one Indian population IE-W-LP4, which also happens to reside in the same geographical region. We also observe a significant enrichment of molecular processes related to ion binding and receptor activity in genes encompassing OG-specific CNVRs. Our results suggest that retention of CNVRs from ancestral natives and de novo acquisition of CNVRs could accelerate the process of adaptation especially in an extreme environment. Additionally, this population would be enormously useful for dissecting genes and delineating the involvement of CNVs in salt adaptation.

Dyslexia is a heritable neurodevelopmental disorder characterized by difficulties in reading and writing. In this study, we describe the identification of a set of 17 polymorphisms located across 1.9Mb region on chromosome 5q31.3, encompassing genes of the PCDHG cluster, TAF7, PCDH1 and ARHGAP26, dominantly inherited with dyslexia in a multi-incident family. Strikingly, the non-risk form of seven variations of the PCDHG cluster, are preponderant in the human lineage, while risk alleles are ancestral and conserved across Neanderthals to non-human primates. Four of these seven ancestral variations (c.460A>C [p.Ile154Leu], c.541G>A [p.Ala181Thr], c.2036G>C [p.Arg679Pro] and c.2059A>G [p.Lys687Glu]) result in amino acid alterations. p.Ile154Leu and p.Ala181Thr are present at EC2: EC3 interacting interface of γA3-PCDH and γA4-PCDH respectively might affect trans-homophilic interaction and hence neuronal connectivity. p.Arg679Pro and p.Lys687Glu are present within the linker region connecting trans-membrane to extracellular domain. Sequence analysis indicated the importance of p.Ile154, p.Arg679 and p.Lys687 in maintaining class specificity. Thus the observed association of PCDHG genes encoding neural adhesion proteins reinforces the hypothesis of aberrant neuronal connectivity in the pathophysiology of dyslexia. Additionally, the striking conservation of the identified variants indicates a role of PCDHG in the evolution of highly specialized cognitive skills critical to reading.

During the course of the COVID-19 pandemic, large-scale genome sequencing of SARS-CoV-2 has been useful in tracking its spread and in identifying variants of concern (VOC). Viral and host factors could contribute to variability within a host that can be captured in next-generation sequencing reads as intra-host single nucleotide variations (iSNVs). Analysing 1347 samples collected till June 2020, we recorded 16 410 iSNV sites throughout the SARS-CoV-2 genome. We found ∼42% of the iSNV sites to be reported as SNVs by 30 September 2020 in consensus sequences submitted to GISAID, which increased to ∼80% by 30th June 2021. Following this, analysis of another set of 1774 samples sequenced in India between November 2020 and May 2021 revealed that majority of the Delta (B.1.617.2) and Kappa (B.1.617.1) lineage-defining variations appeared as iSNVs before getting fixed in the population. Besides, mutations in RdRp as well as RNA-editing by APOBEC and ADAR deaminases seem to contribute to the differential prevalence of iSNVs in hosts. We also observe hyper-variability at functionally critical residues in Spike protein that could alter the antigenicity and may contribute to immune escape. Thus, tracking and functional annotation of iSNVs in ongoing genome surveillance programs could be important for early identification of potential variants of concern and actionable interventions.

Viral infections have a history of abrupt and severe eruptions through the years in the form of pandemics. And yet, definitive therapies or preventive measures are not present. Herbal medicines have been a source of various antiviral compounds such as Oseltamivir, extracted using shikimic acid from star anise (Illicium verum) and Acyclovir from Carissa edulis are FDA (Food and Drug Administration) approved antiviral drugs. In this study, we dissect the anti-coronavirus infection activity of Cissampelos pareira L (Cipa) extract using an integrative approach.

Autonomic modulation is critical during various physiological activities, including orthostatic stimuli and primarily evaluated by heart rate variability (HRV). Orthostatic stress affects people differently suggesting the possibility of identification of predisposed groups to autonomic dysfunction-related disorders in a healthy state. One way to understand this kind of variability is by using Ayurvedic approach that classifies healthy individuals into Prakriti types based on clinical phenotypes. To this end, we explored the differential response to orthostatic stress in different Prakriti types using HRV. HRV was measured in 379 subjects(Vata = 97, Pitta = 68, Kapha = 68, and Mixed Prakriti = 146) from two geographical regions(Vadu and Delhi NCR) for 5 min supine (baseline), 3 min head-up-tilt (HUT) at 60°, and 5 min resupine. We observed that Kapha group had lower baseline HRV than other two groups, although not statistically significant. The relative change (%Δ1&2 ) in various HRV parameters in response to HUT was although minimal in Kapha group. Kapha also had significantly lower change in HR, LF (nu), HF (nu), and LF/HF than Pitta in response to HUT. The relative change (%Δ1 ) in HR and parasympathetic parameters (RMSSD, HF, SD1) was significantly greater in the Vata than in the Kapha. Thus, the low baseline and lower response to HUT in Kapha and the maximum drop in parasympathetic activity of Vata may indicate a predisposition to early autonomic dysfunction and associated conditions. It emphasizes the critical role of Prakriti-based phenotyping in stratifying the differential responses of cardiac autonomic modulation in various postures among healthy individuals across different populations.

Spinocerebellar ataxia type-12 (SCA12) is a neurological disorder caused due to triplet (CAG) repeat expansion in 5' UTR of PPP2R2B. It is one of the most prominent SCA-subtype in Indian population and till date no patient specific models have been described. Human-induced-pluripotent-stem cell (HiPSC) based disease modelling has become the next generation tool for studying various human pathologies. In the present study we established three SCA12 patient specific iPSC lines. All the generated lines have shown pluripotency markers, normal karyotype, in-vitro three germ layers differentiation potential, vector clearance, SCA12 mutation, parental genomic identity and contamination free culture.

Primate-specific Alus harbor different regulatory features, including miRNA targets. In this study, we provide evidence for miRNA-mediated modulation of transcript isoform levels during heat-shock response through exaptation of Alu-miRNA sites in mature mRNA. We performed genome-wide expression profiling coupled with functional validation of miRNA target sites within exonized Alus, and analyzed conservation of these targets across primates. We observed that two miRNAs (miR-15a-3p and miR-302d-3p) elevated in stress response, target RAD1, GTSE1, NR2C1, FKBP9 and UBE2I exclusively within Alu. These genes map onto the p53 regulatory network. Ectopic overexpression of miR-15a-3p downregulates GTSE1 and RAD1 at the protein level and enhances cell survival. This Alu-mediated fine-tuning seems to be unique to humans as evident from the absence of orthologous sites in other primate lineages. We further analyzed signatures of selection on Alu-miRNA targets in the genome, using 1000 Genomes Phase-I data. We found that 198 out of 3177 Alu-exonized genes exhibit signatures of selection within Alu-miRNA sites, with 60 of them containing SNPs supported by multiple evidences (global-FST > 0.3, pair-wise-FST > 0.5, Fay-Wu's H < -20, iHS > 2.0, high ΔDAF) and implicated in p53 network. We propose that by affecting multiple genes, Alu-miRNA interactions have the potential to facilitate population-level adaptations in response to environmental challenges.

Identification and study of genetic variation in recently admixed populations not only provides insight into historical population events but also is a powerful approach for mapping disease loci. We studied a population (OG-W-IP) that is of African-Indian origin and has resided in the western part of India for 500 years; members of this population are believed to be descendants of the Bantu-speaking population of Africa. We have carried out this study by using a set of 18,534 autosomal markers common between Indian, CEPH-HGDP, and HapMap populations. Principal-components analysis clearly revealed that the African-Indian population derives its ancestry from Bantu-speaking west-African as well as Indo-European-speaking north and northwest Indian population(s). STRUCTURE and ADMIXTURE analyses show that, overall, the OG-W-IPs derive 58.7% of their genomic ancestry from their African past and have very little inter-individual ancestry variation (8.4%). The extent of linkage disequilibrium also reveals that the admixture event has been recent. Functional annotation of genes encompassing the ancestry-informative markers that are closer in allele frequency to the Indian ancestral population revealed significant enrichment of biological processes, such as ion-channel activity, and cadherins. We briefly examine the implications of determining the genetic diversity of this population, which could provide opportunities for studies involving admixture mapping.

Heterogeneity amidst healthy individuals at genomic level is being widely acknowledged. This, in turn, is modulated by differential response to environmental cues and treatment regimens, necessitating the need for stratified/personalized therapy. We intend to understand the molecular determinants of Ayurvedic way (ancient Indian system of medicine) of endo-phenotyping individuals into distinct constitution types termed "Prakriti," which forms the basis of personalized treatment. In this study, we explored and analyzed the healthy human gut microbiome structure within three predominant Prakriti groups from a genetically homogenous cohort to discover differentially abundant taxa, using 16S rRNA gene based microbial community profiling. We found Bacteroidetes and Firmicutes as major gut microbial components in varying composition, albeit with similar trend across Prakriti. Multiple species of the core microbiome showed differential abundance within Prakriti types, with gender specific signature taxons. Our study reveals that despite overall uniform composition of gut microbial community, healthy individuals belonging to different Prakriti groups have enrichment of specific bacteria. It highlights the importance of Prakriti based endo-phenotypes to explain the variability amongst healthy individuals in gut microbial flora that have important consequences for an individual's health, disease and treatment.

Alu repeats contribute to phylogenetic novelties in conserved regulatory networks in primates. Our study highlights how exonized Alus could nucleate large-scale mRNA-miRNA interactions. Using a functional genomics approach, we characterize a transcript isoform of an orphan gene, CYP20A1 (CYP20A1_Alu-LT) that has exonization of 23 Alus in its 3'UTR. CYP20A1_Alu-LT, confirmed by 3'RACE, is an outlier in length (9 kb 3'UTR) and widely expressed. Using publically available data sets, we demonstrate its expression in higher primates and presence in single nucleus RNA-seq of 15,928 human cortical neurons. miRanda predicts ∼4,700 miRNA recognition elements (MREs) for ∼1,000 miRNAs, primarily originated within these 3'UTR-Alus. CYP20A1_Alu-LT could be a potential multi-miRNA sponge as it harbors ≥10 MREs for 140 miRNAs and has cytosolic localization. We further tested whether expression of CYP20A1_Alu-LT correlates with mRNAs harboring similar MRE targets. RNA-seq with conjoint miRNA-seq analysis was done in primary human neurons where we observed CYP20A1_Alu-LT to be downregulated during heat shock response and upregulated in HIV1-Tat treatment. In total, 380 genes were positively correlated with its expression (significantly downregulated in heat shock and upregulated in Tat) and they harbored MREs for nine expressed miRNAs which were also enriched in CYP20A1_Alu-LT. MREs were significantly enriched in these 380 genes compared with random sets of differentially expressed genes (P = 8.134e-12). Gene ontology suggested involvement of these genes in neuronal development and hemostasis pathways thus proposing a novel component of Alu-miRNA-mediated transcriptional modulation that could govern specific physiological outcomes in higher primates.

Bioactive fractions obtained from medicinal plants which have been used for the treatment of multiple diseases could exert their effects by targeting common pathways. Prior knowledge of their usage could allow us to identify novel molecular links. In this study, we explored the molecular basis of action of one such herbal formulation Cissampelos pareira L. (Cipa), used for the treatment of female hormone disorders and fever. Transcriptomic studies on MCF7 cell lines treated with Cipa extract carried out using Affymetrix arrays revealed a downregulation of signatures of estrogen response potentially modulated through estrogen receptor α (ERα). Molecular docking analysis identified 38 Cipa constituents that potentially bind (ΔG <  - 7.5) with ERα at the same site as estrogen. The expression signatures in the connectivity map ( https://clue.io/; ) revealed high positive scores with translation inhibitors such as emetine (score: 99.61) and knockdown signatures of genes linked to the antiviral response such as ribosomal protein RPL7 (score: 99.92), which is a reported ERα coactivator. Further, gene knockdown experiments revealed that Cipa exhibits antiviral activity in dengue infected MCF7 cells potentially modulated through estrogen receptor 1. This approach reveals a novel pathway involving the ESR1-RPL7 axis which could be a potential target in dengue viral infection.

With the advent of next-generation sequencing, large-scale initiatives for mining whole genomes and exomes have been employed to better understand global or population-level genetic architecture. India encompasses more than 17% of the world population with extensive genetic diversity, but is under-represented in the global sequencing datasets. This gave us the impetus to perform and analyze the whole genome sequencing of 1029 healthy Indian individuals under the pilot phase of the 'IndiGen' program. We generated a compendium of 55,898,122 single allelic genetic variants from geographically distinct Indian genomes and calculated the allele frequency, allele count, allele number, along with the number of heterozygous or homozygous individuals. In the present study, these variants were systematically annotated using publicly available population databases and can be accessed through a browsable online database named as 'IndiGenomes' http://clingen.igib.res.in/indigen/. The IndiGenomes database will help clinicians and researchers in exploring the genetic component underlying medical conditions. Till date, this is the most comprehensive genetic variant resource for the Indian population and is made freely available for academic utility. The resource has also been accessed extensively by the worldwide community since it's launch.

COVID-19 pneumonia has been associated with severe acute hypoxia, sepsis-like states, thrombosis and chronic sequelae including persisting hypoxia and fibrosis. The molecular hypoxia response pathway has been associated with such pathologies and our recent observations on anti-hypoxic and anti-inflammatory effects of whole aqueous extract of Adhatoda Vasica (AV) prompted us to explore its effects on relevant preclinical mouse models.

Within-population genetic diversity is greatest within Africa, while between-population genetic diversity is directly proportional to geographic distance. The most divergent contemporary human populations include the click-speaking forager peoples of southern Africa, broadly defined as Khoesan. Both intra- (Bantu expansion) and inter-continental migration (European-driven colonization) have resulted in complex patterns of admixture between ancient geographically isolated Khoesan and more recently diverged populations. Using gender-specific analysis and almost 1 million autosomal markers, we determine the significance of estimated ancestral contributions that have shaped five contemporary southern African populations in a cohort of 103 individuals. Limited by lack of available data for homogenous Khoesan representation, we identify the Ju/'hoan (n = 19) as a distinct early diverging human lineage with little to no significant non-Khoesan contribution. In contrast to the Ju/'hoan, we identify ancient signatures of Khoesan and Bantu unions resulting in significant Khoesan- and Bantu-derived contributions to the Southern Bantu amaXhosa (n = 15) and Khoesan !Xun (n = 14), respectively. Our data further suggests that contemporary !Xun represent distinct Khoesan prehistories. Khoesan assimilation with European settlement at the most southern tip of Africa resulted in significant ancestral Khoesan contributions to the Coloured (n = 25) and Baster (n = 30) populations. The latter populations were further impacted by 170 years of East Indian slave trade and intra-continental migrations resulting in a complex pattern of genetic variation (admixture). The populations of southern Africa provide a unique opportunity to investigate the genomic variability from some of the oldest human lineages to the implications of complex admixture patterns including ancient and recently diverged human lineages.