Discovery and Functional Characterization of Two Regulatory Variants Underlying Lupus Susceptibility at 2p13.1

Antidepressant Exposure and DNA Methylation: Insights from a Methylome-Wide Association Study

Importance

Understanding antidepressant mechanisms could help design more effective and tolerated treatments.

Objective

Identify DNA methylation (DNAm) changes associated with antidepressant exposure.

Design

Case-control methylome-wide association studies (MWAS) of antidepressant exposure were performed from blood samples collected between 2006-2011 in Generation Scotland (GS). The summary statistics were tested for enrichment in specific tissues, gene ontologies and an independent MWAS in the Netherlands Study of Depression and Anxiety (NESDA). A methylation profile score (MPS) was derived and tested for its association with antidepressant exposure in eight independent cohorts, alongside prospective data from GS.

Setting

Cohorts; GS, NESDA, FTC, SHIP-Trend, FOR2107, LBC1936, MARS-UniDep, ALSPAC, E-Risk, and NTR.

Participants

Participants with DNAm data and self-report/prescription derived antidepressant exposure.

Main Outcomes and Measures

Whole-blood DNAm levels were assayed by the EPIC/450K Illumina array (9 studies, N exposed = 661, N unexposed = 9,575) alongside MBD-Seq in NESDA (N exposed = 398, N unexposed = 414). Antidepressant exposure was measured by self- report and/or antidepressant prescriptions.

Results

The self-report MWAS (N = 16,536, N exposed = 1,508, mean age = 48, 59% female) and the prescription-derived MWAS (N = 7,951, N exposed = 861, mean age = 47, 59% female), found hypermethylation at seven and four DNAm sites (p < 9.42x10 -8 ), respectively. The top locus was cg26277237 ( KANK1, p self-report = 9.3x10 -13 , p prescription = 6.1x10 -3 ). The self-report MWAS found a differentially methylated region, mapping to DGUOK-AS1 ( p adj = 5.0x10 -3 ) alongside significant enrichment for genes expressed in the amygdala, the "synaptic vesicle membrane" gene ontology and the top 1% of CpGs from the NESDA MWAS (OR = 1.39, p < 0.042). The MPS was associated with antidepressant exposure in meta-analysed data from external cohorts (N studies = 9, N = 10,236, N exposed = 661, f3 = 0.196, p < 1x10 -4 ).

Conclusions and Relevance

Antidepressant exposure is associated with changes in DNAm across different cohorts. Further investigation into these changes could inform on new targets for antidepressant treatments.

3 Key Points

Question: Is antidepressant exposure associated with differential whole blood DNA methylation?Findings: In this methylome-wide association study of 16,536 adults across Scotland, antidepressant exposure was significantly associated with hypermethylation at CpGs mapping to KANK1 and DGUOK-AS1. A methylation profile score trained on this sample was significantly associated with antidepressant exposure (pooled f3 [95%CI]=0.196 [0.105, 0.288], p < 1x10 -4 ) in a meta-analysis of external datasets. Meaning: Antidepressant exposure is associated with hypermethylation at KANK1 and DGUOK-AS1 , which have roles in mitochondrial metabolism and neurite outgrowth. If replicated in future studies, targeting these genes could inform the design of more effective and better tolerated treatments for depression.

A Non-Coding Variant in SLC15A4 Modulates Enhancer Activity and Lysosomal Deacidification Linked to Lupus Susceptibility

Systemic lupus erythematosus (SLE) is a complex autoimmune disease with a strong genetic basis. Despite the identification of several single nucleotide polymorphisms (SNPs) near the SLC15A4 gene that are significantly associated with SLE across multiple populations, specific causal SNP(s) and molecular mechanisms responsible for disease susceptibility are unknown. To address this gap, we employed bioinformatics, expression quantitative trait loci (eQTLs), and 3D chromatin interaction analysis to nominate a likely functional variant, rs35907548, in an active intronic enhancer of SLC15A4. Through luciferase reporter assays followed by chromatin immunoprecipitation (ChIP)-qPCR, we observed significant allele-specific enhancer effects of rs35907548 in diverse cell lines. The rs35907548 risk allele T is associated with increased regulatory activity and target gene expression, as shown by eQTLs and chromosome conformation capture (3C)-qPCR. The latter revealed long-range chromatin interactions between the rs35907548 enhancer and the promoters of SLC15A4, GLTLD1, and an uncharacterized lncRNA. The enhancer-promoter interactions and expression effects were validated by CRISPR/Cas9 knock-out (KO) of the locus in HL60 promyeloblast cells. KO cells also displayed dramatically dysregulated endolysosomal pH regulation. Together, our data show that the rs35907548 risk allele affects multiple aspects of cellular physiology and may directly contribute to SLE.

ARID3a from the ARID family: structure, role in autoimmune diseases and drug discovery

The AT-rich interaction domain (ARID) family of DNA-binding proteins is a group of transcription factors and chromatin regulators with a highly conserved ARID domain that recognizes specific AT-rich DNA sequences. Dysfunction of ARID family members has been implicated in various human diseases including cancers and intellectual disability. Among them, ARID3a has gained increasing attention due to its potential involvement in autoimmunity. In this article we provide an overview of the ARID family, focusing on the structure and biological functions of ARID3a. It explores the role of ARID3a in autoreactive B cells and its contribution to autoimmune diseases such as systemic lupus erythematosus and primary biliary cholangitis. Furthermore, we also discuss the potential for drug discovery targeting ARID3a and present a plan for future research in this field.

Leveraging molecular quantitative trait loci to comprehend complex diseases/traits from the omics perspective

Comprehending the molecular basis of quantitative genetic variation is a principal goal for complex diseases or traits. Molecular quantitative trait loci (molQTLs) have made it possible to investigate the effects of genetic variants hiding behind large-scale omics data. A deeper understanding of molQTL is urgently required in light of the multi-dimensionalization of omics data to more fully elucidate the pertinent biological mechanisms. Herein, we reviewed molQTLs with the corresponding resource from the omics perspective and further discussed the integrative strategy of GWAS-molQTL to infer their causal effects. Subsequently, we described the opportunities and challenges encountered by molQTL. The case studies showed that molQTL is essential for complex diseases and traits, whether single- or multi-omics QTLs. Overall, we highlighted the functional significance of genetic variants to employ the discovery of molQTL in complex diseases and traits.

A Non-Coding Variant in SLC15A4 Modulates Enhancer Activity and Lysosomal Deacidification Linked to Lupus Susceptibility

Systemic lupus erythematosus (SLE) is a complex autoimmune disease with a strong genetic basis. Despite the identification of several single nucleotide polymorphisms (SNPs) near the SLC15A4 gene that are significantly associated with SLE across multiple populations, specific causal SNP(s) and molecular mechanisms responsible for disease susceptibility are unknown. To address this gap, we employed bioinformatics, expression quantitative trait loci (eQTLs), and 3D chromatin interaction analysis to nominate a likely functional variant, rs35907548, in an active intronic enhancer of SLC15A4 . Through luciferase reporter assays followed by chromatin immunoprecipitation (ChIP)-qPCR, we observed significant allele-specific enhancer effects of rs35907548 in diverse cell lines. The rs35907548 risk allele T is associated with increased regulatory activity and target gene expression, as shown by eQTLs and chromosome conformation capture (3C)-qPCR. The latter revealed long-range chromatin interactions between the rs35907548 enhancer and the promoters of SLC15A4, GLTLD1 , and an uncharacterized lncRNA. The enhancer-promoter interactions and expression effects were validated by CRISPR/Cas9 knock-out (KO) of the locus in HL60 promyeloblast cells. KO cells also displayed dramatically dysregulated endolysosomal pH regulation. Together, our data show that the rs35907548 risk allele affects multiple aspects of cellular physiology and may directly contribute to SLE.