Genome-wide enhancer RNA profiling adds molecular links between genetic variation and human cancers.

BACKGROUND: Dysregulation of enhancer transcription occurs in multiple cancers. Enhancer RNAs (eRNAs) are transcribed products from enhancers that play critical roles in transcriptional control. Characterizing the genetic basis of eRNA expression may elucidate the molecular mechanisms underlying cancers. METHODS: Initially, a comprehensive analysis of eRNA quantitative trait loci (eRNAQTLs) was performed in The Cancer Genome Atlas (TCGA), and functional features were characterized using multi-omics data. To establish the first eRNAQTL profiles for colorectal cancer (CRC) in China, epigenomic data were used to define active enhancers, which were subsequently integrated with transcription and genotyping data from 154 paired CRC samples. Finally, large-scale case-control studies (34,585 cases and 69,544 controls) were conducted along with multipronged experiments to investigate the potential mechanisms by which candidate eRNAQTLs affect CRC risk. RESULTS: A total of 300,112 eRNAQTLs were identified across 30 different cancer types, which exert their influence on eRNA transcription by modulating chromatin status, binding affinity to transcription factors and RNA-binding proteins. These eRNAQTLs were found to be significantly enriched in cancer risk loci, explaining a substantial proportion of cancer heritability. Additionally, tumor-specific eRNAQTLs exhibited high responsiveness to the development of cancer. Moreover, the target genes of these eRNAs were associated with dysregulated signaling pathways and immune cell infiltration in cancer, highlighting their potential as therapeutic targets. Furthermore, multiple ethnic population studies have confirmed that an eRNAQTL rs3094296-T variant decreases the risk of CRC in populations from China (OR = 0.91, 95%CI 0.88-0.95, P = 2.92 × 10-7) and Europe (OR = 0.92, 95%CI 0.88-0.95, P = 4.61 × 10-6). Mechanistically, rs3094296 had an allele-specific effect on the transcription of the eRNA ENSR00000155786, which functioned as a transcriptional activator promoting the expression of its target gene SENP7. These two genes synergistically suppressed tumor cell proliferation. Our curated list of variants, genes, and drugs has been made available in CancereRNAQTL ( http://canernaqtl.whu.edu.cn/#/ ) to serve as an informative resource for advancing this field. CONCLUSION: Our findings underscore the significance of eRNAQTLs in transcriptional regulation and disease heritability, pinpointing the potential of eRNA-based therapeutic strategies in cancers.

Developing an optimal stratification model for colorectal cancer screening and reducing racial disparities in multi-center population-based studies.

BACKGROUND: Early detection of colorectal neoplasms can reduce the colorectal cancer (CRC) burden by timely intervention for high-risk individuals. However, effective risk prediction models are lacking for personalized CRC early screening in East Asian (EAS) population. We aimed to develop, validate, and optimize a comprehensive risk prediction model across all stages of the dynamic adenoma-carcinoma sequence in EAS population. METHODS: To develop precision risk-stratification and intervention strategies, we developed three trans-ancestry PRSs targeting colorectal neoplasms: (1) using 148 previously identified CRC risk loci (PRS148); (2) SNPs selection from large-scale meta-analysis data by clumping and thresholding (PRS183); (3) PRS-CSx, a Bayesian approach for genome-wide risk prediction (PRSGenomewide). Then, the performance of each PRS was assessed and validated in two independent cross-sectional screening sets, including 4600 patients with advanced colorectal neoplasm, 4495 patients with non-advanced adenoma, and 21,199 normal individuals from the ZJCRC (Zhejiang colorectal cancer set; EAS) and PLCO (the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial; European, EUR) studies. The optimal PRS was further incorporated with lifestyle factors to stratify individual risk and ultimately tested in the PLCO and UK Biobank prospective cohorts, totaling 350,013 participants. RESULTS: Three trans-ancestry PRSs achieved moderately improved predictive performance in EAS compared to EUR populations. Remarkably, the PRSs effectively facilitated a thorough risk assessment across all stages of the dynamic adenoma-carcinoma sequence. Among these models, PRS183 demonstrated the optimal discriminatory ability in both EAS and EUR validation datasets, particularly for individuals at risk of colorectal neoplasms. Using two large-scale and independent prospective cohorts, we further confirmed a significant dose-response effect of PRS183 on incident colorectal neoplasms. Incorporating PRS183 with lifestyle factors into a comprehensive strategy improves risk stratification and discriminatory accuracy compared to using PRS or lifestyle factors separately. This comprehensive risk-stratified model shows potential in addressing missed diagnoses in screening tests (best NPV = 0.93), while moderately reducing unnecessary screening (best PPV = 0.32). CONCLUSIONS: Our comprehensive risk-stratified model in population-based CRC screening trials represents a promising advancement in personalized risk assessment, facilitating tailored CRC screening in the EAS population. This approach enhances the transferability of PRSs across ancestries and thereby helps address health disparity.

The global distribution and the risk prediction of relapsing fever group Borrelia: a data review with modelling analysis.

BACKGROUND: The recent discovery of emerging relapsing fever group Borrelia (RFGB) species, such as Borrelia miyamotoi, poses a growing threat to public health. However, the global distribution and associated risk burden of these species remain uncertain. We aimed to map the diversity, distribution, and potential infection risk of RFGB. METHODS: We searched PubMed, Web of Science, GenBank, CNKI, and eLibrary from Jan 1, 1874, to Dec 31, 2022, for published articles without language restriction to extract distribution data for RFGB detection in vectors, animals, and humans, and clinical information about human patients. Only articles documenting RFGB infection events were included in this study, and data for RFGB detection in vectors, animals, or humans were composed into a dataset. We used three machine learning algorithms (boosted regression trees, random forest, and least absolute shrinkage and selection operator logistic regression) to assess the environmental, ecoclimatic, biological, and socioeconomic factors associated with the occurrence of four major RFGB species: Borrelia miyamotoi, Borrelia lonestari, Borrelia crocidurae, and Borrelia hermsii; and mapped their worldwide risk level. FINDINGS: We retrieved 13 959 unique studies, among which 697 met the selection criteria and were used for data extraction. 29 RFGB species have been recorded worldwide, of which 27 have been identified from 63 tick species, 12 from 61 wild animals, and ten from domestic animals. 16 RFGB species caused human infection, with a cumulative count of 26 583 cases reported from Jan 1, 1874, to Dec 31, 2022. Borrelia recurrentis (17 084 cases) and Borrelia persica (2045 cases) accounted for the highest proportion of human infection. B miyamotoi showed the widest distribution among all RFGB, with a predicted environmentally suitable area of 6·92 million km2, followed by B lonestari (1·69 million km2), B crocidurae (1·67 million km2), and B hermsii (1·48 million km2). The habitat suitability index of vector ticks and climatic factors, such as the annual mean temperature, have the most significant effect among all predictive models for the geographical distribution of the four major RFGB species. INTERPRETATION: The predicted high-risk regions are considerably larger than in previous reports. Identification, surveillance, and diagnosis of RFGB infections should be prioritised in high-risk areas, especially within low-income regions. FUNDING: National Key Research and Development Program of China.

A genetic variant in the immune-related gene ERAP1 affects colorectal cancer prognosis.

BACKGROUND: Findings on the association of genetic factors and colorectal cancer (CRC) survival are limited and inconsistent, and revealing the mechanism underlying their prognostic roles is of great importance. This study aimed to explore the relationship between functional genetic variations and the prognosis of CRC and further reveal the possible mechanism. METHODS: We first systematically performed expression quantitative trait locus (eQTL) analysis using The Cancer Genome Atlas (TCGA) dataset. Then, the Kaplan-Meier analysis was used to filter out the survival-related eQTL target genes of CRC patients in two public datasets (TCGA and GSE39582 dataset from the Gene Expression Omnibus database). The seven most potentially functional eQTL single nucleotide polymorphisms (SNPs) associated with six survival-related eQTL target genes were genotyped in 907 Chinese CRC patients with clinical prognosis data. The regulatory mechanism of the survival-related SNP was further confirmed by functional experiments. RESULTS: The rs71630754 regulating the expression of endoplasmic reticulum aminopeptidase 1 ( ERAP1 ) was significantly associated with the prognosis of CRC (additive model, hazard ratio [HR]: 1.43, 95% confidence interval [CI]: 1.08-1.88, P = 0.012). The results of dual-luciferase reporter assay and electrophoretic mobility shift assay showed that the A allele of the rs71630754 could increase the binding of transcription factor 3 (TCF3) and subsequently reduce the expression of ERAP1 . The results of bioinformatic analysis showed that lower expression of ERAP1 could affect the tumor immune microenvironment and was significantly associated with severe survival outcomes. CONCLUSION: The rs71630754 could influence the prognosis of CRC patients by regulating the expression of the immune-related gene ERAP1 . TRIAL REGISTRATION: No. NCT00454519 ( https://clinicaltrials.gov/ ).

Prioritization of risk genes in colorectal cancer by integrative analysis of multi-omics data and gene networks.

Genome-wide association studies (GWASs) have identified over 140 colorectal cancer (CRC)-associated loci; however, target genes at the majority of loci and underlying molecular mechanisms are poorly understood. Here, we utilized a Bayesian approach, integrative risk gene selector (iRIGS), to prioritize risk genes at CRC GWAS loci by integrating multi-omics data. As a result, a total of 105 high-confidence risk genes (HRGs) were identified, which exhibited strong gene dependencies for CRC and enrichment in the biological processes implicated in CRC. Among the 105 HRGs, CEBPB, located at the 20q13.13 locus, acted as a transcription factor playing critical roles in cancer. Our subsequent assays indicated the tumor promoter function of CEBPB that facilitated CRC cell proliferation by regulating multiple oncogenic pathways such as MAPK, PI3K-Akt, and Ras signaling. Next, by integrating a fine-mapping analysis and three independent case-control studies in Chinese populations consisting of 8,039 cases and 12,775 controls, we elucidated that rs1810503, a putative functional variant regulating CEBPB, was associated with CRC risk (OR=0.90, 95%CI=0.86-0.93, P=1.07×10-7). The association between rs1810503 and CRC risk was further validated in three additional multi-ancestry populations consisting of 24,254 cases and 58,741 controls. Mechanistically, the rs1810503 A to T allele change weakened the enhancer activity in an allele-specific manner to decrease CEBPB expression via long-range promoter-enhancer interactions, mediated by the transcription factor, REST, and thus decreased CRC risk. In summary, our study provides a genetic resource and a generalizable strategy for CRC etiology investigation, and highlights the biological implications of CEBPB in CRC tumorigenesis, shedding new light on the etiology of CRC.

Healthy dietary patterns, genetic risk, and gastrointestinal cancer incident risk: a large-scale prospective cohort study.

BACKGROUND: Dietary patterns have been associated with several cancers, especially gastrointestinal cancer (GIC). However, whether a healthy dietary pattern could modify the risk of GIC among people with different genetic backgrounds is not clear. OBJECTIVE: The objective of the study was to investigate how dietary patterns and genetic susceptibility contribute to the risk of GIC independently and jointly. METHODS: This large-scale prospective cohort study included 105,463 participants in UK Biobank who were aged 40-72 y and cancer-free at baseline. Dietary intake (Oxford WebQ) was used to calculate dietary pattern scores including dietary approach to stop hypertension (DASH) score and healthful plant-based diet index (hPDI). Genetic risk was quantified by a polygenic risk score (PRS) comprising 129 known GIC-associated loci. Cox proportional hazards regression was performed to estimate the associations of dietary patterns and PRS with GIC incidence after adjusting for potential confounders. RESULTS: Over a median follow-up of 11.70 y, 1,661 participants were diagnosed with GIC. DASH and hPDI were associated with 20% and 36% reductions, respectively, in GIC risk. Low PRS was associated with a 30 % decrease in GIC risk (HR: 0.70; 95% CI: 0.62, 0.79). Participants with healthy dietary scores at high-genetic risk had a lower GIC risk with HR of 0.77 (95% CI: 0.60, 0.98) for DASH and 0.66 (95% CI: 0.52, 0.84) for hPDI than those with unhealthy dietary score. Participants with both high-dietary score and low-genetic risk showed the lowest risk of GIC, with HR of 0.58 (95% CI: 0.45, 0.75) for DASH and 0.45 (95% CI: 0.34, 0.58) for hPDI. CONCLUSIONS: Adherence to DASH and hPDI were associated with a lower risk of some gastrointestinal cancers, and these 2 dietary patterns may partly compensate for genetic predispositions to cancer. Our results advance the development of precision medicine strategies that consider both dietary patterns and genetics to improve gastrointestinal health.

Genome-wide association analysis of left ventricular imaging-derived phenotypes identifies 72 risk loci and yields genetic insights into hypertrophic cardiomyopathy.

Left ventricular regional wall thickness (LVRWT) is an independent predictor of morbidity and mortality in cardiovascular diseases (CVDs). To identify specific genetic influences on individual LVRWT, we established a novel deep learning algorithm to calculate 12 LVRWTs accurately in 42,194 individuals from the UK Biobank with cardiac magnetic resonance (CMR) imaging. Genome-wide association studies of CMR-derived 12 LVRWTs identified 72 significant genetic loci associated with at least one LVRWT phenotype (P < 5 × 10-8), which were revealed to actively participate in heart development and contraction pathways. Significant causal relationships were observed between the LVRWT traits and hypertrophic cardiomyopathy (HCM) using genetic correlation and Mendelian randomization analyses (P < 0.01). The polygenic risk score of inferoseptal LVRWT at end systole exhibited a notable association with incident HCM, facilitating the identification of high-risk individuals. The findings yield insights into the genetic determinants of LVRWT phenotypes and shed light on the biological basis for HCM etiology.

Genetic Variants That Impact Alternative Polyadenylation in Cancer Represent Candidate Causal Risk Loci.

Alternative polyadenylation (APA) is emerging as a major mechanism of posttranscriptional regulation. APA can impact the development and progression of cancer, suggesting that the genetic determinants of APA might play an important role in regulating cancer risk. Here, we depicted a pan-cancer atlas of human APA quantitative trait loci (apaQTL), containing approximately 0.7 million apaQTLs across 32 cancer types. Systematic multiomics analyses indicated that cancer apaQTLs could contribute to APA regulation by altering poly(A) motifs, RNA-binding proteins (RBP), and chromatin regulatory elements and were preferentially enriched in genome-wide association studies (GWAS)-identified cancer susceptibility loci. Moreover, apaQTL-related genes (aGene) were broadly related to cancer signaling pathways, high mutational burden, immune infiltration, and drug response, implicating their potential as therapeutic targets. Furthermore, apaQTLs were mapped in Chinese colorectal cancer tumor tissues and then screened for functional apaQTLs associated with colorectal cancer risk in 17,789 cases and 19,951 controls using GWAS-ChIP data, with independent validation in a large-scale population consisting of 6,024 cases and 10,022 controls. A multi-ancestry-associated apaQTL variant rs1020670 with a C>G change in DNM1L was identified, and the G allele contributed to an increased risk of colorectal cancer. Mechanistically, the risk variant promoted aberrant APA and facilitated higher usage of DNM1L proximal poly(A) sites mediated by the RBP CSTF2T, which led to higher expression of DNM1L with a short 3'UTR. This stabilized DNM1L to upregulate its expression, provoking colorectal cancer cell proliferation. Collectively, these findings generate a resource for understanding APA regulation and the genetic basis of human cancers, providing insights into cancer etiology. SIGNIFICANCE: Cancer risk is mediated by alternative polyadenylation quantitative trait loci, including the rs1020670-G variant that promotes alternative polyadenylation of DNM1L and increases colorectal cancer risk.

Genome-wide enhancer-gene regulatory maps link causal variants to target genes underlying human cancer risk.

Genome-wide association studies have identified numerous variants associated with human complex traits, most of which reside in the non-coding regions, but biological mechanisms remain unclear. However, assigning function to the non-coding elements is still challenging. Here we apply Activity-by-Contact (ABC) model to evaluate enhancer-gene regulation effect by integrating multi-omics data and identified 544,849 connections across 20 cancer types. ABC model outperforms previous approaches in linking regulatory variants to target genes. Furthermore, we identify over 30,000 enhancer-gene connections in colorectal cancer (CRC) tissues. By integrating large-scale population cohorts (23,813 cases and 29,973 controls) and multipronged functional assays, we demonstrate an ABC regulatory variant rs4810856 associated with CRC risk (Odds Ratio = 1.11, 95%CI = 1.05-1.16, P = 4.02 × 10-5) by acting as an allele-specific enhancer to distally facilitate PREX1, CSE1L and STAU1 expression, which synergistically activate p-AKT signaling. Our study provides comprehensive regulation maps and illuminates a single variant regulating multiple genes, providing insights into cancer etiology.

Risk SNP in a transcript of RP11-638I2.4 increases lncRNA-YY1 interaction and pancreatic cancer susceptibility.

Tens of thousands of long non-coding RNAs (lncRNAs) have been identified through RNA-seq analysis, but the biological and pathological significance remains unclear. By integrating the genome-wide lncRNA data with a cross-ancestry meta-analysis of PDAC GWASs, we depicted a comprehensive atlas of pancreatic ductal adenocarcinoma (PDAC)-associated lncRNAs, containing 1,204 lncRNA (445 novel lncRNAs and 759 GENCODE annotated lncRNAs) and 4,368 variants. Furthermore, we found that PDAC-associated lncRNAs could function by altering chromatin activity, transcription factors, and RNA-binding proteins binding affinity. Importantly, genetic variants linked to PDAC are preferentially found at PDAC-associated lncRNA regions, supporting the biological and clinical relevance of PDAC-associated lncRNAs. Finally, we prioritized a novel transcript (MICT00000110172.1) of RP11-638I2.4 as a potential tumor promoter. MICT00000110172.1 is able to reinforce the interaction with YY1, which could reverse the effect of YY1 on pancreatic cancer cell cycle arrest to promote the pancreatic cancer growth. G > A change at rs2757535 in the second exon of MICT00000110172.1 induces a spatial structural change and creates a target region for YY1 binding, which enforces the effect of MICT00000110172.1 in an allele-specific manner, and thus confers susceptibility to tumorigenesis. In summary, our results extend the repertoire of PDAC-associated lncRNAs that could act as a starting point for future functional explorations, and the identification of lncRNA-based target therapy.

Genetic Control of Alternative Splicing and its Distinct Role in Colorectal Cancer Mechanisms.

BACKGROUND & AIMS: Dysregulation of alternative splicing is implicated in many human diseases, and understanding the genetic variation underlying transcript splicing is essential to dissect the molecular mechanisms of cancers. We aimed to provide a comprehensive functional dissection of splicing quantitative trait loci (sQTLs) in cancer and focus on elucidating its distinct role in colorectal cancer (CRC) mechanisms. METHODS: We performed a comprehensive sQTL analysis to identify genetic variants that control messenger RNA splicing across 33 cancer types from The Cancer Genome Atlas and independently validated in our 154 CRC tissues. Then, large-scale, multicenter, multi-ethnic case-control studies (34,585 cases and 76,023 controls) were conducted to examine the association of these sQTLs with CRC risk. A series of biological experiments in vitro and in vivo were performed to investigate the potential mechanisms of the candidate sQTLs and target genes. RESULTS: The molecular characterization of sQTL revealed its distinct role in cancer susceptibility. Tumor-specific sQTL further showed better response to cancer development. In addition, functionally informed polygenic risk score highlighted the potentiality of sQTLs in the CRC prediction. Complemented by large-scale population studies, we identified that the risk allele (T) of a multi-ancestry-associated sQTL rs61746794 significantly increased the risk of CRC in Chinese (odds ratio, 1.20; 95% CI, 1.12-1.29; P = 8.82 × 10-7) and European (odds ratio, 1.11; 95% CI, 1.07-1.16; P = 1.13 × 10-7) populations. rs61746794-T facilitated PRMT7 exon 16 splicing mediated by the RNA-binding protein PRPF8, thus increasing the level of canonical PRMT7 isoform (PRMT7-V2). Overexpression of PRMT7-V2 significantly enhanced the growth of CRC cells and xenograft tumors compared with PRMT7-V1. Mechanistically, PRMT7-V2 functions as an epigenetic writer that catalyzes the arginine methylation of H4R3 and H3R2, subsequently regulating diverse biological processes, including YAP, AKT, and KRAS pathway. A selective PRMT7 inhibitor, SGC3027, exhibited antitumor effects on human CRC cells. CONCLUSIONS: Our study provides an informative sQTLs resource and insights into the regulatory mechanisms linking splicing variants to cancer risk and serving as biomarkers and therapeutic targets.

In Utero and Childhood/Adolescence Exposure to Tobacco Smoke, Genetic Risk, and Cancer Incidence in Adulthood: A Prospective Cohort Study.

OBJECTIVE: To evaluate the association of early-life tobacco smoke exposure, especially interacting with cancer genetic variants, with adult cancer. PARTICIPANTS AND METHODS: We examined the associations of in utero tobacco smoke exposure, age of smoking initiation, and their interaction with genetic risk levels with cancer incidence in 393,081 participants from the UK Biobank. Information on tobacco exposure was obtained by self-reported questionnaires. A cancer polygenic risk score was constructed by weighting and integrating 702 genome-wide association studies-identified risk variants. Cox proportional hazards regression models were used to calculate hazard ratios (HRs) for overall cancer and organ-specific cancer incidence. RESULTS: During 11.8 years of follow-up, 23,450 (5.97%) and 23,413 (6.03%) incident cancers were included in the analyses of in utero exposure and age of smoking initiation, respectively. The HR (95% CI) for incident cancer in participants with in utero exposure to tobacco smoke was 1.04 (1.01-1.07) for overall cancer, 1.59 (1.44-1.75) for respiratory cancer, and 1.09 (1.03-1.17) for gastrointestinal cancer. The relative risk of incident cancer increased with earlier smoking initiation (Ptrend<.001), with the HR (95% CI) of 1.44 (1.36-1.51) for overall cancer, 13.28 (11.39-15.48) for respiratory cancer, and 1.72 (1.54-1.91) for gastrointestinal cancer in smokers with initiation in childhood compared with never smokers. Importantly, a positive additive interaction between age of smoking initiation and genetic risk was observed for overall cancer (Padditive=.04) and respiratory cancer (Padditive=.003) incidence. CONCLUSION: In utero exposure and earlier smoking initiation are associated with overall and organ-specific cancer, and age of smoking initiation interaction with genetic risk is associated with respiratory cancer.

Interaction of metabolism-related pathway gene variants with bisphenol A exposure on serum lipid profiles.

Bisphenol A (BPA) can be metabolized by metabolic enzymes and may induce abnormal lipid metabolism. We hypothesized that BPA exposure and its interaction with metabolism-related genes might be associated with serum lipid profiles. We performed a two-stage study among 955 middle-aged and elderly participants in Wuhan, China. Urinary BPA level was estimated without (BPA, µg/L) or with (BPA/Cr, µg/g) adjustments for urinary creatinine and ln-transformed values (ln-BPA or ln-BPA/Cr) were used to normalize the asymmetrical distributions. A total of 412 metabolism-related gene variants were selected and used for gene-BPA interaction analysis. Multiple linear regression was used to analyze the interactions between BPA exposure and metabolism-related genes on serum lipid profiles. In the discovery stage, both ln-BPA and ln-BPA/Cr was associated with decreased high-density lipoprotein cholesterol (HDL-C). Gene-urinary BPA interaction for IGFBP7 rs9992658 was observed to associate with HDL-C levels in both discovery and validation stages, with Pinteraction equal to 9.87 × 10-4 (ln-BPA) and 1.22 × 10-3 (ln-BPA/Cr) in combined analyses. In addition, the inverse association of urinary BPA with HDL-C levels was only observed among individuals carrying rs9992658 AA genotype, but not in individuals carrying rs9992658 AC or CC genotypes. The interaction between BPA exposure and metabolism-related gene IGFBP7 (rs9992658) was associated with HDL-C levels.

HSPA12A was identified as a key driver in colorectal cancer GWAS loci 10q26.12 and modulated by an enhancer-promoter interaction.

Although genome-wide association studies (GWASs) have identified over 100 colorectal cancer (CRC) risk loci, an understanding of causal genes or risk variants and their biological functions in these loci remain unclear. Recently, genomic loci 10q26.12 with lead SNP rs1665650 was identified as an essential CRC risk loci of Asian populations. However, the functional mechanism of this region has not been fully clarified. Here, we applied an RNA interfering-based on-chip approach to screen for the genes essential for cell proliferation in the CRC risk loci 10q26.12. Notably, HSPA12A had the most significant effect among the identified genes and functioned as a crucial oncogene facilitating cell proliferation. Moreover, we conducted an integrative fine-mapping analysis to identify putative casual variants and further explored their association with CRC risk in a large-scale Chinese population consisting of 4054 cases and 4054 controls and also independently validated in 5208 cases and 20,832 controls from the UK biobank cohort. We identified a risk SNP rs7093835 in the intron of HSPA12A that was significantly associated with an increased risk of CRC (OR 1.23, 95% CI 1.08-1.41, P = 1.92 × 10-3). Mechanistically, the risk variant could facilitate an enhancer-promoter interaction mediated by the transcriptional factor (TF) GRHL1 and ultimately upregulate HSPA12A expression, which provides functional evidence to support our population findings. Collectively, our study reveals the important role of HSPA12A in CRC development and illustrates a novel enhancer-promoter interaction module between HSPA12A and its regulatory elements rs7093835, providing new insights into the etiology of CRC.

Triclosan is associated with breast cancer via oxidative stress and relative telomere length.

Introduction: Triclosan (TCS), a widely prescribed broad-spectrum antibacterial agent, is an endocrine-disrupting chemical. The relationship and biological mechanisms between TCS exposure and breast cancer (BC) are disputed. We aimed to examine the correlation between urinary TCS exposure and BC risk and estimated the mediating effects of oxidative stress and relative telomere length (RTL) in the above association. Methods: This case-control study included 302 BC patients and 302 healthy individuals in Wuhan, China. We detected urinary TCS, three common oxidative stress biomarkers [8-hydroxy-2-deoxyguanosine (8-OHdG), 8-iso-prostaglandin F2α (8-isoPGF2α), 4-hydroxy-2-nonenal-mercapturic acid (HNE-MA)], and RTL in peripheral blood mononuclear cells. Results: Significant associations were observed between log-transformed urinary concentrations of TCS, 8-OHdG, HNE-MA, 8-isoPGF2α, RTL, and BC risk, with the odds ratios (95% confidence intervals) being 1.58 (1.32-1.91), 3.08 (1.55-6.23), 3.39 (2.45-4.77), 3.99 (2.48-6.54), and 1.67 (1.35-2.09), respectively. Continuous TCS exposure was significantly positively correlated with RTL, HNE-MA, and 8-isoPGF2α (all p<0.05) but not with 8-OHdG (p = 0.060) after adjusting for covariates. The mediated proportions of 8-isoPGF22α and RTL in the relationship between TCS and BC risk were 12.84% and 8.95%, respectively (all p<0.001). Discussion: In conclusion, our study provides epidemiological evidence to confirmed the deleterious effects of TCS on BC and indicated the mediating effect of oxidative stress and RTL on the correlation between TCS and BC risk. Moreover, exploring the contribution of TCS to BC can clarify the biological mechanisms of TCS exposure, provide new clues for the pathogenesis of BC, which is of great significance to improving public health systems.

A variant in CYP26B1 associated with esophageal squamous cell carcinoma risk by affecting retinoic acid metabolism.

All-trans retinoic acid (ATRA) is the natural and synthetic analogue of vitamin A, playing an essential tumor suppressive role in multiple cancers including the esophageal squamous cell carcinoma (ESCC). Cytochrome P450 family 26 subfamily B member 1 (CYP26B1) exerts a critical regulator of ATRA levels through specific inactivation of ATRA to hydroxylated forms. Our previous exome-wide analyses revealed a rare missense variant in CYP26B1 significantly associated with ESCC risk in the Chinese population. However, it is still unclear whether there are common variants in CYP26B1 affect the susceptibility of ESCC and the tumor promotion role of CYP26B1 in vivo. In this research, we conducted a two-stage case-control study comprised of 5057 ESCC cases and 5397 controls, followed by a series of biochemical experiments to explore the function of CYP26B1 and its common variants in the tumorigenesis of ESCC. Intriguingly, we identified a missense variant rs2241057[A>G] in the fourth exon of CYP26B1 significantly associated with the ESCC risk (combined odds ratio = 1.28; 95% confidence interval = 1.15-1.42; p = 2.96 × 10-6 ). Through further functional analysis, we demonstrated that ESCC cells with the overexpression of rs2241057[G] had a significant lower level of retinoic acid, compared with the overexpression of rs2241057[A] or the control vector. In addition, the CYP26B1 overexpression and knock-out ESCC cells affected cell proliferation rate both in vitro and in vivo. These results highlighted the carcinogenicity of CYP26B1 related to the ATRA metabolism in ESCC risk.

Saturated fatty acid intake, genetic risk and colorectal cancer incidence: A large-scale prospective cohort study.

Previous investigations mainly focused on the associations of dietary fatty acids with colorectal cancer (CRC) risk, which ignored gene-environment interaction and mechanisms interpretation. We conducted a case-control study (751 cases and 3058 controls) and a prospective cohort study (125 021 participants) to explore the associations between dietary fatty acids, genetic risks, and CRC. Results showed that high intake of saturated fatty acid (SFA) was associated with a higher risk of CRC than low SFA intake (HR =1.22, 95% CI:1.02-1.46). Participants at high genetic risk had a greater risk of CRC with the HR of 2.48 (2.11-2.91) than those at low genetic risk. A multiplicative interaction of genetic risk and SFA intake with incident CRC risk was found (PInteraction = 7.59 × 10-20 ), demonstrating that participants with high genetic risk and high SFA intake had a 3.75-fold greater risk of CRC than those with low genetic risk and low SFA intake. Furthermore, incorporating PRS and SFA into traditional clinical risk factors improved the discriminatory accuracy for CRC risk stratification (AUC from 0.706 to 0.731). Multi-omics data showed that exposure to SFA-rich high-fat dietary (HFD) can responsively induce epigenome reprogramming of some oncogenes and pathological activation of fatty acid metabolism pathway, which may contribute to CRC development through changes in gut microbiomes, metabolites, and tumor-infiltrating immune cells. These findings suggest that individuals with high genetic risk of CRC may benefit from reducing SFA intake. The incorporation of SFA intake and PRS into traditional clinical risk factors will help improve high-risk sub-populations in individualized CRC prevention.

Identification of specific susceptibility loci for the early-onset colorectal cancer.

BACKGROUND: The incidence of early-onset colorectal cancer (EOCRC; patients < 50 years old) has been rising rapidly, whereas the EOCRC genetic susceptibility remains incompletely investigated. Here, we aimed to systematically identify specific susceptible genetic variants for EOCRC. METHODS: Two parallel GWASs were conducted in 17,789 CRC cases (including 1490 EOCRC cases) and 19,951 healthy controls. A polygenic risk score (PRS) model was built based on identified EOCRC-specific susceptibility variants by using the UK Biobank cohort. We also interpreted the potential biological mechanisms of the prioritized risk variant. RESULTS: We identified 49 independent susceptibility loci that were significantly associated with the susceptibility to EOCRC and the diagnosed age of CRC (both P < 5.0×10-4), replicating 3 previous CRC GWAS loci. There are 88 assigned susceptibility genes involved in chromatin assembly and DNA replication pathways, mainly associating with precancerous polyps. Additionally, we assessed the genetic effect of the identified variants by developing a PRS model. Compared to the individuals in the low genetic risk group, the individuals in the high genetic risk group have increased EOCRC risk, and these results were replicated in the UKB cohort with a 1.63-fold risk (95% CI: 1.32-2.02, P = 7.67×10-6). The addition of the identified EOCRC risk loci significantly increased the prediction accuracy of the PRS model, compared to the PRS model derived from the previous GWAS-identified loci. Mechanistically, we also elucidated that rs12794623 may contribute to the early stage of CRC carcinogenesis via allele-specific regulating the expression of POLA2. CONCLUSIONS: These findings will broaden the understanding of the etiology of EOCRC and may facilitate the early screening and individualized prevention.

Genome-Wide Splicing Quantitative Expression Locus Analysis Identifies Causal Risk Variants for Non-Small Cell Lung Cancer.

Alternative RNA splicing is an essential mechanism linking genetic variation to human diseases. While the signals from genome-wide association studies (GWAS) have been linked to expression quantitative trait loci (eQTL) in previous studies, further work is needed to better elucidate the relationship to other genetic regulatory mechanisms, such as splicing QTLs (sQTL). Here, we performed a genome-wide sQTL analysis to identify variants that might affect RNA splicing in 1,010 non-small cell lung cancer (NSCLC) samples from The Cancer Genome Atlas. The identified sQTLs were largely independent of eQTLs and were predominantly enriched in exonic regions, genetic regulatory elements, RNA-binding protein (RBP) binding sites, and known NSCLC risk loci. In addition, target genes affected by sQTLs (sGenes) were involved in multiple processes in cancer, including cell growth, apoptosis, metabolism, immune infiltration, and drug responses, and sGenes were frequently altered genetically in NSCLC. Systematic screening of sQTLs associated with NSCLC risk using GWAS data from 15,474 cases and 12,375 controls identified an sQTL variant rs156697-G allele that was significantly associated with an increased risk of NSCLC. The association between the rs156697-G variant and NSCLC risk was further validated in two additional large population cohorts. The risk variant promoted inclusion of GSTO2 alternative exon 5 and led to higher expression of the GSTO2 full-length isoform (GSTO2-V1) and lower expression of the truncated GSTO2 isoform (GSTO2-V2), which was induced by RBP quaking (QKI). Mechanistically, compared with GSTO2-V1, GSTO2-V2 inhibited NSCLC cells proliferation by increasing S-glutathionylation of AKT1 and thereby functionally blocking AKT1 phosphorylation and activation. Overall, this study provides a comprehensive view of splicing variants linked to NSCLC risk and provides a set of genetic targets with therapeutic potential. SIGNIFICANCE: Analysis of sQTL reveals the role of genetically driven mRNA splicing alterations in NSCLC risk and elucidates that rs156697 variant impacts risk by altering GSTO2 splicing.

Birthweight, genetic risk, and gastrointestinal cancer incidence: a prospective cohort study.

BACKGROUND: The epidemiologic studies investigating the association of birthweight and genetic factors with gastrointestinal cancer remain scarce. The study aimed to prospectively assess the interactions and joint effects of birthweight and genetic risk levels on gastrointestinal cancer incidence in adulthood. METHODS: A total of 254,997 participants were included in the UK Biobank study. We used multivariate restricted cubic splines and Cox regression models to estimate the hazard ratios (HRs) and 95% confidential intervals (CI) for the association between birthweight and gastrointestinal cancer risk, then constructed a polygenic risk score (PRS) to assess its interaction and joint effect with birthweight on the development of gastrointestinal cancer. RESULTS: We documented 2512 incident cases during a median follow-up of 8.88 years. Compare with participants reporting a normal birthweight (2.5-4.5 kg), multivariable-adjusted HR of gastrointestinal cancer incidence for participants with high birthweight (≥4.5 kg) was 1.17 (95%CI: 1.01-1.36). Such association was remarkably observed in pancreatic cancer, with an HR of 1.82 (95%CI: 1.26-2.64). No statistically significant association was observed between low birth weight and gastrointestinal cancers. Participants with high birthweight and high PRS had the highest risk of gastrointestinal cancer (HR: 2.95, 95%CI: 2.19-3.96). CONCLUSION: Our findings highlight that high birthweight is associated with a higher incidence of gastrointestinal cancer, especially for pancreatic cancer. Benefits would be obtained from birthweight control, particularly for individuals with a high genetic risk.KEY MESSAGESThe epidemiologic studies investigating the association of birthweight and genetic factors with gastrointestinal cancer remain scarce.This cohort study of 254,997 adults in the United Kingdom found an association of high birthweight with the incidence of gastrointestinal cancer, especially for pancreatic cancer, and also found that participants with high birthweight and high polygenic risk score had the highest risk of gastrointestinal cancer.Our data suggests a possible effect of in utero or early life exposures on adulthood gastrointestinal cancer, especially for those with a high genetic risk.