An atlas of genetic effects on cellular composition of the tumor microenvironment.

Deciphering the composition of the tumor microenvironment (TME) is critical for understanding tumorigenesis and to design immunotherapies. In the present study, we mapped genetic effects on cell-type proportions using single-cell and bulk RNA sequencing data, identifying 3,494 immunity quantitative trait loci (immunQTLs) across 23 cancer types from The Cancer Genome Atlas. Functional annotation revealed regulatory potential and we further assigned 1,668 genes that regulate TME composition. We constructed a combined immunQTL map by integrating data from European and Chinese colorectal cancer (CRC) samples. A polygenic risk score that incorporates these immunQTLs and hits on a genome-wide association study outperformed in CRC risk stratification within 447,495 multiethnic individuals. Using large-scale population cohorts, we identified that the immunQTL rs1360948 is associated with CRC risk and prognosis. Mechanistically, the rs1360948-G-allele increases CCL2 expression, recruiting regulatory T cells that can exert immunosuppressive effects on CRC progression. Blocking the CCL2-CCR2 axis enhanced anti-programmed cell death protein 1 ligand therapy. Finally, we have established a database (CancerlmmunityQTL2) to serve the research community and advance our understanding of immunogenomic interactions in cancer pathogenesis.

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.

Genetic architectures of the human hippocampus and those involved in neuropsychiatric traits.

BACKGROUND: The hippocampus, with its complex subfields, is linked to numerous neuropsychiatric traits. While most research has focused on its global structure or a few specific subfields, a comprehensive analysis of hippocampal substructures and their genetic correlations across a wide range of neuropsychiatric traits remains underexplored. Given the hippocampus's high heritability, considering hippocampal and subfield volumes (HASV) as endophenotypes for neuropsychiatric conditions is essential. METHODS: We analyzed MRI-derived volumetric data of hippocampal and subfield structures from 41,525 UK Biobank participants. Genome-wide association studies (GWAS) on 24 HASV traits were conducted, followed by genetic correlation, overlap, and Mendelian randomization (MR) analyses with 10 common neuropsychiatric traits. Polygenic risk scores (PRS) based on HASV traits were also evaluated for predicting these traits. RESULTS: Our analysis identified 352 independent genetic variants surpassing a significance threshold of 2.1 × 10-9 within the 24 HASV traits, located across 93 chromosomal regions. Notably, the regions 12q14.3, 17q21.31, 12q24.22, 6q21, 9q33.1, 6q25.1, and 2q24.2 were found to influence multiple HASVs. Gene set analysis revealed enrichment of neural differentiation and signaling pathways, as well as protein binding and degradation. Of 240 HASV-neuropsychiatric trait pairs, 75 demonstrated significant genetic correlations (P < 0.05/240), revealing 433 pleiotropic loci. Particularly, genes like ACBD4, ARHGAP27, KANSL1, MAPT, ARL17A, and ARL17B were involved in over 50 HASV-neuropsychiatric pairs. Leveraging Mendelian randomization analysis, we further confirmed that atrophy in the left hippocampus, right hippocampus, right hippocampal body, and right CA1-3 region were associated with an increased risk of developing Parkinson's disease (PD). Furthermore, PRS for all four HASVs were significantly linked to a higher risk of Parkinson's disease (PD), with the highest hazard ratio (HR) of 1.30 (95% CI 1.18-1.43, P = 6.15 × 10⁻8) for right hippocampal volume. CONCLUSIONS: These findings highlight the extensive distribution of pleiotropic genetic determinants between HASVs and neuropsychiatric traits. Moreover, they suggest a significant potential for effectively managing and intervening in these diseases during their early stages.

Integrated systematic functional screen and fine-mapping decipher the role and genetic regulation of RPS19 in colorectal cancer development.

Despite genome-wide association studies (GWAS) have identified more than 200 risk loci associated with colorectal cancer (CRC), the causal genes or risk variants within these loci and their biological functions remain not fully revealed. Recently, the genomic locus 19q13.2, with the lead SNP rs1800469 was identified as a crucial CRC risk locus in Asian populations. However, the functional mechanism of this region has not been fully elucidated. Here we employed an RNA interfering-based on-chip approach to screen for the genes essential for cell proliferation in the CRC risk locus 19q13.2. Notably, we found that RPS19 exhibited the most significant effect among the identified genes and acted as a critical oncogene facilitating CRC cell proliferation. Subsequently, combining integrative fine-mapping analysis and a large-scale population study consisting of 6027 cases and 6099 controls, we prioritized rs1025497 as a potential causal candidate for CRC risk, demonstrating that rs1025497[A] allele significantly reduced the risk of CRC (OR 0.70, 95% confidence interval = 0.56-0.83, P = 1.12 × 10-6), which was further validated in UK Biobank cohort comprising 5,313 cases and 21,252 controls. Mechanistically, we experimentally elucidated that variant rs1025497 might acted as an allele-specific silencer, inhibiting the expression level of oncogene RPS19 mediated by the transcription suppressive factor HBP1. Taken together, our sturdy unveils the significant role of RPS19 during CRC pathogenesis and delineates its distal regulatory mechanism mediated by rs1025497, advancing our understanding of the etiology of CRC and provided new insights into the personalized medicine of human cancer.

In Utero and Childhood/Adolescence Exposure to Tobacco Smoke, Genetic Risk, and Lung Cancer Incidence and Mortality in Adulthood.

Rationale: The individual effects of early-life tobacco smoke exposure and its interactions with genetic factors on lung cancer in adulthood remain unclear. Objectives: To investigate the associations of early-life tobacco exposures as well as their interactions with polygenic risk scores (PRSs) with lung cancer incidence and mortality. Methods: A total of 432,831 participants from the UK Biobank study were included. We estimated the associations of in utero exposure to tobacco smoke, the age of smoking initiation and their interactions with PRSs with lung cancer incidence and mortality in adulthood using Cox proportional hazard models. Measurements and Main Results: Lung cancer incidence (hazard ratio [HR]: 1.59, 95% confidence interval [CI], 1.44-1.76) increased among participants with in utero tobacco exposure. Multivariable-adjusted HRs (with 95% CIs) of lung cancer incidence for smoking initiation in adulthood, adolescence, and childhood (versus never-smokers) were 6.10 (5.25-7.09), 9.56 (8.31-11.00), and 15.15 (12.90-17.79) (Ptrend < 0.001). Similar findings were observed in lung cancer mortality. Participants with high PRSs and in utero tobacco exposure (versus low PRSs participants without in utero exposure) had an HR of 2.35 for lung cancer incidence (95% CI, 1.97-2.80, Pinteraction = 0.089) and 2.43 for mortality (95% CI, 2.05-2.88, Pinteraction = 0.032). High PRSs with smoking initiation in childhood (versus never-smokers with low PRSs) had HRs of 18.71 for incidence (95% CI, 14.21-24.63, Pinteraction = 0.004) and 19.74 for mortality (95% CI, 14.98-26.01, Pinteraction = 0.033). Conclusions: In utero and childhood/adolescence exposure to tobacco smoke and its interaction with genetic factors may substantially increase the risks of lung cancer incidence and mortality in adulthood.

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.

Recombinant expression of IGF-1 and LR3 IGF-1 fused with xylanase in Pichia pastoris.

Insulin-like growth factor-1 (IGF-1) is a pleiotropic protein hormone and has become an attractive therapeutic target because of its multiple roles in various physiological processes, including growth, development, and metabolism. However, its production is hindered by low heterogenous protein expression levels in various expression systems and hard to meet the needs of clinical and scientific research. Here, we report that human IGF-1 and its analog Long R3 IGF-1 (LR3 IGF-1) are recombinant expressed and produced in the Pichia pastoris (P. pastoris) expression system through being fused with highly expressed xylanase XynCDBFV. Furthermore, purified IGF-1 and LR3 IGF-1 display excellent bioactivity of cell proliferation compared to the standard IGF-1. Moreover, higher heterologous expression levels of the fusion proteins XynCDBFV-IGF-1 and XynCDBFV-LR3 IGF-1 are achieved by fermentation in a 15-L bioreactor, reaching up to about 0.5 g/L XynCDBFV-IGF-1 and 1 g/L XynCDBFV-TEV-LR3 IGF-1. Taken together, high recombinant expression of bioactive IGF-1 and LR3 IGF-1 is acquired with the assistance of xylanase as a fusion partner in P. pastoris, which could be used for both clinical and scientific applications. KEY POINTS: • Human IGF-1 and LR3 IGF-1 are produced in the P. pastoris expression system. • Purified IGF-1 and LR3 IGF-1 show bioactivity comparable to the standard IGF-1. • High heterologous expression of IGF-1 and LR3 IGF-1 is achieved by fermentation in a bioreactor.

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.

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.

CancerImmunityQTL: a database to systematically evaluate the impact of genetic variants on immune infiltration in human cancer.

Tumor-infiltrating immune cells as integral component of the tumor microenvironment are associated with tumor progress, prognosis and responses to immunotherapy. Genetic variants have been demonstrated to impact tumor-infiltrating, underscoring the heritable character of immune landscape. Therefore, identification of immunity quantitative trait loci (immunQTLs), which evaluate the effect of genetic variants on immune cells infiltration, might present a critical step toward fully understanding the contribution of genetic variants in tumor development. Although emerging studies have demonstrated the determinants of germline variants on immune infiltration, no database has yet been developed to systematically analyze immunQTLs across multiple cancer types. Using genotype data from TCGA database and immune cell fractions estimated by CIBERSORT, we developed a computational pipeline to identify immunQTLs in 33 cancer types. A total of 913 immunQTLs across different cancer types were identified. Among them, 5 immunQTLs are associated with patient overall survival. Furthermore, by integrating immunQTLs with GWAS data, we identified 527 immunQTLs overlapping with known GWAS linkage disequilibrium regions. Finally, we constructed a user-friendly database, CancerImmunityQTL (http://www.cancerimmunityqtl-hust.com/) for users to browse, search and download data of interest. This database provides an informative resource to understand the germline determinants of immune infiltration in human cancer and benefit from personalized cancer immunotherapy.

Hepatocellular carcinoma risk variant modulates lncRNA HLA-DQB1-AS1 expression via a long-range enhancer-promoter interaction.

Substantial evidence highlighted the critical role of long non-coding RNAs (lncRNA) in driving hepatocarcinogenesis. We hypothesized that functional variants in genome-wide association studies (GWASs) associated loci might alter the expression levels of lncRNAs and contribute to the development of hepatocellular carcinoma (HCC). Here, we prioritized potentially cis-expression quantitative trait loci-based single nucleotide polymorphism (SNP)-lncRNA association together with the physical interaction by the analyses from Hi-C data in GWAS loci of chronic hepatitis B and HCC. Subsequently, by leveraging two-stage case-control study (1738 hepatitis B [HBV]) related HCC cases and 1988 HBV persistent carriers) and biological assays, we identified that rs2647046 was significantly associated with HCC risk (odds ratio = 1.26, 95% CI = 1.11 to 1.43, P = 4.14 × 10-4). Luciferase reporter assays and electrophoretic mobility shift assays showed that rs2647046 A allele significantly increased transcriptional activity via influencing transcript factor binding affinity. Allele-specific chromosome conformation capture assays revealed that enhancer with rs2647046 interacted with the HLA-DQB1-AS1 promoter to allele-specifically influence its expression by CTCF-mediated long-range loop. Cell proliferation assays indicated that HLA-DQB1-AS1 is a potential oncogene in HCC. Our study showed HLA-DQB1-AS1 regulated by a causal SNP in a long-range interaction manner conferred the susceptibility to HCC, suggesting an important mechanism of modulating lncRNA expression for risk-associated SNPs in the etiology of HCC.

Clinical characteristics and risk factors of fatal patients with COVID-19: a retrospective cohort study in Wuhan, China.

BACKGROUND: The coronavirus disease 2019 (COVID-19) has caused a global pandemic, resulting in considerable mortality. The risk factors, clinical treatments, especially comprehensive risk models for COVID-19 death are urgently warranted. METHODS: In this retrospective study, 281 non-survivors and 712 survivors with propensity score matching by age, sex, and comorbidities were enrolled from January 13, 2020 to March 31, 2020. RESULTS: Higher SOFA, qSOFA, APACHE II and SIRS scores, hypoxia, elevated inflammatory cytokines, multi-organ dysfunction, decreased immune cell subsets, and complications were significantly associated with the higher COVID-19 death risk. In addition to traditional predictors for death risk, including APACHE II (AUC = 0.83), SIRS (AUC = 0.75), SOFA (AUC = 0.70) and qSOFA scores (AUC = 0.61), another four prediction models that included immune cells subsets (AUC = 0.90), multiple organ damage biomarkers (AUC = 0.89), complications (AUC = 0.88) and inflammatory-related indexes (AUC = 0.75) were established. Additionally, the predictive accuracy of combining these risk factors (AUC = 0.950) was also significantly higher than that of each risk group alone, which was significant for early clinical management for COVID-19. CONCLUSIONS: The potential risk factors could help to predict the clinical prognosis of COVID-19 patients at an early stage. The combined model might be more suitable for the death risk evaluation of COVID-19.

N6-methyladenosine mRNA methylation of PIK3CB regulates AKT signalling to promote PTEN-deficient pancreatic cancer progression.

OBJECTIVE: Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal cancers worldwide. Thus far, most drugs have failed to significantly improve patient survival. N6-methyladenosine (m6A) plays an important role in the progression of PDAC, but its aberrant regulation driven by germline variants in human diseases remains unclear. DESIGN: We first performed an exome-wide association analysis in 518 PDAC patients with overall survival and replicated in an independent population containing 552 PDAC patients. Then, a series of biochemical experiments in vitro and in vivo were conducted to investigate potential mechanisms of the candidate variant and its target gene PIK3CB underlying the PDAC progression. Moreover, the PIK3CB-selective inhibitor KIN-193 was used to block PDAC tumour growth. RESULTS: We identified a missense variant rs142933486 in PIK3CB that is significantly associated with the overall survival of PDAC by reducing the PIK3CB m6A level, which facilitated its mRNA and protein expression levels mediated by the m6A 'writer' complex (METTL13/METTL14/WTAP) and the m6A 'reader' YTHDF2. The upregulation of PIK3CB is widely found in PDAC tumour tissues and significantly correlated with the poor prognosis of PDAC, especially in PTEN-deficient patients. We further demonstrated that PIK3CB overexpression substantially enhanced the proliferation and migration abilities of PTEN-deficient PDAC cells and activated AKT signalling pathway. Remarkably, KIN-193, a PIK3CB-selective inhibitor, is shown to serve as an effective anticancer agent for blocking PTEN-deficient PDAC. CONCLUSIONS: These findings demonstrate aberrant m6A homoeostasis as an oncogenic mechanism in PDAC and highlight the potential of PIK3CB as a therapeutic target for this disease.

Clinical characteristics and risk factors associated with COVID-19 disease severity in patients with cancer in Wuhan, China: a multicentre, retrospective, cohort study.

BACKGROUND: COVID-19 has spread globally. Epidemiological susceptibility to COVID-19 has been reported in patients with cancer. We aimed to systematically characterise clinical features and determine risk factors of COVID-19 disease severity for patients with cancer and COVID-19. METHODS: In this multicentre, retrospective, cohort study, we included all adult patients (aged ≥18 years) with any type of malignant solid tumours and haematological malignancy who were admitted to nine hospitals in Wuhan, China, with laboratory-confirmed COVID-19 between Jan 13 and March 18, 2020. Enrolled patients were statistically matched (2:1) with patients admitted with COVID-19 who did not have cancer with propensity score on the basis of age, sex, and comorbidities. Demographic characteristics, laboratory examinations, illness severity, and clinical interventions were compared between patients with COVID-19 with or without cancer as well as between patients with cancer with non-severe or severe COVID-19. COVID-19 disease severity was defined on admission on the basis of the WHO guidelines. Univariable and multivariable logistic regression, adjusted for age, sex, comorbidities, cancer type, tumour stage, and antitumour treatments, were used to explore risk factors associated with COVID-19 disease severity. This study was registered in the Chinese Clinical Trial Register, ChiCTR2000030807. FINDINGS: Between Jan 13 and March 18, 2020, 13 077 patients with COVID-19 were admitted to the nine hospitals in Wuhan and 232 patients with cancer and 519 statistically matched patients without cancer were enrolled. Median follow-up was 29 days (IQR 22-38) in patients with cancer and 27 days (20-35) in patients without cancer. Patients with cancer were more likely to have severe COVID-19 than patients without cancer (148 [64%] of 232 vs 166 [32%] of 519; odds ratio [OR] 3·61 [95% CI 2·59-5·04]; p<0·0001). Risk factors previously reported in patients without cancer, such as older age; elevated interleukin 6, procalcitonin, and D-dimer; and reduced lymphocytes were validated in patients with cancer. We also identified advanced tumour stage (OR 2·60, 95% CI 1·05-6·43; p=0·039), elevated tumour necrosis factor α (1·22, 1·01-1·47; p=0·037), elevated N-terminal pro-B-type natriuretic peptide (1·65, 1·03-2·78; p=0·032), reduced CD4+ T cells (0·84, 0·71-0·98; p=0·031), and reduced albumin-globulin ratio (0·12, 0·02-0·77; p=0·024) as risk factors of COVID-19 severity in patients with cancer. INTERPRETATION: Patients with cancer and COVID-19 were more likely to deteriorate into severe illness than those without cancer. The risk factors identified here could be helpful for early clinical surveillance of disease progression in patients with cancer who present with COVID-19. FUNDING: China National Natural Science Foundation.

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.

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.

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.

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.

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.

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.