An integrative multi-context Mendelian randomization method for identifying risk genes across human tissues.

Mendelian randomization (MR) provides valuable assessments of the causal effect of exposure on outcome, yet the application of conventional MR methods for mapping risk genes encounters new challenges. One of the issues is the limited availability of expression quantitative trait loci (eQTLs) as instrumental variables (IVs), hampering the estimation of sparse causal effects. Additionally, the often context- or tissue-specific eQTL effects challenge the MR assumption of consistent IV effects across eQTL and GWAS data. To address these challenges, we propose a multi-context multivariable integrative MR framework, mintMR, for mapping expression and molecular traits as joint exposures. It models the effects of molecular exposures across multiple tissues in each gene region, while simultaneously estimating across multiple gene regions. It uses eQTLs with consistent effects across more than one tissue type as IVs, improving IV consistency. A major innovation of mintMR involves employing multi-view learning methods to collectively model latent indicators of disease relevance across multiple tissues, molecular traits, and gene regions. The multi-view learning captures the major patterns of disease relevance and uses these patterns to update the estimated tissue relevance probabilities. The proposed mintMR iterates between performing a multi-tissue MR for each gene region and joint learning the disease-relevant tissue probabilities across gene regions, improving the estimation of sparse effects across genes. We apply mintMR to evaluate the causal effects of gene expression and DNA methylation for 35 complex traits using multi-tissue QTLs as IVs. The proposed mintMR controls genome-wide inflation and offers insights into disease mechanisms.

The association of cigarette smoking with DNA methylation and gene expression in human tissue samples.

Cigarette smoking adversely affects many aspects of human health, and epigenetic responses to smoking may reflect mechanisms that mediate or defend against these effects. Prior studies of smoking and DNA methylation (DNAm), typically measured in leukocytes, have identified numerous smoking-associated regions (e.g., AHRR). To identify smoking-associated DNAm features in typically inaccessible tissues, we generated array-based DNAm data for 916 tissue samples from the GTEx (Genotype-Tissue Expression) project representing 9 tissue types (lung, colon, ovary, prostate, blood, breast, testis, kidney, and muscle). We identified 6,350 smoking-associated CpGs in lung tissue (n = 212) and 2,735 in colon tissue (n = 210), most not reported previously. For all 7 other tissue types (sample sizes 38-153), no clear associations were observed (false discovery rate 0.05), but some tissues showed enrichment for smoking-associated CpGs reported previously. For 1,646 loci (in lung) and 22 (in colon), smoking was associated with both DNAm and local gene expression. For loci detected in both lung and colon (e.g., AHRR, CYP1B1, CYP1A1), top CpGs often differed between tissues, but similar clusters of hyper- or hypomethylated CpGs were observed, with hypomethylation at regulatory elements corresponding to increased expression. For lung tissue, 17 hallmark gene sets were enriched for smoking-associated CpGs, including xenobiotic- and cancer-related gene sets. At least four smoking-associated regions in lung were impacted by lung methylation quantitative trait loci (QTLs) that co-localize with genome-wide association study (GWAS) signals for lung function (FEV1/FVC), suggesting epigenetic alterations can mediate the effects of smoking on lung health. Our multi-tissue approach has identified smoking-associated regions in disease-relevant tissues, including effects that are shared across tissue types.

Mapping potential pathways from polygenic liability through brain structure to psychological problems across the transition to adolescence.

BACKGROUND: We used a polygenic score for externalizing behavior (extPGS) and structural MRI to examine potential pathways from genetic liability to conduct problems via the brain across the adolescent transition. METHODS: Three annual assessments of child conduct problems, attention-deficit/hyperactivity problems, and internalizing problems were conducted across across 9-13 years of age among 4,475 children of European ancestry in the Adolescent Brain Cognitive DevelopmentSM Study (ABCD Study®). RESULTS: The extPGS predicted conduct problems in each wave (R2 = 2.0%-2.9%). Bifactor models revealed that the extPRS predicted variance specific to conduct problems (R2 = 1.7%-2.1%), but also variance that conduct problems shared with other measured problems (R2 = .8%-1.4%). Longitudinally, extPGS predicted levels of specific conduct problems (R2 = 2.0%), but not their slope of change across age. The extPGS was associated with total gray matter volume (TGMV; R2 = .4%) and lower TGMV predicted both specific conduct problems (R2 = 1.7%-2.1%) and the variance common to all problems in each wave (R2 = 1.6%-3.1%). A modest proportion of the polygenic liability specific to conduct problems in each wave was statistically mediated by TGMV. CONCLUSIONS: Across the adolescent transition, the extPGS predicted both variance specific to conduct problems and variance shared by all measured problems. The extPGS also was associated with TGMV, which robustly predicted conduct problems. Statistical mediation analyses suggested the hypothesis that polygenic variation influences individual differences in brain development that are related to the likelihood of conduct problems during the adolescent transition, justifying new research to test this causal hypothesis.

Returning personal genetic information on susceptibility to arsenic toxicity to research participants in Bangladesh.

BACKGROUND: There is growing consensus that researchers should offer to return genetic results to participants, but returning results in lower-resource countries has received little attention. In this study, we return results on genetic susceptibility to arsenic toxicity to participants in a Bangladeshi cohort exposed to arsenic through naturally-contaminated drinking water. We examine the impact on behavioral changes related to exposure reduction. METHODS: We enrolled participants from the Health Effects of Arsenic Longitudinal Study who had (1) high arsenic (≥150 µg/g creatinine) in a recent urine sample and (2) existing data on genetic variants impacting arsenic metabolism efficiency (AS3MT and FTCD). We used genetic data to recruit three study groups, each with n = 103: (1) efficient metabolizers (low-risk), (2) inefficient metabolizers (high-risk), and (3) a randomly-selected control group (NCT05072132). At baseline, all participants received information on the effects of arsenic and how to reduce exposure by switching to a low arsenic well. The two intervention groups also received their arsenic metabolism efficiency status (based on their genetic results). Changes in behavior and arsenic exposure were assessed using questionnaires and urine arsenic measures after six months. RESULTS: Clear decreases in urine arsenic after six months were observed for all three groups. The inefficient group self-reported higher levels of attempted switching to lower arsenic wells than the other groups; however, there was no detectable difference in urine arsenic reduction among the three groups. Participants showed strong interest in receiving genetic results and found them useful. The inefficient group experienced higher levels of anxiety than the other groups. Among the efficient group, that receiving genetic results did not appear to hinder behavioral change. CONCLUSION: Returning genetic results increased self-reported exposure-reducing behaviors but did not have a detectable impact on reducing urine arsenic over and above a one-on-one educational intervention.

Integrative cross-omics and cross-context analysis elucidates molecular links underlying genetic effects on complex traits.

Genetic effects on functionally related 'omic' traits often co-occur in relevant cellular contexts, such as tissues. Motivated by the multi-tissue methylation quantitative trait loci (mQTLs) and expression QTLs (eQTLs) analysis, we propose X-ING (Cross-INtegrative Genomics) for cross-omics and cross-context integrative analysis. X-ING takes as input multiple matrices of association statistics, each obtained from different omics data types across multiple cellular contexts. It models the latent binary association status of each statistic, captures the major association patterns among omics data types and contexts, and outputs the posterior mean and probability for each input statistic. X-ING enables the integration of effects from different omics data with varying effect distributions. In the multi-tissue cis-association analysis, X-ING shows improved detection and replication of mQTLs by integrating eQTL maps. In the trans-association analysis, X-ING reveals an enrichment of trans-associations in many disease/trait-relevant tissues.

An integrative multi-context Mendelian randomization method for identifying risk genes across human tissues.

Mendelian randomization (MR) provides valuable assessments of the causal effect of exposure on outcome, yet the application of conventional MR methods for mapping risk genes encounters new challenges. One of the issues is the limited availability of expression quantitative trait loci (eQTLs) as instrumental variables (IVs), hampering the estimation of sparse causal effects. Additionally, the often context/tissue-specific eQTL effects challenge the MR assumption of consistent IV effects across eQTL and GWAS data. To address these challenges, we propose a multi-context multivariable integrative MR framework, mintMR, for mapping expression and molecular traits as joint exposures. It models the effects of molecular exposures across multiple tissues in each gene region, while simultaneously estimating across multiple gene regions. It uses eQTLs with consistent effects across more than one tissue type as IVs, improving IV consistency. A major innovation of mintMR involves employing multi-view learning methods to collectively model latent indicators of disease relevance across multiple tissues, molecular traits, and gene regions. The multi-view learning captures the major patterns of disease-relevance and uses these patterns to update the estimated tissue relevance probabilities. The proposed mintMR iterates between performing a multi-tissue MR for each gene region and joint learning the disease-relevant tissue probabilities across gene regions, improving the estimation of sparse effects across genes. We apply mintMR to evaluate the causal effects of gene expression and DNA methylation for 35 complex traits using multi-tissue QTLs as IVs. The proposed mintMR controls genome-wide inflation and offers new insights into disease mechanisms.

DNA methylation correlates of chronological age in diverse human tissue types.

BACKGROUND: While the association of chronological age with DNA methylation (DNAm) in whole blood has been extensively studied, the tissue-specificity of age-related DNAm changes remains an active area of research. Studies investigating the association of age with DNAm in tissues such as brain, skin, immune cells, fat, and liver have identified tissue-specific and non-specific effects, thus, motivating additional studies of diverse human tissue and cell types. RESULTS: Here, we performed an epigenome-wide association study, leveraging DNAm data (Illumina EPIC array) from 961 tissue samples representing 9 tissue types (breast, lung, colon, ovary, prostate, skeletal muscle, testis, whole blood, and kidney) from the Genotype-Tissue Expression (GTEx) project. We identified age-associated CpG sites (false discovery rate < 0.05) in 8 tissues (all except skeletal muscle, n = 47). This included 162,002 unique hypermethylated and 90,626 hypomethylated CpG sites across all tissue types, with 130,137 (80%) hypermethylated CpGs and 74,703 (82%) hypomethylated CpG sites observed in a single tissue type. While the majority of age-associated CpG sites appeared tissue-specific, the patterns of enrichment among genomic features, such as chromatin states and CpG islands, were similar across most tissues, suggesting common mechanisms underlying cellular aging. Consistent with previous findings, we observed that hypermethylated CpG sites are enriched in regions with repressed polycomb signatures and CpG islands, while hypomethylated CpG sites preferentially occurred in non-CpG islands and enhancers. To gain insights into the functional effects of age-related DNAm changes, we assessed the correlation between DNAm and local gene expression changes to identify age-related expression quantitative trait methylation (age-eQTMs). We identified several age-eQTMs present in multiple tissue-types, including in the CDKN2A, HENMT1, and VCWE regions. CONCLUSION: Overall, our findings will aid future efforts to develop biomarkers of aging and understand mechanisms of aging in diverse human tissue types.

The Impact of Inherited Genetic Variation on DNA Methylation in Prostate Cancer and Benign Tissues of African American and European American Men.

BACKGROUND: American men of African ancestry (AA) have higher prostate cancer incidence and mortality rates compared with American men of European ancestry (EA). Differences in genetic susceptibility mechanisms may contribute to this disparity. METHODS: To gain insights into the regulatory mechanisms of prostate cancer susceptibility variants, we tested the association between SNPs and DNA methylation (DNAm) at nearby CpG sites across the genome in benign and cancer prostate tissue from 74 AA and 74 EA men. Genome-wide SNP data (from benign tissue) and DNAm were generated using Illumina arrays. RESULTS: Among AA men, we identified 6,298 and 2,641 cis-methylation QTLs (meQTL; FDR of 0.05) in benign and tumor tissue, respectively, with 6,960 and 1,700 detected in EA men. We leveraged genome-wide association study (GWAS) summary statistics to identify previously reported prostate cancer GWAS signals likely to share a common causal variant with a detected meQTL. We identified nine GWAS-meQTL pairs with strong evidence of colocalization (four in EA benign, three in EA tumor, two in AA benign, and three in AA tumor). Among these colocalized GWAS-meQTL pairs, we identified colocalizing expression quantitative trait loci (eQTL) impacting four eGenes with known roles in tumorigenesis. CONCLUSIONS: These findings highlight epigenetic regulatory mechanisms by which prostate cancer-risk SNPs can modify local DNAm and/or gene expression in prostate tissue. IMPACT: Overall, our findings showed general consistency in the meQTL landscape of AA and EA men, but meQTLs often differ by tissue type (normal vs. cancer). Ancestry-based linkage disequilibrium differences and lack of AA representation in GWAS decrease statistical power to detect colocalization for some regions.