Integrating muti-omics data to identify tissue-specific DNA methylation biomarkers for cancer risk.

The relationship between tissue-specific DNA methylation and cancer risk remains inadequately elucidated. Leveraging resources from the Genotype-Tissue Expression consortium, here we develop genetic models to predict DNA methylation at CpG sites across the genome for seven tissues and apply these models to genome-wide association study data of corresponding cancers, namely breast, colorectal, renal cell, lung, ovarian, prostate, and testicular germ cell cancers. At Bonferroni-corrected P < 0.05, we identify 4248 CpGs that are significantly associated with cancer risk, of which 95.4% (4052) are specific to a particular cancer type. Notably, 92 CpGs within 55 putative novel loci retain significant associations with cancer risk after conditioning on proximal signals identified by genome-wide association studies. Integrative multi-omics analyses reveal 854 CpG-gene-cancer trios, suggesting that DNA methylation at 309 distinct CpGs might influence cancer risk through regulating the expression of 205 unique cis-genes. These findings substantially advance our understanding of the interplay between genetics, epigenetics, and gene expression in cancer etiology.

Case-Case Genome-Wide Analyses Identify Subtype-Informative Variants That Confer Risk for Breast Cancer.

Breast cancer includes several subtypes with distinct characteristic biological, pathologic, and clinical features. Elucidating subtype-specific genetic etiology could provide insights into the heterogeneity of breast cancer to facilitate the development of improved prevention and treatment approaches. In this study, we conducted pairwise case-case comparisons among five breast cancer subtypes by applying a case-case genome-wide association study (CC-GWAS) approach to summary statistics data of the Breast Cancer Association Consortium. The approach identified 13 statistically significant loci and eight suggestive loci, the majority of which were identified from comparisons between triple-negative breast cancer (TNBC) and luminal A breast cancer. Associations of lead variants in 12 loci remained statistically significant after accounting for previously reported breast cancer susceptibility variants, among which, two were genome-wide significant. Fine mapping implicated putative functional/causal variants and risk genes at several loci, e.g., 3q26.31/TNFSF10, 8q22.3/NACAP1/GRHL2, and 8q23.3/LINC00536/TRPS1, for TNBC as compared with luminal cancer. Functional investigation further identified rs16867605 at 8q22.3 as a SNP that modulates the enhancer activity of GRHL2. Subtype-informative polygenic risk scores (PRS) were derived, and patients with a high subtype-informative PRS had an up to two-fold increased risk of being diagnosed with TNBC instead of luminal cancers. The CC-GWAS PRS remained statistically significant after adjusting for TNBC PRS derived from traditional case-control GWAS in The Cancer Genome Atlas and the African Ancestry Breast Cancer Genetic Consortium. The CC-GWAS PRS was also associated with overall survival and disease-specific survival among patients with breast cancer. Overall, these findings have advanced our understanding of the genetic etiology of breast cancer subtypes, particularly for TNBC. Significance: The discovery of subtype-informative genetic risk variants for breast cancer advances our understanding of the etiologic heterogeneity of breast cancer, which could accelerate the identification of targets and personalized strategies for prevention and treatment.

Large-scale integration of omics and electronic health records to identify potential risk protein biomarkers and therapeutic drugs for cancer prevention and intervention.

Identifying risk protein targets and their therapeutic drugs is crucial for effective cancer prevention. Here, we conduct integrative and fine-mapping analyses of large genome-wide association studies data for breast, colorectal, lung, ovarian, pancreatic, and prostate cancers, and characterize 710 lead variants independently associated with cancer risk. Through mapping protein quantitative trait loci (pQTL) for these variants using plasma proteomics data from over 75,000 participants, we identify 365 proteins associated with cancer risk. Subsequent colocalization analysis identifies 101 proteins, including 74 not reported in previous studies. We further characterize 36 potential druggable proteins for cancers or other disease indications. Analyzing >3.5 million electronic health records, we uncover five drugs (Haloperidol, Trazodone, Tranexamic Acid, Haloperidol, and Captopril) associated with increased cancer risk and two drugs (Caffeine and Acetazolamide) linked to reduced colorectal cancer risk. This study offers novel insights into therapeutic drugs targeting risk proteins for cancer prevention and intervention.

Using genome and transcriptome data from African-ancestry female participants to identify putative breast cancer susceptibility genes.

African-ancestry (AA) participants are underrepresented in genetics research. Here, we conducted a transcriptome-wide association study (TWAS) in AA female participants to identify putative breast cancer susceptibility genes. We built genetic models to predict levels of gene expression, exon junction, and 3' UTR alternative polyadenylation using genomic and transcriptomic data generated in normal breast tissues from 150 AA participants and then used these models to perform association analyses using genomic data from 18,034 cases and 22,104 controls. At Bonferroni-corrected P < 0.05, we identified six genes associated with breast cancer risk, including four genes not previously reported (CTD-3080P12.3, EN1, LINC01956 and NUP210L). Most of these genes showed a stronger association with risk of estrogen-receptor (ER) negative or triple-negative than ER-positive breast cancer. We also replicated the associations with 29 genes reported in previous TWAS at P < 0.05 (one-sided), providing further support for an association of these genes with breast cancer risk. Our study sheds new light on the genetic basis of breast cancer and highlights the value of conducting research in AA populations.

Genome-wide association analyses of breast cancer in women of African ancestry identify new susceptibility loci and improve risk prediction.

We performed genome-wide association studies of breast cancer including 18,034 cases and 22,104 controls of African ancestry. Genetic variants at 12 loci were associated with breast cancer risk (P < 5 × 10-8), including associations of a low-frequency missense variant rs61751053 in ARHGEF38 with overall breast cancer (odds ratio (OR) = 1.48) and a common variant rs76664032 at chromosome 2q14.2 with triple-negative breast cancer (TNBC) (OR = 1.30). Approximately 15.4% of cases with TNBC carried six risk alleles in three genome-wide association study-identified TNBC risk variants, with an OR of 4.21 (95% confidence interval = 2.66-7.03) compared with those carrying fewer than two risk alleles. A polygenic risk score (PRS) showed an area under the receiver operating characteristic curve of 0.60 for the prediction of breast cancer risk, which outperformed PRS derived using data from females of European ancestry. Our study markedly increases the population diversity in genetic studies for breast cancer and demonstrates the utility of PRS for risk prediction in females of African ancestry.

Epigenome-wide association study of lung cancer among never smokers in two prospective cohorts in Shanghai, China.

BACKGROUND: The aetiology of lung cancer among individuals who never smoked remains elusive, despite 15% of lung cancer cases in men and 53% in women worldwide being unrelated to smoking. Epigenetic alterations, particularly DNA methylation (DNAm) changes, have emerged as potential drivers. Yet, few prospective epigenome-wide association studies (EWAS), primarily focusing on peripheral blood DNAm with limited representation of never smokers, have been conducted. METHODS: We conducted a nested case-control study of 80 never-smoking incident lung cancer cases and 83 never-smoking controls within the Shanghai Women's Health Study and Shanghai Men's Health Study. DNAm was measured in prediagnostic oral rinse samples using Illumina MethylationEPIC array. Initially, we conducted an EWAS to identify differentially methylated positions (DMPs) associated with lung cancer in the discovery sample of 101 subjects. The top 50 DMPs were further evaluated in a replication sample of 62 subjects, and results were pooled using fixed-effect meta-analysis. RESULTS: Our study identified three DMPs significantly associated with lung cancer at the epigenome-wide significance level of p<8.22×10-8. These DMPs were identified as cg09198866 (MYH9; TXN2), cg01411366 (SLC9A10) and cg12787323. Furthermore, examination of the top 1000 DMPs indicated significant enrichment in epithelial regulatory regions and their involvement in small GTPase-mediated signal transduction pathways. Additionally, GrimAge acceleration was identified as a risk factor for lung cancer (OR=1.19 per year; 95% CI 1.06 to 1.34). CONCLUSIONS: While replication in a larger sample size is necessary, our findings suggest that DNAm patterns in prediagnostic oral rinse samples could provide novel insights into the underlying mechanisms of lung cancer in never smokers.

Large-Scale Alternative Polyadenylation-Wide Association Studies to Identify Putative Cancer Susceptibility Genes.

Alternative polyadenylation (APA) modulates mRNA processing in the 3'-untranslated regions (3' UTR), affecting mRNA stability and translation efficiency. Research into genetically regulated APA has the potential to provide insights into cancer risk. In this study, we conducted large APA-wide association studies to investigate associations between APA levels and cancer risk. Genetic models were built to predict APA levels in multiple tissues using genotype and RNA sequencing data from 1,337 samples from the Genotype-Tissue Expression project. Associations of genetically predicted APA levels with cancer risk were assessed by applying the prediction models to data from large genome-wide association studies of six common cancers among European ancestry populations: breast, ovarian, prostate, colorectal, lung, and pancreatic cancers. A total of 58 risk genes (corresponding to 76 APA sites) were associated with at least one type of cancer, including 25 genes previously not linked to cancer susceptibility. Of the identified risk APAs, 97.4% and 26.3% were supported by 3'-UTR APA quantitative trait loci and colocalization analyses, respectively. Luciferase reporter assays for four selected putative regulatory 3'-UTR variants demonstrated that the risk alleles of 3'-UTR variants, rs324015 (STAT6), rs2280503 (DIP2B), rs1128450 (FBXO38), and rs145220637 (LDHA), significantly increased the posttranscriptional activities of their target genes compared with reference alleles. Furthermore, knockdown of the target genes confirmed their ability to promote proliferation and migration. Overall, this study provides insights into the role of APA in the genetic susceptibility to common cancers. Significance: Systematic evaluation of associations of alternative polyadenylation with cancer risk reveals 58 putative susceptibility genes, highlighting the contribution of genetically regulated alternative polyadenylation of 3'UTRs to genetic susceptibility to cancer.

Genetic correlation and causal associations between psychiatric disorders and lung cancer risk.

BACKGROUND: Patients with certain psychiatric disorders have increased lung cancer incidence. However, establishing a causal relationship through traditional epidemiological methods poses challenges. METHODS: Available summary statistics of genome-wide association studies of cigarette smoking, lung cancer, and eight psychiatric disorders, including attention deficit/hyperactivity disorder (ADHD), autism, depression, major depressive disorder, bipolar disorder, insomnia, neuroticism, and schizophrenia (range N: 46,350-1,331,010) were leveraged to estimate genetic correlations using Linkage Disequilibrium Score Regression and assess causal effect of each psychiatric disorder on lung cancer using two-sample Mendelian randomization (MR) models, comprising inverse-variance weighted (IVW), weighted median, MR-Egger, pleiotropy residual sum and outlier testing (MR-PRESSO), and a constrained maximum likelihood approach (cML-MR). RESULTS: Significant positive correlations were observed between each psychiatric disorder and both smoking and lung cancer (all FDR < 0.05), except for the correlation between autism and lung cancer. Both univariable and the cML-MA MR analyses demonstrated that liability to schizophrenia, depression, ADHD, or insomnia was associated with an increased risk of overall lung cancer. Genetic liability to insomnia was linked specifically to squamous cell carcinoma (SCC), while genetic liability to ADHD was associated with an elevated risk of both SCC and small cell lung cancer (all P < 0.05). The later was further supported by multivariable MR analyses, which accounted for smoking. LIMITATIONS: Participants were constrained to European ancestry populations. Causal estimates from binary psychiatric disorders may be biased. CONCLUSION: Our findings suggest appropriate management of several psychiatric disorders, particularly ADHD, may potentially reduce the risk of developing lung cancer.

Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes.

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

Drug-target Mendelian randomization revealed a significant association of genetically proxied metformin effects with increased prostate cancer risk.

The association between metformin use and risk of prostate cancer remains controversial, while data from randomized trials is lacking. We aim to evaluate the association of genetically proxied metformin effects with prostate cancer risk using a drug-target Mendelian randomization (MR) approach. Summary statistics for prostate cancer were obtained from the Prostate Cancer Association Group to Investigate Cancer Associated Alterations in the Genome Consortium (79,148 cases and 61,106 controls). Cis-expression quantitative trait loci (cis-eQTL) variants in the gene targets of metformin were identified in the GTEx project and eQTLGen consortium. We also obtained male-specific genome-wide association study data for type 2 diabetes, body mass index (BMI), total testosterone, bioavailable testosterone, estradiol, and sex hormone binding globulin for mediation analysis. Inverse-variance weighted (IVW) regression, weighted median, MR-Egger regression, and MR-PRESSO were performed in the main MR analysis. Multivariable MR was used to identify potential mediators and genetic colocalization analysis was performed to assess any shared genetic basis between two traits of interest. We found that genetically proxied metformin effects (1-SD HbA1c reduction, equivalent to 6.75 mmol/mol) were associated with higher risk of prostate cancer (odds ratioIVW [ORIVW]: 1.55, 95% confidence interval, CI: 1.23-1.96, p = 3.0 × 10-3). Two metformin targets, mitochondrial complex I (ORIVW: 1.48, 95% CI: 1.07-2.03, p = 0.016) and gamma-secretase complex (ORIVW: 2.58, 95%CI :1.47-4.55, p = 0.001), showed robust associations with prostate cancer risk, and their effects were partly mediated through BMI (16.4%) and total testosterone levels (34.3%), respectively. These results were further supported by colocalization analysis that expressions of NDUFA13 and BMI, APH1A, and total testosterone may be influenced by shared genetic factors, respectively. In summary, our study indicated that genetically proxied metformin effects may be associated with an increased risk of prostate cancer. Repurposing metformin for prostate cancer prevention in general populations is not supported by our findings.

Gut Microbiome of Patients With Breast Cancer in Vietnam.

PURPOSE: Gut microbiota play an important role in human health, including cancer. Cancer and its treatment, in turn, may alter the gut microbiome. To understand this complex relationship, we profiled the gut microbiome of 356 Vietnamese patients with breast cancer. MATERIALS AND METHODS: Stool samples were collected before chemotherapy, with 162 pre- and 194 postsurgery. The gut microbiome was measured by shotgun metagenomic sequencing. Associations of gut microbial diversity, taxa abundance, and gut microbiome health index (GMHI) with sociodemographic, clinical factors, and tumor characteristics were evaluated. RESULTS: Postsurgery samples were associated with significantly lower α- and ß-diversities (P < .001) and showed significant differences in the abundance of 15% of 2,864 investigated taxa (false discovery rate [FDR] <0.1) compared with presurgery samples. An unhealthy gut microbiome was prevalent among patients with breast cancer, with a mean GMHI of -0.79 and -2.81 in pre- and postsurgery stool samples, respectively. In an analysis of 162 presurgery stool samples, diagnosis delay was significantly associated with lower α-diversity, variation in ß-diversity, an increased abundance of species Enorma massiliensis, and a decreased abundance of Faecalicoccus pleomorphus. High intake of fiber was significantly associated with lower α-diversity and a higher abundance of species belonging to Bifidobacterium, Prevotella, and Bacteroides gena (FDR < 0.1). We did not find that cancer stage and subtype, menopausal status, comorbidity, antibiotic use during 3 months before stool collection, or physical activity was significantly associated with α- and ß-diversities or GMHI although a few significant differences were observed in taxa abundance. CONCLUSION: Our study revealed that diagnosis delay, high fiber intake, and breast cancer surgery, which is always followed by antibiotic prophylaxis in Vietnam, led to a less diverse and unhealthy gut microbiome among patients with breast cancer.

Gut microbiome in association with chemotherapy-induced toxicities among patients with breast cancer.

BACKGROUND: Little research has focused on the relationship between gut microbiome and chemotherapy-induced toxicity. METHODS: This prospective study involves 301 patients with breast cancer who had prechemotherapy stool samples collected. Gut microbiome was sequenced by shotgun metagenomics; associations with chemotherapy-induced toxicities during first-line treatment by gut microbial diversity, composition, and metabolic pathways with severe (i.e., grade ≥3) hematological and gastrointestinal toxicities were evaluated via multivariable logistic regression. RESULTS: High prechemotherapy α-diversity was associated with a significantly reduced risk of both severe hematological toxicity (odds ratio [OR] = 0.94; 95% CI, 0.89-0.99; p = .048) and neutropenia (OR = 0.94; 95% CI, 0.89-0.99; p = .016). A high abundance of phylum Synergistota, class Synergistia, and order Synergistales were significantly associated with a reduced risk of severe neutropenia; conversely, enrichment of phylum Firmicutes C, class Negativicutes, phylum Firmicutes I, and class Bacilli A, order Paenibacillales were significantly associated with an increased risk of severe neutropenia (p range: 0.012-2.32 × 10-3; false discovery rate <0.1). Significant positive associations were also observed between severe nausea/vomiting and high Chao1 indexes, ß-diversity (p < .05), 20 species belonging to the family Lachnospiraceae, Oscillospiraceae, and Ruminococcaceae (p value range: 6.14 × 10-3 to 1.33 × 10-5; false discovery rate <0.1), and three metabolic pathways involved in reductive tricarboxylic acid cycle I and cycle II, and an incomplete reductive tricarboxylic acid cycle (p < .01). Conversely, a high abundance of species Odoribacter laneus and the pathway related to the L-proline biosynthesis II were inversely associated with severe nausea/vomiting. CONCLUSIONS: Our study suggests that gut microbiota may be a potential preventive target to reduce chemotherapy-induced toxicity.

Genetic correlation and causal associations between circulating C-reactive protein levels and lung cancer risk.

PURPOSE: We aimed to characterize genetic correlations and causal associations between circulating C-reactive protein (CRP) levels and the risk of lung cancer (LC). METHODS: Leveraging summary statistics from genome-wide association studies of circulating CRP levels among 575,531 individuals of European ancestry, and LC risk among 29,266 cases and 56,450 controls, we investigated genetic associations of circulating CRP levels with the risk of overall lung cancer and its histological subtypes, by using linkage disequilibrium score (LDSC) regression and Mendelian randomization (MR) analyses. RESULTS: Significant positive genetic correlations between circulating CRP levels and the risk of LC and its histological subtypes were identified from LDSC regression, with correlation coefficients ranging from 0.12 to 0.26, and all false discovery adjusted p < 0.05. Univariable MR demonstrated a nominal association between CRP levels and an increased risk of lung squamous cell carcinoma (SCC) (inverse variance-weighted OR = 1.15, 95% CI 1.01-1.30). However, this association disappeared when multivariable MR included cigarettes per day and/or body mass index. By using our recently developed constrained maximum likelihood-based MR method, we identified significant associations of CRP levels with the risk of overall LC (OR 1.06, 95% CI 1.03-1.09), SCC (OR 1.06, 95% CI 1.02-1.09), and small cell lung cancer (SCLC, OR 1.09, 95% CI 1.03-1.15). Moreover, most univariable and multivariable MR analyses also revealed consistent CRP-SCLC associations. CONCLUSION: There may be a genetic and causal association between circulating CRP levels and the risk of SCLC, which is in line with previous population-based observational studies.

Novel breast cancer susceptibility loci under linkage peaks identified in African ancestry consortia.

BACKGROUND: Expansion of genome-wide association studies across population groups is needed to improve our understanding of shared and unique genetic contributions to breast cancer. We performed association and replication studies guided by a priori linkage findings from African ancestry (AA) relative pairs. METHODS: We performed fixed-effect inverse-variance weighted meta-analysis under three significant AA breast cancer linkage peaks (3q26-27, 12q22-23, and 16q21-22) in 9241 AA cases and 10 193 AA controls. We examined associations with overall breast cancer as well as estrogen receptor (ER)-positive and negative subtypes (193,132 SNPs). We replicated associations in the African-ancestry Breast Cancer Genetic Consortium (AABCG). RESULTS: In AA women, we identified two associations on chr12q for overall breast cancer (rs1420647, OR = 1.15, p = 2.50×10-6; rs12322371, OR = 1.14, p = 3.15×10-6), and one for ER-negative breast cancer (rs77006600, OR = 1.67, p = 3.51×10-6). On chr3, we identified two associations with ER-negative disease (rs184090918, OR = 3.70, p = 1.23×10-5; rs76959804, OR = 3.57, p = 1.77×10-5) and on chr16q we identified an association with ER-negative disease (rs34147411, OR = 1.62, p = 8.82×10-6). In the replication study, the chr3 associations were significant and effect sizes were larger (rs184090918, OR: 6.66, 95% CI: 1.43, 31.01; rs76959804, OR: 5.24, 95% CI: 1.70, 16.16). CONCLUSION: The two chr3 SNPs are upstream to open chromatin ENSR00000710716, a regulatory feature that is actively regulated in mammary tissues, providing evidence that variants in this chr3 region may have a regulatory role in our target organ. Our study provides support for breast cancer variant discovery using prioritization based on linkage evidence.

Novel insights into genetic susceptibility for colorectal cancer from transcriptome-wide association and functional investigation.

BACKGROUND: Transcriptome-wide association studies have been successful in identifying candidate susceptibility genes for colorectal cancer (CRC). To strengthen susceptibility gene discovery, we conducted a large transcriptome-wide association study and an alternative splicing transcriptome-wide association study in CRC using improved genetic prediction models and performed in-depth functional investigations. METHODS: We analyzed RNA-sequencing data from normal colon tissues and genotype data from 423 European descendants to build genetic prediction models of gene expression and alternative splicing and evaluated model performance using independent RNA-sequencing data from normal colon tissues of the Genotype-Tissue Expression Project. We applied the verified models to genome-wide association studies (GWAS) summary statistics among 58 131 CRC cases and 67 347 controls of European ancestry to evaluate associations of genetically predicted gene expression and alternative splicing with CRC risk. We performed in vitro functional assays for 3 selected genes in multiple CRC cell lines. RESULTS: We identified 57 putative CRC susceptibility genes, which included the 48 genes from transcriptome-wide association studies and 15 genes from splicing transcriptome-wide association studies, at a Bonferroni-corrected P value less than .05. Of these, 16 genes were not previously implicated in CRC susceptibility, including a gene PDE7B (6q23.3) at locus previously not reported by CRC GWAS. Gene knockdown experiments confirmed the oncogenic roles for 2 unreported genes, TRPS1 and METRNL, and a recently reported gene, C14orf166. CONCLUSION: This study discovered new putative susceptibility genes of CRC and provided novel insights into the biological mechanisms underlying CRC development.

A large-scale microRNA transcriptome-wide association study identifies two susceptibility microRNAs, miR-1307-5p and miR-192-3p, for colorectal cancer risk.

Transcriptome-wide association studies (TWAS) have identified many putative susceptibility genes for colorectal cancer (CRC) risk. However, susceptibility miRNAs, critical dysregulators of gene expression, remain unexplored. We genotyped DNA samples from 313 CRC East Asian patients and performed small RNA sequencing in their normal colon tissues distant from tumors to build genetic models for predicting miRNA expression. We applied these models and data from genome-wide association studies (GWAS) including 23 942 cases and 217 267 controls of East Asian ancestry to investigate associations of predicted miRNA expression with CRC risk. Perturbation experiments separately by promoting and inhibiting miRNAs expressions and further in vitro assays in both SW480 and HCT116 cells were conducted. At a Bonferroni-corrected threshold of P < 4.5 × 10-4, we identified two putative susceptibility miRNAs, miR-1307-5p and miR-192-3p, located in regions more than 500 kb away from any GWAS-identified risk variants in CRC. We observed that a high predicted expression of miR-1307-5p was associated with increased CRC risk, while a low predicted expression of miR-192-3p was associated with increased CRC risk. Our experimental results further provide strong evidence of their susceptible roles by showing that miR-1307-5p and miR-192-3p play a regulatory role, respectively, in promoting and inhibiting CRC cell proliferation, migration, and invasion, which was consistently observed in both SW480 and HCT116 cells. Our study provides additional insights into the biological mechanisms underlying CRC development.

Identification of candidate DNA methylation biomarkers related to Alzheimer's disease risk by integrating genome and blood methylome data.

Alzheimer disease (AD) is a common neurodegenerative disease with a late onset. It is critical to identify novel blood-based DNA methylation biomarkers to better understand the extent of the molecular pathways affected in AD. Two sets of blood DNA methylation genetic prediction models developed using different reference panels and modelling strategies were leveraged to evaluate associations of genetically predicted DNA methylation levels with AD risk in 111,326 (46,828 proxy) cases and 677,663 controls. A total of 1,168 cytosine-phosphate-guanine (CpG) sites showed a significant association with AD risk at a false discovery rate (FDR) < 0.05. Methylation levels of 196 CpG sites were correlated with expression levels of 130 adjacent genes in blood. Overall, 52 CpG sites of 32 genes showed consistent association directions for the methylation-gene expression-AD risk, including nine genes (CNIH4, THUMPD3, SERPINB9, MTUS1, CISD1, FRAT2, CCDC88B, FES, and SSH2) firstly reported as AD risk genes. Nine of 32 genes were enriched in dementia and AD disease categories (P values ranged from 1.85 × 10-4 to 7.46 × 10-6), and 19 genes in a neurological disease network (score = 54) were also observed. Our findings improve the understanding of genetics and etiology for AD.

Large-scale alternative polyadenylation (APA)-wide association studies to identify putative susceptibility genes in human common cancers.

Alternative polyadenylation (APA) modulates mRNA processing in the 3' untranslated regions (3'UTR), which affect mRNA stability and translation efficiency. Here, we build genetic models to predict APA levels in multiple tissues using sequencing data of 1,337 samples from the Genotype-Tissue Expression, and apply these models to assess associations between genetically predicted APA levels and cancer risk with data from large genome-wide association studies of six common cancers, including breast, ovary, prostate, colorectum, lung, and pancreas among European-ancestry populations. At a Bonferroni-corrected P □<□0.05, we identify 58 risk genes, including seven in newly identified loci. Using luciferase reporter assays, we demonstrate that risk alleles of 3'UTR variants, rs324015 ( STAT6 ), rs2280503 ( DIP2B ), rs1128450 ( FBXO38 ) and rs145220637 ( LDAH ), could significantly increase post-transcriptional activities of their target genes compared to reference alleles. Further gene knockdown experiments confirm their oncogenic roles. Our study provides additional insight into the genetic susceptibility of these common cancers.

Genome-wide association study of prostate-specific antigen levels in 392,522 men identifies new loci and improves cross-ancestry prediction.

We conducted a multi-ancestry genome-wide association study of prostate-specific antigen (PSA) levels in 296,754 men (211,342 European ancestry; 58,236 African ancestry; 23,546 Hispanic/Latino; 3,630 Asian ancestry; 96.5% of participants were from the Million Veteran Program). We identified 318 independent genome-wide significant (p≤5e-8) variants, 184 of which were novel. Most demonstrated evidence of replication in an independent cohort (n=95,768). Meta-analyzing discovery and replication (n=392,522) identified 447 variants, of which a further 111 were novel. Out-of-sample variance in PSA explained by our new polygenic risk score reached 16.9% (95% CI=16.1%-17.8%) in European ancestry, 9.5% (95% CI=7.0%-12.2%) in African ancestry, 18.6% (95% CI=15.8%-21.4%) in Hispanic/Latino, and 15.3% (95% CI=12.7%-18.1%) in Asian ancestry, and lower for higher age. Our study highlights how including proportionally more participants from underrepresented populations improves genetic prediction of PSA levels, with potential to personalize prostate cancer screening.

Enhancing Disease Risk Gene Discovery by Integrating Transcription Factor-Linked Trans-located Variants into Transcriptome-Wide Association Analyses.

Transcriptome-wide association studies (TWAS) have been successful in identifying putative disease susceptibility genes by integrating gene expression predictions with genome-wide association studies (GWAS) data. However, current TWAS models only consider cis-located variants to predict gene expression. Here, we introduce transTF-TWAS, which includes transcription factor (TF)-linked trans-located variants for model building. Using data from the Genotype-Tissue Expression project, we predict alternative splicing and gene expression and applied these models to large GWAS datasets for breast, prostate, and lung cancers. Our analysis revealed 887 putative cancer susceptibility genes, including 465 in regions not yet reported by previous GWAS and 137 in known GWAS loci but not yet reported previously, at Bonferroni-corrected P < 0.05. We demonstrate that transTF-TWAS surpasses other approaches in both building gene prediction models and identifying disease-associated genes. These results have shed new light on several genetically driven key regulators and their associated regulatory networks underlying disease susceptibility.