Multiple allelic associations from genes involved in energy metabolism were identified in celiac disease.

Energy metabolism is a critical factor that influences disease pathogenesis. Recent high-throughput genomic studies have enabled us to look into disease biology with greater details. Celiac disease (CD) is an inflammatory autoimmune disease where ~60 non-HLA genes were identified which in conjunction with HLA genes explain ~55% of the disease heritability. In this study we aimed to identify susceptibility energy metabolism genes and investigate their role in CD. We re-analysed published Immunochip genotyping data, which were originally analysed for CD association studies in north Indian and Dutch population. 269 energy metabolism genes were tested. Meta-analysis was done for the identified SNPs. To validate the functional implications of identified markers and/or genes, in silico functional annotation was performed. Six SNPs were identified in north Indians, of which three markers from two loci were replicated in Dutch. rs2071592 (PMeta=5.01e-75) and rs2251824 (PMeta=1.87e-14) from ATP6V1G2-NFKBIL1-DDX39B locus and rs4947331 (PMeta= 9.85e-13) from NEU1 locus were found significantly associated. Identified genes are key regulators of cellular energy metabolism and associated with several immune mediated diseases. In silico functional annotation showed significant biological relevance of these novel markers and genes. FDI approved therapeutics against ATP6V1G2 and NEU1 are currently in use to treat chronic and inflammatory diseases. This study identified two pathogenic loci, originally involved in energy metabolism. Extensive investigation showed their synergistic role in CD pathogenesis by promoting immune mediated enteric inflammation. Proposed CD pathogenesis model in this study needs to be tested through tissue-on-chip and in vivo methods to ensure its translational application.

Contributions of human ACE2 and TMPRSS2 in determining host-pathogen interaction of COVID-19.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection is at present an emerging global public health crisis. Angiotensin converting enzyme 2 (ACE2) and trans-membrane protease serine 2 (TMPRSS2) are the two major host factors that contribute to the virulence of SARS-CoV-2 and pathogenesis of coronavirus disease-19 (COVID-19). Transmission of SARS-CoV-2 from animal to human is considered a rare event that necessarily requires strong evolutionary adaptations. Till date no other human cellular receptors are identified beside ACE2 for SARS-CoV-2 entry inside the human cell. Proteolytic cleavage of viral spike (S)-protein and ACE2 by TMPRSS2 began the entire host-pathogen interaction initiated with the physical binding of ACE2 to S-protein. SARS-CoV-2 S-protein binds to ACE2 with much higher affinity and stability than that of SARS-CoVs. Molecular interactions between ACE2-S and TMPRSS2-S are crucial and preciously mediated by specific residues. Structural stability, binding affinity and level of expression of these three interacting proteins are key susceptibility factors for COVID-19. Specific protein-protein interactions (PPI) are being identified that explains uniqueness of SARS-CoV-2 infection. Amino acid substitutions due to naturally occurring genetic polymorphisms potentially alter these PPIs and poses further clinical heterogeneity of COVID-19. Repurposing of several phytochemicals and approved drugs against ACE2, TMPRSS2 and S-protein have been proposed that could inhibit PPI between them. We have also identified some novel lead phytochemicals present in Azadirachta indica and Aloe barbadensis which could be utilized for further in vitro and in vivo anti-COVID-19 drug discovery. Uncovering details of ACE2-S and TMPRSS2-S interactions would further contribute to future research on COVID-19.

Assessment of risk conferred by coding and regulatory variations of TMPRSS2 and CD26 in susceptibility to SARS-CoV-2 infection in human.

At present, more than 200 countries and territories are directly affected by the coronavirus disease-19 (COVID-19) pandemic. Incidence and case fatality rate are significantly higher among elderly individuals (age>60 years), type 2 diabetes and hypertension patients. Cellular receptor ACE2, serine protease TMPRSS2 and exopeptidase CD26 (also known as DPP4) are the three membrane bound proteins potentially implicated in SARS-CoV-2 infection. We hypothesised that common variants from TMPRSS2 and CD26 may play critical role in infection susceptibility of predisposed population or group of individuals. Coding (missense) and regulatory variants from TMPRSS2 and CD26 were studied across 26 global populations. Two missense and five regulatory SNPs were identified to have differential allelic frequency. Significant linkage disequilibrium (LD) signature was observed in different populations. Modelled protein-protein interaction (PPI) predicted strong molecular interaction between these two receptors and SARS-CoV-2 spike protein (S1 domain). However, two missense SNPs, rs12329760 (TMPRSS2) and rs1129599 (CD26), were not found to be involved physically in the said interaction. Four regulatory variants (rs112657409, rs11910678, rs77675406 and rs713400) from TMPRSS2 were found to influence the expression of TMPRSS2 and pathologically relevant MX1. rs13015258 a 50 UTR variant from CD26 have significant role in regulation of expression of key regulatory genes that could be involved in SARS-CoV-2 internalization. Overexpression of CD26 through epigenetic modification at rs13015258-C allele was found critical and could explain the higher SARS-CoV-2 infected fatality rate among type 2 diabetes.

Association study identified biologically relevant receptor genes with synergistic functions in celiac disease.

Receptors are essential mediators of cellular physiology, which facilitate molecular and cellular cross-talk with the environment. Nearly 20% of the all known celiac disease (CD) genes are receptors by function. We hypothesized that novel biologically relevant susceptibility receptor genes act in synergy in CD pathogenesis. We attempted to identify novel receptor genes in CD by re-analyzing published Illumina Immunochip dense genotype data for a north Indian and  a European (Dutch) cohort. North Indian dataset was screened for 269 known receptor genes. Association statistics for SNPs were considered with minor allele frequency >15% and association P ≤ 0.005 to attend desired study power. Identified markers were tested for cross-ethnic replication in a European CD dataset. Markers were analyzed in-silico to explain their functional significance in CD. Six novel SNPs from MOG (rs29231, p = 1.21e-11), GABBR1 (rs3025643, p = 1.60e-7), OR2H2 (rs1233388, p = 0.0002), ABCF1 (rs9262119, p = 0.0005), ADRA1A (rs10102024, p = 0.003), and ACVR2A (rs7560426, p = 0.004) were identified in north Indians, of which three genes namely, GABBR1 (rs3025643, p = 5.38e-8), OR2H2 (rs1233388, p = 3.29e-5) and ABCF1 (rs9262119, p = 0.0002) were replicated in Dutch. Tissue specific functional annotation, potential epigenetic regulation, co-expression, protein-protein interaction and pathway enrichment analyses indicated differential expression and synergistic function of key genes that could alter cellular homeostasis, ubiquitination mediated phagosome pathway and cellular protein processing to contribute for CD. At present multiple therapeutic compounds/drugs are available targeting GABBR1 and ADRA1A, which could be tested for their effectiveness against CD in controlled drug trials.