An autoimmune disease risk variant: A trans master regulatory effect mediated by IRF1 under immune stimulation?

Trans-eQTL hotspots shape complex traits by modulating cellular states

Regulatory genetic variation shapes gene expression, providing an important mechanism connecting DNA variation and complex traits. The causal relationships between gene expression and complex traits remain poorly understood. Here, we integrated transcriptomes and 46 genetically complex growth traits in a large cross between two strains of the yeast Saccharomyces cerevisiae. We discovered thousands of genetic correlations between gene expression and growth, suggesting potential functional connections. Local regulatory variation was a minor source of these genetic correlations. Instead, genetic correlations tended to arise from multiple independent trans-acting regulatory loci. Trans-acting hotspots that affect the expression of numerous genes accounted for particularly large fractions of genetic growth variation and of genetic correlations between gene expression and growth. Genes with genetic correlations were enriched for similar biological processes across traits but with heterogeneous direction of effect. Our results reveal how trans-acting regulatory hotspots shape complex traits by altering cellular states.

Lack of evidence for GWAS signals of exfoliation glaucoma working via monogenic loss-of-function mutation in the nearest gene

PURPOSE

Exfoliation syndrome (XFS) is a systemic disease of elastin-rich tissues involving a deposition of fibrillar exfoliative material (XFM) in the anterior chamber of the eye, which can promote glaucoma. The purpose of this study was to create mice with CRISPR/Cas9-induced variations in candidate genes identified from human genome-wide association studies (GWAS) and screen them for indices of XFS.

METHODS

Variants predicted to be deleterious were sought in the Agpat1, Cacna1a, Loxl1, Pomp, Rbms3, Sema6a, and Tlcd5 genes of C57BL/6J mice using CRISPR/Cas9-based gene editing. Strains were phenotyped by slit-lamp, SD-OCT imaging, and fundus exams at 1-5 mos of age. Smaller cohorts of 12-mos-old mice were also studied.

RESULTS

Deleterious variants were identified in six targets; Pomp was recalcitrant to targeting. Multiple alleles of some targets were isolated, yielding 12 strains. Across all genotypes and ages, 277 mice were assessed by 902 slit-lamp exams, 928 SD-OCT exams, and 358 fundus exams. Homozygosity for Agpat1 or Cacna1a mutations led to early lethality; homozygosity for Loxl1 mutations led to pelvic organ prolapse, preventing aging. Loxl1 homozygotes exhibited a conjunctival phenotype of potential relevance to XFS. Multiple other genotype-specific phenotypes were variously identified. XFM was not observed in any mice.

CONCLUSIONS

This study did not detect XFM in any of the strains. This may have been due to species-specific differences, background dependence, or insufficient aging. Alternatively, it is possible that the current candidates, selected based on proximity to GWAS signals, are not effectors acting via monogenic loss-of-function mechanisms.

The relationship between extreme inter-individual variation in macrophage gene expression and genetic susceptibility to inflammatory bowel disease

The differentiation of resident intestinal macrophages from blood monocytes depends upon signals from the macrophage colony-stimulating factor receptor (CSF1R). Analysis of genome-wide association studies (GWAS) indicates that dysregulation of macrophage differentiation and response to microorganisms contributes to susceptibility to chronic inflammatory bowel disease (IBD). Here, we analyzed transcriptomic variation in monocyte-derived macrophages (MDM) from affected and unaffected sib pairs/trios from 22 IBD families and 6 healthy controls. Transcriptional network analysis of the data revealed no overall or inter-sib distinction between affected and unaffected individuals in basal gene expression or the temporal response to lipopolysaccharide (LPS). However, the basal or LPS-inducible expression of individual genes varied independently by as much as 100-fold between subjects. Extreme independent variation in the expression of pairs of HLA-associated transcripts (HLA-B/C, HLA-A/F and HLA-DRB1/DRB5) in macrophages was associated with HLA genotype. Correlation analysis indicated the downstream impacts of variation in the immediate early response to LPS. For example, variation in early expression of IL1B was significantly associated with local SNV genotype and with subsequent peak expression of target genes including IL23A, CXCL1, CXCL3, CXCL8 and NLRP3. Similarly, variation in early IFNB1 expression was correlated with subsequent expression of IFN target genes. Our results support the view that gene-specific dysregulation in macrophage adaptation to the intestinal milieu is associated with genetic susceptibility to IBD.

Gene Signature of Regulatory T Cells Isolated from Children with Selective IgA Deficiency and Common Variable Immunodeficiency

Selective IgA deficiency (SIgAD) is the most common form and common variable immunodeficiency (CVID) is the most symptomatic form of predominant antibody deficiency. Despite differences in the clinical picture, a similar genetic background is suggested. A common feature of both disorders is the occurrence of autoimmune conditions. Regulatory T cells (Tregs) are the major immune cell type that maintains autoimmune tolerance. As the different types of abnormalities of Treg cells have been associated with autoimmune disorders in primary immunodeficiency (PID) patients, in our study we aimed to analyze the gene expression profiles of Treg cells in CVID and SIgAD patients compared to age-matched healthy controls. The transcriptome-wide gene profiling was performed by microarray technology. As a result, we analyzed and visualized gene expression patterns of isolated population of Treg cells. We showed the differences at the gene level between patients with and without autoimmunizations. Our findings suggest that the gene signatures of Treg cells isolated from SIgAD and CVID patients differ from age-matched healthy controls and from each other, presenting transcriptional profiles enriched in innate immune or Th response, respectively. The occurrence of autoimmunity in both types of PID is associated with down-regulation of class I IFNs signaling pathways. In summary, our findings improve our understanding of Treg dysfunctions in patients with common PIDs and associated autoimmunity.

Haplotype-aware modeling of cis-regulatory effects highlights the gaps remaining in eQTL data

Expression Quantitative Trait Loci (eQTLs) are critical to understanding the mechanisms underlying disease-associated genomic loci. Nearly all protein-coding genes in the human genome have been associated with one or more eQTLs. Here we introduce a multi-variant generalization of allelic Fold Change (aFC), aFC-n, to enable quantification of the cis-regulatory effects in multi-eQTL genes under the assumption that all eQTLs are known and conditionally independent. Applying aFC-n to 458,465 eQTLs in the Genotype-Tissue Expression (GTEx) project data, we demonstrate significant improvements in accuracy over the original model in estimating the eQTL effect sizes and in predicting genetically regulated gene expression over the current tools. We characterize some of the empirical properties of the eQTL data and use this framework to assess the current state of eQTL data in terms of characterizing cis-regulatory landscape in individual genomes. Notably, we show that 77.4% of the genes with an allelic imbalance in a sample show 0.5 log2 fold or more of residual imbalance after accounting for the eQTL data underlining the remaining gap in characterizing regulatory landscape in individual genomes. We further contrast this gap across tissue types, and ancestry backgrounds to identify its correlates and guide future studies.

multi-GPA-Tree: Statistical approach for pleiotropy informed and functional annotation tree guided prioritization of GWAS results

Genome-wide association studies (GWAS) have successfully identified over two hundred thousand genotype-trait associations. Yet some challenges remain. First, complex traits are often associated with many single nucleotide polymorphisms (SNPs), most with small or moderate effect sizes, making them difficult to detect. Second, many complex traits share a common genetic basis due to 'pleiotropy' and and though few methods consider it, leveraging pleiotropy can improve statistical power to detect genotype-trait associations with weaker effect sizes. Third, currently available statistical methods are limited in explaining the functional mechanisms through which genetic variants are associated with specific or multiple traits. We propose multi-GPA-Tree to address these challenges. The multi-GPA-Tree approach can identify risk SNPs associated with single as well as multiple traits while also identifying the combinations of functional annotations that can explain the mechanisms through which risk-associated SNPs are linked with the traits. First, we implemented simulation studies to evaluate the proposed multi-GPA-Tree method and compared its performance with existing statistical approaches. The results indicate that multi-GPA-Tree outperforms existing statistical approaches in detecting risk-associated SNPs for multiple traits. Second, we applied multi-GPA-Tree to a systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA), and to a Crohn's disease (CD) and ulcertive colitis (UC) GWAS, and functional annotation data including GenoSkyline and GenoSkylinePlus. Our results demonstrate that multi-GPA-Tree can be a powerful tool that improves association mapping while facilitating understanding of the underlying genetic architecture of complex traits and potential mechanisms linking risk-associated SNPs with complex traits.

Integrated analysis of Mendelian Randomization and Bayesian colocalization reveals bidirectional causal association between inflammatory bowel disease and psoriasis

BACKGROUND

Observational studies have suggested an association between inflammatory bowel disease [IBD] and psoriasis. However, the detailed genetic basis, causality, and direction of this association remain unclear.

METHODS

Bidirectional two-sample Mendelian Randomization [MR] analysis was conducted using summary statistics from published genome-wide association studies. Bayesian Colocalization and multivariable MR [MVMR] analyses were performed to identify candidate variants and risk genes involved in the shared genetic basis between IBD, psoriasis, and their subtypes.

RESULTS

Genetically predicted IBD and Crohn's disease [CD] were associated with an increased risk of psoriasis, psoriasis vulgaris [PsV], and psoriatic arthritis [PsA] (IBD on psoriasis: pooled odds ratio [OR] 1.09, 95% confidence interval [CI] 1.04-1.14, p = .0001; CD on psoriasis: pooled OR 1.10, 95% CI 1.06-1.15, p < .0001) and vice versa (psoriasis on IBD: pooled OR 1.11, 95%CI 1.02-1.21), whereas CD only exhibited a unidirectional association with psoriasis. Colocalization analysis revealed eight candidate genetic variants and risk genes (including LINC00824, CDKAL1, IL10, IL23R, DNAJC27, LPP, RUNX3, and RGS14) associated with a shared genetic basis. Among these, IL23R, DNAJC27, LPP, and RGS14 were further validated by MVMR analysis.

CONCLUSION

Our findings indicated bidirectional causal associations between IBD and psoriasis (including PsV and PsA), which were attributed primarily to CD rather than Ulcerative colitis [UC]. Furthermore, we identified several candidate variants and risk genes involved in the shared genetic basis of IBD and psoriasis. Acquiring a better understanding of the shared genetic architecture underlying IBD and psoriasis would help improve clinical strategies.

Inferring cell-type-specific causal gene regulatory networks during human neurogenesis

BACKGROUND

Genetic variation influences both chromatin accessibility, assessed in chromatin accessibility quantitative trait loci (caQTL) studies, and gene expression, assessed in expression QTL (eQTL) studies. Genetic variants can impact either nearby genes (cis-eQTLs) or distal genes (trans-eQTLs). Colocalization between caQTL and eQTL, or cis- and trans-eQTLs suggests that they share causal variants. However, pairwise colocalization between these molecular QTLs does not guarantee a causal relationship. Mediation analysis can be applied to assess the evidence supporting causality versus independence between molecular QTLs. Given that the function of QTLs can be cell-type-specific, we performed mediation analyses to find epigenetic and distal regulatory causal pathways for genes within two major cell types of the developing human cortex, progenitors and neurons.

RESULTS

We find that the expression of 168 and 38 genes is mediated by chromatin accessibility in progenitors and neurons, respectively. We also find that the expression of 11 and 12 downstream genes is mediated by upstream genes in progenitors and neurons. Moreover, we discover that a genetic locus associated with inter-individual differences in brain structure shows evidence for mediation of SLC26A7 through chromatin accessibility, identifying molecular mechanisms of a common variant association to a brain trait.

CONCLUSIONS

In this study, we identify cell-type-specific causal gene regulatory networks whereby the impacts of variants on gene expression were mediated by chromatin accessibility or distal gene expression. Identification of these causal paths will enable identifying and prioritizing actionable regulatory targets perturbing these key processes during neurodevelopment.

A lncRNA from an inflammatory bowel disease risk locus maintains intestinal host-commensal homeostasis

Inflammatory bowel diseases (IBD) are known to have complex, genetically influenced etiologies, involving dysfunctional interactions between the intestinal immune system and the microbiome. Here, we characterized how the RNA transcript from an IBD-associated long non-coding RNA locus ("CARINH-Colitis Associated IRF1 antisense Regulator of Intestinal Homeostasis") protects against IBD. We show that CARINH and its neighboring gene coding for the transcription factor IRF1 together form a feedforward loop in host myeloid cells. The loop activation is sustained by microbial factors, and functions to maintain the intestinal host-commensal homeostasis via the induction of the anti-inflammatory factor IL-18BP and anti-microbial factors called guanylate-binding proteins (GBPs). Extending these mechanistic insights back to humans, we demonstrate that the function of the CARINH/IRF1 loop is conserved between mice and humans. Genetically, the T allele of rs2188962, the most probable causal variant of IBD within the CARINH locus from the human genetics study, impairs the inducible expression of the CARINH/IRF1 loop and thus increases genetic predisposition to IBD. Our study thus illustrates how an IBD-associated lncRNA maintains intestinal homeostasis and protects the host against colitis.

Functional Characterization of Genetic Variant Effects on Expression

Thousands of common genetic variants in the human population have been associated with disease risk and phenotypic variation by genome-wide association studies (GWAS). However, the majority of GWAS variants fall into noncoding regions of the genome, complicating our understanding of their regulatory functions, and few molecular mechanisms of GWAS variant effects have been clearly elucidated. Here, we set out to review genetic variant effects, focusing on expression quantitative trait loci (eQTLs), including their utility in interpreting GWAS variant mechanisms. We discuss the interrelated challenges and opportunities for eQTL analysis, covering determining causal variants, elucidating molecular mechanisms of action, and understanding context variability. Addressing these questions can enable better functional characterization of disease-associated loci and provide insights into fundamental biological questions of the noncoding genetic regulatory code and its control of gene expression. Expected final online publication date for the Annual Review of Biomedical Data Science, Volume 5 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

GREEN-DB: a framework for the annotation and prioritization of non-coding regulatory variants from whole-genome sequencing data

Non-coding variants have long been recognized as important contributors to common disease risks, but with the expansion of clinical whole genome sequencing, examples of rare, high-impact non-coding variants are also accumulating. Despite recent advances in the study of regulatory elements and the availability of specialized data collections, the systematic annotation of non-coding variants from genome sequencing remains challenging. Here, we propose a new framework for the prioritization of non-coding regulatory variants that integrates information about regulatory regions with prediction scores and HPO-based prioritization. Firstly, we created a comprehensive collection of annotations for regulatory regions including a database of 2.4 million regulatory elements (GREEN-DB) annotated with controlled gene(s), tissue(s) and associated phenotype(s) where available. Secondly, we calculated a variation constraint metric and showed that constrained regulatory regions associate with disease-associated genes and essential genes from mouse knock-outs. Thirdly, we compared 19 non-coding impact prediction scores providing suggestions for variant prioritization. Finally, we developed a VCF annotation tool (GREEN-VARAN) that can integrate all these elements to annotate variants for their potential regulatory impact. In our evaluation, we show that GREEN-DB can capture previously published disease-associated non-coding variants as well as identify additional candidate disease genes in trio analyses.

Transcription factor regulation of eQTL activity across individuals and tissues

Tens of thousands of genetic variants associated with gene expression (cis-eQTLs) have been discovered in the human population. These eQTLs are active in various tissues and contexts, but the molecular mechanisms of eQTL variability are poorly understood, hindering our understanding of genetic regulation across biological contexts. Since many eQTLs are believed to act by altering transcription factor (TF) binding affinity, we hypothesized that analyzing eQTL effect size as a function of TF level may allow discovery of mechanisms of eQTL variability. Using GTEx Consortium eQTL data from 49 tissues, we analyzed the interaction between eQTL effect size and TF level across tissues and across individuals within specific tissues and generated a list of 10,098 TF-eQTL interactions across 2,136 genes that are supported by at least two lines of evidence. These TF-eQTLs were enriched for various TF binding measures, supporting with orthogonal evidence that these eQTLs are regulated by the implicated TFs. We also found that our TF-eQTLs tend to overlap genes with gene-by-environment regulatory effects and to colocalize with GWAS loci, implying that our approach can help to elucidate mechanisms of context-specificity and trait associations. Finally, we highlight an interesting example of IKZF1 TF regulation of an APBB1IP gene eQTL that colocalizes with a GWAS signal for blood cell traits. Together, our findings provide candidate TF mechanisms for a large number of eQTLs and offer a generalizable approach for researchers to discover TF regulators of genetic variant effects in additional QTL datasets.

Interferon Lambda in the Pathogenesis of Inflammatory Bowel Diseases

Interferon λ (IFN-λ) is critical for host viral defense at mucosal surfaces and stimulates immunomodulatory signals, acting on epithelial cells and few other cell types due to restricted IFN-λ receptor expression. Epithelial cells of the intestine play a critical role in the pathogenesis of Inflammatory Bowel Disease (IBD), and the related type II interferons (IFN-γ) have been extensively studied in the context of IBD. However, a role for IFN-λ in IBD onset and progression remains unclear. Recent investigations of IFN-λ in IBD are beginning to uncover complex and sometimes opposing actions, including pro-healing roles in colonic epithelial tissues and potentiation of epithelial cell death in the small intestine. Additionally, IFN-λ has been shown to act through non-epithelial cell types, such as neutrophils, to protect against excessive inflammation. In most cases IFN-λ demonstrates an ability to coordinate the host antiviral response without inducing collateral hyperinflammation, suggesting that IFN-λ signaling pathways could be a therapeutic target in IBD. This mini review discusses existing data on the role of IFN-λ in the pathogenesis of inflammatory bowel disease, current gaps in the research, and therapeutic potential of modulating the IFN-λ-stimulated response.

TF-RBP-AS Triplet Analysis Reveals the Mechanisms of Aberrant Alternative Splicing Events in Kidney Cancer: Implications for Their Possible Clinical Use as Prognostic and Therapeutic Biomarkers

Aberrant alternative splicing (AS) is increasingly linked to cancer; however, how AS contributes to cancer development still remains largely unknown. AS events (ASEs) are largely regulated by RNA-binding proteins (RBPs) whose ability can be modulated by a variety of genetic and epigenetic mechanisms. In this study, we used a computational framework to investigate the roles of transcription factors (TFs) on regulating RBP-AS interactions. A total of 6519 TF–RBP–AS triplets were identified, including 290 TFs, 175 RBPs, and 16 ASEs from TCGA–KIRC RNA sequencing data. TF function categories were defined according to correlation changes between RBP expression and their targeted ASEs. The results suggested that most TFs affected multiple targets, and six different classes of TF-mediated transcriptional dysregulations were identified. Then, regulatory networks were constructed for TF–RBP–AS triplets. Further pathway-enrichment analysis showed that these TFs and RBPs involved in triplets were enriched in a variety of pathways that were associated with cancer development and progression. Survival analysis showed that some triplets were highly associated with survival rates. These findings demonstrated that the integration of TFs into alternative splicing regulatory networks can help us in understanding the roles of alternative splicing in cancer.