Meta-analysis of genome-wide association data identifies four new susceptibility loci for colorectal cancer

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.

Synchronous gastric and colon cancers: Important to consider hereditary syndromes and chronic inflammatory disease associations

In this editorial we comment on the manuscript, describing management and surveillance strategies in synchronous and metachronous, gastric and colon cancers. Synchronous or metachronous primary malignancies at different sites of the gastrointestinal tract pose a unique diagnostic and therapeutic challenge. Multidisciplinary services and strategies are required for the management of multiple site primary malignancies, to provide the best oncological outcomes. Although this study highlights the dual cancers in 76 sporadic cases, the authors excluded 55 patients due to combination of factors which includes; incomplete clinical data, genetic syndrome, gastric stump cancers. In addition, the authors did not elaborate if any patients presented with signet ring cell morphology, E-cadherin mutations or presence of inflammatory bowel disease. Genetic and mutational errors and epithelial field defects from chronic inflammatory diseases of the gastrointestinal tract are important when considering synchronous gastric and colonic cancers. We will briefly discuss these in this editorial.

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.

BMP signaling in cancer stemness and differentiation

The BMP (Bone morphogenetic protein) signaling pathway plays a central role in metazoan biology, intricately shaping embryonic development, maintaining tissue homeostasis, and influencing disease progression. In the context of cancer, BMP signaling exhibits context-dependent dynamics, spanning from tumor suppression to promotion. Cancer stem cells (CSCs), a modest subset of neoplastic cells with stem-like attributes, exert substantial influence by steering tumor growth, orchestrating therapy resistance, and contributing to relapse. A comprehensive grasp of the intricate interplay between CSCs and their microenvironment is pivotal for effective therapeutic strategies. Among the web of signaling pathways orchestrating cellular dynamics within CSCs, BMP signaling emerges as a vital conductor, overseeing CSC self-renewal, differentiation dynamics, and the intricate symphony within the tumor microenvironment. Moreover, BMP signaling's influence in cancer extends beyond CSCs, intricately regulating cellular migration, invasion, and metastasis. This multifaceted role underscores the imperative of comprehending BMP signaling's contributions to cancer, serving as the foundation for crafting precise therapies to navigate multifaceted challenges posed not only by CSCs but also by various dimensions of cancer progression. This article succinctly encapsulates the diverse roles of the BMP signaling pathway across different cancers, spanning glioblastoma multiforme (GBM), diffuse intrinsic pontine glioma (DIPG), colorectal cancer, acute myeloid leukemia (AML), lung cancer, prostate cancer, and osteosarcoma. It underscores the necessity of unraveling underlying mechanisms and molecular interactions. By delving into the intricate tapestry of BMP signaling's engagement in cancers, researchers pave the way for meticulously tailored therapies, adroitly leveraging its dualistic aspects-whether as a suppressor or promoter-to effectively counter the relentless march of tumor progression.

Long-Term Simulation of Microgravity Induces Changes in Gene Expression in Breast Cancer Cells

Microgravity changes the gene expression pattern in various cell types. This study focuses on the breast cancer cell lines MCF-7 (less invasive) and MDA-MB-231 (triple-negative, highly invasive). The cells were cultured for 14 days under simulated microgravity (s-µg) conditions using a random positioning machine (RPM). We investigated cytoskeletal and extracellular matrix (ECM) factors as well as focal adhesion (FA) and the transmembrane proteins involved in different cellular signaling pathways (MAPK, PAM and VEGF). The mRNA expressions of 24 genes of interest (TUBB, ACTB, COL1A1, COL4A5, LAMA3, ITGB1, CD44, VEGF, FLK1, EGFR, SRC, FAK1, RAF1, AKT1, ERK1, MAPK14, MAP2K1, MTOR, RICTOR, VCL, PXN, CDKN1, CTNNA1 and CTNNB1) were determined by quantitative real-time PCR (qPCR) and studied using STRING interaction analysis. Histochemical staining was carried out to investigate the morphology of the adherent cells (ADs) and the multicellular spheroids (MCSs) after RPM exposure. To better understand this experimental model in the context of breast cancer patients, a weighted gene co-expression network analysis (WGCNA) was conducted to obtain the expression profiles of 35 breast cell lines from the HMS LINCS Database. The qPCR-verified genes were searched in the mammalian phenotype database and the human genome-wide association studies (GWAS) Catalog. The results demonstrated the positive association between the real metastatic microtumor environment and MCSs with respect to the extracellular matrix, cytoskeleton, morphology, different cellular signaling pathway key proteins and several other components. In summary, the microgravity-engineered three-dimensional MCS model can be utilized to study breast cancer cell behavior and to assess the therapeutic efficacies of drugs against breast cancer in the future.

An integrated analysis of the cancer genome atlas data discovers a hierarchical association structure across thirty three cancer types

Cancer cells harbor molecular alterations at all levels of information processing. Genomic/epigenomic and transcriptomic alterations are inter-related between genes, within and across cancer types and may affect clinical phenotypes. Despite the abundant prior studies of integrating cancer multi-omics data, none of them organizes these associations in a hierarchical structure and validates the discoveries in extensive external data. We infer this Integrated Hierarchical Association Structure (IHAS) from the complete data of The Cancer Genome Atlas (TCGA) and compile a compendium of cancer multi-omics associations. Intriguingly, diverse alterations on genomes/epigenomes from multiple cancer types impact transcriptions of 18 Gene Groups. Half of them are further reduced to three Meta Gene Groups enriched with (1) immune and inflammatory responses, (2) embryonic development and neurogenesis, (3) cell cycle process and DNA repair. Over 80% of the clinical/molecular phenotypes reported in TCGA are aligned with the combinatorial expressions of Meta Gene Groups, Gene Groups, and other IHAS subunits. Furthermore, IHAS derived from TCGA is validated in more than 300 external datasets including multi-omics measurements and cellular responses upon drug treatments and gene perturbations in tumors, cancer cell lines, and normal tissues. To sum up, IHAS stratifies patients in terms of molecular signatures of its subunits, selects targeted genes or drugs for precision cancer therapy, and demonstrates that associations between survival times and transcriptional biomarkers may vary with cancer types. These rich information is critical for diagnosis and treatments of cancers.

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.

Role of the E3 ubiquitin-ligase Hakai in intestinal inflammation and cancer bowel disease

The E3 ubiquitin-ligases are important for cellular protein homeostasis and their deregulation is implicated in cancer. The E3 ubiquitin-ligase Hakai is involved in tumour progression and metastasis, through the regulation of the tumour suppressor E-cadherin. Hakai is overexpressed in colon cancer, however, the implication in colitis-associated cancer is unknown. Here, we investigated the potential role of Hakai in intestinal inflammation and cancer bowel disease. Several mouse models of colitis and associated cancer were used to analyse Hakai expression by immunohistochemistry. We also analysed Hakai expression in patients with inflamed colon biopsies from ulcerative colitis and Crohn's disease. By Hakai interactome analysis, it was identified Fatty Acid Synthase (FASN) as a novel Hakai-interacting protein. Moreover, we show that Hakai induces FASN ubiquitination and degradation via lysosome, thus regulating FASN-mediated lipid accumulation. An inverse expression of FASN and Hakai was detected in inflammatory AOM/DSS mouse model. In conclusion, Hakai regulates FASN ubiquitination and degradation, resulting in the regulation of FASN-mediated lipid accumulation, which is associated to the development of inflammatory bowel disease. The interaction between Hakai and FASN may be an important mechanism for the homeostasis of intestinal barrier function and in the pathogenesis of this disease.

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.

Identification and functional validation of HLA-C as a potential gene involved in colorectal cancer in the Korean population

BACKGROUND

Colorectal cancer (CRC) is the third most common cancer worldwide and is influenced by environmental and genetic factors. Although numerous genetic loci for CRC have been identified, the overall understanding of the genetic factors is yet to be elucidated. We sought to discover new genes involved in CRC applying genetic association analysis and functional study.

RESULTS

We conducted exome array analysis on 194 CRC and 600 control subjects for discovering new candidate CRC genes. Fisher's exact test detected one exome-wide significant functional locus for CRC on SMCO1 (P < 10^-6) and two suggestive functional loci on HLA-C and NUTM1 (10^-6 ≤ P < 10^-4). To evaluate the biological role of three candidate CRC genes, the differential expression of these genes between CRC and non-cancer colorectal cells was analyzed using qRT-PCR and publicly available gene expression data. Of three genes, HLA-C consistently revealed the significant down-regulation in CRC cells. In addition, we detected a reduction in cell viability in the HLA-C overexpression CRC cell line, implying the functional relevance of HLA-C in CRC. To understand the underlying mechanism exerted by HLA-C in CRC development, we conducted RNA sequencing analyses of HLA-C overexpression CRC cells and non-cancer colorectal cells. Pathway analysis detected that significantly down-regulated genes in HLA-C overexpression CRC cells were highly enriched in cancer-related signaling pathways such as JAK/STAT, ErbB, and Hedgehog signaling pathways.

CONCLUSIONS

Exome array CRC case-control analysis followed by functional validation demonstrated that HLA-C likely exerts its influence on CRC development via cancer-related signaling pathways.

A Bibliometric Analysis Based on Web of Science: Current Perspectives and Potential Trends of SMAD7 in Oncology

Background: The number of publications on SMAD7 in the field of oncology is increasing rapidly with an upward tendency. In most cases, the mechanisms of carcinogenesis usually relate to disorders of signaling activity. Considering the crucial role of SMAD7 in the crosstalk of multiple signaling pathways, it is necessary to clarify and define the dominant research topics, core authors, and their cumulative research contributions, as well as the cooperative relationships among documents or researchers.

Methods: Altogether, 3477 documents were retrieved from the Web of Science Core Collection with the following criteria: TS= (SMAD7 OR SMAD7-protein OR Small-Mothers-Against-Decapentaplegic-7) refined by WEB OF SCIENCE CATEGORY (ONCOLOGY) AND [excluding] PUBLICATION YEARS (2021) AND DOCUMENT TYPES (ARTICLE OR REVIEW) AND LANGUAGES (ENGLISH) AND WEB OF SCIENCE INDEX (Web of Science Core Collection, SCI), and the timespan of 2011–2020. Bibliometric visualization analysis was conducted with CiteSpace and VOSviewer.

Results: The number of documents grew each year. A total of 2703 articles and 774 reviews were identified from 86 countries/regions, 3524 organizations, 928 journals, and 19,745 authors. China was the most prolific country, with 1881 documents. Contributions from China, the United States, and Germany were the most substantial. The most influential author was Lan Huiyao at The Chinese University of Hong Kong, with 24 publications and 2348 total citations. The bibliometric analysis showed that multilateral cooperation among diverse institutions or investigators was beneficial to high-quality outputs. The keyword “PPAR-gamma” exhibited the strongest burst in recent years, suggesting a potent research focus in the future.

Conclusion: Research on SMAD7 in oncology is continuously developing. Bibliometrics is an interesting tool to present the characteristics of publication years, main authors, and productive organizations in a visualized way. It is worth mentioning that a prospective focus might be the specific mechanism of the interaction of PPAR-gamma with SMAD7 in oncology. In all, bibliometric analysis provides an overview and identifies potential research trends for further studies in this academic field.

F-box only protein 9 and its role in cancer

The F-box proteins (FBP), substrate recognition subunit of the SCF (Skp1-Cullin1-F-box protein complex) E3 ligase, play important roles in the ubiquitylation and subsequent degradation of the target proteins from several cellular processes. Disorders of F-box protein-mediated proteolysis lead to human malignancies. FBP plays an important role in many cellular processes, including cell proliferation, cell cycle, apoptosis, migration, invasion, and metastasis, suggesting that it can be associated with tumorigenesis, cancer development and progression. However, the expression and function of FBXO9 (F-box only protein 9) differ in various types of human cancer. Due to the ability to regulate the stability and activity of oncogenes and tumor-suppressor genes, and the physiological functions of many of the F-box proteins remain subtle, further genetic and mechanistic studies will elaborate and help define FBXO9's role. Targeting F-box protein or F-box protein signaling pathways could be an effective strategy for preventing or treating human cancer. This review is presented to summarize the part of FBXO9 in different types of human cancer and its regulation mechanism, and to pave the way to design FBXO9-targeting anticancer therapies.

An overview of pharmacological activities of baicalin and its aglycone baicalein: New insights into molecular mechanisms and signaling pathways

The flavonoids, baicalin, and its aglycone baicalein possess multi-fold therapeutic properties and are mainly found in the roots of Oroxylum indicum (L.) Kurz and Scutellaria baicalensis Georgi. These flavonoids have been reported to possess various pharmacological properties, including antibacterial, antiviral, anticancer, anticonvulsant, anti-oxidant, hepatoprotective, and neuroprotective effects. The pharmacological properties of baicalin and baicalein are due to their abilities to scavenge reactive oxygen species (ROS) and interaction with various signaling molecules associated with apoptosis, inflammation, autophagy, cell cycle, mitochondrial dynamics, and cytoprotection. In this review, we summarized the molecular mechanisms underlying the chemopreventive and chemotherapeutic applications of baicalin and baicalein in the treatment of cancer and inflammatory diseases. In addition, the preventive effects of baicalin and baicalein on mitochondrial dynamics and functions were highlighted with a particular emphasis on their anti-oxidative and cytoprotective properties. The current review highlights could be useful for future prospective studies to further improve the pharmacological applications of baicalein and baicalin. These studies should define the threshold for optimal drug exposure, dose optimization and focus on therapeutic drug monitoring, objective disease markers, and baicalin/baicalein drug levels.

Vitamin D treatment induces in vitro and ex vivo transcriptomic changes indicating anti-tumor effects

Vitamin D deficiency is associated with risk of several common cancers, including colorectal cancer (CRC). Here we have utilized patient derived epithelial organoids (ex vivo) and CRC cell lines (in vitro) to show that calcitriol (1,25OHD) increased the expression of the CRC tumor suppressor gene, CDH1, at both the transcript and protein level. Whole genome expression analysis demonstrated significant differential expression of a further six genes after 1,25OHD treatment, including genes with established links to carcinogenesis GADD45, EFTUD1 and KIAA1199. Furthermore, gene ontologies relevant to carcinogenesis were enriched by 1,25OHD treatment (e.g., 'regulation of Wnt signaling pathway', 'regulation of cell death'), with common enriched processes across in vitro and ex vivo cultures including 'negative regulation of cell proliferation', 'regulation of cell migration' and 'regulation of cell differentiation'. Our results identify genes and pathways that are modifiable by calcitriol that have links to CRC tumorigenesis. Hence the findings provide potential mechanism to the epidemiological and clinical trial data indicating a causal association between vitamin D and CRC. We suggest there is strong rationale for further well-designed trials of vitamin D supplementation as a novel CRC chemopreventive and chemotherapeutic agent.

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.

Evidence for shared genetic risk factors between lymphangioleiomyomatosis and pulmonary function

Introduction

Lymphangioleiomyomatosis (LAM) is a rare low-grade metastasising disease characterised by cystic lung destruction. The genetic basis of LAM remains incompletely determined, and the disease cell-of-origin is uncertain. We analysed the possibility of a shared genetic basis between LAM and cancer, and LAM and pulmonary function.

Methods

The results of genome-wide association studies of LAM, 17 cancer types and spirometry measures (forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC), FEV1/FVC ratio and peak expiratory flow (PEF)) were analysed for genetic correlations, shared genetic variants and causality. Genomic and transcriptomic data were examined, and immunodetection assays were performed to evaluate pleiotropic genes.

Results

There were no significant overall genetic correlations between LAM and cancer, but LAM correlated negatively with FVC and PEF, and a trend in the same direction was observed for FEV1. 22 shared genetic variants were uncovered between LAM and pulmonary function, while seven shared variants were identified between LAM and cancer. The LAM-pulmonary function shared genetics identified four pleiotropic genes previously recognised in LAM single-cell transcriptomes: ADAM12, BNC2, NR2F2 and SP5. We had previously associated NR2F2 variants with LAM, and we identified its functional partner NR3C1 as another pleotropic factor. NR3C1 expression was confirmed in LAM lung lesions. Another candidate pleiotropic factor, CNTN2, was found more abundant in plasma of LAM patients than that of healthy women.

Conclusions

This study suggests the existence of a common genetic aetiology between LAM and pulmonary function.

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.

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.

Genetic mapping of novel modifiers for ApcMin induced intestinal polyps' development using the genetic architecture power of the collaborative cross mice

BACKGROUND

Familial adenomatous polyposis is an inherited genetic disease, characterized by colorectal polyps. It is caused by inactivating mutations in the Adenomatous polyposis coli (Apc) gene. Mice carrying a nonsense mutation in the Apc gene at R850, which is designated ApcMin/+ (Multiple intestinal neoplasia), develop intestinal adenomas. Several genetic modifier loci of Min (Mom) were previously mapped, but so far, most of the underlying genes have not been identified. To identify novel modifier loci associated with ApcMin/+, we performed quantitative trait loci (QTL) analysis for polyp development using 49 F1 crosses between different Collaborative Cross (CC) lines and C57BL/6 J-ApcMin/+mice. The CC population is a genetic reference panel of recombinant inbred lines, each line independently descended from eight genetically diverse founder strains. C57BL/6 J-ApcMin/+ males were mated with females from 49 CC lines. F1 offspring were terminated at 23 weeks and polyp counts from three sub-regions (SB1-3) of small intestinal and colon were recorded.

RESULTS

The number of polyps in all these sub-regions and colon varied significantly between the different CC lines. At 95% genome-wide significance, we mapped nine novel QTL for variation in polyp number, with distinct QTL associated with each intestinal sub-region. QTL confidence intervals varied in width between 2.63-17.79 Mb. We extracted all genes in the mapped QTL at 90 and 95% CI levels using the BioInfoMiner online platform to extract, significantly enriched pathways and key linker genes, that act as regulatory and orchestrators of the phenotypic landscape associated with the ApcMin/+ mutation.

CONCLUSIONS

Genomic structure of the CC lines has allowed us to identify novel modifiers and confirmed some of the previously mapped modifiers. Key genes involved mainly in metabolic and immunological processes were identified. Future steps in this analysis will be to identify regulatory elements - and possible epistatic effects - located in the mapped QTL.

Meta-Analyzing Multiple Omics Data With Robust Variable Selection

High-throughput omics data are becoming more and more popular in various areas of science. Given that many publicly available datasets address the same questions, researchers have applied meta-analysis to synthesize multiple datasets to achieve more reliable results for model estimation and prediction. Due to the high dimensionality of omics data, it is also desirable to incorporate variable selection into meta-analysis. Existing meta-analyzing variable selection methods are often sensitive to the presence of outliers, and may lead to missed detections of relevant covariates, especially for lasso-type penalties. In this paper, we develop a robust variable selection algorithm for meta-analyzing high-dimensional datasets based on logistic regression. We first search an outlier-free subset from each dataset by borrowing information across the datasets with repeatedly use of the least trimmed squared estimates for the logistic model and together with a hierarchical bi-level variable selection technique. We then refine a reweighting step to further improve the efficiency after obtaining a reliable non-outlier subset. Simulation studies and real data analysis show that our new method can provide more reliable results than the existing meta-analysis methods in the presence of outliers.

Individual and Joint Associations of Genetic Risk and Healthy Lifestyle Score with Colorectal Neoplasms Among Participants of Screening Colonoscopy

Genetic and lifestyle factors contribute to colorectal cancer risk. We investigated their individual and joint associations with various stages of colorectal carcinogenesis. We assessed associations of a polygenic risk score (PRS) and a healthy lifestyle score (HLS) with presence of nonadvanced adenomas and advanced neoplasms among 2,585 participants of screening colonoscopy from Germany. The PRS and HLS individually showed only weak associations with presence of nonadvanced adenomas; stronger associations were observed with advanced neoplasms (ORs, 95% CI, for highest vs. lowest risk tertile: PRS 2.27, 1.78-2.88; HLS 1.96, 1.53-2.51). The PRS was associated with higher odds of advanced neoplasms among carriers of any neoplasms (1.65, 1.23-2.22). Subjects in the highest risk tertile (vs. lowest tertile) of both scores had higher risks for nonadvanced adenomas (1.77, 1.09-2.86), for advanced neoplasms (3.95, 2.53-6.16) and, among carriers of any neoplasms, for advanced versus nonadvanced neoplasms (2.26, 1.31-3.92). Both scores were individually associated with increased risk of nonadvanced adenomas and, much more pronounced, advanced neoplasms. The similarly strong association in relative terms across all levels of genetic risk implies that a healthy lifestyle may be particularly beneficial in those at highest genetic risk, given that the same relative risk reduction in this group would imply a stronger absolute risk reduction. Genetic factors may be of particular relevance for the transition of nonadvanced to advanced adenomas. PREVENTION RELEVANCE: Genetic factors have strong impact on the risk of colorectal neoplasms, which may be reduced by healthy lifestyle. Similarly strong associations in relative terms across all levels of genetic risk imply that a healthy lifestyle may be beneficial due to higher absolute risk reduction in those at highest genetic risk.

Variation rs9929218 and risk of the colorectal Cancer and adenomas: A meta-analysis

BACKGROUNDS

Genome-wide association studies (GWAS) have identified multiple common CRC-related (colorectal cancer) SNPs (single nucleotide polymorphisms) including the Cadherin 1(CDH1) rs9929218 may act by increasing the risk of colorectal cancer, colorectal adenoma, or both. These studies, however, reported inconsistent associations.

METHODS

To derive a more accurate approximation of the connection, we carried out a meta-analysis of 12 published pieces of research including 11,590 controls and 8192 cases. We used odds ratios (ORs) and 95% confidence intervals (CIs) to evaluate the associations' strength.

RESULTS

Meta-analysis implied considerable association between CRC and rs9929218 (OR = 1.21, 95%CI 1.04-1.42 for GG versus AA; OR = 1.22, 95%CI 1.05-1.42 for GG/AG versus AA). In the subgroup analyses, significantly increased risks were found among Europeans.

CONCLUSIONS

In summary, our meta-analysis studies in different populations confirmed that SNP rs9929218 is significantly associated with CRC risk and that this variant may have a greater impact on Europeans.

Polygenic risk score is a predictor of adenomatous polyps at screening colonoscopy

BACKGROUND

Single nucleotide polymorphism (SNP)-based polygenic risk scoring is predictive of colorectal cancer (CRC) risk. However, few studies have investigated the association of genetic risk score (GRS) with detection of adenomatous polyps at screening colonoscopy.

METHODS

We randomly selected 1769 Caucasian subjects who underwent screening colonoscopy from the Genomic Health Initiative (GHI), a biobank of NorthShore University HealthSystem. Outcomes from initial screening colonoscopy were recorded. Twenty-two CRC risk-associated SNPs were obtained from the Affymetrix™ SNP array and used to calculate an odds ratio (OR)-weighted and population-standardized GRS. Subjects with GRS of < 0.5, 0.5-1.5, and > 1.5 were categorized as low, average and elevated risk.

RESULTS

Among 1,769 subjects, 520 (29%) had 1 or more adenomatous polyps. GRS was significantly higher in subjects with adenomatous polyps than those without; mean (95% confidence interval) was 1.02 (1.00-1.05) and 0.97 (0.95-0.99), respectively, p < 0.001. The association remained significant after adjusting for age, gender, body mass index, and family history, p < 0.001. The detection rate of adenomatous polyps was 10.8%, 29.0% and 39.7% in subjects with low, average and elevated GRS, respectively, p-trend < 0.001. Higher GRS was also associated with early age diagnosis of adenomatous polyps, p < 0.001. In contrast, positive family history was not associated with risk and age of adenomatous polyps.

CONCLUSIONS

GRS was significantly associated with adenomatous polyps in subjects undergoing screening colonoscopy. This result may help in stratifying average risk patients and facilitating personalized colonoscopy screening strategies.

Familial Colorectal Cancer and Genetic Susceptibility: Colorectal Risk Variants in First-Degree Relatives of Patients With Colorectal Cancer

INTRODUCTION

Epidemiological studies estimate that having a first-degree relative (FDR) with colorectal cancer (CRC) increases 2-fold to 3-fold the risk of developing the disease. Because FDRs of CRC patients are more likely to co-inherit CRC risk variants, we aimed to evaluate potential differences in genotype distribution of single nucleotide polymorphisms (SNPs) related to CRC risk between FDRs of patients with nonsyndromic CRC (cases) and individuals with no family history of CRC (controls).

METHODS

We designed a case-control study comprising 750 cases and 750 Spanish Caucasian controls matched by sex, age, and histological findings after colonoscopy. Genomic DNA from all participants was genotyped for 88 SNPs associated with CRC risk using the MassArray (Sequenom) platform.

RESULTS

Ten of the 88 SNPs analyzed revealed significant associations (P < 0.05) with a family history of CRC in our population. The most robust associations were found for the rs17094983G>A SNP in the long noncoding RNA LINC01500 (odds ratio = 0.72; 95% confidence interval: 0.58-0.88, log-additive model), and the rs11255841T>A SNP in the long noncoding RNA LINC00709 (odds ratio = 2.04; 95% confidence interval: 1.19-3.51, dominant model). Of interest, the observed associations were in the same direction than those reported for CRC risk.

DISCUSSION

FDRs of CRC patients show significant differences in genotype distribution of SNPs related to CRC risk as compared to individuals with no family history of CRC. Genotyping of CRC risk variants in FDRs of CRC patients may help to identify subjects at risk that would benefit from stricter surveillance and CRC screening programs.

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.

Rhophilin-2 Upregulates Glutamine Synthetase by Stabilizing c-Myc Protein and Confers Resistance to Glutamine Deprivation in Lung Cancer

Introduction

RHPN2, a member of rhophilin family of rho-binding proteins, regulates actin cytoskeleton and vesicular trafficking, and promotes mesenchymal transformation in cancer. We have found that RHPN2 was significantly mutated in lung adenocarcinoma (LUAD). However, the role of RHPN2 in lung cancer is not fully understood.

Methods

In the present study, we investigated the expression of RHPN2 in 125 patients with LUAD by qRT-PCR and correlated its expression with clinical characteristics. The effects of RHPN2 on the proliferation and invasion of lung cancer cells were determined by CCK-8 and in vitro transwell assays, clonogenic assay, and xenograft mouse model. The RhoA pull down assay and Western blotting were performed to elucidate the mechanism of RNPN2 in tumorigenesis of lung cancer.

Results

RHPN2 was overexpressed in tumors from LUAD, and high levels of RHPN2 were associated with poor prognosis of LUAD patients. RHPN2 was required for proliferation and invasion of lung cancer cells. Intriguingly, overexpression of RHPN2 conferred the resistance to glutamine depletion in lung cancer cells. Mechanistic studies revealed that ectopic overexpression of RHPN2 promoted the stability of c-Myc protein via phosphorylation at Ser62 and increased c-Myc target glutamine synthetase (GS). Analysis of GS expression in clinical sample showed that the expression of GS was elevated in tumor cells. Kaplan-Meier analysis revealed that high levels of GS were significantly associated with worse overall survival time of the patients with LUAD.

Conclusions

Taken together, this study suggested that RHPN2 was involved in tumorigenesis of lung cancer via modulating c-Myc stability and the expression of its target GS in lung adenocarcinoma, which links RHPN2 and glutamine metabolism.

Unique roles of rare variants in the genetics of complex diseases in humans

Genome-wide association studies have identified >10,000 genetic variants associated with various phenotypes and diseases. Although the majority are common variants, rare variants with >0.1% of minor allele frequency have been investigated by imputation and using disease-specific custom SNP arrays. Rare variants sequencing analysis mainly revealed have played unique roles in the genetics of complex diseases in humans due to their distinctive features, in contrast to common variants. Unique roles are hypothesis-free evidence for gene causality, a precise target of functional analysis for understanding disease mechanisms, a new favorable target for drug development, and a genetic marker with high disease risk for personalized medicine. As whole-genome sequencing continues to identify more rare variants, the roles associated with rare variants will also increase. However, a better estimation of the functional impact of rare variants across whole genome is needed to enhance their contribution to improvements in human health.

Use of Polygenic Risk Scores to Select Screening Intervals After Negative Findings From Colonoscopy

BACKGROUND & AIMS

Polygenic risk scores (PRSs) could help to define starting ages for colorectal cancer (CRC) screening. However, the role of PRS in determining the length of screening interval after negative findings from colonoscopies is unclear. We aimed to evaluate CRC risk according to PRS and time since last negative colonoscopy.

METHODS

We collected data from 3827 cases and 2641 CRC-free controls in a population-based case-control study in Germany. We constructed a polygenic risk scoring system, based on 90 single-nucleotide polymorphisms, associated with risk of CRC in people of European descent. Participants were classified as having low, medium, or high genetic risk according to tertiles of PRSs among controls. Multiple logistic regression models were used to assess CRC risk according to PRS and time since last negative colonoscopy.

RESULTS

Compared to individuals without colonoscopy in the low PRS category, a 42%-85% lower risk of CRC was observed for individuals who had a negative finding from colonoscopy within 10 years. Beyond 10 years after a negative finding from colonoscopy, significantly lower risk only persisted for the low and medium PRS groups, but not for the high PRS group. Adjusted odds ratios were 0.44 (95% CI, 0.29-0.68), 0.51 (95% CI, 0.34-0.77), and 0.85 (95% CI, 0.58-1.23) in the low, medium, and high PRS group, respectively. Within any time interval, risks were lower for distal than for proximal CRCs.

CONCLUSIONS

Based on findings from a population-based case-control study, the recommended 10-year screening interval for colonoscopy may not need to be shortened among people with high PRSs, but could potentially be prolonged for people with low and medium PRSs. Studies are needed to address personalized time intervals for repeat colonoscopies in average-risk screening cohorts.

The Impact of a Comprehensive Risk Prediction Model for Colorectal Cancer on a Population Screening Program

BACKGROUND

In many countries, population colorectal cancer (CRC) screening is based on age and family history, though more precise risk prediction could better target screening. We examined the impact of a CRC risk prediction model (incorporating age, sex, lifestyle, genomic, and family history factors) to target screening under several feasible screening scenarios.

METHODS

We estimated the model's predicted CRC risk distribution in the Australian population. Predicted CRC risks were categorized into screening recommendations under 3 proposed scenarios to compare with current recommendations: 1) highly tailored, 2) 3 risk categories, and 3) 4 sex-specific risk categories. Under each scenario, for 35- to 74-year-olds, we calculated the number of CRC screens by immunochemical fecal occult blood testing (iFOBT) and colonoscopy and the proportion of predicted CRCs over 10 years in each screening group.

RESULTS

Currently, 1.1% of 35- to 74-year-olds are recommended screening colonoscopy and 56.2% iFOBT, and 5.7% and 83.2% of CRCs over 10 years were predicted to occur in these groups, respectively. For the scenarios, 1) colonoscopy was recommended to 8.1% and iFOBT to 37.5%, with 36.1% and 50.1% of CRCs in each group; 2) colonoscopy was recommended to 2.4% and iFOBT to 56.0%, with 13.2% and 76.9% of cancers in each group; and 3) colonoscopy was recommended to 5.0% and iFOBT to 54.2%, with 24.5% and 66.5% of cancers in each group.

CONCLUSIONS

A highly tailored CRC screening scenario results in many fewer screens but more cancers in those unscreened. Category-based scenarios may provide a good balance between number of screens and cancers detected and are simpler to implement.

Estimation of Absolute Risk of Colorectal Cancer Based on Healthy Lifestyle, Genetic Risk, and Colonoscopy Status in a Population-Based Study

BACKGROUND & AIMS

Estimates of absolute risk of colorectal cancer (CRC) are needed to facilitate communication and better inform the public about the potentials and limits of cancer prevention.

METHODS

Using data from a large population-based case-control study in Germany (Darmkrebs: Chancen der Verhütung durch Screening [DACHS] study, which began in 2003) and population registry data, we calculated 30-year absolute risk estimates for development of CRC based on a healthy lifestyle score (derived from 5 modifiable lifestyle factors: smoking, alcohol consumption, diet, physical activity, and body fatness), a polygenic risk score (based on 90 single-nucleotide polymorphisms), and colonoscopy history.

RESULTS

We analyzed data from 4220 patients with CRC and 3338 individuals without CRC. Adherence to a healthy lifestyle and colonoscopy in the preceding 10 years were associated with a reduced relative risk of CRC in men and women. We observed a higher CRC risk in participants with high or intermediate genetic risk scores. For 50-year-old men and women without a colonoscopy, the absolute risk of CRC varied according to the polygenic risk score and the healthy lifestyle score (men, 3.5%-13.4%; women, 2.5%-10.6%). For 50-year-old men and women with a colonoscopy, the absolute risk of developing CRC was much lower but still varied according to the polygenic risk score and the healthy lifestyle score (men, 1.2%-4.8%; women, 0.9%-4.2%). Among all risk factor profiles, the 30-year absolute risk estimates consistently decreased with adherence to a healthy lifestyle.

CONCLUSIONS

In a population-based study, we found that a colonoscopy can drastically reduce the absolute risk of CRC and that the genetically predetermined risk of CRC can be further reduced by adherence to a healthy lifestyle. Our results show the magnitude of CRC prevention possible through colonoscopy and lifestyle at a predefined genetic risk. This observational study has been registered in the German Clinical Trials Register (DRKS00011793), which is a primary registry in the World Health Organization Registry Network.

A Systematic Analysis of Interactions between Environmental Risk Factors and Genetic Variation in Susceptibility to Colorectal Cancer

BACKGROUND

The underlying etiology of colorectal cancer includes both genetic variation and environmental exposures. The main aim of this study was to search for interaction effects between well-established environmental risk factors and published common genetic variants exerting main effects on colorectal cancer risk.

METHODS

We used a two-phase approach: (i) discovery phase (2,652 incident colorectal cancer cases and 10,608 controls from UK Biobank) and (ii) validation phase (1,656 cases and 2,497 controls from the Study of Colorectal Cancer in Scotland). Interactions with nominal P < 0.05 in phase I were taken forward for validation in phase II. Furthermore, we constructed a weighted genetic risk score (GRS) of colorectal cancer risk for each individual and studied interactions between the GRS and the environmental risk factors.

RESULTS

Seventy of the 1,500 tested interactions were nominally significant in phase I. After testing these 70 interactions in phase II, an interaction between rs11903757 (2q32.3) and body mass index (BMI) was nominally significant (P = 0.02) with the same direction of effect. The rs11903757*BMI interaction was also significant (ratio of ORs = 1.26; 95% confidence interval, 1.10-1.44; P interaction = 6.03 × 10-4; P heterogeneity = 0.63) in a meta-analysis combining results from both phases. No interactions were significant in phase II after accounting for multiple testing. No interactions involving the GRS were found with statistical significance.

CONCLUSIONS

Limited evidence of gene-environment interactions in colorectal cancer risk was observed. There are potential modifications of the rs11903757 effect by BMI on colorectal cancer risk.

IMPACT

Our findings might contribute to identifying subpopulations with different susceptibility to the effect of BMI on colorectal cancer risk.

Detecting Shared Genetic Architecture Among Multiple Phenotypes by Hierarchical Clustering of Gene-Level Association Statistics

Emerging large-scale biobanks pairing genotype data with phenotype data present new opportunities to prioritize shared genetic associations across multiple phenotypes for molecular validation. Past research, by our group and others, has shown gene-level tests of association produce biologically interpretable characterization of the genetic architecture of a given phenotype. Here, we present a new method, Ward clustering to identify Internal Node branch length outliers using Gene Scores (WINGS), for identifying shared genetic architecture among multiple phenotypes. The objective of WINGS is to identify groups of phenotypes, or "clusters," sharing a core set of genes enriched for mutations in cases. We validate WINGS using extensive simulation studies and then combine gene-level association tests with WINGS to identify shared genetic architecture among 81 case-control and seven quantitative phenotypes in 349,468 European-ancestry individuals from the UK Biobank. We identify eight prioritized phenotype clusters and recover multiple published gene-level associations within prioritized clusters.

Genome-wide association meta-analyses combining multiple risk phenotypes provide insights into the genetic architecture of cutaneous melanoma susceptibility

Most genetic susceptibility to cutaneous melanoma remains to be discovered. Meta-analysis genome-wide association study (GWAS) of 36,760 cases of melanoma (67% newly genotyped) and 375,188 controls identified 54 significant (P < 5 × 10-8) loci with 68 independent single nucleotide polymorphisms. Analysis of risk estimates across geographical regions and host factors suggests the acral melanoma subtype is uniquely unrelated to pigmentation. Combining this meta-analysis with GWAS of nevus count and hair color, and transcriptome association approaches, uncovered 31 potential secondary loci for a total of 85 cutaneous melanoma susceptibility loci. These findings provide insights into cutaneous melanoma genetic architecture, reinforcing the importance of nevogenesis, pigmentation and telomere maintenance, together with identifying potential new pathways for cutaneous melanoma pathogenesis.

Risk SNP-Mediated Enhancer-Promoter Interaction Drives Colorectal Cancer through Both FADS2 and AP002754.2

Although genome-wide association studies (GWAS) have identified more than 100 colorectal cancer risk loci, most of the biological mechanisms associated with these loci remain unclear. Here we first performed a comprehensive expression quantitative trait loci analysis in colorectal cancer tissues adjusted for multiple confounders to test the determinants of germline variants in established GWAS susceptibility loci on mRNA and long noncoding RNA (lncRNA) expression. Combining integrative functional genomic/epigenomic analyses and a large-scale population study consisting of 6,024 cases and 10,022 controls, we then prioritized rs174575 with a C>G change as a potential causal candidate for colorectal cancer at 11q12.2, as its G allele was associated with an increased risk of colorectal cancer (OR = 1.26; 95% confidence interval = 1.17-1.36; P = 2.57 × 10^-9). rs174575 acted as an allele-specific enhancer to distally facilitate expression of both FADS2 and lncRNA AP002754.2 via long-range enhancer-promoter interaction loops, which were mediated by E2F1. AP002754.2 further activated a transcriptional activator that upregulated FADS2 expression. FADS2, in turn, was overexpressed in colorectal cancer tumor tissues and functioned as a potential oncogene that facilitated colorectal cancer cell proliferation and xenograft growth in vitro and in vivo by increasing the metabolism of PGE2, an oncogenic molecule involved in colorectal cancer tumorigenesis. Our findings represent a novel mechanism by which a noncoding variant can facilitate long-range genome interactions to modulate the expression of multiple genes including not only mRNA, but also lncRNA, which provides new insights into the understanding of colorectal cancer etiology. SIGNIFICANCE: This study provides an oncogenic regulatory circuit among several oncogenes including E2F1, FADS2, and AP002754.2 underlying the association of rs174575 with colorectal cancer risk, which is driven by long-range enhancer-promoter interaction loops. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/9/1804/F1.large.jpg.

A New Polygenic Model for Nonfamilial Colorectal Cancer Inheritance Based on the Genetic Architecture of the Azoxymethane-Induced Mouse Model

The azoxymethane model of colorectal cancer (CRC) was used to gain insights into the genetic heterogeneity of nonfamilial CRC. We observed significant differences in susceptibility parameters across 40 mouse inbred strains, with 6 new and 18 of 24 previously identified mouse CRC modifier alleles detected using genome-wide association analysis. Tumor incidence varied in F1 as well as intercrosses and backcrosses between resistant and susceptible strains. Analysis of inheritance patterns indicates that resistance to CRC development is inherited as a dominant characteristic genome-wide, and that susceptibility appears to occur in individuals lacking a large-effect, or sufficient numbers of small-effect, polygenic resistance alleles. Our results suggest a new polygenic model for inheritance of nonfamilial CRC, and that genetic studies in humans aimed at identifying individuals with elevated susceptibility should be pursued through the lens of absence of dominant resistance alleles rather than for the presence of susceptibility alleles.

Identification of Novel Loci and New Risk Variant in Known Loci for Colorectal Cancer Risk in East Asians

BACKGROUND

Risk variants identified so far for colorectal cancer explain only a small proportion of familial risk of this cancer, particularly in Asians.

METHODS

We performed a genome-wide association study (GWAS) of colorectal cancer in East Asians, including 23,572 colorectal cancer cases and 48,700 controls. To identify novel risk loci, we selected 60 promising risk variants for replication using data from 58,131 colorectal cancer cases and 67,347 controls of European descent. To identify additional risk variants in known colorectal cancer loci, we performed conditional analyses in East Asians.

RESULTS

An indel variant, rs67052019 at 1p13.3, was found to be associated with colorectal cancer risk at P = 3.9 × 10-8 in Asians (OR per allele deletion = 1.13, 95% confidence interval = 1.08-1.18). This association was replicated in European descendants using a variant (rs2938616) in complete linkage disequilibrium with rs67052019 (P = 7.7 × 10-3). Of the remaining 59 variants, 12 showed an association at P < 0.05 in the European-ancestry study, including rs11108175 and rs9634162 at P < 5 × 10-8 and two variants with an association near the genome-wide significance level (rs60911071, P = 5.8 × 10-8; rs62558833, P = 7.5 × 10-8) in the combined analyses of Asian- and European-ancestry data. In addition, using data from East Asians, we identified 13 new risk variants at 11 loci reported from previous GWAS.

CONCLUSIONS

In this large GWAS, we identified three novel risk loci and two highly suggestive loci for colorectal cancer risk and provided evidence for potential roles of multiple genes and pathways in the etiology of colorectal cancer. In addition, we showed that additional risk variants exist in many colorectal cancer risk loci identified previously.

IMPACT

Our study provides novel data to improve the understanding of the genetic basis for colorectal cancer risk.

Colorectal cancer susceptibility variants and risk of conventional adenomas and serrated polyps: results from three cohort studies

BACKGROUND

Increasing evidence suggests that conventional adenomas (CAs) and serrated polyps (SPs) represent two distinct groups of precursor lesions for colorectal cancer (CRC). The influence of common genetic variants on risk of CAs and SPs remain largely unknown.

METHODS

Among 27 426 participants within three prospective cohort studies, we created a weighted genetic risk score (GRS) based on 40 CRC-related single nucleotide polymorphisms (SNPs) identified in previous genome-wide association studies; and we examined the association of GRS (per one standard deviation increment) with risk of CAs, SPs and synchronous CAs and SPs, by multivariable logistic regression. We also analysed individual variants in the secondary analysis.

RESULTS

During 18-20 years of follow-up, we documented 2952 CAs, 1585 SPs and 794 synchronous CAs and SPs. Higher GRS was associated with increased risk of CAs [odds ratio (OR) = 1.17, 95% confidence interval (CI): 1.12-1.21] and SPs (OR = 1.09, 95% CI: 1.03-1.14), with a stronger association for CAs than SPs (Pheterogeneity=0.01). An even stronger association was found for patients with synchronous CAs and SPs (OR = 1.32), advanced CAs (OR = 1.22) and multiple CAs (OR = 1.25). Different sets of variants were associated with CAs and SPs, with a Spearman correlation coefficient of 0.02 between the ORs associating the 40 SNPs with the two lesions. After correcting for multiple testing, three variants were associated with CAs (rs3802842, rs6983267 and rs7136702) and two with SPs (rs16892766 and rs4779584).

CONCLUSIONS

Common genetic variants play a potential role in the conventional and serrated pathways of CRC. Different sets of variants are identified for the two pathways, further supporting the aetiological heterogeneity of CRC.

Glycosylation of immunoglobulin G is regulated by a large network of genes pleiotropic with inflammatory diseases

Effector functions of immunoglobulin G (IgG) are regulated by the composition of a glycan moiety, thus affecting activity of the immune system. Aberrant glycosylation of IgG has been observed in many diseases, but little is understood about the underlying mechanisms. We performed a genome-wide association study of IgG N-glycosylation (N = 8090) and, using a data-driven network approach, suggested how associated loci form a functional network. We confirmed in vitro that knockdown of IKZF1 decreases the expression of fucosyltransferase FUT8, resulting in increased levels of fucosylated glycans, and suggest that RUNX1 and RUNX3, together with SMARCB1, regulate expression of glycosyltransferase MGAT3. We also show that variants affecting the expression of genes involved in the regulation of glycoenzymes colocalize with variants affecting risk for inflammatory diseases. This study provides new evidence that variation in key transcription factors coupled with regulatory variation in glycogenes modifies IgG glycosylation and has influence on inflammatory diseases.

Establishing a valid approach for estimating familial risk of cancer explained by common genetic variants

We critically examined existing approaches for the estimation of the excess familial risk of cancer that can be attributed to identified common genetic risk variants and propose an alternative, more straightforward approach for calculating this proportion using well-established epidemiological methodology. We applied the underlying equations of the traditional approaches and the new epidemiological approach for colorectal cancer (CRC) in a large population-based case-control study in Germany with 4,447 cases and 3,480 controls, who were recruited from 2003 to 2016 and for whom interview, medical and genomic data were available. Having a family history of CRC (FH) was associated with a 1.77-fold risk increase in our study population (95% CI 1.52-2.07). Traditional approaches yielded estimates of the FH-associated risk explained by 97 common genetics variants from 9.6% to 23.1%, depending on various assumptions. Our alternative approach resulted in smaller and more consistent estimates of this proportion, ranging from 5.4% to 14.3%. Commonly employed methods may lead to strongly divergent and possibly exaggerated estimates of excess familial risk of cancer explained by associated known common genetic variants. Our results suggest that familial risk and risk associated with known common genetic variants might reflect two complementary major sources of risk.

Evaluation of gene-environment interactions for colorectal cancer susceptibility loci using case-only and case-control designs

BACKGROUND

Genome-wide association studies (GWAS) have identified more than 40 colorectal cancer susceptibility loci, but only a small fraction of heritability was explained. To account for missing heritability, we investigated gene-environment interactions (G × Es) between GWAS-identified single-nucleotide polymorphisms (SNPs) and established risk or protective factors for colorectal cancer using both case-only and case-control study designs.

METHODS

Data on 703 colorectal cancer cases and 1406 healthy controls from the National Cancer Center in Korea were used. We tested interactions between 31 GWAS-identified SNPs and 13 established risk or protective factors for colorectal cancer (family history, body mass index, history of colorectal polyps, inflammatory bowel disease, and diabetes mellitus, alcohol drinking, smoking, regular exercise, regular aspirin use, postmenopausal hormone replace therapy, red meat and processed meat intake, and dairy consumption). Logistic regression models were used to assess G × Es for colorectal cancer risk.

RESULTS

The SNP rs4444235 at 14q22.2 interacted with regular exercise in colorectal cancer (p_case-only = 2.4 × 10^- 3, p_case-control = 1.5 × 10^- 3). The risk allele (C) of rs4444235 increased the risk of colorectal cancer in regularly exercising individuals (OR = 1.47, 95% CI = 1.02-2.10) but decreased the risk in non-exercising individuals (OR = 0.76, 95% CI = 0.62-0.94). Furthermore, the G × E between the SNP rs2423279 at 20p12.3 and regular aspirin use was statistically significant (p_case-only = 7.7 × 10^- 3, p_case-control = 1.6 × 10^- 3). The additive effect of the risk allele (T) of rs2423279 on colorectal cancer risk was increased among regular aspirin users (OR = 4.62, 95% CI = 1.97-10.80).

CONCLUSION

Our results suggest that SNP rs4444235 at 14q22.2 and SNP rs2423279 at 20p12.3 may interact with regular exercise and aspirin use in colorectal carcinogenesis.

Single nucleotide polymorphisms associated with susceptibility for development of colorectal cancer: Case-control study in a Basque population

Given the significant population diversity in genetic variation, we aimed to investigate whether single nucleotide polymorphisms (SNPs) previously identified in studies of colorectal cancer (CRC) susceptibility were also relevant to the population of the Basque Country (North of Spain). We genotyped 230 CRC cases and 230 healthy controls for 48 previously reported CRC-susceptibility SNPs. Only the rs6687758 in DUPS10 exhibited a statistically significant association with CRC risk based on the crude analysis. The rs6687758 AG genotype conferred about 2.13-fold increased risk for CRC compared to the AA genotype. Moreover, we found significant associations in cases between smoking status, physical activity, and the rs6687758 SNP. The results of a Genetic Risk Score (GRS) showed that the risk alleles were more frequent in cases than controls and the score was associated with CRC in crude analysis. In conclusion, we have confirmed a CRC susceptibility locus and the existence of associations between modifiable factors and the rs6687758 SNP; moreover, the GRS was associated with CRC. However, further experimental validations are needed to establish the role of this SNP, the function of the gene identified, as well as the contribution of the interaction between environmental factors and this locusto the risk of CRC.

Cost-utility analysis of genetic polymorphism universal screening in colorectal cancer prevention by detection of high-risk individuals

BACKGROUND

In the past 15 years numerous studies have been published on the involvement of low-penetrance susceptibility genes on the risk for developing colorectal cancer (CRC).

AIM

To perform an economic analysis of blood genetic testing in CRC screening in a population-based nationwide setting using polymorphisms in prostaglandin E₂ pathway genes as proof of concept.

METHODS

A cost-utility analysis was performed from a societal perspective in Portugal comparing two strategies: blood genetic testing by the age of 40 versus no genetic screening under different assumptions of the cost of genetic testing (€10 and €30) and expected risk (1.5 to 5-fold). The adopted threshold was set at €44,870 (USD 50,000). The primary outcome was the incremental cost-effectiveness ratio (ICER) for a base case scenario.

RESULTS

Polymorphism genotyping provided cost-utility only under the assumption of a 5-fold increased risk in the general population, providing ICERs of €44,356 and €30,389 for €30 and €10 tests, respectively.

CONCLUSION

Blood genetic screening for colorectal cancer has cost-utility only under specific assumptions of increased CRC risk and conservative cost estimates. Future studies should focus on defining genetic profiles because single-gene approaches are very unlikely to be cost-effective considering their modest predictive value.

Head-to-Head Comparison of Family History of Colorectal Cancer and a Genetic Risk Score for Colorectal Cancer Risk Stratification

OBJECTIVES

Family history (FH) is associated with increased risk of colorectal cancer (CRC). We aimed to examine the potential for CRC risk stratification by known common genetic variants beyond FH in a large population-based case-control study from Germany.

METHODS

Four thousand four hundred forty-seven cases and 3,480 controls recruited in 2003-2016 were included for whom comprehensive interview, medical, and genomic data were available. Associations with CRC risk were estimated from multiple logistic regression models for FH and a genetic risk score (GRS) based on 90 previously identified risk variants.

RESULTS

CRC in a first-degree relative was associated with a 1.71-fold (95% confidence interval 1.47-2.00) increase in CRC risk. A higher risk increase (odds ratio 2.06, 95% confidence interval 1.78-2.39) was estimated for the GRS when it was dichotomized at a cutoff yielding the same positivity rate as FH among controls. Furthermore, the GRS provides substantial additional risk stratification in both people with and especially without FH. Equal or even slightly higher risks were observed for participants without FH with a GRS in the upper 20% compared with participants with FH with a GRS below median. The observed patterns were confirmed in a replication study.

DISCUSSION

In contrast to common perception, known genetic variants do not primarily reflect some minor share of the familial excess risk of CRC, but rather reflect a substantial share of risk independent of FH.

A 16q22.1 variant confers susceptibility to colorectal cancer as a distal regulator of ZFP90

Genome-wide association studies (GWASs) implicate 16q22.1 locus in risk for colorectal cancer (CRC). However, the underlying oncogenic mechanisms remain unknown. Here, through comprehensive filtration, we prioritized rs7198799, a common SNP in the second intron of the CDH1, as the putative causal variant. In addition, we found an association of CRC-risk allele C of rs7198799 with elevated transcript level of biological plausible candidate gene ZFP90 via expression quantitative trait loci analysis. Mechanistically, causal variant rs7198799 resides in an enhancer element and remotely regulate ZFP90 expression by targeting the transcription factor NFATC2. Remarkably, CRISPR/Cas9-guided single-nucleotide editing demonstrated the direct effect of rs7198799 on ZFP90 expression and CRC cellular malignant phenotype. Furthermore, ZFP90 affects several oncogenic pathways, including BMP4, and promotes carcinogenesis in patients and in animal models with ZFP90 specific genetic manipulation. Taken together, these findings reveal a risk SNP-mediated long-range regulation on the NFATC2-ZFP90-BMP4 pathway underlying the initiation of CRC.

Ability of known susceptibility SNPs to predict colorectal cancer risk for persons with and without a family history

Before SNP-based risk can be incorporated in colorectal cancer (CRC) screening, the ability of these SNPs to estimate CRC risk for persons with and without a family history of CRC, and the screening implications need to be determined. We estimated the association with CRC of a 45 SNP-based risk using 1181 cases and 999 controls, and its correlation with CRC risk predicted from detailed family history. We estimated the predicted change in the distribution across predefined risk categories, and implications for recommended screening commencement age, from adding SNP-based risk to family history. The inter-quintile risk ratio for colorectal cancer risk of the SNP-based risk was 3.28 (95% CI 2.54-4.22). SNP-based and family history-based risks were not correlated (r = 0.02). For persons with no first-degree relatives with CRC, screening could commence 4 years earlier for women (5 years for men) in the highest quintile of SNP-based risk. For persons with two first-degree relatives with CRC, screening could commence 16 years earlier for men and women in the highest quintile, and 7 years earlier for the lowest quintile. This 45 SNP panel in conjunction with family history, can identify people who could benefit from earlier screening. Risk reclassification by 45 SNPs could inform targeted screening for CRC prevention, particularly in clinical genetics settings when mutations in high-risk genes cannot be identified. Yet to be determined is cost-effectiveness, resources requirements, community, patient and clinician acceptance, and feasibility with potentially ethical, legal and insurance implications.

Mendelian randomisation: A powerful and inexpensive method for identifying and excluding non-genetic risk factors for colorectal cancer

Colorectal cancer (CRC) is the third most common cancer in economically developed countries and a major cause of cancer-related mortality. The importance of lifestyle and diet as major determinants of CRC risk is suggested by differences in CRC incidence between countries and in migration studies. Previous observational epidemiological studies have identified associations between modifiable environmental risk factors and CRC, but these studies can be susceptible to reverse causation and confounding, and their results can therefore conflict. Mendelian randomisation (MR) analysis represents an approach complementary to conventional observational studies examining associations between exposures and disease. The MR strategy employs allelic variants as instrumental variables (IVs), which act as proxies for non-genetic exposures. These allelic variants are randomly assigned during meiosis and can therefore inform on life-long exposure, whilst not being subject to reverse causation. In previous studies MR frameworks have associated several modifiable factors with CRC risk, including adiposity, hyperlipidaemia, fatty acid profile and alcohol consumption. In this review we detail the use of MR to investigate and discover CRC risk factors, and its future applications.

HSPB1 rs2070804 polymorphism is associated with the depth of primary tumor

BACKGROUND

Colorectal cancer (CRC) is the third most common cancer in the world. Genome-wide association studies are a powerful method to analyze the status of single-nucleotide polymorphisms (SNPs) in specific genes. Heat shock proteins (HSPs) were found to be involved in the cancer progression and chemoresistance. However, there is still no further study about polymorphisms of HSP beta-1 (HSPB1) in colorectal cancer. We proposed the SNP of HSPB1 may be correlated with the progression and metastasis in colon cancer.

METHODS

We recruited 379 colorectal cancer patients and categorized as four stages following the UICC TNM system. Then, we selected tagging SNPs of HSPB1 by 10% minimum allelic frequency in Han Chinese population from the HapMap database and analyze with the Chi-square test.

RESULTS

We demonstrated the association of HSPB1 genetic polymorphisms rs2070804 with tumor depth with colorectal cancer. But, there is a lack of association between HSPB1 genetic polymorphisms and colorectal cancer invasion, recurrence or metastasis.

CONCLUSIONS

The polymorphisms of HSPB1 seemed to change the tumor behavior of colorectal cancer. HSPB1 rs2070804 polymorphism is associated with the depth of the primary tumor. But, there is no further correlation with other to the clinical parameters such as cancer invasiveness, local recurrence, or distant metastasis.

Systematic Functional Interrogation of Genes in GWAS Loci Identified ATF1 as a Key Driver in Colorectal Cancer Modulated by a Promoter-Enhancer Interaction

Genome-wide association studies (GWASs) have identified approximately 100 colorectal cancer (CRC) risk loci. However, the causal genes in these loci have not been systematically interrogated. We conducted a high-throughput RNA-interference functional screen to identify the genes essential for proliferation in the CRC risk loci of Asian populations. We found that ATF1, located in the 12q13.12 region, functions as an oncogene that facilitates cell proliferation; ATF1 has the most significant effect of the identified genes and promotes CRC xenograft growth by affecting cell apoptosis. Next, by integrating a fine-mapping analysis, a two-stage affected-control study consisting of 6,213 affected individuals and 10,388 controls, and multipronged experiments, we elucidated that two risk variants, dbSNP: rs61926301 and dbSNP: rs7959129, that located in the ATF1 promoter and first intron, respectively, facilitate a promoter-enhancer interaction, mediated by the synergy of SP1 and GATA3, to upregulate ATF1 expression, thus synergistically predisposing to CRC risk (OR = 1.77, 95% CI = 1.42-2.21, p = 3.16 × 10^-7; P_{multiplicative-interaction} = 1.20 × 10^-22; P_{additive-interaction} = 6.50 × 10^-3). Finally, we performed RNA-seq and ChIP-seq assays in CRC cells treated with ATF1 overexpression in order to dissect the target programs of ATF1. Results showed that ATF1 activates a subset of genes, including BRAF, NRAS, MYC, BIRC2, DAAM1, MAML2, STAT1, ID1, and NKD2, related to apoptosis, Wnt, TGF-β, and MAPK pathways, and these effects could cooperatively increase the risk of CRC. These findings reveal the clinical potential of ATF1 in CRC development and illuminate a promoter-enhancer interaction module between the ATF1 regulatory elements dbSNP: rs61926301 and dbSNP: rs7959129, and they bring us closer to understanding the molecular drivers of cancer.

A germline mutation in Rab43 gene identified from a cancer family predisposes to a hereditary liver-colon cancer syndrome

Background

Hereditary cancer syndromes have inherited germline mutations which predispose to benign and malignant tumors. Understanding of the molecular causes in hereditary cancer syndromes has advanced cancer treatment and prevention. However, the causal genes of many hereditary cancer syndromes remain unknown due to their rare frequency of mutation.

Methods

A large Chinese family with a history of hereditary liver-colon cancer syndrome was studied. The genomic DNA was extracted from the blood samples of involved family members, whole-exome sequencing was performed to identify genetic variants. Functional validation of a candidate variant was carried out using gene expression, gene knockout and immunohistochemistry.

Results

The whole-exome of the proband diagnosed with colon cancer was sequenced in comparison with his mother. A total of 13 SNVs and 16 InDels were identified. Among these variants, we focused on a mutation of Rab43 gene, a GTPase family member involving in protein trafficking, for further validation. Sanger DNA sequencing confirmed a mutation (c: 128810106C > T, p: A158T) occurred in one allele of Rab43 gene from the proband, that heterozygous mutation also was verified in the genome of the proband’s deceased father with liver cancer, but not in his healthy mother and sister. Ectopic expression of the Rab43 A158T mutant in Huh7 cells led to more enhanced cell growth, proliferation and migration compared to the expression of wild type Rab43. Conversely, knockout of Rab43 in HepG2 cells resulted in slow cell growth and multiple nuclei formation and impaired activation of Akt. Finally, a positive correlation between the expression levels of Rab43 protein and cancer development in that family was confirmed.

Conclusions

A germline mutation of Rab43 gene is identified to be associated with the onset of a familial liver-colon cancer syndrome. Our finding points to a potential role of protein trafficking in the tumorigenesis of the familial cancer syndrome, and helps the genetic counseling to the affected family members.

Electronic supplementary material

The online version of this article (10.1186/s12885-019-5845-4) contains supplementary material, which is available to authorized users.