Large-scale genetic study in East Asians identifies six new loci associated with colorectal cancer risk

A pleiotropy scan to discover new susceptibility loci for pancreatic ductal adenocarcinoma

Pleiotropic variants (i.e. genetic polymorphisms influencing more than one phenotype) are often associated with cancer risk. A scan of pleiotropic variants was successfully conducted 10 years ago in relation to pancreatic ductal adenocarcinoma susceptibility. However, in the last decade, genetic association studies performed on several human traits have greatly increased the number of known pleiotropic variants. Based on the hypothesis that variants already associated with a least one trait have a higher probability of association with other traits, 61 052 variants reported to be associated by at least one genome-wide association study with at least one human trait were tested in the present study consisting of two phases (discovery and validation), comprising a total of 16 055 pancreatic ductal adenocarcinoma (PDAC) cases and 212 149 controls. The meta-analysis of the two phases showed two loci (10q21.1-rs4948550 (P = 6.52 × 10-5) and 7q36.3-rs288762 (P = 3.03 × 10-5) potentially associated with PDAC risk. 10q21.1-rs4948550 shows a high degree of pleiotropy and it is also associated with colorectal cancer risk while 7q36.3-rs288762 is situated 28,558 base pairs upstream of the Sonic Hedgehog (SHH) gene, which is involved in the cell-differentiation process and PDAC etiopathogenesis. In conclusion, none of the single nucleotide polymorphisms (SNPs) showed a formally statistically significant association after correction for multiple testing. However, given their pleiotropic nature and association with various human traits including colorectal cancer, the two SNPs showing the best associations with PDAC risk merit further investigation through fine mapping and ad hoc functional studies.

Population genomics advances in frontier ethnic minorities in China

China, with its large geographic span, possesses rich genetic diversity across vast frontier regions in addition to the Han Chinese majority. Importantly, demographic events and various natural and cultural environments in Chinese frontier regions have shaped the genomic diversity of ethnic minorities via local adaptations. Thus, insights into the genetic diversity and adaptive evolution of these under-represented ethnic groups are crucial for understanding evolutionary scenarios and biomedical implications in East Asian populations. Here, we focus on ethnic minorities in Chinese frontier regions and review research advances regarding genomic diversity, genetic structure, population history, genetic admixture, and local adaptation. We first provide an overview of the extensive genetic diversity across populations in different Chinese frontier regions. Next, we summarize research progress regarding genetic ancestry, demographic history, the adaptive process, and the archaic identification of multiple ethnic minorities in different Chinese frontier regions. Finally, we discuss the gaps and opportunities in genomic studies of Chinese populations and the need for a more comprehensive understanding of genomic diversity and the evolution of populations of East Asian ancestry in the post-genomic era.

Systematic functional interrogation of genome-wide association studies locus 17p13.3 deciphered role and genetic control of FAM57A in colorectal cancer development

Objective

Genome-wide association studies (GWAS) have identified over 150 risk loci linked to colorectal cancer (CRC), including the 17p13.3 locus with the tag single nucleotide polymorphism (SNP) rs12603526 in the Asian population. However, the specific causal gene and the functional regulatory mechanisms in this region remain unresolved, necessitating further investigation to elucidate the underlying mechanisms of CRC.

Methods

We employed an RNA interference-based functional approach to identify genes critical for CRC cell proliferation at the GWAS locus 17p13.3. Bioinformatic fine-mapping analysis was conducted to prioritize causal variants. A large-scale study involving 7,013 cases and 7,329 controls from a Chinese population, along with another cohort of 5,158 cases and 20,632 controls from the UK Biobank, was performed to validate the association between the candidate variant and the gene. A series of biological experiments was conducted to explore the function of the candidate gene and its regulatory mechanisms.

Results

We identified FAM57A as a key oncogene that promotes CRC cell proliferation, and confirmed its carcinogenic role through in vitro proliferation assays. The variant rs526835 was prioritized as a causal candidate for CRC risk, located in a functional region with enhancer properties, and showed a significant quantitative association with FAM57A expression. The rs526835 [T] variant was associated with a 1.17-fold increase in CRC risk [95% confidence interval (95% CI): 1.11-1.23, P=1.23×10-9] in the large-scale Chinese cohort, which was further corroborated in the UK Biobank cohort. Mechanistically, we demonstrated that rs526835 enhances a promoter-enhancer interaction mediated by the transcription factor JUN, leading to increased expression of FAM57A.

Conclusions

We reveal the underlying mechanisms of CRC predisposition at the GWAS locus 17p13.3. Additionally, our findings highlight the critical role of FAM57A in CRC pathogenesis and introduce a novel enhancer-promoter interaction between FAM57A and rs526835, which could inform future precision prevention and personalized cancer therapies.

A literature review of genetics and epigenetics of HCV-related hepatocellular carcinoma: translational impact

Background and Objective

Hepatocellular carcinoma (HCC) poses a significant global health burden and ranks as the fifth most prevalent cancer on a global scale. Hepatitis C virus (HCV) infection remains one of the major risk factors for HCC development. HCC is a heterogeneous disease, and the development of HCC caused by HCV is intricate and involves various factors, including genetic susceptibility, viral factors, immune response due to chronic inflammation, alcohol abuse, and metabolic dysfunction associated with fatty liver disease. In this review, we provide a comprehensive and updated review of research on the genetics and epigenetic mechanisms implicated in developing HCC associated with HCV infection. We also discuss the potential translational implications, including novel biomarkers and drugs for treatment.

Methods

A comprehensive literature search was conducted in June 2023 in PubMed and Embase databases.

Key Content and Findings

Recent findings indicate that a variety of genetic and epigenetic processes are involved in the pathogenesis and continue to exist even after the complete elimination of HCV. The deregulation of the epigenome has been identified as a significant factor in the deletrious effects of liver disease, especially during the initial stages when genetic alterations are uncommon. The enduring "epigenetic memory" of gene expression is believed to be regulated by epigenetic mechanisms, indicating that alterations caused by HCV infection continue to exist and are linked to the risk of development of liver cancer even after successful treatment. Systems biology analytical methods will be required to delineate the magnitude and significance of both genetic and epigenomic alterations in tumor evolution.

Conclusions

By facilitating a more profound understanding of these aspects, this will ultimately foster the advancement of novel therapies and ultimately improve outcomes for patients.

ZMIZ1 Regulates Proliferation, Autophagy and Apoptosis of Colon Cancer Cells by Mediating Ubiquitin-Proteasome Degradation of SIRT1

There are nearly 1.15 million new cases of colon cancer, as well as 586,858 deaths from colon cancer worldwide in 2020. The aim of this study is to reveal whether ZMIZ1 can control the fate of colon cancer cells and the mechanism by which it functions. Specific shRNA transfection was used to knock down the expression of ZMIZ1 in colon cancer cell lines (HCT116 and HT29), and cell proliferation was detected using EdU and CCK-8 reagents, apoptosis by flow cytometry, and autophagy by western blot. The interaction of ZMIZ1 and SIRT1 was analyzed. Knockdown of ZMIZ1 significantly inhibited autophagy and proliferation, and induced apoptosis of HCT116 and HT29 cells. The mRNA level of SIRT1 was not affected by ZMIZ1 knockdown, but the protein level of SIRT1 was significantly decreased and the protein level of the SIRT1-specific substrate, acetylated FOXO3a, was reduced. Immunoprecipitation assays identified the interaction between SIRT1 and ZMIZ1 in HCT116 and HT29 cells. ZMIZ1 increased intracellular ubiquitination of SIRT1. Knockdown or pharmacological inhibition of SIRT1 neutralized the effects of ZMIZ knockdown on proliferation, autophagy and apoptosis in HCT116 and HT29 cells. ZMIZ1 may control the fate of colon cancer cells through the SIRT1/FOXO3a axis. Targeting ZMIZ1 would be beneficial for the treatment of colon cancer.

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.

Genetic risk impacts the association of menopausal hormone therapy with colorectal cancer risk

BACKGROUND

Menopausal hormone therapy (MHT), a common treatment to relieve symptoms of menopause, is associated with a lower risk of colorectal cancer (CRC). To inform CRC risk prediction and MHT risk-benefit assessment, we aimed to evaluate the joint association of a polygenic risk score (PRS) for CRC and MHT on CRC risk.

METHODS

We used data from 28,486 postmenopausal women (11,519 cases and 16,967 controls) of European descent. A PRS based on 141 CRC-associated genetic variants was modeled as a categorical variable in quartiles. Multiplicative interaction between PRS and MHT use was evaluated using logistic regression. Additive interaction was measured using the relative excess risk due to interaction (RERI). 30-year cumulative risks of CRC for 50-year-old women according to MHT use and PRS were calculated.

RESULTS

The reduction in odds ratios by MHT use was larger in women within the highest quartile of PRS compared to that in women within the lowest quartile of PRS (p-value = 2.7 × 10-8). At the highest quartile of PRS, the 30-year CRC risk was statistically significantly lower for women taking any MHT than for women not taking any MHT, 3.7% (3.3%-4.0%) vs 6.1% (5.7%-6.5%) (difference 2.4%, P-value = 1.83 × 10-14); these differences were also statistically significant but smaller in magnitude in the lowest PRS quartile, 1.6% (1.4%-1.8%) vs 2.2% (1.9%-2.4%) (difference 0.6%, P-value = 1.01 × 10-3), indicating 4 times greater reduction in absolute risk associated with any MHT use in the highest compared to the lowest quartile of genetic CRC risk.

CONCLUSIONS

MHT use has a greater impact on the reduction of CRC risk for women at higher genetic risk. These findings have implications for the development of risk prediction models for CRC and potentially for the consideration of genetic information in the risk-benefit assessment of MHT use.

Mitochondrial genome transfer drives metabolic reprogramming in adjacent colonic epithelial cells promoting TGFβ1-mediated tumor progression

Although nontumor components play an essential role in colon cancer (CC) progression, the intercellular communication between CC cells and adjacent colonic epithelial cells (CECs) remains poorly understood. Here, we show that intact mitochondrial genome (mitochondrial DNA, mtDNA) is enriched in serum extracellular vesicles (EVs) from CC patients and positively correlated with tumor stage. Intriguingly, circular mtDNA transferred via tumor cell-derived EVs (EV-mtDNA) enhances mitochondrial respiration and reactive oxygen species (ROS) production in CECs. Moreover, the EV-mtDNA increases TGFβ1 expression in CECs, which in turn promotes tumor progression. Mechanistically, the intercellular mtDNA transfer activates the mitochondrial respiratory chain to induce the ROS-driven RelA nuclear translocation in CECs, thereby transcriptionally regulating TGFβ1 expression and promoting tumor progression via the TGFβ/Smad pathway. Hence, this study highlights EV-mtDNA as a major driver of paracrine metabolic crosstalk between CC cells and adjacent CECs, possibly identifying it as a potential biomarker and therapeutic target for CC.

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.

Germline genetic regulation of the colorectal tumor immune microenvironment

OBJECTIVE

To evaluate the contribution of germline genetics to regulating the briskness and diversity of T cell responses in CRC, we conducted a genome-wide association study to examine the associations between germline genetic variation and quantitative measures of T cell landscapes in 2,876 colorectal tumors from participants in the Molecular Epidemiology of Colorectal Cancer Study (MECC).

METHODS

Germline DNA samples were genotyped and imputed using genome-wide arrays. Tumor DNA samples were extracted from paraffin blocks, and T cell receptor clonality and abundance were quantified by immunoSEQ (Adaptive Biotechnologies, Seattle, WA). Tumor infiltrating lymphocytes per high powered field (TILs/hpf) were scored by a gastrointestinal pathologist. Regression models were used to evaluate the associations between each variant and the three T-cell features, adjusting for sex, age, genotyping platform, and global ancestry. Three independent datasets were used for replication.

RESULTS

We identified a SNP (rs4918567) near RBM20 associated with clonality at a genome-wide significant threshold of 5 × 10- 8, with a consistent direction of association in both discovery and replication datasets. Expression quantitative trait (eQTL) analyses and in silico functional annotation for these loci provided insights into potential functional roles, including a statistically significant eQTL between the T allele at rs4918567 and higher expression of ADRA2A (P = 0.012) in healthy colon mucosa.

CONCLUSIONS

Our study suggests that germline genetic variation is associated with the quantity and diversity of adaptive immune responses in CRC. Further studies are warranted to replicate these findings in additional samples and to investigate functional genomic mechanisms.

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.

Circulating white blood cell traits and colorectal cancer risk: A Mendelian randomisation study

Observational studies have suggested a protective role for eosinophils in colorectal cancer (CRC) development and implicated neutrophils, but the causal relationships remain unclear. Here, we aimed to estimate the causal effect of circulating white blood cell (WBC) counts (N = ~550 000) for basophils, eosinophils, monocytes, lymphocytes and neutrophils on CRC risk (N = 52 775 cases and 45 940 controls) using Mendelian randomisation (MR). For comparison, we also examined this relationship using individual-level data from UK Biobank (4043 incident CRC cases and 332 773 controls) in a longitudinal cohort analysis. The inverse-variance weighted (IVW) MR analysis suggested a protective effect of increased basophil count and eosinophil count on CRC risk [OR per 1-SD increase: 0.88, 95% CI: 0.78-0.99, P = .04; OR: 0.93, 95% CI: 0.88-0.98, P = .01]. The protective effect of eosinophils remained [OR per 1-SD increase: 0.88, 95% CI: 0.80-0.97, P = .01] following adjustments for all other WBC subtypes, to account for genetic correlation between the traits, using multivariable MR. A protective effect of increased lymphocyte count on CRC risk was also found [OR: 0.84, 95% CI: 0.76-0.93, P = 6.70e-4] following adjustment. Consistent with MR results, a protective effect for eosinophils in the cohort analysis in the fully adjusted model [RR per 1-SD increase: 0.96, 95% CI: 0.93-0.99, P = .02] and following adjustment for the other WBC subtypes [RR: 0.96, 95% CI: 0.93-0.99, P = .001] was observed. Our study implicates peripheral blood immune cells, in particular eosinophils and lymphocytes, in CRC development, highlighting a need for mechanistic studies to interrogate these relationships.

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.

FADS2 function at the major cancer hotspot 11q13 locus alters fatty acid metabolism in cancer

Dysregulation of fatty acid metabolism and de novo lipogenesis is a key driver of several cancer types through highly unsaturated fatty acid (HUFA) signaling precursors such as arachidonic acid. The human chromosome 11q13 locus has long been established as the most frequently amplified in a variety of human cancers. The fatty acid desaturase genes (FADS1, FADS2 and FADS3) responsible for HUFA biosynthesis localize to the 11q12-13.1 region. FADS2 activity is promiscuous, catalyzing biosynthesis of several unsaturated fatty acids by Δ6, Δ8, and Δ4 desaturation. Our main aim here is to review known and putative consequences of FADS2 dysregulation due to effects on the 11q13 locus potentially driving various cancer types. FADS2 silencing causes synthesis of sciadonic acid (5Z,11Z,14Z-20:3) in MCF7 cells and breast cancer in vivo. 5Z,11Z,14Z-20:3 is structurally identical to arachidonic acid (5Z,8Z,11Z,14Z-20:4) except it lacks the internal Δ8 double bond required for prostaglandin and leukotriene synthesis, among other eicosanoids. Palmitic acid has substrate specificity for both SCD and FADS2. Melanoma, prostate, liver and lung cancer cells insensitive to SCD inhibition show increased FADS2 activity and sapienic acid biosynthesis. Elevated serum mead acid levels found in hepatocellular carcinoma patients suggest an unsatisfied demand for arachidonic acid. FADS2 circular RNAs are at high levels in colorectal and lung cancer tissues. FADS2 circular RNAs are associated with shorter overall survival in colorectal cancer patients. The evidence thusfar supports an effort for future research on the role of FADS2 as a tumor suppressor in a range of neoplastic disorders.

Unlocking the potential of milk whey protein components in colorectal cancer prevention and therapy

Extensive research from large prospective cohort studies and meta-analytical investigations over recent decades have consistently indicated that dairy foods have protective effects, reducing the risk of colorectal cancer. Most of the literature has explored the potential role of milk minerals and vitamins in managing colorectal cancer. Yet, there is a paucity of a comprehensive summary of the anticancer attributes of milk protein components and their underlying mechanisms of action. Recent advancements have spotlighted the potential of whey proteins, including β-lactoglobulin, α-lactalbumin, serum albumin, and lactoferrin, as promising candidates for both the prevention and treatment of colorectal cancer. Notably, whey proteins have demonstrated a more pronounced capacity for suppressing carcinogen-induced tumors when compared to casein. Their strong binding affinity enables them to serve as effective carriers for small molecules or drugs targeting colon cancer therapy. Furthermore, numerous studies have underscored the anti-inflammatory and antioxidant prowess of whey proteins in cancer prevention. Additionally, whey proteins have been shown to trigger apoptosis, hinder tumor cell proliferation, and impede metastasis. This comprehensive review, therefore, not only substantiates the significance of incorporating whey protein components into a balanced daily diet but also underscores their potential in safeguarding against the onset and progression of colorectal cancer.

Dissecting the pathogenic effects of smoking and its hallmarks in blood DNA methylation on colorectal cancer risk

BACKGROUND

Tobacco smoking is suggested as a risk factor for colorectal cancer (CRC), but the complex relationship and the potential pathway are not fully understood.

METHODS

We performed two-sample Mendelian randomisation (MR) analyses with genetic instruments for smoking behaviours and related DNA methylation in blood and summary-level GWAS data of colorectal cancer to disentangle the relationship. Colocalization analyses and prospective gene-environment interaction analyses were also conducted as replication.

RESULTS

Convincing evidence was identified for the pathogenic effect of smoking initiation on CRC risk and suggestive evidence was observed for the protective effect of smoking cessation in the univariable MR analyses. Multivariable MR analysis revealed that these associations were independent of other smoking phenotypes and alcohol drinking. Genetically predicted methylation at CpG site cg17823346 [ZMIZ1] were identified to decrease CRC risk; while genetically predicted methylation at cg02149899 would increase CRC risk. Colocalization and gene-environment interaction analyses added further evidence to the relationship between epigenetic modification at cg17823346 [ZMIZ1] as well as cg02149899 and CRC risk.

DISCUSSION

Our study confirms the significant association between tobacco smoking, DNA methylation and CRC risk and yields a novel insight into the pathogenic effect of tobacco smoking on CRC risk.

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.

Changing Colorectal Cancer Trends in Asians: Epidemiology and Risk Factors

Once an infrequent disease in parts of Asia, the rate of colorectal cancer in recent decades appears to be steadily increasing. Colorectal cancer represents one of the most important causes of cancer mortality worldwide, including in many regions in Asia. Rapid changes in socioeconomic and lifestyle habits have been attributed to the notable increase in the incidence of colorectal cancers in many Asian countries. Through published data from the International Agency for Cancer Research (IARC), we utilized available continuous data to determine which Asian nations had a rise in colorectal cancer rates. We found that East and South East Asian countries had a significant rise in colorectal cancer rates. Subsequently, we summarized here the known genetics and environmental risk factors for colorectal cancer among populations in this region as well as approaches to screening and early detection that have been considered across various countries in the region.

FADS1-arachidonic acid axis enhances arachidonic acid metabolism by altering intestinal microecology in colorectal cancer

Colonocyte metabolism shapes the microbiome. Metabolites are the main mediators of information exchange between intestine and microbial communities. Arachidonic acid (AA) is an essential polyunsaturated fatty acid and its role in colorectal cancer (CRC) remains unexplored. In this study, we show that AA feeding promotes tumor growth in AOM/DSS and intestinal specific Apc^-/- mice via modulating the intestinal microecology of increased gram-negative bacteria. Delta-5 desaturase (FADS1), a rate-limiting enzyme, is upregulated in CRC and effectively mediates AA synthesis. Functionally, FADS1 regulates CRC tumor growth via high AA microenvironment-induced enriched gram-negative microbes. Elimination of gram-negative microbe abolishes FADS1 effect. Mechanistically, gram-negative microbes activate TLR4/MYD88 pathway in CRC cells that contributes FADS1-AA axis to metabolize to prostaglandin E2 (PGE2). Cumulatively, we report a potential cancer-promoting mechanism of FADS1-AA axis in CRC that converts raising synthesized AA to PGE2 via modulating the intestinal microecology of gram-negative.

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.

Developing and validating polygenic risk scores for colorectal cancer risk prediction in East Asians

Several polygenic risk scores (PRSs) have been developed to predict the risk of colorectal cancer (CRC) in European descendants. We used genome-wide association study (GWAS) data from 22 702 cases and 212 486 controls of Asian ancestry to develop PRSs and validated them in two case-control studies (1454 Korean and 1736 Chinese). Eleven PRSs were derived using three approaches: GWAS-identified CRC risk SNPs, CRC risk variants identified through fine-mapping of known risk loci and genome-wide risk prediction algorithms. Logistic regression was used to estimate odds ratios (ORs) and area under the curve (AUC). PRS115-EAS , a PRS with 115 GWAS-reported risk variants derived from East-Asian data, validated significantly better than PRS115-EUR derived from European descendants. In the Korea validation set, OR per SD increase of PRS115-EAS was 1.63 (95% CI = 1.46-1.82; AUC = 0.63), compared with OR of 1.44 (95% CI = 1.29-1.60, AUC = 0.60) for PRS115-EUR . PRS115-EAS/EUR derived using meta-analysis results of both populations slightly improved the AUC to 0.64. Similar but weaker associations were found in the China validation set. Individuals among the highest 5% of PRS115-EAS/EUR have a 2.52-fold elevated CRC risk compared with the medium (41-60th) risk group and have a 12% to 20% risk of developing CRC by age 85. PRSs constructed using results from fine-mapping and genome-wide algorithms did not perform as well as PRS115-EAS and PRS115-EAS/EUR in risk prediction, possibly due to a small sample size. Our results indicate that CRC PRSs are promising in predicting CRC risk in East Asians and highlights the importance of using population-specific data to build CRC risk prediction models.

Genetic architecture of tuberculosis susceptibility: A comprehensive research synopsis, meta-analyses, and epidemiological evidence

To date, many studies have been conducted to investigate associations between variants and tuberculosis risk; however, the results have been inconclusive. Here, we systematically provide a summary of the understanding of the genetic architecture of tuberculosis susceptibility. We searched PubMed, Embase and Web of Science to identify genetic association studies of tuberculosis published through October 31, 2021. We conducted meta-analyses for the genetic association with tuberculosis risk. We graded levels of cumulative epidemiological evidence of significant associations with risk of tuberculosis and false-positive report probability tests. We performed functional annotations for these variants using data from the Encyclopedia of DNA Elements (ENCODE) Project and other databases. We identified 703 eligible articles comprising 298,074 cases and 879,593 controls through screening a total of 24,398 citations. Meta-analyses were conducted for 614 genetic variants in 469 genes or loci. We found 39 variants that were nominally significantly associated with tuberculosis risk. Cumulative epidemiological evidence for a significant association was graded strong for 9 variants in or near 9 genes. Among them, 5 variants were associated with tuberculosis risk in at least three main ethnicity (African, Asian and White) which together explained approximately 9.59% of the familial relative risk of tuberculosis. Data from ENCODE and other databases suggested that 8 of these 9 genetic variants with strong evidence might fall within putative functional regions. Our study summarizes the current literature on the genetic architecture of tuberculosis susceptibility and provides useful data for designing future studies to investigate the genetic association with tuberculosis risk.

Fatty acid desaturase 1 (FADS1) is a cancer marker for patient survival and a potential novel target for precision cancer treatment

Fatty Acid Desaturase-1 (FADS1) or delta 5 desaturase (D5D) is a rate-limiting enzyme involved in the biosynthesis of long-chain polyunsaturated fatty acids (LC-PUFAs), i.e., arachidonic acid (ARA) and eicosapentaenoic (EPA). These LC-PUFAs and their metabolites play essential and broad roles in cancer cell proliferation, metastasis, and tumor microenvironment. However, the role of FADS1 in cancers remains incompletely understood. Utilizing The Cancer Genome Atlas (TCGA) database, we explored the role of FADS1 across different cancer types using multiple bioinformatics and statistical tools. Moreover, we studied the impact of a FADS1 inhibitor (D5D-IN-326) on proliferation of multiple cancer cell lines. We identified that FADS1 gene is a predictor for cancer survival in multiple cancer types. Compared to normal tissue, the mRNA expression of FADS1 is significantly increased in primary tumors while even higher in metastatic and recurrent tumors. Mechanistically, pathway analysis demonstrated that FADS1 is associated with cholesterol biosynthesis and cell cycle control genes. Interestingly, FADS1 expression is higher when TP53 is mutated. Tumors with increased FADS1 expression also demonstrated an increased signatures of fibroblasts and macrophages infiltration among most cancer types. Our in vitro assays showed that D5D-IN-326 significantly inhibited cell proliferation of kidney, colon, breast, and lung cancer cell lines in a dose-dependent manner. Lastly, single nucleotide polymorphisms (SNPs) which are well-established expression quantitative trait loci (eQTLs) for FADS1 in normal human tissues are also significantly correlated with FADS1 expression in tumors of multiple tissue types, potentially serving as a marker to stratify cancer patients with high/low FADS1 expression in their tumor tissue. Our study suggests that FADS1 plays multiple roles in cancer biology and is potentially a novel target for precision cancer treatment.

NXN suppresses metastasis of hepatocellular carcinoma by promoting degradation of Snail through binding to DUB3

The poor prognosis of hepatocellular carcinoma (HCC) could be attributed to its high metastasis rate. Here, we report the role of nucleoredoxin (NXN), a multifunctional redox-active protein, in HCC metastasis. The expression of NXN in HCC tissues was measured by immunohistochemistry. The role of NXN on HCC proliferation was determined by CCK-8, EdU and colony formation assays in vitro and subcutaneous tumor formation model in vivo. Transwell and wound healing assays and tail vein injection model were performed to assess the function of NXN on HCC metastasis. Co-immunoprecipitation assay was performed to examine the interaction among NXN, Snail and DUB3. Our results showed that NXN was downregulated in HCC tissues compared to adjacent liver tissues. Patients with low NXN expression had shorter overall survival (OS) time (P < 0.001) than those with high NXN expression. Biologically, ectopic expression of NXN significantly inhibited the proliferation and metastasis of HCC cells both in vitro and in vivo by suppressing epithelial-mesenchymal transition (EMT). Mechanistically, NXN promoted the ubiquitin-proteasome-mediated degradation of Snail through interaction with DUB3. Further, depletion of Snail abolished NXN-inhibited cell proliferation and metastasis. In summary, NXN suppressed the proliferation and metastasis of HCC by inhibiting DUB3-mediated deubiquitylation of Snail protein. Our study demonstrates that NXN, DUB3 and Snail complex functioned as an important regulatory mechanism of HCC progression and indicates a potential therapeutic approach for the treatment of HCC metastasis.

ZMIZ proteins: partners in transcriptional regulation and risk factors for human disease

Coregulator proteins interact with signal-dependent transcription factors to modify their transcriptional activity. ZMIZ1 and ZMIZ2 (zinc finger MIZ-type containing 1 and 2) are coregulators with nonredundant functions that share unique structural characteristics. Among other interacting domains, they possess a MIZ (Msx-interacting zinc finger) that relates them to members of the protein inhibitor of activated STAT (PIAS) family and provides them the capacity to function as SUMO E3 ligases. The ZMIZ proteins stimulate the activity of various signaling pathways, including the androgen receptor (AR), P53, SMAD3/4, WNT/β-catenin, and NOTCH1 pathways, and interact with the BAF chromatin remodeling complex. Due to their molecular versatility, ZMIZ proteins have pleiotropic effects and thus are important for embryonic development and for human diseases. Both have been widely associated with cancer, and ZMIZ1 has been very frequently identified as a risk allele for several autoimmune conditions and other disorders. Moreover, mutations in the coding region of the ZMIZ1 gene are responsible for a severe syndromic neurodevelopmental disability. Because the actions of coregulators are highly gene-specific, a better knowledge of the associations that exist between the function of the ZMIZ coregulators and human pathologies is expected to potentiate the use of ZMIZ1 and ZMIZ2 as new drug targets for diseases such as hormone-dependent cancers.

Lipid peroxidation biomarkers associated with height and obesity measures in the opposite direction in women

OBJECTIVE

This study aimed to investigate whether and how lipid peroxidation markers are associated with height and obesity measures.

METHODS

In two independent samples of women (Study 1: n = 1,005; Study 2: n = 1,158), systemic levels of lipid peroxidation were assessed by urinary markers F₂ -isoprostanes (F₂ -IsoPs) and its major metabolite (F₂ -IsoP-M), with gas chromatography/negative ion chemical ionization mass spectrometry assays. Anthropometric parameters were directly measured and genetically estimated, and they were used in the primary analysis and in a Mendelian randomization analysis in relation to lipid peroxidation, respectively, with general linear models.

RESULTS

After adjusting for potential confounders, it was found that measured adult height was inversely associated with levels of F₂ -IsoPs (β = -0.89, p < 0.001) and F₂ -IsoP-M (β = -0.71, p = 0.003), whereas obesity measures were positively associated with F₂ -IsoP-M (β = 1.81, p < 0.001 for BMI; and β = 0.77, p < 0.001 for waist circumference). Results were consistent between the two study samples. The opposite associations were further replicated when using genetically determined measures of height and obesity in the Mendelian randomization analysis. Moreover, analyses mutually adjusted for height and obesity measures suggested that these associations were independent of one another.

CONCLUSIONS

This study, for the first time, to our knowledge, reveals that a shared biological process (lipid peroxidation) is associated with both height and obesity measures but in the opposite direction.

Large-scale Integrated Analysis of Genetics and Metabolomic Data Reveals Potential Links Between Lipids and Colorectal Cancer Risk

BACKGROUND

The etiology of colorectal cancer is not fully understood.

METHODS

Using genetic variants and metabolomics data including 217 metabolites from the Framingham Heart Study (n = 1,357), we built genetic prediction models for circulating metabolites. Models with prediction R2 > 0.01 (Nmetabolite = 58) were applied to predict levels of metabolites in two large consortia with a combined sample size of approximately 46,300 cases and 59,200 controls of European and approximately 21,700 cases and 47,400 controls of East Asian (EA) descent. Genetically predicted levels of metabolites were evaluated for their associations with colorectal cancer risk in logistic regressions within each racial group, after which the results were combined by meta-analysis.

RESULTS

Of the 58 metabolites tested, 24 metabolites were significantly associated with colorectal cancer risk [Benjamini-Hochberg FDR (BH-FDR) < 0.05] in the European population (ORs ranged from 0.91 to 1.06; P values ranged from 0.02 to 6.4 × 10-8). Twenty one of the 24 associations were replicated in the EA population (ORs ranged from 0.26 to 1.69, BH-FDR < 0.05). In addition, the genetically predicted levels of C16:0 cholesteryl ester was significantly associated with colorectal cancer risk in the EA population only (OREA: 1.94, 95% CI, 1.60-2.36, P = 2.6 × 10-11; OREUR: 1.01, 95% CI, 0.99-1.04, P = 0.3). Nineteen of the 25 metabolites were glycerophospholipids and triacylglycerols (TAG). Eighteen associations exhibited significant heterogeneity between the two racial groups (PEUR-EA-Het < 0.005), which were more strongly associated in the EA population. This integrative study suggested a potential role of lipids, especially certain glycerophospholipids and TAGs, in the etiology of colorectal cancer.

CONCLUSIONS

This study identified potential novel risk biomarkers for colorectal cancer by integrating genetics and circulating metabolomics data.

IMPACT

The identified metabolites could be developed into new tools for risk assessment of colorectal cancer in both European and EA populations.

Risk Stratification for Early-Onset Colorectal Cancer Using a Combination of Genetic and Environmental Risk Scores: An International Multi-Center Study

BACKGROUND

The incidence of colorectal cancer (CRC) among individuals aged younger than 50 years has been increasing. As screening guidelines lower the recommended age of screening initiation, concerns including the burden on screening capacity and costs have been recognized, suggesting that an individualized approach may be warranted. We developed risk prediction models for early-onset CRC that incorporate an environmental risk score (ERS), including 16 lifestyle and environmental factors, and a polygenic risk score (PRS) of 141 variants.

METHODS

Relying on risk score weights for ERS and PRS derived from studies of CRC at all ages, we evaluated risks for early-onset CRC in 3486 cases and 3890 controls aged younger than 50 years. Relative and absolute risks for early-onset CRC were assessed according to values of the ERS and PRS. The discriminatory performance of these scores was estimated using the covariate-adjusted area under the receiver operating characteristic curve.

RESULTS

Increasing values of ERS and PRS were associated with increasing relative risks for early-onset CRC (odds ratio per SD of ERS = 1.14, 95% confidence interval [CI] = 1.08 to 1.20; odds ratio per SD of PRS = 1.59, 95% CI = 1.51 to 1.68), both contributing to case-control discrimination (area under the curve = 0.631, 95% CI = 0.615 to 0.647). Based on absolute risks, we can expect 26 excess cases per 10 000 men and 21 per 10 000 women among those scoring at the 90th percentile for both risk scores.

CONCLUSIONS

Personal risk scores have the potential to identify individuals at differential relative and absolute risk for early-onset CRC. Improved discrimination may aid in targeted CRC screening of younger, high-risk individuals, potentially improving outcomes.

Prostate Cancer Susceptibility Loci Identified in GATA2 and ZMIZ1 in Chinese Population

Background

Common genetic risk variants for prostate cancer (PCa) have been identified at approximately 170 loci using genome-wide association studies (GWAS), most of which were identified in European populations. Recently, GWAS were applied to a large Japanese cohort and identified 12 novel susceptibility loci associated with PCa risk. In this study, we aim to investigate PCa susceptibility loci in the Chinese population. The study data will be used to promote PCa risk control in China.

Methods

A total of 235 PCa patients and 252 control subjects (all unrelated) were enrolled in this case-control PCa study. Nine single nucleotide polymorphisms (SNPs) were genotyped in GATA2 (rs73862213, rs2335052, and rs10934857), ZMIZ1 (rs704017, rs77911174, and rs3740259), and SUN2 (rs78397383, rs5750680, and rs138705) genes. The associations between the candidate SNPs and PCa were analyzed using multiple-factor logistic regression and haplotype analysis.

Results

The allele frequency distributions of rs73862213 and rs2335052 in the GATA2 gene and rs704017 and rs77911174 in the ZMIZ1 gene were found to be significantly different between PCa cases and controls. Haplotype analysis revealed that the G-C-A haplotype of the GATA2 gene (order of SNPs: rs73862213-rs2335052-rs10934857) and the G-G-G haplotype of the ZMIZ1 gene (order of SNPs: rs704017-rs77911174-rs3740259) were associated with increased PCa risk. None of the SUN2 haplotypes were associated with PCa.

Conclusions

Our study data indicates that the minor alleles of rs73862213 and rs2335052 in the GATA2 gene and rs704017 and rs77911174 in the ZMIZ1 gene were associated with increased PCa risk. These findings greatly extended our knowledge of the etiology of PCa.

Metabolomics and the Multi-Omics View of Cancer

Cancer is widely regarded to be a genetic disease. Indeed, over the past five decades, the genomic perspective on cancer has come to almost completely dominate the field. However, this genome-only view is incomplete and tends to portray cancer as a disease that is highly heritable, driven by hundreds of complex genetic interactions and, consequently, difficult to prevent or treat. New evidence suggests that cancer is not as heritable or purely genetic as once thought and that it really is a multi-omics disease. As highlighted in this review, the genome, the exposome, and the metabolome all play roles in cancer’s development and manifestation. The data presented here show that >90% of cancers are initiated by environmental exposures (the exposome) which lead to cancer-inducing genetic changes. The resulting genetic changes are, then, propagated through the altered DNA of the proliferating cancer cells (the genome). Finally, the dividing cancer cells are nourished and sustained by genetically reprogrammed, cancer-specific metabolism (the metabolome). As shown in this review, all three “omes” play roles in initiating cancer. Likewise, all three “omes” interact closely, often providing feedback to each other to sustain or enhance tumor development. Thanks to metabolomics, these multi-omics feedback loops are now much more evident and their roles in explaining the hallmarks of cancer are much better understood. Importantly, this more holistic, multi-omics view portrays cancer as a disease that is much more preventable, easier to understand, and potentially, far more treatable.

Identification of potential functional variants and genes at 18q21.1 associated with the carcinogenesis of colorectal cancer

Genome-wide association studies (GWAS) have identified more than 160 susceptibility loci for colorectal cancer (CRC). The effects of these variants, particularly their mechanisms, however, remain unclear. In this study, a comprehensive functional annotation of CRC-related GWAS signals was firstly conducted to identify the potential causal variants. We found that the SNP rs7229639 in intron 3 of SMAD7 at 18q21.1 might serve as a putative functional variant in CRC. The SNP rs7229639 is located in a region with evidence of regulatory potential. Dual-luciferase reporter assays revealed that three other SNPs (rs77544449, rs60385309 and rs72917785), in strong linkage disequilibrium (LD) with rs7229639, exhibited allele-specific enhancer activity, of which one of the target genes may conceivably be LIPG, as suggested by eQTL association data and Hi-C data. We also verified that LIPG promoted malignancy of CRC cells in vitro, with supporting clinical data indicating that LIPG is upregulated and correlated with a poor prognosis in CRC. Finally, pitavastatin was observed to exhibit an anti-CRC activity and modest inhibition of LIPG mRNA levels. Collectively, our data suggest that these functional variants at 18q21.1 are involved in the pathogenesis of CRC by modulating enhancer activity, and possibly LIPG expression, thus indicating a promising therapeutic target for CRC. The results of functional annotation in our investigation could also serve as an inventory for CRC susceptibility SNPs and offer guides for post-GWAS downstream functional studies.

In the latest statistics, the incidence and mortality rate of colorectal cancer (CRC) remains high. Genome-wide association studies (GWAS) have become a powerful tool for identifying genetic susceptibility loci that confer significant risk on disease, and have identified more than 160 risk loci associated with CRC. However, it has proven quite difficult to identify the regulatory variants and target genes involved behind these GWAS signals. Here, we take advantage of multi-omics data and multiple biological experiments to reveal new biological pathways affecting susceptibility to CRC. We show that a specific genetic variant, rs7229639, and three other high linked functional variants (rs77544449, rs60385309 and rs72917785) at 18q21.1 might regulate the expression of LIPG, a gene that was shown to exhibit an oncogenic function by our in-vitro experiments and clinical data analysis. The link between genetic variants, gene expression and CRC phenotype established by us could provide references for follow-up basic and clinical studies.

Coding and regulatory variants are associated with serum protein levels and disease

Circulating proteins can be used to diagnose and predict disease-related outcomes. A deep serum proteome survey recently revealed close associations between serum protein networks and common disease. In the current study, 54,469 low-frequency and common exome-array variants were compared to 4782 protein measurements in the serum of 5343 individuals from the AGES Reykjavik cohort. This analysis identifies a large number of serum proteins with genetic signatures overlapping those of many diseases. More specifically, using a study-wide significance threshold, we find that 2021 independent exome array variants are associated with serum levels of 1942 proteins. These variants reside in genetic loci shared by hundreds of complex disease traits, highlighting serum proteins' emerging role as biomarkers and potential causative agents of a wide range of diseases.

Colorectal cancer risk susceptibility loci in a Swedish population

To search for colorectal cancer (CRC) risk loci, Swedish samples were used for a genome-wide haplotype analysis. A logistic regression model was employed in 2663 CRC cases and 1642 controls in the discovery analysis. Three analyses were done, on all, familial-, and nonfamilial CRC samples and only results with odds ratio (OR) > 1 were analyzed. single nucleotide polymorphism (SNP) analysis did not generate any statistically significant results. Haplotype analysis suggested novel loci, on chromosome 2q36.1 (OR = 1.71, p value = 5.6924 × 10^-8 ) in all CRC samples, chromosome 1q43 (OR = 4.04 p value = 3.24 × 10^-8 ) in familial CRC samples, and two hits in nonfamilial CRC samples, chromosomes 2q36.1 (OR = 1.71 p value = 5.69 × 10^-8 ) and 3p24.3 (OR = 1.62 p value = 6.21 × 10^-9 ). Moreover, one locus on chromosome 20q13.33 was suggested in analyses of all samples, and five more novel loci were suggested on chromosomes 10q25.3, 15q,22.31, 17p11.2, 1p34.2, and 3q24. The haplotypes from the analysis of all samples were replicated in a second study of CRC cases and controls from the same part of Sweden. In summary, using haplotype analysis in Swedish CRC samples, the best hits were novel loci and the locus on chromosomes 2q36.1 and 20q13.33 suggested in the analysis of all samples were confirmed in a second cohort. The ORs were often higher than ORs from published genome-wide association study (GWAS). The study suggested it was possible that a risk locus could involve more than one gene, and that haplotypes could give information on the gene or genes possibly involved in the risk at specific locus.

zmiz1a zebrafish mutants have defective erythropoiesis, altered expression of autophagy genes, and a deficient response to vitamin D

ZMIZ1 is a transcriptional coactivator that is related to members of the protein inhibitor of activated STAT (PIAS) family. ZMIZ1 regulates the activity of various transcription factors including the androgen receptor, p53, and Smad3. ZMIZ1 also interacts with Notch1 and selectively regulates Notch1 target genes relevant for T cell development and leukemogenesis in mammals. Human ZMIZ1 is additionally characterized as a latitude-dependent autoimmune disease (LDAD) risk gene, as it is responsive to vitamin D and has been associated with at least eleven blood cell traits. To address the function of ZMIZ1 in fish, we introduced CRISPR/Cas9 mutations in the zmiz1a gene in zebrafish. We observed that inactivation of zmiz1a in developing zebrafish larvae results in lethality at 15 days post fertilization (dpf) and delayed erythroid maturation. Differential gene expression analysis indicated that 15 dpf zmiz1a-null larvae had altered expression of autophagy genes, and erythrocytes that lacked Zmiz1a function exhibited an accumulation of mitochondrial DNA. Furthermore, we observed that autophagy gene expression was dysregulated at earlier stages of development, which suggests the involvement of Zmiz1a in the regulation of autophagy genes beyond the process of red blood cell differentiation. Finally, we showed that the loss of Zmiz1a decreased the capacity of the embryos to respond to vitamin D, indicating additional participation of Zmiz1a as a mediator of vitamin D activity.

Association of four genetic variants with colorectal cancer in Kazakhstan population

The study was conducted to search for polymorphisms located in the 10th chromosome associated with colorectal adenocarcinoma in representatives of the Kazakhstan population. Study was performed with 282 colorectal cancer (CRC) patients and 159 controls. Genotyping of SNPs was performed by QuantStudio 12K Flex PCR. For four significant SNPs inheritance model analysis was performed. Increasing risk of CRC was noted for rs10795668 in log-additive model (OR = 1.45, 95% CI: 1.05–1.99, p = 0.023); for rs1035209 in log-additive model (OR = 1.79, 95% CI: 1.18–2.72, p = 0.003); for rs11190164 in log-additive model (OR = 1.67, 95% CI: 1.17–2.38, p = 0.004). Decreasing risk of CRC was noted for rs10506868 in log-additive model (OR = 0.56, 95% CI: 0.37–0.85, p = 0.006). We detected SNPs that are associated with CRC risk in the Kazakhstan population.

Repurposing non-invasive prenatal testing data: Population study of single nucleotide variants associated with colorectal cancer and Lynch syndrome

In our previous work, genomic data generated through non-invasive prenatal testing (NIPT) based on low-coverage massively parallel whole-genome sequencing of total plasma DNA of pregnant women in Slovakia was described as a valuable source of population specific data. In the present study, these data were used to determine the population allele frequency of common risk variants located in genes associated with colorectal cancer (CRC) and Lynch syndrome (LS). Allele frequencies of identified variants were compared with six world populations to detect significant differences between populations. Finally, variants were interpreted, functional consequences were searched for and clinical significance of variants was investigated using publicly available databases. Although the present study did not identify any pathogenic variants associated with CRC or LS in the Slovak population using NIPT data, significant differences were observed in the allelic frequency of risk CRC variants previously reported in genome-wide association studies and common variants located in genes associated with LS. As Slovakia is one of the leading countries with the highest incidence of CRC among male patients in the world, there is a need for studies dedicated to investigating the cause of such a high incidence of CRC in Slovakia. The present study also assumed that extensive cross-country data aggregation of NIPT results would represent an unprecedented source of information concerning human genome variation in cancer research.

Delta-5-desaturase: A novel therapeutic target for cancer management

•

D5D is an independent prognostic factor in cancer.

•

D5D aggravates cancer progression via mediating AA/PGE₂ production from DGLA.

•

AA/PGE₂ promotes cancer progression via regulating the tumor microenvironment.

•

Inhibition of D5D redirects COX-2 catalyzed DGLA peroxidation, producing 8-HOA.

•

8-HOA suppress cancer by regulating proliferation, apoptosis, and metastasis.

Delta-5 desaturase (D5D) is a rate-limiting enzyme that introduces double-bonds to the delta-5 position of the n-3 and n-6 polyunsaturated fatty acid chain. Since fatty acid metabolism is a vital factor in cancer development, several recent studies have revealed that D5D activity and expression could be an independent prognostic factor in cancers. However, the mechanistic basis of D5D in cancer progression is still controversial. The classical concept believes that D5D could aggravate cancer progression via mediating arachidonic acid (AA)/prostaglandin E₂ production from dihomo-γ-linolenic acid (DGLA), resulting in activation of EP receptors, inflammatory pathways, and immunosuppression. On the contrary, D5D may prevent cancer progression through activating ferroptosis, which is iron-dependent cell death. Suppression of D5D by RNA interference and small-molecule inhibitor has been identified as a promising anti-cancer strategy. Inhibition of D5D could shift DGLA peroxidation pattern from generating AA to a distinct anti-cancer free radical byproduct, 8-hydroxyoctanoic acid, resulting in activation of apoptosis pathway and simultaneously suppression of cancer cell survival, proliferation, migration, and invasion. Hence, understanding the molecular mechanisms of D5D on cancer may therefore facilitate the development of novel therapeutical applications. Given that D5D may serve as a promising target in cancer, in this review, we provide an updated summary of current knowledge on the role of D5D in cancer development and potentially useful therapeutic strategies.

Non-Lynch Familial and Early-Onset Colorectal Cancer Explained by Accumulation of Low-Risk Genetic Variants

A large proportion of familial and/or early-onset cancer patients do not carry pathogenic variants in known cancer predisposing genes. We aimed to assess the contribution of previously validated low-risk colorectal cancer (CRC) alleles to familial/early-onset CRC (fCRC) and to serrated polyposis. We estimated the association of CRC with a 92-variant-based weighted polygenic risk score (wPRS) using 417 fCRC patients, 80 serrated polyposis patients, 1077 hospital-based incident CRC patients, and 1642 controls. The mean wPRS was significantly higher in fCRC than in controls or sporadic CRC patients. fCRC patients in the highest (20th) wPRS quantile were at four-fold greater CRC risk than those in the middle quantile (10th). Compared to low-wPRS fCRC, a higher number of high-wPRS fCRC patients had developed multiple primary CRCs, had CRC family history, and were diagnosed at age ≥50. No association with wPRS was observed for serrated polyposis. In conclusion, a relevant proportion of mismatch repair (MMR)-proficient fCRC cases might be explained by the accumulation of low-risk CRC alleles. Validation in independent cohorts and development of predictive models that include polygenic risk score (PRS) data and other CRC predisposing factors will determine the implementation of PRS into genetic testing and counselling in familial and early-onset CRC.

Impact of Amerind ancestry and FADS genetic variation on omega-3 deficiency and cardiometabolic traits in Hispanic populations

Long chain polyunsaturated fatty acids (LC-PUFAs) have critical signaling roles that regulate dyslipidemia and inflammation. Genetic variation in the FADS gene cluster accounts for a large portion of interindividual differences in circulating and tissue levels of LC-PUFAs, with the genotypes most strongly predictive of low LC-PUFA levels at strikingly higher frequencies in Amerind ancestry populations. In this study, we examined relationships between genetic ancestry and FADS variation in 1102 Hispanic American participants from the Multi-Ethnic Study of Atherosclerosis. We demonstrate strong negative associations between Amerind genetic ancestry and LC-PUFA levels. The FADS rs174537 single nucleotide polymorphism (SNP) accounted for much of the AI ancestry effect on LC-PUFAs, especially for low levels of n-3 LC-PUFAs. Rs174537 was also strongly associated with several metabolic, inflammatory and anthropomorphic traits including circulating triglycerides (TGs) and E-selectin in MESA Hispanics. Our study demonstrates that Amerind ancestry provides a useful and readily available tool to identify individuals most likely to have FADS-related n-3 LC-PUFA deficiencies and associated cardiovascular risk.

Implications of Lifestyle Factors and Polygenic Risk Score for Absolute Risk Prediction of Colorectal Neoplasm and Risk-Adapted Screening

Background: Estimation of absolute risk of developing colorectal neoplasm is essential for personalized colorectal cancer (CRC) screening. We developed models to determine relative and absolute risks of colorectal neoplasm based on lifestyle and genetic variants and to validate their application in risk-adapted screening.

Methods: We prospectively collected data from 203 advanced neoplasms, 464 non-advanced adenomas, and 1,213 healthy controls from a CRC screening trial in China in 2018–2019. The risk prediction model based on four lifestyle factors and a polygenic risk score (PRS) consisted of 19 CRC-associated single-nucleotide polymorphisms. We assessed the relative and 10-year absolute risks of developing colorectal neoplasm and the yield of a risk-adapted screening approach incorporating risk models, fecal immunochemical test, and colonoscopy.

Results: Compared to the participants with favorable lifestyle and lower PRS, those with unfavorable lifestyle and higher PRS had 2.87- and 3.79-fold higher risk of colorectal neoplasm in males and females, respectively. For a 50-year-old man or a 50-year-old woman with the highest risk profile, the estimated 10-year absolute risk of developing colorectal neoplasm was 6.59% (95% CI: 6.53–6.65%) and 4.19% (95% CI: 4.11–4.28%), respectively, compared to 2.80% (95% CI: 2.78–2.81%) for men and 2.24% (95% CI: 2.21–2.27%) for women with the lowest risk profile. The positive predictive value for advanced neoplasm was 31.7%, and the number of colonoscopies needed to detect one advanced neoplasm was 3.2.

Conclusion: The risk models, absolute risk estimates, and risk-adapted screening presented in our study would contribute to developing effective personalized CRC prevention and screening strategies.

Genetic risk factors for colorectal cancer in multiethnic Indonesians

Colorectal cancer is a common cancer in Indonesia, yet it has been understudied in this resource-constrained setting. We conducted a genome-wide association study focused on evaluation and preliminary discovery of colorectal cancer risk factors in Indonesians. We administered detailed questionnaires and collecting blood samples from 162 colorectal cancer cases throughout Makassar, Indonesia. We also established a control set of 193 healthy individuals frequency matched by age, sex, and ethnicity. A genome-wide association analysis was performed on 84 cases and 89 controls passing quality control. We evaluated known colorectal cancer genetic variants using logistic regression and established a genome-wide polygenic risk model using a Bayesian variable selection technique. We replicate associations for rs9497673, rs6936461 and rs7758229 on chromosome 6; rs11255841 on chromosome 10; and rs4779584, rs11632715, and rs73376930 on chromosome 15. Polygenic modeling identified 10 SNP associated with colorectal cancer risk. This work helps characterize the relationship between variants in the SCL22A3, SCG5, GREM1, and STXBP5-AS1 genes and colorectal cancer in a diverse Indonesian population. With further biobanking and international research collaborations, variants specific to colorectal cancer risk in Indonesians will be identified.

Identifying Novel Susceptibility Genes for Colorectal Cancer Risk From a Transcriptome-Wide Association Study of 125,478 Subjects

BACKGROUND AND AIMS

Susceptibility genes and the underlying mechanisms for the majority of risk loci identified by genome-wide association studies (GWAS) for colorectal cancer (CRC) risk remain largely unknown. We conducted a transcriptome-wide association study (TWAS) to identify putative susceptibility genes.

METHODS

Gene-expression prediction models were built using transcriptome and genetic data from the 284 normal transverse colon tissues of European descendants from the Genotype-Tissue Expression (GTEx), and model performance was evaluated using data from The Cancer Genome Atlas (n = 355). We applied the gene-expression prediction models and GWAS data to evaluate associations of genetically predicted gene-expression with CRC risk in 58,131 CRC cases and 67,347 controls of European ancestry. Dual-luciferase reporter assays and knockdown experiments in CRC cells and tumor xenografts were conducted.

RESULTS

We identified 25 genes associated with CRC risk at a Bonferroni-corrected threshold of P < 9.1 × 10-6, including genes in 4 novel loci, PYGL (14q22.1), RPL28 (19q13.42), CAPN12 (19q13.2), MYH7B (20q11.22), and MAP1L3CA (20q11.22). In 9 known GWAS-identified loci, we uncovered 9 genes that have not been reported previously, whereas 4 genes remained statistically significant after adjusting for the lead risk variant of the locus. Through colocalization analysis in GWAS loci, we additionally identified 12 putative susceptibility genes that were supported by TWAS analysis at P < .01. We showed that risk allele of the lead risk variant rs1741640 affected the promoter activity of CABLES2. Knockdown experiments confirmed that CABLES2 plays a vital role in colorectal carcinogenesis.

CONCLUSIONS

Our study reveals new putative susceptibility genes and provides new insight into the biological mechanisms underlying CRC development.

Performance of common genetic variants in risk prediction for colorectal cancer in Chinese: A two-stage and multicenter study

The efficacy of susceptible variants derived from genome-wide association studies (GWAs) optimizing discriminatory accuracy of colorectal cancer (CRC) in Chinese remains unclear. In the present validation study, we assessed 75 recently identified variants from GWAs. A risk predictive model combining 19 variants using the least absolute shrinkage and selection operator (LASSO) statistics offered certain clinical advantages. This model demonstrated an area under the receiver operating characteristic (AUC) of 0.61 during training analysis and yielded robust AUCs from 0.59 to 0.61 during validation analysis in three independent centers. The individuals carrying the highest quartile of risk score revealed over 2-fold risks of CRC (ranging from 2.12 to 2.90) compared with those who presented the lowest quartile of risk score. This genetic model offered the possibility of partitioning risk within the average risk population, which might serve as a first step toward developing individualized CRC prevention strategies in China.

SNPs associated with colorectal cancer at 15q13.3 affect risk enhancers that modulate GREM1 gene expression

Several genome wide association studies of colorectal cancer (CRC) have identified single nucleotide polymorphisms (SNPs) on chromosome 15q13.3 associated with CRC risk. To identify functional variant(s) underlying this association, we investigated SNPs in linkage disequilibrium with the risk‐associated SNP rs4779584 that overlapped regulatory regions/enhancer elements characterized in colon‐related tissues and cells. We identified several SNP‐containing regulatory regions that exhibited enhancer activity in vitro, including one SNP (rs1406389) that correlated with allele‐specific effects on enhancer activity. Deletion of either this enhancer or another enhancer that had previously been reported in this region correlated with decreased expression of GREM1 following CRISPR/Cas9 genome editing. That GREM1 is one target of these enhancers was further supported by an expression quantitative trait loci correlation between rs1406389 and GREM1 expression in the transverse but not sigmoid colon in the Genotype‐Tissue Expression dataset. Taken together, we conclude that the 15q13.3 region contains at least two functional variants that map to distinct enhancers and impact CRC risk through modulation of GREM1 expression.

Multi-omics analysis to identify susceptibility genes for colorectal cancer

Most genetic variants for colorectal cancer (CRC) identified in genome-wide association studies (GWAS) are located in intergenic regions, implying pathogenic dysregulations of gene expression. However, comprehensive assessments of target genes in CRC remain to be explored. We conducted a multi-omics analysis using transcriptome and/or DNA methylation data from the Genotype-Tissue Expression, The Cancer Genome Atlas and the Colonomics projects. We identified 116 putative target genes for 45 GWAS-identified variants. Using summary-data-based Mendelian randomization approach (SMR), we demonstrated that the CRC susceptibility for 29 out of the 45 CRC variants may be mediated by cis-effects on gene regulation. At a cutoff of the Bonferroni-corrected PSMR < 0.05, we determined 66 putative susceptibility genes, including 39 genes that have not been previously reported. We further performed in vitro assays for two selected genes, DIP2B and SFMBT1, and provide functional evidence that they play a vital role in colorectal carcinogenesis via disrupting cell behavior, including migration, invasion and epithelial-mesenchymal transition. Our study reveals a large number of putative novel susceptibility genes and provides additional insight into the underlying mechanisms for CRC genetic risk loci.

Association between ALDH2 and ADH1B Polymorphisms and the Risk for Colorectal Cancer in Koreans

Purpose

Excessive alcohol consumption has been linked to an increased risk of colorectal cancer (CRC). We evaluated the association between alcohol-related genetic variants and CRC risk.

Materials and Methods

The study cohort consisted of 5,435 CRC cases and 3,553 population-based cancer-free controls. Genotype data were generated from germline DNA using the Infinium OncoArray-500K BeadChip in 2,535 cases and 2,287 controls and the Infinium Multi-Ethnic Global BeadChip in 2,900 cases and 1,266 controls. The associations between aldehyde dehydrogenase 2 (ALDH2) rs671 and alcohol dehydrogenase 1B (ADH1B) rs1229984 polymorphisms and CRC risk were assessed using multivariate logistic regression analyses.

Results

Compared with the major homozygous ALDH2 genotype (GG), heterozygous or minor homozygous ALDH2 genotype (GA or AA, related to a low alcohol consumption) was significantly associated with a reduced risk for CRC in men (odds ratio [OR], 0.78; 95% confidence interval [CI], 0.68 to 0.90), but not in women (OR, 0.70; 95% CI, 0.47 to 1.05). A stronger association was found among regular drinkers (OR, 0.58; 95% CI, 0.47 to 0.71 in men and OR, 0.33; 95% CI, 0.18 to 0.58 in women). No association of CRC risk with ADH1B rs1229984 genotype was found. The association between alcohol-related combined genotypes and risk of CRC was significant (p for linear=0.001). The combined genotype with the highest genetically predicted alcohol consumption (ALDH2 rs671 GG and ADH1B rs1229984 AG/GG) was associated with a high risk for CRC (OR, 1.35; 95% CI, 1.11 to 1.63).

Conclusion

Our study provides strong evidence for a possible causal association between alcohol consumption and CRC risk.

Association of MDM4 Gene rs4245739 Polymorphism with the Risk and Clinical Characteristics of Colorectal Cancer in a Chinese Han Population

Background

Studies show that MDM4 may play a pivotal role in colorectal cancer (CRC). Recently, a host of studies suggest that MDM4 gene rs4245739 polymorphism may modify the risk of different cancers.

Methods

In this study, we were interested whether MDM4 gene rs4245739 polymorphism correlated with the risk and clinical characteristics of CRC. Logistic regression was adopted to estimate the association of rs4245739 polymorphism and CRC risk.

Results

We enrolled 444 CRC patients and 530 controls and found MDM4 gene rs4245739 polymorphism may decrease the risk of CRC. Stratified analyses uncovered that this variant was connected to a less risk of CRC in females, non-drinkers, non-smokers, and people under 60 years old. Additionally, rs4245739 polymorphism was related to TNM staging, pathological type, tumor size, and location of CRC. Furthermore, this polymorphism was significantly linked with the survival of CRC.

Conclusion

Totally, this study suggests that MDM4 rs4245739 polymorphism is linked with the risk and clinical characteristics of CRC.

High dimensional model representation of log likelihood ratio: binary classification with SNP data

BACKGROUND

Developing binary classification rules based on SNP observations has been a major challenge for many modern bioinformatics applications, e.g., predicting risk of future disease events in complex conditions such as cancer. Small-sample, high-dimensional nature of SNP data, weak effect of each SNP on the outcome, and highly non-linear SNP interactions are several key factors complicating the analysis. Additionally, SNPs take a finite number of values which may be best understood as ordinal or categorical variables, but are treated as continuous ones by many algorithms.

METHODS

We use the theory of high dimensional model representation (HDMR) to build appropriate low dimensional glass-box models, allowing us to account for the effects of feature interactions. We compute the second order HDMR expansion of the log-likelihood ratio to account for the effects of single SNPs and their pairwise interactions. We propose a regression based approach, called linear approximation for block second order HDMR expansion of categorical observations (LABS-HDMR-CO), to approximate the HDMR coefficients. We show how HDMR can be used to detect pairwise SNP interactions, and propose the fixed pattern test (FPT) to identify statistically significant pairwise interactions.

RESULTS

We apply LABS-HDMR-CO and FPT to synthetically generated HAPGEN2 data as well as to two GWAS cancer datasets. In these examples LABS-HDMR-CO enjoys superior accuracy compared with several algorithms used for SNP classification, while also taking pairwise interactions into account. FPT declares very few significant interactions in the small sample GWAS datasets when bounding false discovery rate (FDR) by 5%, due to the large number of tests performed. On the other hand, LABS-HDMR-CO utilizes a large number of SNP pairs to improve its prediction accuracy. In the larger HAPGEN2 dataset FTP declares a larger portion of SNP pairs used by LABS-HDMR-CO as significant.

CONCLUSION

LABS-HDMR-CO and FPT are interesting methods to design prediction rules and detect pairwise feature interactions for SNP data. Reliably detecting pairwise SNP interactions and taking advantage of potential interactions to improve prediction accuracy are two different objectives addressed by these methods. While the large number of potential SNP interactions may result in low power of detection, potentially interacting SNP pairs, of which many might be false alarms, can still be used to improve prediction accuracy.

Using genetics to decipher the link between type 2 diabetes and cancer: shared aetiology or downstream consequence?

Recent developments in the field of genetics have accelerated our understanding of the aetiology of complex diseases. Type 2 diabetes mellitus and cancer are no exception, with large-scale genome-wide association studies (GWAS) facilitating exploration of the underlying pathology. Here, we discuss how genetics studies can be used to investigate the relationship between these complex diseases. Observational epidemiological studies consistently report that people with type 2 diabetes have a higher risk of several types of cancer. Indeed, type 2 diabetes and cancer share many common risk factors, such as obesity, ageing, poor diet and low levels of physical activity. However, questions remain regarding the biological mechanisms that link these two diseases. Large-scale GWAS of type 2 diabetes and cancer allow us to consider the evidence for shared genetic architecture. Several shared susceptibility genes have been identified, yet tissue specificity and direction of effect must be taken into account when considering common genetic aetiology. We also consider how GWAS, and associated techniques such as Mendelian randomisation, allow us to dissect the link between the two diseases and address questions such as 'Does type 2 diabetes cause cancer or is the increased risk observed driven by higher adiposity or another associated metabolic feature?' Graphical abstract.

Genome-wide association meta-analysis identifies GP2 gene risk variants for pancreatic cancer

Pancreatic cancer is the fourth leading cause of cancer-related deaths in Japan. To identify risk loci, we perform a meta-analysis of three genome-wide association studies comprising 2,039 pancreatic cancer patients and 32,592 controls in the Japanese population. Here, we identify 3 (13q12.2, 13q22.1, and 16p12.3) genome-wide significant loci (P < 5.0 × 10⁻⁸), of which 16p12.3 has not been reported in the Western population. The lead single nucleotide polymorphism (SNP) at 16p12.3 is rs78193826 (odds ratio = 1.46, 95% confidence interval = 1.29-1.66, P = 4.28 × 10⁻⁹), an Asian-specific, nonsynonymous glycoprotein 2 (GP2) gene variant. Associations between selected GP2 gene variants and pancreatic cancer are replicated in 10,822 additional cases and controls of East Asian origin. Functional analyses using cell lines provide supporting evidence of the effect of rs78193826 on KRAS activity. These findings suggest that GP2 gene variants are probably associated with pancreatic cancer susceptibility in populations of East Asian ancestry.

Previous genome-wide association studies have identified risk loci for pancreatic cancer but were centered on individuals of European ancestry. Here the authors identify GP2 gene variants associated with pancreatic cancer susceptibility in populations of East Asian ancestry.