ILRUN, a Human Plasma Lipid GWAS Locus, Regulates Lipoprotein Metabolism in Mice

Liver eQTL meta-analysis illuminates potential molecular mechanisms of cardiometabolic traits

Understanding the molecular mechanisms of complex traits is essential for developing targeted interventions. We analyzed liver expression quantitative-trait locus (eQTL) meta-analysis data on 1,183 participants to identify conditionally distinct signals. We found 9,013 eQTL signals for 6,564 genes; 23% of eGenes had two signals, and 6% had three or more signals. We then integrated the eQTL results with data from 29 cardiometabolic genome-wide association study (GWAS) traits and identified 1,582 GWAS-eQTL colocalizations for 747 eGenes. Non-primary eQTL signals accounted for 17% of all colocalizations. Isolating signals by conditional analysis prior to coloc resulted in 37% more colocalizations than using marginal eQTL and GWAS data, highlighting the importance of signal isolation. Isolating signals also led to stronger evidence of colocalization: among 343 eQTL-GWAS signal pairs in multi-signal regions, analyses that isolated the signals of interest resulted in higher posterior probability of colocalization for 41% of tests. Leveraging allelic heterogeneity, we predicted causal effects of gene expression on liver traits for four genes. To predict functional variants and regulatory elements, we colocalized eQTL with liver chromatin accessibility QTL (caQTL) and found 391 colocalizations, including 73 with non-primary eQTL signals and 60 eQTL signals that colocalized with both a caQTL and a GWAS signal. Finally, we used publicly available massively parallel reporter assays in HepG2 to highlight 14 eQTL signals that include at least one expression-modulating variant. This multi-faceted approach to unraveling the genetic underpinnings of liver-related traits could lead to therapeutic development.

Genome-wide search identified DNA methylation sites that regulate the metabolome

Background: Identifying DNA methylation sites that regulate the metabolome is important for several purposes. In this study, publicly available GWAS data were integrated to find methylation sites that impact metabolome through a discovery and replication scheme and by using Mendelian randomization. Results: The outcome of analyses revealed 107 methylation sites associated with 84 metabolites at the genome-wide significance level (p<5e^-8) at both the discovery and replication stages. A large percentage of the observed associations (85%) were with lipids, significantly higher than expected (p = 0.0003). A number of CpG (methylation) sites showed specificity e.g., cg20133200 within PFKP was associated with glucose only and cg10760299 within GATM impacted the level of creatinine; in contrast, there were sites associated with numerous metabolites e.g., cg20102877 on the 2p23.3 region was associated with 39 metabolites. Integrating transcriptome data enabled identifying genes (N = 82) mediating the impact of methylation sites on the metabolome and cardiometabolic traits. For example, PABPC4 mediated the impact of cg15123755-HDL on type-2 diabetes. KCNK7 mediated the impact of cg21033440-lipids on hypertension. POC5, ILRUN, FDFT1, and NEIL2 mediated the impact of CpG sites on obesity through metabolic pathways. Conclusion: This study provides a catalog of DNA methylation sites that regulate the metabolome for downstream applications.

Identification of a novel role for the immunomodulator ILRUN in the development of several T cell subsets in mice

Inflammation and lipid regulator with UBA-like and NBR1-like domains (ILRUN) is a protein-encoding gene associated with innate immune signaling, lipid metabolism and cancer. In the context of innate immunity, ILRUN inhibits IRF3-mediated transcription of antimicrobial and proinflammatory cytokines by inducing degradation of the transcriptional coactivators CBP and p300. There remains a paucity of information, however, regarding the innate immune roles of ILRUN beyond in vitro analyses. To address this, we utilize a knockout mouse model to investigate the effect of ILRUN on cytokine expression in splenocytes and on the development of immune cell populations in the spleen and thymus. We show elevated production of tumor necrosis factor and interleukin-6 cytokines in ILRUN-deficient splenocytes following stimulation with the innate immune ligands polyinosinic:polycytidylic acid or lipopolysaccharide. Differences were also observed in the populations of several T cell subsets, including regulatory, mucosal-associated invariant and natural killer. These data identify novel functions for ILRUN in the development of certain immune cell populations and support previous in vitro findings that ILRUN negatively regulates the synthesis of pathogen-stimulated cytokines. This establishes the ILRUN knockout mouse model as a valuable resource for further study of the functions of ILRUN in health and disease.

Genetic Mapping of Multiple Traits Identifies Novel Genes for Adiposity, Lipids, and Insulin Secretory Capacity in Outbred Rats

Despite the successes of human genome-wide association studies, the causal genes underlying most metabolic traits remain unclear. We used outbred heterogeneous stock (HS) rats, coupled with expression data and mediation analysis, to identify quantitative trait loci (QTLs) and candidate gene mediators for adiposity, glucose tolerance, serum lipids, and other metabolic traits. Physiological traits were measured in 1,519 male HS rats, with liver and adipose transcriptomes measured in >410 rats. Genotypes were imputed from low-coverage whole-genome sequencing. Linear mixed models were used to detect physiological and expression QTLs (pQTLs and eQTLs, respectively), using both single nucleotide polymorphism (SNP)- and haplotype-based models for pQTL mapping. Genes with cis-eQTLs that overlapped pQTLs were assessed as causal candidates through mediation analysis. We identified 14 SNP-based pQTLs and 19 haplotype-based pQTLs, of which 10 were in common. Using mediation, we identified the following genes as candidate mediators of pQTLs: Grk5 for fat pad weight and serum triglyceride pQTLs on Chr1, Krtcap3 for fat pad weight and serum triglyceride pQTLs on Chr6, Ilrun for a fat pad weight pQTL on Chr20, and Rfx6 for a whole pancreatic insulin content pQTL on Chr20. Furthermore, we verified Grk5 and Ktrcap3 using gene knockdown/out models, thereby shedding light on novel regulators of obesity.

ILRUN Promotes Atherosclerosis Through Lipid-Dependent and Lipid-Independent Factors

BACKGROUND

Common genetic variation in close proximity to the ILRUN gene are significantly associated with coronary artery disease as well as with plasma lipid traits. We recently demonstrated that hepatic inflammation and lipid regulator with ubiquitin-associated domain-like and NBR1-like domains (ILRUN) regulates lipoprotein metabolism in vivo in mice. However, whether ILRUN, which is expressed in vascular cells, directly impacts atherogenesis remains unclear. We sought to determine the role of ILRUN in atherosclerosis development in mice.

METHODS

For our study, we generated global Ilrun-deficient (IlrunKO) male and female mice on 2 hyperlipidemic backgrounds: low density lipoprotein receptor knockout (LdlrKO) and apolipoprotein E knockout (ApoeKO; double knockout [DKO]).

RESULTS

Compared with littermate control mice (single LdlrKO or ApoeKO), deletion of Ilrun in DKO mice resulted in significantly attenuated both early and advanced atherosclerotic lesion development, as well as reduced necrotic area. DKO mice also had significantly decreased plasma cholesterol levels, primarily attributable to non-HDL (high-density lipoprotein) cholesterol. Hepatic-specific reconstitution of ILRUN in DKO mice on the ApoeKO background normalized plasma lipids, but atherosclerotic lesion area and necrotic area remained reduced in DKO mice. Further analysis showed that loss of Ilrun increased efferocytosis receptor MerTK expression in macrophages, enhanced in vitro efferocytosis, and significantly improved in situ efferocytosis in advanced lesions.

CONCLUSIONS

Our results support ILRUN as an important novel regulator of atherogenesis that promotes lesion progression and necrosis. It influences atherosclerosis through both plasma lipid-dependent and lipid-independent mechanisms. These findings support ILRUN as the likely causal gene responsible for genetic association of variants with coronary artery disease at this locus and suggest that suppression of ILRUN activity might be expected to reduce atherosclerosis.

ILRUN Downregulates ACE2 Expression and Blocks Infection of Human Cells by SARS-CoV-2

The human protein-coding gene ILRUN (inflammation and lipid regulator with UBA-like and NBR1-like domains; previously C6orf106) was identified as a proviral factor for Hendra virus infection and was recently characterized to function as an inhibitor of type I interferon expression. Here, we have utilized transcriptome sequencing (RNA-seq) to define cellular pathways regulated by ILRUN in the context of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection of Caco-2 cells. We find that inhibition of ILRUN expression by RNA interference alters transcription profiles of numerous cellular pathways, including upregulation of the SARS-CoV-2 entry receptor ACE2 and several other members of the renin-angiotensin aldosterone system. In addition, transcripts of the SARS-CoV-2 coreceptors TMPRSS2 and CTSL were also upregulated. Inhibition of ILRUN also resulted in increased SARS-CoV-2 replication, while overexpression of ILRUN had the opposite effect, identifying ILRUN as a novel antiviral factor for SARS-CoV-2 replication. This represents, to our knowledge, the first report of ILRUN as a regulator of the renin-angiotensin-aldosterone system (RAAS). IMPORTANCE There is no doubt that the current rapid global spread of COVID-19 has had significant and far-reaching impacts on our health and economy and will continue to do so. Research in emerging infectious diseases, such as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is growing rapidly, with new breakthroughs in the understanding of host-virus interactions to assist with the development of innovative and exciting therapeutic strategies. Here, we present the first evidence that modulation of the human protein-coding gene ILRUN functions as an antiviral factor for SARS-CoV-2 infection, likely through its newly identified role in regulating the expression of SARS-CoV-2 entry receptors ACE2, TMPRSS2, and CTSL. These data improve our understanding of biological pathways that regulate host factors critical to SARS-CoV-2 infection, contributing to the development of antiviral strategies to deal with the current SARS-CoV-2 pandemic.

Gene-based association analysis identified a novel gene associated with systemic lupus erythematosus

OBJECTIVE

Systemic lupus erythematosus (SLE) is a complex autoimmune disease with strong genetic predisposition. Genome-wide association studies (GWAS) of SLE have identified more than 50 robust susceptibility loci. However, traditional individual SNP-based GWAS have made it difficult to identify variants with small effects. Moreover, variants revealed by GWAS only explain a limited disease heritability, suggesting that many susceptibility genes remain uncovered.

METHODS

We first curated the published SLE GWAS data from 1047 SLE patients and 1205 healthy controls of Chinese ancestry and performed a gene-based association study. Then quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was conducted to verify novel genes identified above.

RESULTS

Gene-based association study identified 10 SLE-associated genes, nine of which were reported by previous GWAS, the other one, ILRUN, is a newly identified gene and was further validated by qRT-PCR. Gene expression analysis of Gene Expression Omnibus (GEO) datasets also showed that the expression of ILRUN in patients with SLE was lower than that in normal subjects.

CONCLUSION

In this study, gene-based association study and qRT-PCR identified that ILRUN is a novel susceptibility gene of SLE. ILRUN may regulate inflammation and antiviral response through its effect on the transcription of type I interferons )I-IFN, and participate in the pathogenesis of SLE.

Learning a genome-wide score of human-mouse conservation at the functional genomics level

Identifying genomic regions with functional genomic properties that are conserved between human and mouse is an important challenge in the context of mouse model studies. To address this, we develop a method to learn a score of evidence of conservation at the functional genomics level by integrating information from a compendium of epigenomic, transcription factor binding, and transcriptomic data from human and mouse. The method, Learning Evidence of Conservation from Integrated Functional genomic annotations (LECIF), trains neural networks to generate this score for the human and mouse genomes. The resulting LECIF score highlights human and mouse regions with shared functional genomic properties and captures correspondence of biologically similar human and mouse annotations. Analysis with independent datasets shows the score also highlights loci associated with similar phenotypes in both species. LECIF will be a resource for mouse model studies by identifying loci whose functional genomic properties are likely conserved.