Methylation Quantitative Trait Loci Research Articles

An Epigenome-Wide Association Study of DNA Methylation and Proliferative Retinopathy over 28 Years in Type 1 Diabetes

To perform a prospective epigenome-wide association study of DNA methylation (DNAm) and 28-year proliferative diabetic retinopathy (PDR) incidence in type 1 diabetes (T1D). Prospective observational cohort study. The Pittsburgh Epidemiology of Diabetes Complications (EDC) study of childhood-onset (< 17 years) T1D. Stereoscopic fundus photographs were taken in fields 1, 2, and 4 at baseline, 2, 4, 6, 8, 16, 23, and 28 years after DNAm measurements. The photos were graded using the modified Airlie House System. In those free of PDR at baseline (n=265; mean T1D duration of 18 years at baseline), whole blood DNAm (EPIC array) at 683 597 CpGs was analyzed in Cox models for time to event. Associations between significant CpGs and clinical risk factors were assessed; genetic variants associated with DNAm were identified (methylation quantitative trait loci [meQTLs]). Mendelian randomization was used to examine evidence of causal associations between DNAm and PDR. Post hoc regional and functional analyses were performed. Proliferative diabetic retinopathy was defined as the first instance of a grade of ≥ 60 in at least 1 eye or pan-retinal photocoagulation for PDR. Follow-up time was calculated from the study visit at which DNAm data were available (baseline) until PDR incidence or censoring (December 31, 2018 or last follow-up). PDR incidence was 53% over 28-years' follow-up. Greater DNAm of cg27512687 (KIF16B) was associated with reduced PDR incidence (P=6.3×10-9; false discovery rate [FDR]: < 0.01); 113 cis-meQTLs (P<5×10-8) were identified. Mendelianrandomization analysis using the sentinel meQTL as the instrumental variable supported a potentially causal association between cg27512687 and PDR. Cg27512687 was also associated with lower pulse rate and albumin excretion rate and higher estimated glomerular filtration rate, but its association with PDR remained independently significant after adjustment for those factors. In regional analyses, DNAm of FUT4, FKBP1A, and RIN2 was also associated with PDR incidence. DNA methylation of KIF16B, FUT4, FKBP1A, and RIN2 was associated with PDR incidence, supporting roles for epigenetic regulation of iron clearance, developmental pathways, and autophagy in PDR pathogenesis. Further study of those loci may provide insight into novel targets for interventions to prevent or delay PDR in T1D. Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Circadian clock-related genome-wide mendelian randomization identifies putatively genes for ulcerative colitis and its comorbidity

BackgroundCircadian rhythm is crucial to the function of the immune system. Disorders of the circadian rhythm can contribute to inflammatory diseases such as Ulcerative colitis (UC). This Mendelian Randomization (MR) analysis applies genetic tools to represent the aggregated statistical results of exposure to circadian rhythm disorders and UC and its comorbidities, allowing for causal inferences.MethodsSummary statistics of protein, DNA methylation and gene expression quantitative trait loci in individuals of European ancestry (pQTL, mQTL, and eQTL, respectively) were used. Genetic variants located within or near 152 circadian clock-related genes and closely related to circadian rhythm disorders were selected as instrumental variables. Causal relationships with UC and its comorbidities were then estimated through employed Summary data-based Mendelian Randomization (SMR) and Inverse-Variance-Weighted MR (IVW-MR).ResultsThrough preliminary SMR analysis, we identified a potential causal relationship between circadian clock-related genes and UC along with its comorbidities, which was further confirmed by IVW-MR analysis. Our study identified strong evidence of positive correlation involving seven overlapping genes (CSNK1E, OPRL1, PIWIL2, RORC, MAX, PPP5C, and AANAT) through MWAS and TWAS in UC, four overlapping genes (OPRL1, CHRNB2, FBXL17, and SIRT1) in UC with PSC, and three overlapping genes (ARNTL, USP7, and KRAS) in UC with arthropathy.ConclusionsThis SMR study demonstrates the causal effect of circadian rhythm disorders in UC and its comorbidities. Furthermore, our investigation pinpointed candidate genes that could potentially serve as drug targets.

Sex-dependent placental methylation quantitative trait loci provide insight into the prenatal origins of childhood onset traits and conditions

Molecular quantitative trait loci (QTLs) allow us to understand the biology captured in genome-wide association studies (GWASs). The placenta regulates fetal development and shows sex differences in DNA methylation. We therefore hypothesized that placental methylation QTL (mQTL) explain variation in genetic risk for childhood onset traits, and that effects differ by sex. We analyzed 411 term placentas from two studies and found 49,252 methylation (CpG) sites with mQTL and 2,489 CpG sites with sex-dependent mQTL. All mQTL were enriched in regions that typically affect gene expression in prenatal tissues. All mQTL were also enriched in GWAS results for growth- and immune-related traits, but male- and female-specific mQTL were more enriched than cross-sex mQTL. mQTL colocalized with trait loci at 777 CpG sites, with 216 (28%) specific to males or females. Overall, mQTL specific to male and female placenta capture otherwise overlooked variation in childhood traits.

Multi-omics and pathway analyses of genome-wide associations implicate regulation and immunity in verbal declarative memory performance

BackgroundUncovering the functional relevance underlying verbal declarative memory (VDM) genome-wide association study (GWAS) results may facilitate the development of interventions to reduce age-related memory decline and dementia.MethodsWe performed multi-omics and pathway enrichment analyses of paragraph (PAR-dr) and word list (WL-dr) delayed recall GWAS from 29,076 older non-demented individuals of European descent. We assessed the relationship between single-variant associations and expression quantitative trait loci (eQTLs) in 44 tissues and methylation quantitative trait loci (meQTLs) in the hippocampus. We determined the relationship between gene associations and transcript levels in 53 tissues, annotation as immune genes, and regulation by transcription factors (TFs) and microRNAs. To identify significant pathways, gene set enrichment was tested in each cohort and meta-analyzed across cohorts. Analyses of differential expression in brain tissues were conducted for pathway component genes.ResultsThe single-variant associations of VDM showed significant linkage disequilibrium (LD) with eQTLs across all tissues and meQTLs within the hippocampus. Stronger WL-dr gene associations correlated with reduced expression in four brain tissues, including the hippocampus. More robust PAR-dr and/or WL-dr gene associations were intricately linked with immunity and were influenced by 31 TFs and 2 microRNAs. Six pathways, including type I diabetes, exhibited significant associations with both PAR-dr and WL-dr. These pathways included fifteen MHC genes intricately linked to VDM performance, showing diverse expression patterns based on cognitive status in brain tissues.ConclusionsVDM genetic associations influence expression regulation via eQTLs and meQTLs. The involvement of TFs, microRNAs, MHC genes, and immune-related pathways contributes to VDM performance in older individuals.

Genome-Wide Interaction Analysis With DASH Diet Score Identified Novel Loci for Systolic Blood Pressure.

The Dietary Approaches to Stop Hypertension (DASH) diet score lowers blood pressure (BP). We examined interactions between genotype and the DASH diet score in relation to systolic BP. We analyzed up to 9 420 585 single nucleotide polymorphisms in up to 127 282 individuals of 6 population groups (91% of European population) from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (n=35 660) and UK Biobank (n=91 622) and performed European population-specific and cross-population meta-analyses. We identified 3 loci in European-specific analyses and an additional 4 loci in cross-population analyses at Pinteraction<5e-8. We observed a consistent interaction between rs117878928 at 15q25.1 (minor allele frequency, 0.03) and the DASH diet score (Pinteraction=4e-8; P for heterogeneity, 0.35) in European population, where the interaction effect size was 0.42±0.09 mm Hg (Pinteraction=9.4e-7) and 0.20±0.06 mm Hg (Pinteraction=0.001) in Cohorts for Heart and Aging Research in Genomic Epidemiology and the UK Biobank, respectively. The 1 Mb region surrounding rs117878928 was enriched with cis-expression quantitative trait loci (eQTL) variants (P=4e-273) and cis-DNA methylation quantitative trait loci variants (P=1e-300). Although the closest gene for rs117878928 is MTHFS, the highest narrow sense heritability accounted by single nucleotide polymorphisms potentially interacting with the DASH diet score in this locus was for gene ST20 at 15q25.1. We demonstrated gene-DASH diet score interaction effects on systolic BP in several loci. Studies with larger diverse populations are needed to validate our findings.

Epigenome-wide association study of DNA methylation in maternal blood leukocytes with BMI in pregnancy and gestational weight gain

ObjectivesWe aimed to discover CpG sites with differential DNA methylation in peripheral blood leukocytes associated with body mass index (BMI) in pregnancy and gestational weight gain (GWG) in women of European and South Asian ancestry. Furthermore, we aimed to investigate how the identified sites were associated with methylation quantitative trait loci, gene ontology, and cardiometabolic parameters.MethodsIn the Epigenetics in pregnancy (EPIPREG) sample we quantified maternal DNA methylation in peripheral blood leukocytes in gestational week 28 with Illumina’s MethylationEPIC BeadChip. In women with European (n = 303) and South Asian (n = 164) ancestry, we performed an epigenome-wide association study of BMI in gestational week 28 and GWG between gestational weeks 15 and 28 using a meta-analysis approach. Replication was performed in the Norwegian Mother, Father, and Child Cohort Study, the Study of Assisted Reproductive Technologies (MoBa-START) (n = 877, mainly European/Norwegian).ResultsWe identified one CpG site significantly associated with GWG (p 5.8 × 10−8) and five CpG sites associated with BMI at gestational week 28 (p from 4.0 × 10–8 to 2.1 × 10–10). Of these, we were able to replicate three in MoBa-START; cg02786370, cg19758958 and cg10472537. Two sites are located in genes previously associated with blood pressure and BMI. DNA methylation at the three replicated CpG sites were associated with levels of blood pressure, lipids and glucose in EPIPREG (p from 1.2 × 10−8 to 0.04).ConclusionsWe identified five CpG sites associated with BMI at gestational week 28, and one with GWG. Three of the sites were replicated in an independent cohort. Several genetic variants were associated with DNA methylation at cg02786379 and cg16733643 suggesting a genetic component influencing differential methylation. The identified CpG sites were associated with cardiometabolic traits.ClinicalTrials.gov registration noNot applicable

Epigenetic link between Agent Orange exposure and type 2 diabetes in Korean veterans.

Conflicting findings have been reported regarding the association between Agent Orange (AO) exposure and type 2 diabetes. This study aimed to examine whether AO exposure is associated with the development of type 2 diabetes and to verify the causal relationship between AO exposure and type 2 diabetes by combining DNA methylation with DNA genotype analyses. An epigenome-wide association study and DNA genotype analyses of the blood of AO-exposed and AO-unexposed individuals with type 2 diabetes and that of healthy controls were performed. Methylation quantitative trait locus and Mendelian randomisation analyses were performed to evaluate the causal effect of AO-exposure-identified CpGs on type 2 diabetes. AO-exposed individuals with type 2 diabetes were associated with six hypermethylated CpG sites (cg20075319, cg21757266, cg05203217, cg20102280, cg26081717, and cg21878650) and one hypo-methylated CpG site (cg07553761). Methylation quantitative trait locus analysis showed the methylation levels of some CpG sites (cg20075319, cg20102280, and cg26081717) to be significantly different. Mendelian randomisation analysis showed that CpG sites that were differentially methylated in AO-exposed individuals were causally associated with type 2 diabetes; the reverse causal effect was not significant. These findings reflect the need for further epigenetic studies on the causal relationship between AO exposure and type 2 diabetes.

PTGES2 and RNASET2 identified as novel potential biomarkers and therapeutic targets for basal cell carcinoma: insights from proteome-wide mendelian randomization, colocalization, and MR-PheWAS analyses.

Basal cell carcinoma (BCC) is the most common skin cancer, lacking reliable biomarkers or therapeutic targets for effective treatment. Genome-wide association studies (GWAS) can aid in identifying drug targets, repurposing existing drugs, predicting clinical trial side effects, and reclassifying patients in clinical utility. Hence, the present study investigates the association between plasma proteins and skin cancer to identify effective biomarkers and therapeutic targets for BCC. Proteome-wide mendelian randomization was performed using inverse-variance-weight and Wald Ratio methods, leveraging 1Mb cis protein quantitative trait loci (cis-pQTLs) in the UK Biobank Pharma Proteomics Project (UKB-PPP) and the deCODE Health Study, to determine the causal relationship between plasma proteins and skin cancer and its subtypes in the FinnGen R10 study and the SAIGE database of Lee lab. Significant association with skin cancer and its subtypes was defined as a false discovery rate (FDR) < 0.05. pQTL to GWAS colocalization analysis was executed using a Bayesian model to evaluate five exclusive hypotheses. Strong colocalization evidence was defined as a posterior probability for shared causal variants (PP.H4) of ≥0.85. Mendelian randomization-Phenome-wide association studies (MR-PheWAS) were used to evaluate potential biomarkers and therapeutic targets for skin cancer and its subtypes within a phenome-wide human disease category. PTGES2, RNASET2, SF3B4, STX8, ENO2, and HS3ST3B1 (besides RNASET2, five other plasma proteins were previously unknown in expression quantitative trait loci (eQTL) and methylation quantitative trait loci (mQTL)) were significantly associated with BCC after FDR correction in the UKB-PPP and deCODE studies. Reverse MR showed no association between BCC and these proteins. PTGES2 and RNASET2 exhibited strong evidence of colocalization with BCC based on a posterior probability PP.H4 >0.92. Furthermore, MR-PheWAS analysis showed that BCC was the most significant phenotype associated with PTGES2 and RNASET2 among 2,408 phenotypes in the FinnGen R10 study. Therefore, PTGES2 and RNASET2 are highlighted as effective biomarkers and therapeutic targets for BCC within the phenome-wide human disease category. The study identifies PTGES2 and RNASET2 plasma proteins as novel, reliable biomarkers and therapeutic targets for BCC, suggesting more effective clinical application strategies for patients.

Multi-omic insight into the molecular networks of mitochondrial dysfunction in the pathogenesis of inflammatory bowel disease

Mitochondrial dysfunction has been linked to the development of inflammatory bowel disease (IBD), but the genetic pathophysiology was not fully elucidated. We employed Mendelian randomization and colocalization analyses to investigate the associations between mitochondrial-related genes and IBD via integrating multi-omics. Summary-level data of mitochondrial gene methylation, expression and protein abundance levels were obtained from corresponding methylation, expression and protein quantitative trait loci studies, respectively. We obtained genetic associations with IBD and its two subtypes from the Inflammatory Bowel Disease Genetics Consortium (discovery), the UK Biobank (replication), and the FinnGen study (replication). We performed summary-data-based Mendelian randomization analysis to assess the associations of mitochondrial gene-related molecular features with IBD. Colocalization analysis was further conducted to assess whether the identified signal pairs shared a causal genetic variant. After integrating the multi-omics data between mQTL-eQTL and eQTL-pQTL, we identified two mitochondrial genes, i.e., PARK7 and ACADM, with tier 1 evidence for their associations with IBD and ulcerative colitis (UC). PDK1 and FISI genes were associated with UC risk with tier 2 and tier 3 evidence, respectively. The methylation of cg05467918 in ACADM was associated with lower expression of ACADM, which fits with the positive effect of cg05467918 methylation on UC risk. Consistently, the inverse associations between gene methylation and gene expression were also observed in PARK7 (cg10385390) and PDK1 (cg17679246), which were corroborated with the protective role in UC. At circulating protein level, genetically predicted higher levels of PARK7 (OR 0.36, 95% CI 0.25-0.52) and HINT1 (OR 0.47, 95% CI 0.30-0.74) were inversely associated with IBD risk; genetically predicted higher level of HINT1 was associated with a decreased risk of Crohn's disease (CD) (OR 0.26, 95% CI 0.14-0.49) and a higher level of ACADM (OR 0.67, 95% CI 0.55-0.83), PDK1 (OR 0.63, 95% CI 0.49-0.81), FIS1 (OR 0.63, 95% CI 0.47-0.83) was associated with a decreased risk of UC. We found that the mitochondrial PARK7 gene was putatively associated with IBD risk, and mitochondrial FIS1, PDK1, and ACADM genes were associated with UC risk with evidence from multi-omics levels. This study identified mitochondrial genes in relation to IBD, which may enhance the understanding of the pathogenic mechanisms of IBD development. XL is supported by the Natural Science Fund for Distinguished Young Scholars of Zhejiang Province (LR22H260001) and Healthy Zhejiang One Million People Cohort (K-20230085).

High-resolution DNA methylation screening of the major histocompatibility complex in multiple sclerosis.

The HLA-DRB1 gene in the major histocompatibility complex (MHC) region in chromosome 6p21 is the strongest genetic factor identified as influencing multiple sclerosis (MS) susceptibility. DNA methylation changes associated with MS have been consistently detected at the MHC region. However, understanding the full scope of epigenetic regulations of the MHC remains incomplete, due in part to the limited coverage of this region by standard whole genome bisulfite sequencing or array-based methods. We developed and validated an MHC capture protocol coupled with bisulfite sequencing and conducted a comprehensive analysis of the MHC methylation landscape in blood samples from 147 treatment naïve MS study participants and 129 healthy controls. We identified 132 differentially methylated region (DMRs) within MHC region associated with disease status. The DMRs overlapped with established MS risk loci. Integration of the MHC methylome with human leukocyte antigen (HLA) genetic data indicate that the methylation changes are significantly associated with HLA genotypes. Using DNA methylation quantitative trait loci (mQTL) mapping and the causal inference test (CIT), we identified 643 cis-mQTL-DMRs paired associations, including 71 DMRs possibly mediating causal relationships between 55 single nucleotide polymorphisms (SNPs) and MS risk. The results describe MS-associated methylation changes in MHC region and highlight the association between HLA genotypes and methylation changes. Results from the mQTL and CIT analyses provide evidence linking MHC region variations, methylation changes, and disease risk for MS.

Abstract C003: Analysis of genetic control of differential DNA methylation between breast cancer patients of African and European ancestry

Abstract Black women with significant African ancestry experience a disproportionately higher incidence of the most aggressive and invasive breast cancer subtypes, along with worse clinical outcomes across all population groups in the United States. These health disparities arise from complex interactions between genetic, sociocultural, and environmental factors. Although there is a critical need to understand how these non-genetic factors contribute to disease risk or etiology, they are difficult to measure reliably, and their long-term impacts are not well understood. Epigenetic mechanisms that regulate gene expression, such as DNA methylation (DNAm), have been proposed to capture the molecular impact of adverse environments on tumor biology. While DNAm levels can change in response to environmental cues, they are also regulated by individual genotypes at methylation quantitative trait loci (meQTLs), presenting a challenge for modeling environmental effects on cancer molecular phenotypes. We hypothesize that ancestry differential DNAm in breast cancer can be partially explained by ancestry differential genotype frequencies of SNPs at meQTL loci. We leveraged multi-omic data from The Cancer Genome Atlas (TCGA) breast cancer (BRCA) cohort to perform integrative analysis for patients of European and African ancestry (n=578), including ancestry-based differential methylation analysis, meQTL analysis and association testing of genotypes and ancestry. We identified 804 CpGs differentially methylated between European and African ancestry BRCA patients. meQTL analysis at these sites identified 130 cis-meQTLs and 1,019 trans-meQTLs. ~84% of the identified cis- and ~80% of the identified trans-meQTLs show evidence of differential genotype frequencies by ancestry. Gene set enrichment analysis of the meQTL target genes implicated genes involved in regulation of cell differentiation and developmental processes, phosphatase activity and phosphorylation processes. Our findings show that ancestry specific DNAm differences display strong genetic associations via ancestry associated meQTLs (a-meQTLs) affecting fundamental biological functions in breast cancer. This result points to the potential for confounding due to genetic ancestry in studies that use DNA methylation as a molecular readout of environmental exposures. Citation Format: Kyriaki Founta, Nyasha Chambwe. Analysis of genetic control of differential DNA methylation between breast cancer patients of African and European ancestry [abstract]. In: Proceedings of the 16th AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2023 Sep 29-Oct 2;Orlando, FL. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2023;32(12 Suppl):Abstract nr C003.

Integration of Alzheimer’s Disease GWAS with Molecular QTLs Reveals Prioritized Risk Genes and Disease‐associated Molecular Alterations

AbstractBackgroundDespite the success of genome‐wide association studies (GWAS) to identify genetic variants associated with complex diseases, as the majority of these variants reside in non‐coding genome, one important task is to link these variants to disease‐associated molecular alterations in risk genes. To this end, we developed a systematic gene prioritization methodology that integrated our recent Alzheimer’s disease (AD) GWAS meta‐analysis with molecular quantitative trait loci (QTL) in disease‐relevant tissues and identified highly likely candidate risk genes and molecular alterations within the majority of 42 novel loci identified. In this subsequent study, we updated our pipeline with new molecular QTL catalogues and analyzed the complete landscape of AD risk.MethodWe performed GWAS‐QTL integration in AD‐relevant tissues and cells by querying the lead variants in each locus in the molecular QTL catalogues, (ii) investigating genetic colocalization between GWAS and molecular QTL signals using coloc, and (iii) conducting transcriptome, methylome, and proteome‐wide association studies (TWAS/MWAS/PWAS) using the EADB stage I meta‐analysis GWAS on 85,934 cases and 401,577 controls. The updated gene prioritization pipeline included new analyses based on new brain protein expression QTL (pQTL), cell‐type‐specific expression QTL (ct‐eQTL), expression QTL (eQTL), histone acetylation QTL (haQTL), and methylation QTL (mQTL) catalogues. Moreover, we used Nanopore RNA sequencing to further investigate risk‐associated splicing events revealed by splicing QTL (sQTL) integration.ResultOur brain ct‐eQTL coloc results showed that while the most risk‐associated gene expression alterations are specific to microglia, the proportion of astrocyte‐specific eQTL colocalizations were doubled compared to the previous studies. These included EGFR, a previously prioritized gene due to bulk brain eQTL coloc hits, whose downregulation was regulated by the protective signal specifically in astrocytes. Importantly, in a previously unresolved locus we could prioritize HS3ST5 through astrocyte‐specific eQTLs as well. Furthermore, pQTL coloc and PWAS prioritized SNX1, a previously prioritized gene at a lower confidence together with neighboring FAM96A gene, by implicating its genetic downregulation with increased AD risk.ConclusionOur updated pipeline systematically prioritized AD risk genes and disease‐associated molecular alterations in the majority of risk loci, and contributed to better interpretation of AD GWAS and understanding of genetically‐associated molecular mechanisms underlying AD risk.

Establishing the contribution of genetic risk loci in the Lewy dody diseases

AbstractBackgroundDementia with Lewy bodies (DLB), Parkinson’s disease (PD) and Parkinson’s disease dementia (PDD) are termed Lewy body diseases (LBDs). These neurodegenerative diseases are classified by the accumulation of alpha‐synuclein in neurons, forming Lewy bodies (LB). Whilst pathologically similar, symptomatic staging is different. Previous genome wide association studies (GWAS) have identified classic Alzheimer’s disease (AD) and PD loci as risk loci in DLB. Therefore, we have used polygenic risk score (PRS) analysis to investigate the contribution of AD and PD variants in the LBDs.MethodDNA was extracted from bulk brain tissue from 481 individuals (162 controls, 95 DLB, 156 PDD and 68 PD) and genotyped on the Illumina global screening array. After a strict quality control and imputation pipeline, PRSs for multiple relevant traits using reference GWAS summary statistics were calculated using PRSice2.ResultStudy groups have been sourced consisting of cases with PD, PDD and DLB based on LB deposition and clinical symptom staging. All groups, including controls are stratified for age and concomitant AD pathology. Assessing multiple PRSs for association with pathological development and clinical outcomes we have identified the contribution of genetic risk from related neurodegenerative and neuropsychiatric disorders in our cohort of PD, PDD and DLB samples.ConclusionWe have collated a well powered study cohort to interrogate the genetic basis of the Lewy body dementias and the relative contribution of AD and PD risk loci in the LBDs. We are currently generating methylation quantitative trait loci (mQTLs) for this dataset using previously generated matched DNA methylation data in the brain.

An Epigenome‐wide association study of psychosis in Alzheimer’s disease dorsolateral prefrontal cortex

AbstractBackgroundPsychosis is a debilitating syndrome occurring in 40‐60% of people with Alzheimer’s disease (AD) and corresponds with a more severe disease course. Evidence suggests that psychosis in AD (AD+P) is associated with a distinct profile of neurobiological changes, but little is known about the molecular processes driving aetiology. In this study, we performed an epigenome‐wide association study (EWAS) to investigate DNA methylation associated with AD+P in the dorsolateral prefrontal cortex of post‐mortem brain samples.MethodBrain samples were obtained from the University of Pittsburgh Alzheimer’s disease Research Centre (PITT‐ADRC). AD pathology was assessed and classified using CERAD neurotic plaque density score, Braak neurofibrillary tangle stages and NIA‐RI criteria. We used the presence or absence of delusions and hallucinations, Mini‐Mental State Exam (MMSE) and Clinical Dementia Rating (CDR) to categorise our samples into AD with Psychosis (AD+P) and AD without Psychosis (AD‐P) groups. DNA was extracted, bisulfite‐treated and profiled on the Illumina Methylation EPIC Array and genotyping was carried out on the Illumina Global Screening Array (GSA).ResultAfter data quality control, 153 samples (40 AD‐P, 113 AD+P) remained for further analysis. The normalized beta values and the effect estimate values were used to extract Differentially Methylated positions (DMP) and differentially methylated regions (DMRs), respectively. To identify methylation quantitative trait loci (mQTLs) in both cis and trans, SNP‐CpG pairs with a p‐value threshold of &lt;1e‐5 were calculated. We performed a colocalization analysis for the most significant DMPs and mQTLs, using publicly available summary statistics of relevant genome‐wide association studies (GWAS).ConclusionThe development of effective therapies for AD+P is an urgent priority and requires a better understanding of the molecular mechanisms underlying this syndrome. To address this, we have collated a well‐powered study cohort to interrogate the epigenetic basis of AD+P in the brain and have identified several loci of interest that require further study.

Epigenomic profiling of the infrapatellar fat pad in osteoarthritis.

Osteoarthritis is a prevalent, complex disease of the joints, and affects multiple intra-articular tissues. Here, we have examined genome-wide DNA methylation profiles of primary infrapatellar fat pad and matched blood samples from 70 osteoarthritis patients undergoing total knee replacement surgery. Comparing the DNA methylation profiles between these tissues reveal widespread epigenetic differences. We produce the first genome-wide methylation quantitative trait locus (mQTL) map of fat pad, and make the resource available to the wider community. Using two-sample Mendelian randomization and colocalization analyses, we resolve osteoarthritis GWAS signals and provide insights into the molecular mechanisms underpinning disease aetiopathology. Our findings provide the first view of the epigenetic landscape of infrapatellar fat pad primary tissue in osteoarthritis.

DNA methylation sites in early adulthood characterised by pubertal timing and development: a twin study

BackgroundPuberty is a highly heritable and variable trait, with environmental factors having a role in its eventual timing and development. Early and late pubertal onset are both associated with various diseases developing later in life, and epigenetic characterisation of pubertal timing and development could lead to important insights. Blood DNA methylation, reacting to both genotype and environment, has been associated with puberty; however, such studies are relatively scarce. We investigated peripheral blood DNA methylation profiles (using Illumina 450 K and EPIC platforms) of 1539 young adult Finnish twins associated with pubertal development scale (PDS) at ages 12 and 14 as well as pubertal age (PA).ResultsFixed effect meta-analysis of the two platforms on 347,521 CpGs in common identified 58 CpG sites associated (p < 1 × 10−5) with either PDS or PA. All four CpGs associated with PA and 45 CpGs associated with PDS were sex-specific. Thirteen CpGs had a high heritability (h2: 0.51–0.98), while one CpG site (mapped to GET4) had a high shared environmental component accounting for 68% of the overall variance in methylation at the site. Utilising twin discordance analysis, we found 6 CpG sites (5 associated with PDS and 1 with PA) that had an environmentally driven association with puberty. Furthermore, genes with PDS- or PA-associated CpGs were consistently linked to various developmental processes and diseases such as breast, prostate and ovarian cancer, while methylation quantitative trait loci of associated CpG sites were enriched in immune pathways developing during puberty.ConclusionsBy identifying puberty-associated DNA methylation sites and examining the effects of sex, environment and genetics, we shed light on the intricate interplay between environment and genetics in the context of puberty. Through our comprehensive analysis, we not only deepen the understanding of the significance of both genetic and environmental factors in the complex processes of puberty and its timing, but also gain insights into potential links with disease risks.

DNA methylation is not associated with sensitization to or dietary introduction of highly allergenic foods in a subset of the CHILD cohort at age 1 year.

In the first year of life, DNA methylation (DNAm) patterns are established and are particularly susceptible to exposure-induced changes. Some of these changes may leave lasting effects by persistently altering gene expression or cell type composition or function, contributing to disease. In this discovery study, we investigated DNAm associations with sensitization to peanut, egg, or cow's milk and hypothesized that genes demonstrating DNAm differences in immune cells may play a role in the development of food sensitization. Infant sensitization (a skin prick test wheal size that is at least 2 mm greater than the negative control) was measured to peanut, egg, and cow's milk at age 1 year, and ages of food introduction were reported prospectively. PBMC DNAm was measured in blood samples at 1 year in 144 infants, oversampled for atopy or wheeze. Statistical analysis of Illumina 450k array DNAm data was conducted in R with adjustment for clinical and genetic covariables and a minimum effect size of 1%, false discovery rate of 5%, and medium-confidence false discovery rate threshold of 20%. There were no DNAm differences between infants with and without peanut, egg, or cow's milk sensitization. Borderline significant sites with high effect sizes were enriched for methylation quantitative trait loci, hinting at genetic factors influencing DNAm at these sites. DNAm patterns did not differ by peanut or egg introduction before or after 12 months. This small pilot study did not show differences inmethylation by food sensitization or introduction, but it diddemonstrate DNAm patterns linked to genetic variants.

Identifying actionable druggable targets for breast cancer: Mendelian randomization and population-based analyses

Drug repurposing provides a cost-effective approach to address the need for breast cancer prevention and therapeutics. We aimed to identify actionable druggable targets using Mendelian randomization (MR) and then validate the candidate drugs using population-based analyses. We identified genetic instruments for 1406 actionable targets of approved non-oncological drugs based on gene expression, DNA methylation, and protein expression quantitative trait loci (eQTL, mQTL, and pQTL, respectively). Genome-wide association study (GWAS) summary statistics were obtained from the Breast Cancer Association Consortium (122,977 cases, 105,974 controls). We further conducted a nested case-control study using data retrieved from Swedish registers to validate the candidate drugs that were identified from MR analyses. We identified six significant MR associations with gene expression levels (TUBB, MDM2, CSK, ULK3, MC1R and KCNN4) and two significant associations with gene methylation levels across 21 CpG islands (RPS23 and MAPT). Results from the nested case-control study showed that the use of raloxifene (targeting MAPT) was associated with 35% reduced breast cancer risk (odds ratio, OR, 0.65; 95% confidence interval, CI, 0.51-0.83). However, usage of estradiol, tolterodine, and nitrofurantoin (also targeting MAPT) was associated with increased breast cancer risk, with adjusted ORs and 95% CI of 1.10 (1.07-1.13), 1.16 (1.09-1.24), and 1.09 (1.05-1.13), respectively. The effect of raloxifene and nitrofurantoin lost significance in further validation analyses using active-comparator and new-user design. This large-scale MR analysis, combined with population-based validation, identified eight druggable target genes for breast cancer and suggested that raloxifene is an effective chemoprevention against breast cancer. Swedish Research Council, Cancerfonden, Crafoordska Stiftelsen, Allmänna Sjukhusets i Malmö Stiftelsen för bekämpande av cancer, 111 Project and MAS cancer.

Integrated Analyses of Single-Cell Transcriptome and Mendelian Randomization Reveal the Protective Role of Resistin in Sepsis Survival in Intensive Care Unit.

The high morbidity and mortality rates associated with sepsis highlight the challenges of finding specific remedies for this condition in the intensive care unit (ICU). This study aimed to explore the differentially expressed genes (DEGs) specific to cell types in sepsis and investigate the role of resistin in the survival of sepsis patients through Mendelian randomization (MR) analyses. We used single-cell and bulk transcriptome data to identify cell type-specific DEGs between sepsis and healthy controls. MR analyses were then conducted to investigate the causal relationships between resistin (one of the identified DEGs) levels and the survival of sepsis patients. Additionally, we utilized meQTL (methylation quantitative trait loci) to identify cytosine-phosphate-guanine (CpG) sites that may directly affect sepsis. We identified 560 cell type-specific DEGs between sepsis and healthy controls. Notably, we observed the upregulation of resistin levels in macrophages during sepsis. In bulk transcriptome, RETN is also upregulated in sepsis samples compared with healthy controls. MR analyses revealed a negative association existed between the expression of resistin, at both gene and protein levels, and the mortality or severity of sepsis patients in ICU. Moreover, there were no associations observed between resistin levels and death or organ failure due to other causes. We also identified three methylation CpG sites, located in RETN or its promoter region-cg06633066, cg22322184, and cg02346997-that directly affected both resistin protein levels and sepsis death in the ICU. Our findings suggest that resistin may provide feasible protection for sepsis patients, particularly those with severe cases, without serious side effects. Therefore, resistin could be a potential drug candidate for sepsis treatment. Additionally, we identified two CpG sites, cg06633066 and cg22322184, that were associated with RETN protein levels and sepsis death, providing novel insights into the underlying mechanisms of sepsis.

Endoplasmic reticulum stress related genome-wide Mendelian randomization identifies therapeutic genes for ulcerative colitis and Crohn's disease.

Background: Endoplasmic reticulum stress (ERS) is an important pathophysiological mechanism in ulcerative colitis (UC) and Crohn's disease (CD). ERS-related genes may be influenced by genetic factors and intestinal inflammation. However, the role of ERS as a trigger or potential etiological factor for UC and CD is unclear, as the expression of ERS-related genes in UC and CD may be the cause or subsequent changes in intestinal inflammation. Here, we used a three-step summary data-based Mendelian randomization (SMR) approach integrating multi-omics data to identify putative causal effects of ERS-related genes in UC and CD. Methods: Genome-wide association study (GWAS) summary data for UC (6,968 cases and 20,464 controls) and CD (5,956 cases and 14,927 controls) were extracted as outcome, and DNA methylation quantitative trait loci (mQTL, 1,980 participants) data and expression QTL data (eQTL, 31,684 participants) from the blood were obtained as exposure. The ERS-related genes were extracted from the GeneCards database, and then the GWAS summary data were integrated with the mQTL and eQTL data associated with ERS genes by SMR. Sensitivity analysis included two-sample MR analysis, power calculations, Bayesian co-localization analysis, and phenotype scanning were performed to evaluate the robustness of the results. Results: A total of 1,193 ERS-related genes were obtained. The three-step SMR analysis showed that cg24011261 CpG site regulating GPX1 expression was associated with a low risk of UC, whereas GPX1 expression regulated by a combination of cg05055782, cg24011261, and cg05551922 CpG sites was associated with a low risk of CD. Sensitivity analysis further supports these findings. Conclusion: This multi-omics integration study identifies a causal relationship between the role of ERS in UC and CD and suggests potential new therapeutic targets for clinical practice.