Identification of Genetic Susceptibility Loci for Colorectal Tumors in a Genome-Wide Meta-analysis

The Role of HOX Transcription Factors in Cancer Predisposition and Progression

Homeobox (HOX) transcription factors, encoded by a subset of homeodomain superfamily genes, play pivotal roles in many aspects of cellular physiology, embryonic development, and tissue homeostasis. Findings over the past decade have revealed that mutations in HOX genes can lead to increased cancer predisposition, and HOX genes might mediate the effect of many other cancer susceptibility factors by recognizing or executing altered genetic information. Remarkably, several lines of evidence highlight the interplays between HOX transcription factors and cancer risk loci discovered by genome-wide association studies, thereby gaining molecular and biological insight into cancer etiology. In addition, deregulated HOX gene expression impacts various aspects of cancer progression, including tumor angiogenesis, cell autophagy, proliferation, apoptosis, tumor cell migration, and metabolism. In this review, we will discuss the fundamental roles of HOX genes in cancer susceptibility and progression, highlighting multiple molecular mechanisms of HOX involved gene misregulation, as well as their potential implications in clinical practice.

Systematic evaluation of cancer-specific genetic risk score for 11 types of cancer in The Cancer Genome Atlas and Electronic Medical Records and Genomics cohorts

BACKGROUND

Genetic risk score (GRS) is an odds ratio (OR)-weighted and population-standardized method for measuring cumulative effect of multiple risk-associated single nucleotide polymorphisms (SNPs). We hypothesize that GRS is a valid tool for risk assessment of most common cancers.

METHODS

Utilizing genotype and phenotype data from The Cancer Genome Atlas (TCGA) and Electronic Medical Records and Genomics (eMERGE), we tested 11 cancer-specific GRSs (bladder, breast, colorectal, glioma, lung, melanoma, ovarian, pancreatic, prostate, renal, and thyroid cancer) for association with the respective cancer type. Cancer-specific GRSs were calculated, for the first time in these cohorts, based on previously published risk-associated SNPs using the Caucasian subjects in these two cohorts.

RESULTS

Mean cancer-specific GRS in the population controls of eMERGE approximated the expected value of 1.00 (between 0.98 and 1.02) for all 11 types of cancer. Mean cancer-specific GRS was consistently higher in respective cancer patients than controls for all 11 types of cancer (P < 0.05). When subjects were categorized into low-, average-, and high-risk groups based on cancer-specific GRS (<0.5, 0.5-1.5, and >1.5, respectively), significant dose-response associations of higher cancer-specific GRS with higher OR of respective type of cancer were found for nine types of cancer (P_-trend  < 0.05). More than 64% subjects in the population controls of eMERGE can be classified as high risk for at least one type of these cancers.

CONCLUSION

Validity of GRS for predicting cancer risk is demonstrated for most types of cancer. If confirmed in larger studies, cancer-specific GRS may have the potential for developing personalized cancer screening strategy.

Genetic variant predictors of gene expression provide new insight into risk of colorectal cancer

Genome-wide association studies have reported 56 independently associated colorectal cancer (CRC) risk variants, most of which are non-coding and believed to exert their effects by modulating gene expression. The computational method PrediXcan uses cis-regulatory variant predictors to impute expression and perform gene-level association tests in GWAS without directly measured transcriptomes. In this study, we used reference datasets from colon (n = 169) and whole blood (n = 922) transcriptomes to test CRC association with genetically determined expression levels in a genome-wide analysis of 12,186 cases and 14,718 controls. Three novel associations were discovered from colon transverse models at FDR ≤ 0.2 and further evaluated in an independent replication including 32,825 cases and 39,933 controls. After adjusting for multiple comparisons, we found statistically significant associations using colon transcriptome models with TRIM4 (discovery P = 2.2 × 10- 4, replication P = 0.01), and PYGL (discovery P = 2.3 × 10- 4, replication P = 6.7 × 10- 4). Interestingly, both genes encode proteins that influence redox homeostasis and are related to cellular metabolic reprogramming in tumors, implicating a novel CRC pathway linked to cell growth and proliferation. Defining CRC risk regions as one megabase up- and downstream of one of the 56 independent risk variants, we defined 44 non-overlapping CRC-risk regions. Among these risk regions, we identified genes associated with CRC (P < 0.05) in 34/44 CRC-risk regions. Importantly, CRC association was found for two genes in the previously reported 2q25 locus, CXCR1 and CXCR2, which are potential cancer therapeutic targets. These findings provide strong candidate genes to prioritize for subsequent laboratory follow-up of GWAS loci. This study is the first to implement PrediXcan in a large colorectal cancer study and findings highlight the utility of integrating transcriptome data in GWAS for discovery of, and biological insight into, risk loci.

Meeting report from the joint IARC-NCI international cancer seminar series: a focus on colorectal cancer

Despite significant progress in our understanding of the etiology, biology and genetics of colorectal cancer, as well as important clinical advances, it remains the third most frequently diagnosed cancer worldwide and is the second leading cause of cancer death. Based on demographic projections, the global burden of colorectal cancer would be expected to rise by 72% from 1.8 million new cases in 2018 to over 3 million in 2040 with substantial increases anticipated in low- and middle-income countries. In this meeting report, we summarize the content of a joint workshop led by the National Cancer Institute and the International Agency for Research on Cancer, which was held to summarize the important achievements that have been made in our understanding of colorectal cancer etiology, genetics, early detection and treatment and to identify key research questions that remain to be addressed.

Combinations of single nucleotide polymorphisms identified in genome-wide association studies determine risk for colorectal cancer

Genome-wide association studies (GWASs) have identified single nucleotide polymorphisms (SNPs) associated with colorectal cancer (CRC) risk, but whether these SNPs have additive effects on the risk of CRC remains unclear. We performed a systematic analysis of GWAS-identified SNPs using GWAS datasets from China (2,248 patients and 3,173 controls) and Europe (4,461 patients and 4,140 controls). We analyzed 58 independent variants from DNA samples from Chinese populations and found 19 SNPs that were significantly associated with CRC risk. We identified two genetic risk scores (GRSs) based on 58 and 19 SNPs, which were significantly associated with an increased risk of CRC. A decision curve analysis showed higher predictive power for the 58 SNPs. Using all the 58 SNPs to assess 5-year absolute risk (AR), we found that, at a cutoff of 0.4% (two times the median AR) and 0.6% (three times the median AR), approximately 32.76 and 16.45% of Chinese individuals were grouped as high risk for developing CRC, respectively. Risk stratification analysis further indicated that the population in the top 30% risk group accounted for 46.71% of the CRC cases. In addition, the 58 SNPs could explain approximately 1.13% of the phenotypic variance in Chinese populations. Similar findings were found in European populations. Combinations of SNPs identified in GWASs may therefore be useful for identifying individuals at high risk for CRC.

Association of an miR-502-binding site polymorphism in the 3'-untranslated region of SET8 with colorectal cancer

The histone methyltransferase SET8 is regulated by microRNA-502 through the binding site in its 3'-untranslated region, and the rs16917496 polymorphism at the miR-502-binding site in the SET8 gene has been implicated in a number of cancer types. The rs16917496 polymorphism including CC, CT and TT genotypes was analyzed in patients with colorectal cancer; the CC genotype was identified to be independently associated with longer post-operative survival times using multivariate analysis (relative risk, 2.406; 95% confidence interval, 1.017-5.691; P=0.046). In addition, decreased SET8 expression was associated with the SET8 CC genotype and longer survival times for patients with colorectal cancer. The results of the present study indicated that miR-502 mediates SET8 expression at least partly by altering the binding affinity between miR-502 and SET8 so as to modify the colorectal cancer outcome. The results indicate that SET8 may be a novel target for colorectal cancer therapy.

Habitual consumption of long-chain n-3 PUFAs and fish attenuates genetically associated long-term weight gain

BACKGROUND

A growing amount of data suggests that n-3 (ω-3) polyunsaturated fatty acid (PUFA) intake may modify the genetic association with weight change.

OBJECTIVES

We aimed to prospectively test interactions of habitual consumption of n-3 PUFAs or fish, the major food source, with overall genetic susceptibility on long-term weight change.

DESIGN

Gene-diet interactions were examined in 11,330 women from the Nurses' Health Study (NHS), 6773 men from the Health Professionals Follow-Up Study (HPFS), and 6254 women from the Women's Health Initiative (WHI).

RESULTS

In the NHS and HPFS cohorts, food-sourced long-chain n-3 PUFA intake showed directionally consistent interactions with genetic risk score on long-term changes in BMI (P-interaction = 0.01 in the HPFS, 0.15 in the NHS, and 0.01 in both cohorts combined). Such interactions were successfully replicated in the WHI, an independent cohort (P-interaction = 0.02 in the WHI and 0.01 in the combined 3 cohorts). The genetic associations with changes in BMI (in kg/m2) consistently decreased (0.15, 0.10, 0.07, and -0.14 per 10 BMI-increasing alleles) across the quartiles of long-chain n-3 PUFAs in the combined cohorts. In addition, high fish intake also attenuated the genetic associations with long-term changes in BMI in the HPFS (P-interaction = 0.01), NHS (P-interaction = 0.03), WHI (P-interaction = 0.10), and the combined cohorts (P-interaction = 0.01); and the differences in BMI changes per 10 BMI-increasing alleles were 0.16, 0.06, -0.08, and -0.18, respectively, across the categories (≤1, 1∼4, 4∼6, and ≥7 servings/wk) of total fish intake. Similar interactions on body weight were observed for fish intake (P-interaction = 0.003) and long-chain n-3 PUFA intake (P-interaction = 0.12).

CONCLUSION

Our study provides replicable evidence to show that high intakes of fish and long-chain n-3 PUFAs are associated with an attenuation of the genetic association with long-term weight gain based on results from 3 prospective cohorts of Caucasians.

KATZLGO: Large-Scale Prediction of LncRNA Functions by Using the KATZ Measure Based on Multiple Networks

Aggregating evidences have shown that long non-coding RNAs (lncRNAs) generally play key roles in cellular biological processes such as epigenetic regulation, gene expression regulation at transcriptional and post-transcriptional levels, cell differentiation, and others. However, most lncRNAs have not been functionally characterized. There is an urgent need to develop computational approaches for function annotation of increasing available lncRNAs. In this article, we propose a global network-based method, KATZLGO, to predict the functions of human lncRNAs at large scale. A global network is constructed by integrating three heterogeneous networks: lncRNA-lncRNA similarity network, lncRNA-protein association network, and protein-protein interaction network. The KATZ measure is then employed to calculate similarities between lncRNAs and proteins in the global network. We annotate lncRNAs with Gene Ontology (GO) terms of their neighboring protein-coding genes based on the KATZ similarity scores. The performance of KATZLGO is evaluated on a manually annotated lncRNA benchmark and a protein-coding gene benchmark with known function annotations. KATZLGO significantly outperforms state-of-the-art computational method both in maximum F-measure and coverage. Furthermore, we apply KATZLGO to predict functions of human lncRNAs and successfully map 12,318 human lncRNA genes to GO terms.

GWAS identifies two novel colorectal cancer loci at 16q24.1 and 20q13.12

Colorectal cancer (CRC) is the fourth leading cause of cancer mortality worldwide. Genome-wide association studies (GWAS) identified more than 50 CRC loci. However, most of the previous studies were conducted in European population, and host genetic factors among Japanese population are largely remained to be identified. To identify novel loci in the Japanese population, here, we performed a large-scale GWAS using 6692 cases and 27 178 controls followed by a replication analysis using more than 11 000 case-control samples. We found the significant association of 10 loci (P < 5 × 10-8), including 2 novel loci on 16q24.1 (IRF8-FOXF1, rs847208, P = 3.15 × 10-9 and odds ratio = 1.107 with 95% confidence interval (CI) of 1.071-1.145) and 20q13.12 (TOX2, rs6065668, P = 4.47 × 10-11 and odds ratio = 0.897 with 95% CI of 0.868-0.926). Moreover, 35 previously reported single nucleotide polymorphisms (SNPs) in 24 regions were validated in the Japanese population (P < 0.05) with the same risk allele as in the previous studies. SNP rs6065668 was significantly associated with TOX2 expression in the sigmoid colon. In addition, nucleotide substitutions in the regulatory region of TOX2 were predicted to alter the binding of several transcription factors, including KLF5. Our findings elucidate the important role of genetic variations in the development of CRC in the Japanese population.

Discovery of common and rare genetic risk variants for colorectal cancer

To further dissect the genetic architecture of colorectal cancer (CRC), we performed whole-genome sequencing of 1,439 cases and 720 controls, imputed discovered sequence variants and Haplotype Reference Consortium panel variants into genome-wide association study data, and tested for association in 34,869 cases and 29,051 controls. Findings were followed up in an additional 23,262 cases and 38,296 controls. We discovered a strongly protective 0.3% frequency variant signal at CHD1. In a combined meta-analysis of 125,478 individuals, we identified 40 new independent signals at P < 5 × 10-8, bringing the number of known independent signals for CRC to ~100. New signals implicate lower-frequency variants, Krüppel-like factors, Hedgehog signaling, Hippo-YAP signaling, long noncoding RNAs and somatic drivers, and support a role for immune function. Heritability analyses suggest that CRC risk is highly polygenic, and larger, more comprehensive studies enabling rare variant analysis will improve understanding of biology underlying this risk and influence personalized screening strategies and drug development.

Functional Polymorphisms in DNA Repair Genes Are Associated with Sporadic Colorectal Cancer Susceptibility and Clinical Outcome

DNA repair processes are involved in both the onset and treatment efficacy of colorectal cancer (CRC). A change of a single nucleotide causing an amino acid substitution in the corresponding protein may alter the efficiency of DNA repair, thus modifying the CRC susceptibility and clinical outcome. We performed a candidate gene approach in order to analyze the association of non-synonymous single nucleotide polymorphisms (nsSNPs) in the genes covering the main DNA repair pathways with CRC risk and clinical outcome modifications. Our candidate polymorphisms were selected according to the foremost genomic and functional prediction databases. Sixteen nsSNPs in 12 DNA repair genes were evaluated in cohorts from the Czech Republic and Austria. Apart from the tumor-node-metastasis (TNM) stage, which occurred as the main prognostic factor in all of the performed analyses, we observed several significant associations of different nsSNPs with survival and clinical outcomes in both cohorts. However, only some of the genes (REV3L, POLQ, and NEIL3) were prominently defined as prediction factors in the classification and regression tree analysis; therefore, the study suggests their association for patient survival. In summary, we provide observational and bioinformatics evidence that even subtle alterations in specific proteins of the DNA repair pathways may contribute to CRC susceptibility and clinical outcome.

Diagnostics for Pleiotropy in Mendelian Randomization Studies: Global and Individual Tests for Direct Effects

Diagnosing pleiotropy is critical for assessing the validity of Mendelian randomization (MR) analyses. The popular MR-Egger method evaluates whether there is evidence of bias-generating pleiotropy among a set of candidate genetic instrumental variables. In this article, we propose a statistical method-global and individual tests for direct effects (GLIDE)-for systematically evaluating pleiotropy among the set of genetic variants (e.g., single nucleotide polymorphisms (SNPs)) used for MR. As a global test, simulation experiments suggest that GLIDE is nearly uniformly more powerful than the MR-Egger method. As a sensitivity analysis, GLIDE is capable of detecting outliers in individual variant-level pleiotropy, in order to obtain a refined set of genetic instrumental variables. We used GLIDE to analyze both body mass index and height for associations with colorectal cancer risk in data from the Genetics and Epidemiology of Colorectal Cancer Consortium and the Colon Cancer Family Registry (multiple studies). Among the body mass index-associated SNPs and the height-associated SNPs, several individual variants showed evidence of pleiotropy. Removal of these potentially pleiotropic SNPs resulted in attenuation of respective estimates of the causal effects. In summary, the proposed GLIDE method is useful for sensitivity analyses and improves the validity of MR.

Chemoprevention of Colorectal Cancer by Dietary Compounds

Colorectal cancer is one of the leading causes of death, and the third most diagnosed type of cancer, worldwide. It is most common amongst men and women over 50 years old. Risk factors include smoking, alcohol, diet, physical inactivity, genetics, alterations in gut microbiota, and associated pathologies (diabetes, obesity, chronic inflammatory bowel diseases). This review will discuss, in detail, the chemopreventive properties of some dietary compounds (phenolic compounds, carotenoids, iridoids, nitrogen compounds, organosulfur compounds, phytosterols, essential oil compounds, polyunsaturated fatty acids and dietary fiber) against colorectal cancer. We present recent data, focusing on in vitro, laboratory animals and clinical trials with the previously mentioned compounds. The chemopreventive properties of the dietary compounds involve multiple molecular and biochemical mechanisms of action, such as inhibition of cell growth, inhibition of tumor initiation, inhibition of adhesion, migration and angiogenesis, apoptosis, interaction with gut microbiota, regulation of cellular signal transduction pathways and xenobiotic metabolizing enzymes, etc. Moreover, this review will also focus on the natural dietary compounds' bioavailability, their synergistic protective effect, as well as the association with conventional therapy. Dietary natural compounds play a major role in colorectal chemoprevention and continuous research in this field is needed.

Genetic susceptibility in the development of colorectal adenomas according to family history of colorectal cancer

Our study aimed to evaluate the relevance of genetic susceptibility in the development of colorectal adenomas (CRA) and its relationship with the presence of family history of colorectal cancer (CRC). Genomic DNA from 750 cases (first degree relatives of patients with CRC) and 750 controls (subjects with no family history of CRC) was genotyped for 99 single nucleotide polymorphisms (SNPs) previously associated with CRC/CRA risk by GWAS and candidate gene studies by using the MassArray™ (Sequenom) platform. Cases and controls were matched by gender, age and histological lesion. Eight hundred and fifty‐eight patients showed no neoplastic lesions, whereas 288 patients showed low‐risk adenomas, and 354 patients presented high‐risk adenomas. Two SNPs (rs10505477, rs6983267) in the CASC8 gene were associated with a reduced risk of CRA in controls (log‐additive models, OR: 0.67, 95%CI:0.54–0.83, and OR:0.66, 95%CI:0.54–0.84, respectively). Stratified analysis by histological lesion revealed the association of rs10505477 and rs6983267 variants with reduced risk of low‐ and high‐risk adenomas in controls, being this effect stronger in low‐risk adenomas (log‐additive models, OR:0.63, 95%CI:0.47–0.84 and OR:0.64, 95%CI:0.47–0.86, respectively). Moreover, 2 SNPs (rs10795668, rs11255841) in the noncoding LINC00709 gene were significantly associated with a reduced risk of low‐risk adenomas in cases (recessive models, OR:0.22, 95%CI:0.06–0.72, and OR:0.08, 95%CI:0.03–0.61) and controls (dominant models, OR:0.50, 95%CI:0.34–0.75, and OR:0.52, 95%CI:0.35–0.78, respectively). In conclusion, some variants associated with CRC risk (rs10505477, rs6983267, rs10795668 and rs11255841) are also involved in the susceptibility to CRA and specific subtypes. These associations are influenced by the presence of family history of CRC.

What's new?

While numerous candidate gene variants have been associated with colorectal cancer, little is known about the relevance of genetic susceptibility or influence of family history in the development of precancerous colorectal adenomas. In the present study, certain genetic variants previously associated with colorectal cancer risk, including two variants in the CASC8 gene and two in the lnc‐RNA LINC00709 gene, were found to be also involved in susceptibility to colorectal adenomas. The associations were modified by family history of colorectal cancer. The results could have implications for colorectal cancer screening and the identification of individuals at increased risk of colorectal adenoma.

Leptin gene variants and colorectal cancer risk: Sex-specific associations

BACKGROUND

High levels of serum leptin and low levels of serum adiponectin are strongly correlated with obesity, a well-established risk factor for colorectal cancer (CRC). Growing evidence suggests that dysregulation of leptin and adiponectin levels may play an etiological role in colorectal carcinogenesis. We evaluated 20 candidate variants in 4 genes previously shown to alter serum leptin and adiponectin levels for associations with obesity (BMI>30 kg/m2) and CRC risk.

METHODS

We analyzed 6,246 CRC cases and 7,714 population-based controls from 11 studies within the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). Associations of each variant with obesity or CRC were evaluated using multivariate logistic regression models stratified by sex and adjusted for age, a study variable, and the first three principal components of genetic ancestry. Gene-specific False Discovery Rate (FDR)-adjusted p-values <0.05 denoted statistical significance.

RESULTS

Two variants in the leptin gene showed statistically significant associations with CRC among women: LEP rs2167270 (OR = 1.13, 95% CI: 1.06-1.21) and LEP rs4731426 (OR = 1.09, 95% CI: 1.02-1.17). These associations remained significant after adjustment for obesity, suggesting that leptin SNPs may influence CRC risk independent of obesity. We observed statistically significant interactions of the leptin variants with hormone replacement therapy (HRT) for CRC risk; these variant associations were strengthened when analyses were restricted to post-menopausal women with low estrogen exposure, as estimated by 'never use' of HRT and/or non-obese BMI. No variants were associated with CRC among men.

CONCLUSIONS

Leptin gene variants may exhibit sex-specific associations with CRC risk. Endogenous and exogenous estrogen exposure may modify the association between these variants, leptin levels, and CRC risk.

A genome-wide association study of energy intake and expenditure

Excessive energy intake or insufficient energy expenditure, which result in energy imbalance, contribute to the development of obesity. Obesity-related genes, such as FTO, are associated with energy traits. No genome-wide association studies (GWAS) have been conducted to detect the genetic associations with energy-related traits, including energy intake and energy expenditure, among European-ancestry populations. In this study, we conducted a genome-wide study using pooled GWAS including 12,030 European-ancestry women and 6,743 European-ancestry men to identify genetic variants associated with these two energy traits. We observed a statistically significant genome-wide SNP heritability for energy intake of 6.05% (95%CI = (1.76, 10.34), P = 0.006); the SNP heritability for expenditure was not statistically significantly greater than zero. We discovered three SNPs on chromosome 12q13 near gene ANKRD33 that were genome-wide significantly associated with increased total energy intake among all men. We also identified signals on region 2q22 that were associated with energy expenditure among lean people. Body mass index related SNPs were found to be significantly associated with energy intake and expenditure through SNP set analyses. Larger GWAS studies of total energy traits are warranted to explore the genetic basis of energy intake, including possible differences between men and women, and the association between total energy intake and other downstream phenotypes, such as diabetes and chronic diseases.

Genetic Risk Score Is Associated With Prevalence of Advanced Neoplasms in a Colorectal Cancer Screening Population

BACKGROUND & AIMS

The presence of specific single nucleotide polymorphisms (SNPs) can be used to calculate an individual's risk for colorectal cancer (CRC), called a genetic risk score (GRS). We investigated whether GRS can identify individuals with clinically relevant neoplasms in a screening colonoscopy population.

METHODS

We derived a GRS based on 48 SNPs associated with CRC, identified in a comprehensive literature search. We obtained genetic data from 1043 participants (50-79 years old) in a screening colonoscopy study in Germany, recruited from 2005 through 2013 (294 with advanced neoplasms, 249 with non-advanced adenoma (NAAs), and 500 without neoplasms). Each participant was assigned a GRS by aggregating their risk alleles (0, 1, or 2). Risk of advanced neoplasms and NAA according to GRS was calculated by multiple logistic regression. Risk advancement periods were calculated. We replicated our findings using data from a subset of the Tennessee Colorectal Polyp Study.

RESULTS

An increased GRS was associated with higher prevalence of advanced neoplasms, but not NAAs. Participants in the middle and upper tertiles of GRS had a 2.2-fold and 2.7-fold increase in risk, respectively, of advanced neoplasms compared to those in the lower tertile. Adjusted odds ratios (ORs) were 1.09 (95% confidence interval [CI], 0.76-1.57) for NAA in the middle tertile and 1.05 (95% CI, 0.70-1.55) for NAA in the upper tertile. The ORs were largest for proximal advanced neoplasms for participants in the middle tertile (OR, 3.55; 95% CI 1.85-6.82) and the upper tertile (OR, 3.61; 95% CI 1.84-7.10). The risk advancement period for medium vs low GRS was 13.4 years (95% CI 4.8-22.0) and for high vs low GRS was 17.5 years (95% CI, 7.8-27.3).

CONCLUSIONS

In a genetic analysis of participants in a CRC screening study in Germany, an increased GRS (based on CRC-associated SNPs) was associated with increased prevalence of advanced neoplasms. These findings might be used in defining risk-adapted screening ages.

Coding variants in NOD-like receptors: An association study on risk and survival of colorectal cancer

Nod-like receptors (NLRs) are important innate pattern recognition receptors and regulators of inflammation or play a role during development. We systematically analysed 41 non-synonymous single nucleotide polymorphisms (SNPs) in 21 NLR genes in a Czech discovery cohort of sporadic colorectal cancer (CRC) (1237 cases, 787 controls) for their association with CRC risk and survival. Five SNPs were found to be associated with CRC risk and eight with survival at 5% significance level. In a replication analysis using data of two large genome-wide association studies (GWASs) from Germany (DACHS: 1798 cases and 1810 controls) and Scotland (2210 cases and 9350 controls) the associations found in the Czech discovery set were not confirmed. However, expression analysis in human gut-related tissues and immune cells revealed that the NLRs associated with CRC risk or survival in the discovery set were expressed in primary human colon or rectum cells, CRC tissue and/or cell lines, providing preliminary evidence for a potential involvement of NLRs in general in CRC development and/or progression. Most interesting was the finding that the enigmatic development-related NLRP5 (also known as MATER) was not expressed in normal colon tissue but in colon cancer tissue and cell lines. Future studies may show whether regulatory variants instead of coding variants might affect the expression of NLRs and contribute to CRC risk and survival.

Determining Risk of Colorectal Cancer and Starting Age of Screening Based on Lifestyle, Environmental, and Genetic Factors

BACKGROUND & AIMS

Guidelines for initiating colorectal cancer (CRC) screening are based on family history but do not consider lifestyle, environmental, or genetic risk factors. We developed models to determine risk of CRC, based on lifestyle and environmental factors and genetic variants, and to identify an optimal age to begin screening.

METHODS

We collected data from 9748 CRC cases and 10,590 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium and the Colorectal Transdisciplinary study, from 1992 through 2005. Half of the participants were used to develop the risk determination model and the other half were used to evaluate the discriminatory accuracy (validation set). Models of CRC risk were created based on family history, 19 lifestyle and environmental factors (E-score), and 63 CRC-associated single-nucleotide polymorphisms identified in genome-wide association studies (G-score). We evaluated the discriminatory accuracy of the models by calculating area under the receiver operating characteristic curve values, adjusting for study, age, and endoscopy history for the validation set. We used the models to project the 10-year absolute risk of CRC for a given risk profile and recommend ages to begin screening in comparison to CRC risk for an average individual at 50 years of age, using external population incidence rates for non-Hispanic whites from the Surveillance, Epidemiology, and End Results program registry.

RESULTS

In our models, E-score and G-score each determined risk of CRC with greater accuracy than family history. A model that combined both scores and family history estimated CRC risk with an area under the receiver operating characteristic curve value of 0.63 (95% confidence interval, 0.62-0.64) for men and 0.62 (95% confidence interval, 0.61-0.63) for women; area under the receiver operating characteristic curve values based on only family history ranged from 0.53 to 0.54 and those based only E-score or G-score ranged from 0.59 to 0.60. Although screening is recommended to begin at age 50 years for individuals with no family history of CRC, starting ages calculated based on combined E-score and G-score differed by 12 years for men and 14 for women, for individuals with the highest vs the lowest 10% of risk.

CONCLUSIONS

We used data from 2 large international consortia to develop CRC risk calculation models that included genetic and environmental factors along with family history. These determine risk of CRC and starting ages for screening with greater accuracy than the family history only model, which is based on the current screening guideline. These scoring systems might serve as a first step toward developing individualized CRC prevention strategies.

Mendelian randomisation study of age at menarche and age at menopause and the risk of colorectal cancer

Background

Substantial evidence supports an association between use of menopausal hormone therapy and decreased colorectal cancer (CRC) risk, indicating a role of exogenous sex hormones in CRC development. However, findings on endogenous oestrogen exposure and CRC are inconsistent.

Methods

We used a Mendelian randomisation approach to test for a causal effect of age at menarche and age at menopause as surrogates for endogenous oestrogen exposure on CRC risk. Weighted genetic risk scores based on 358 single-nucleotide polymorphisms associated with age at menarche and 51 single-nucleotide polymorphisms associated with age at menopause were used to estimate the association with CRC risk using logistic regression in 12,944 women diagnosed with CRC and 10,741 women without CRC from three consortia. Sensitivity analyses were conducted to address pleiotropy and possible confounding by body mass index.

Results

Genetic risk scores for age at menarche (odds ratio per year 0.98, 95% confidence interval: 0.95–1.02) and age at menopause (odds ratio 0.98, 95% confidence interval: 0.94–1.01) were not significantly associated with CRC risk. The sensitivity analyses yielded similar results.

Conclusions

Our study does not support a causal relationship between genetic risk scores for age at menarche and age at menopause and CRC risk.

Expression level of TBX3 gene in renal carcinoma and its clinical significance

The aim of the current study was to investigate and discuss the function of T-box 3 (TBX3) gene expression in the pathogenesis of renal carcinoma. The carcinoma, adjacent and normal renal tissues of 210 patients with renal carcinoma who presented to The Central Hospital of Wuhan, Tongji Medical College from March, 2006 to March, 2012 were collected to extract total RNAs. The total RNAs were reverse-transcribed into complementary DNA (cDNA), and quantitative polymerase chain reaction (qPCR) was applied to detect the expression of TBX3 gene in these tissues, followed by its association with the prognosis of renal carcinoma as well as clinical features. A comparison of the renal carcinoma tissues with the adjacent tissues showed that TBX3 gene was obviously highly expressed in renal carcinoma tissues (P<0.05). In addition, compared with normal renal tissues, TBX3 gene was obviously highly expressed in renal carcinoma tissues (P<0.05). There was no significant difference in the expression levels of TBX3 gene in normal renal tissues and adjacent tissues (P=0.15). The expression of TBX3 gene in renal carcinoma tissues was not associated with patient age, sex and tumor size (P>0.05), but it was associated with tumor-node-metastasis (TNM) staging and lymph node metastasis (P<0.05). The Kaplan-Meier survival analysis revealed that the median survival time of patients in the positive TBX3 gene expression group (37.5 months) was shorter than that in the negative TBX3 gene expression group (66 months), and there was a statistical difference (P<0.05). The 3- and 5-year survival rates in the negative TBX3 gene expression group were 74 and 62%, respectively, and the 3- and 5-year survival rates in the positive TBX3 gene expression group were 52 and 32%, respectively, and the differences were significant (P<0.05). The results suggest that TBX3 gene is highly expressed in renal carcinoma tissues, and it is associated with TNM staging, lymph node metastasis and distant metastasis, which may be involved in the occurrence and metastasis of renal carcinoma.

Recurrent, low-frequency coding variants contributing to colorectal cancer in the Swedish population

Genome-wide association studies (GWAS) have identified dozens of common genetic variants associated with risk of colorectal cancer (CRC). However, the majority of CRC heritability remains unclear. In order to discover low-frequency, high-risk CRC susceptibility variants in Swedish population, we genotyped 1 515 CRC patients enriched for familial cases, and 12 108 controls. Case/control association analysis suggested eight novel variants associated with CRC risk (OR 2.0-17.6, p-value < 2.0E-07), comprised of seven coding variants in genes RAB11FIP5, POTEA, COL27A1, MUC5B, PSMA8, MYH7B, and PABPC1L as well as one variant downstream of NEU1 gene. We also confirmed 27 out of 30 risk variants previously reported from GWAS in CRC with a mixed European population background. This study identified rare, coding sequence variants associated with CRC risk through analysis in a relatively homogeneous population. The segregation data suggest a complex mode of inheritance in seemingly dominant pedigrees.

Genetic Mechanisms of Immune Evasion in Colorectal Cancer

To understand the genetic drivers of immune recognition and evasion in colorectal cancer, we analyzed 1,211 colorectal cancer primary tumor samples, including 179 classified as microsatellite instability-high (MSI-high). This set includes The Cancer Genome Atlas colorectal cancer cohort of 592 samples, completed and analyzed here. MSI-high, a hypermutated, immunogenic subtype of colorectal cancer, had a high rate of significantly mutated genes in important immune-modulating pathways and in the antigen presentation machinery, including biallelic losses of B2M and HLA genes due to copy-number alterations and copy-neutral loss of heterozygosity. WNT/β-catenin signaling genes were significantly mutated in all colorectal cancer subtypes, and activated WNT/β-catenin signaling was correlated with the absence of T-cell infiltration. This large-scale genomic analysis of colorectal cancer demonstrates that MSI-high cases frequently undergo an immunoediting process that provides them with genetic events allowing immune escape despite high mutational load and frequent lymphocytic infiltration and, furthermore, that colorectal cancer tumors have genetic and methylation events associated with activated WNT signaling and T-cell exclusion.Significance: This multi-omic analysis of 1,211 colorectal cancer primary tumors reveals that it should be possible to better monitor resistance in the 15% of cases that respond to immune blockade therapy and also to use WNT signaling inhibitors to reverse immune exclusion in the 85% of cases that currently do not. Cancer Discov; 8(6); 730-49. ©2018 AACR.This article is highlighted in the In This Issue feature, p. 663.

Risk of eighteen genome-wide association study-identified genetic variants for colorectal cancer and colorectal adenoma in Han Chinese

Background

Recent genome-wide association studies (GWAS) identified eighteen single-nucleotide polymorphisms (SNPs) to be significantly associated with the risk of colorectal cancer (CRC). However, overall results of the following replications are inconsistent and little is known about whether these associations also exit in colorectal adenomas (CRA).

Methods

The SNP genotyping was performed using a Sequenom MassARRAY to investigate the association of these eighteen SNPs with colorectal neoplasm in a case-control study consisted of 1049 colorectal cancers, 283 adenomas, and 1030 controls.

Results

Two of these SNPs, rs10505477 and rs719725, showed evidence of an association in both CRC and CRA in our study population. Besides, seven SNPs (rs10808555, rs7014346, rs7837328, rs704017, rs11196172, rs4779584, and rs7229639) were significantly associated with CRC, and another one SNP rs11903757 was over-represented in CRA compared with controls. The strongest association was provided by rs11196172 (OR = 2.02, 95% CI = 1.66 - 2.46, P < 0.0001) and rs11903757 (OR = 1.96, 95% CI = 1.28 - 3.00, P = 0.0026).

Conclusion

These results suggest that some previously reported SNP associations also have impact on CRC and CRA predispositions in the Han Chinese population. A part of genetic risk to CRC is possibly mediated by susceptibility to adenomas.

Effects of interactions between common genetic variants and alcohol consumption on colorectal cancer risk

BACKGROUND

Genome-wide association studies (GWAS) have identified approximately 40 common genetic loci associated with colorectal cancer risk. To investigate possible gene-environment interactions (GEIs) between GWAS-identified single-nucleotide polymorphisms (SNPs) and alcohol consumption with respect to colorectal cancer, a hospital-based case-control study was conducted.

RESULTS

Higher levels of alcohol consumption as calculated based on a standardized definition of a drink (1 drink=12.5g of ethanol) were associated with increased risk of colorectal cancer (OR=2.47, 95% CI=1.62-3.76 for heavy drinkers [>50g/day] compared to never drinkers; ptrend<0.01). SNP rs6687758 near the DUSP10 gene at 1q41 had a statistically significant interaction with alcohol consumption in analyses of standardized drinks (p=4.6×10-3), although this did not surpass the corrected threshold for multiple testing. When stratified by alcohol consumption levels, in an additive model the risk of colorectal cancer associated with the G allele of rs6687758 tended to increase among individuals in the heavier alcohol consumption strata. A statistically significant association between rs6687758 and colorectal cancer risk was observed among moderate alcohol drinkers who consumed between >12.5 and ≤50g of alcohol per day (OR=1.46, 95% CI=1.01-2.11).

METHODS

A total of 2,109 subjects (703 colorectal cancer patients and 1,406 healthy controls) were recruited from the Korean National Cancer Center. For genotyping, 30 GWAS-identified SNPs were selected. A logistic regression model was used to evaluate associations of SNPs and alcohol consumption with colorectal cancer risk. We also tested GEIs between SNPs and alcohol consumption using a logistic model with multiplicative interaction terms.

CONCLUSIONS

Our results suggest that SNP rs6687758 at 1q41 may interact with alcohol consumption in the etiology of colorectal cancer.

Strongly enhanced colorectal cancer risk stratification by combining family history and genetic risk score

Background and aim

Family history (FH) and genetic risk scores (GRSs) are increasingly used for risk stratification for colorectal cancer (CRC) screening. However, they were mostly considered alternatively rather than jointly. The aim of this study was to assess the potential of individual and joint risk stratification for CRC by FH and GRS.

Patients and methods

A GRS was built based on the number of risk alleles in 53 previously identified single-nucleotide polymorphisms among 2,363 patients with a first diagnosis of CRC and 2,198 controls in DACHS [colorectal cancer: chances for prevention through screening], a population-based case-control study in Germany. Associations between GRS and FH with CRC risk were quantified by multiple logistic regression.

Results

A total of 316 cases (13.4%) and 214 controls (9.7%) had a first-degree relative (FDR) with CRC (adjusted odds ratio [aOR] 1.86, 95% CI 1.52–2.29). A GRS in the highest decile was associated with a 3.0-fold increased risk of CRC (aOR 3.00, 95% CI 2.24–4.02) compared with the lowest decile. This association was tentatively more pronounced in older age groups. FH and GRS were essentially unrelated, and their joint consideration provided more accurate risk stratification than risk stratification based on each of the variables individually. For example, risk was 6.1-fold increased in the presence of both FH in a FDR and a GRS in the highest decile (aOR 6.14, 95% CI 3.47–10.84) compared to persons without FH and a GRS in the lowest decile.

Conclusion

Both FH and the so far identified genetic variants carry essentially independent risk information and in combination provide great potential for CRC risk stratification.

Identification of a functional variant for colorectal cancer risk mapping to chromosome 5q31.1

Genome-wide association studies (GWASs) have established chromosome 5q31.1 as a risk locus for colorectal cancer (CRC). We previously identified a potentially regulatory single nucleotide polymorphism (SNP) rs17716310 within 5q31.1. Now, we extended our study with another independent Chinese population, functional assays and analyses of TCGA (The Cancer Genome Atlas) data. Significant associations between rs17716310 and CRC risk were found in Present Study including 1075 CRC cases and 1999 controls (additive model: OR = 1.149, 95% CI = 1.027–1.286, P = 0.016), and in Combined Study including 1766 cases and 2708 controls (additive model: OR = 1.145, 95% CI = 1.045–1.254, P = 0.004). Dual luciferase reporter gene assays indicated that the variant C allele obviously increased transcriptional activity. Using TCGA datasets, we indicated rs17716310 as a cis expression quantitative trait locus (eQTL) for the gene SMAD5, whose expression was significantly higher in CRC tissues. These findings suggested that the functional polymorphism rs17716310 A > C might be a genetic modifier for CRC, promoting the expression of SMAD5 that belonged to the transforming growth factor beta (TGF-β) signaling pathway.

Effects of interactions between common genetic variants and smoking on colorectal cancer

Background

Although genome-wide association studies (GWAS) have identified variants in approximately 40 susceptibility loci for colorectal cancer (CRC), there are few studies on the interactions between identified single-nucleotide polymorphisms (SNPs) and lifestyle risk factors. We evaluated whether smoking could modify associations between these genetic variants and CRC risk.

Methods

A total of 703 CRC patients and 1406 healthy controls were included in this case-control study from the National Cancer Center in Korea. Thirty CRC susceptibility SNPs identified in previous GWAS were genotyped. A logistic regression model was used to examine associations between the SNPs and smoking behaviors by sex. The interaction was estimated by including an additional interaction term in the model.

Results

In men, an increased CRC risk was observed for longer durations (OR_{>28 vs. ≤28years} = 1.49 (95% CI = 1.11–1.98)), greater quantities (OR_{≥20 vs. <20cigarettes/day} = 2.12 (1.61–2.79)), and longer pack-years of smoking (OR_{≥21 vs. <21pack-years} = 1.78 (1.35–2.35)). In women, longer pack-years of smoking significantly increased CRC risk (OR_{≥5 vs. <5pack-years} = 6.11 (1.10–34.00)). Moreover, there were significant interactions between smoking status and the polymorphisms rs1957636 at 14q22.3 (P_interaction = 5.5 × 10⁻⁴) and rs4813802 at 20p12.3 (P_interaction = 0.04) in men. Interactions between smoking status and the rs6687758 at 1q41 (P_interaction = 0.03), duration and the rs174537 at 11q12.2 (P_interaction = 0.05), and pack-years and the rs4813802 (P_interaction = 0.04) were also found in women.

Conclusions

Associations between susceptibility SNPs and CRC risk may be modified by smoking behaviors, supporting the existence of gene-smoking interactions.

Electronic supplementary material

The online version of this article (10.1186/s12885-017-3886-0) contains supplementary material, which is available to authorized users.

The Mbd4 DNA glycosylase protects mice from inflammation-driven colon cancer and tissue injury

Much of the global cancer burden is associated with longstanding inflammation accompanied by release of DNA-damaging reactive oxygen and nitrogen species. Here, we report that the Mbd4 DNA glycosylase is protective in the azoxymethane/dextran sodium sulfate (AOM/DSS) mouse model of inflammation-driven colon cancer. Mbd4 excises T and U from T:G and U:G mismatches caused by deamination of 5-methylcytosine and cytosine. Since the rate of deamination is higher in inflamed tissues, we investigated the role of Mbd4 in inflammation-driven tumorigenesis. In the AOM/DSS assay, Mbd4-/- mice displayed more severe clinical symptoms, decreased survival, and a greater tumor burden than wild-type (WT) controls. The increased tumor burden in Mbd4-/- mice did not arise from impairment of AOM-induced apoptosis in the intestinal crypt. Histopathological analysis indicated that the colonic epithelium of Mbd4-/- mice is more vulnerable than WT to DSS-induced tissue damage. We investigated the role of the Mbd4-/- immune system in AOM/DSS-mediated carcinogenesis by repeating the assay on WT and Mbd4-/- mice transplanted with WT bone marrow. Mbd4-/- mice with WT bone marrow behaved similarly to Mbd4-/- mice. Together, our results indicate that the colonic epithelium of Mbd4-/- mice is more vulnerable to DSS-induced injury, which exacerbates inflammation-driven tissue injury and cancer.

SMAD7 loci contribute to risk of hepatocellular carcinoma and clinicopathologic development among Chinese Han population

Genome-wide association studies (GWAS) have identified three loci at 18q21 (rs4939827, rs7240004, and rs7229639), which maps to SMAD7 loci, were associated with risk of diseases of the digestive system. However, their associations with hepatocellular carcinoma (HCC) risk remain unknown. A case-control study was conducted to assess genetic associations with HCC risk and clinicopathologic development among Chinese Han population. Three SNPs were genotyped among 1,000 HCC cases and 1,000 controls using Sequenom Mass-ARRAY technology. We observed statistically significant associations for the three SMAD7 loci and HCC risk. Each copy of minor allele was associated with a 1.24–1.36 fold increased risk of HCC. We also found that significant differences were observed between rs4939827 and clinical TNM stage and vascular invasion, as well as rs7240004 and vascular invasion. We also established a genetic risk score (GRS) by summing the risk alleles. The GRS was significantly associated with increased risk of HCC and vascular invasion. Our data revealed the SMAD7 loci is associated with HCC susceptibility and its clinicopathologic development.

Colorectal cancer susceptibility loci as predictive markers of rectal cancer prognosis after surgery

To understand the molecular mechanism of rectal cancer and develop markers for disease prognostication, we generated and explored a dataset from 243 rectal cancer patients by gene expression microarray analysis of cancer samples and matched controls, and SNP-arrays of germline DNA. We found that two of the loci most strongly linked with colorectal cancer (CRC) risk, 8q24 (upstream of MYC) and 18q21 (in the intron of SMAD7), as well as 20q13 (in the intron of LAMA5), are tightly associated with the prognosis of rectal cancer patients. For SNPs on 18q21 (rs12953717 and rs4464148) and 20q13 (rs4925386), alleles that correlate with higher risk for the development of CRC are associated with shorter disease free survival (DFS). However, for rs6983267 on 8q24, the low risk allele is associated with a higher risk for recurrence and metastasis after surgery, and importantly, is strongly correlated with the resistance of CRC cell lines to chemoradiotherapy (CRT). We also found that although MYC expression is dramatically increased in cancer, patients with higher levels of MYC have a better prognosis. The expression of SMAD7 is weakly correlated with DFS. Notably, the presence of the 8q24 and 18q21 SNP alleles is not correlated with expression levels of MYC and SMAD7. rs4464148, and probably rs6983267 and rs4925386, are linked with overall survival time of patients. In conclusion, we show that several CRC risk SNPs detect subpopulations of rectal cancer patients with poor prognosis, and that rs6983267 probably affects prognosis through interfering with the resistance of cancer cells to CRT.

Enrichment of colorectal cancer associations in functional regions: Insight for using epigenomics data in the analysis of whole genome sequence-imputed GWAS data

BACKGROUND

The evaluation of less frequent genetic variants and their effect on complex disease pose new challenges for genomic research. To investigate whether epigenetic data can be used to inform aggregate rare-variant association methods (RVAM), we assessed whether variants more significantly associated with colorectal cancer (CRC) were preferentially located in non-coding regulatory regions, and whether enrichment was specific to colorectal tissues.

METHODS

Active regulatory elements (ARE) were mapped using data from 127 tissues and cell-types from NIH Roadmap Epigenomics and Encyclopedia of DNA Elements (ENCODE) projects. We investigated whether CRC association p-values were more significant for common variants inside versus outside AREs, or 2) inside colorectal (CR) AREs versus AREs of other tissues and cell-types. We employed an integrative epigenomic RVAM for variants with allele frequency <1%. Gene sets were defined as ARE variants within 200 kilobases of a transcription start site (TSS) using either CR ARE or ARE from non-digestive tissues. CRC-set association p-values were used to evaluate enrichment of less frequent variant associations in CR ARE versus non-digestive ARE.

RESULTS

ARE from 126/127 tissues and cell-types were significantly enriched for stronger CRC-variant associations. Strongest enrichment was observed for digestive tissues and immune cell types. CR-specific ARE were also enriched for stronger CRC-variant associations compared to ARE combined across non-digestive tissues (p-value = 9.6 × 10-4). Additionally, we found enrichment of stronger CRC association p-values for rare variant sets of CR ARE compared to non-digestive ARE (p-value = 0.029).

CONCLUSIONS

Integrative epigenomic RVAM may enable discovery of less frequent variants associated with CRC, and ARE of digestive and immune tissues are most informative. Although distance-based aggregation of less frequent variants in CR ARE surrounding TSS showed modest enrichment, future association studies would likely benefit from joint analysis of transcriptomes and epigenomes to better link regulatory variation with target genes.

Inherited variation in circadian rhythm genes and risks of prostate cancer and three other cancer sites in combined cancer consortia

Circadian disruption has been linked to carcinogenesis in animal models, but the evidence in humans is inconclusive. Genetic variation in circadian rhythm genes provides a tool to investigate such associations. We examined associations of genetic variation in nine core circadian rhythm genes and six melatonin pathway genes with risk of colorectal, lung, ovarian and prostate cancers using data from the Genetic Associations and Mechanisms in Oncology (GAME-ON) network. The major results for prostate cancer were replicated in the Prostate, Lung, Colorectal and Ovarian (PLCO) cancer screening trial, and for colorectal cancer in the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). The total number of cancer cases and controls was 15,838/18,159 for colorectal, 14,818/14,227 for prostate, 12,537/17,285 for lung and 4,369/9,123 for ovary. For each cancer site, we conducted gene-based and pathway-based analyses by applying the summary-based Adaptive Rank Truncated Product method (sARTP) on the summary association statistics for each SNP within the candidate gene regions. Aggregate genetic variation in circadian rhythm and melatonin pathways were significantly associated with the risk of prostate cancer in data combining GAME-ON and PLCO, after Bonferroni correction (ppathway  < 0.00625). The two most significant genes were NPAS2 (pgene  = 0.0062) and AANAT (pgene  = 0.00078); the latter being significant after Bonferroni correction. For colorectal cancer, we observed a suggestive association with the circadian rhythm pathway in GAME-ON (ppathway  = 0.021); this association was not confirmed in GECCO (ppathway  = 0.76) or the combined data (ppathway  = 0.17). No significant association was observed for ovarian and lung cancer. These findings support a potential role for circadian rhythm and melatonin pathways in prostate carcinogenesis. Further functional studies are needed to better understand the underlying biologic mechanisms.

A Novel Approach for Pathway Analysis of GWAS Data Highlights Role of BMP Signaling and Muscle Cell Differentiation in Colorectal Cancer Susceptibility

Genome-wide association studies (GWAS) have revolutionized the field of gene mapping. As the GWAS field matures, it is becoming clear that for many complex traits, a proportion of the missing heritability is attributable to common variants of individually small effect. Detecting these small effects individually can be difficult, and statistical power would be increased if relevant variants could be grouped together for testing. Here, we propose a VEGAS2Pathway approach that aggregates association strength of individual markers into pre-specified biological pathways. It accounts for gene size and linkage disequilibrium between markers using simulations from the multivariate normal distribution. Pathway size is taken into account via a resampling approach. Importantly, since the approach only requires summary data, the method can easily be applied in all GWASs, including meta-analysis, singleton-based, family-based, and DNA-pooling-based designs. This approach is implemented in a user-friendly web page https://vegas2.qimrberghofer.edu.au and a command line tool. The web implementation uses gene-sets from the gene ontology (GO), curated gene-sets from MSigDB (containing canonical pathways and gene-sets from BIOCARTA, REACTOME, KEGG databases), PANTHER, and pathway commons databases, enabling analysis of a wide range of complex traits. We applied this method on a colorectal cancer GWAS meta-analysis data set (10,934 cases, 12,328 controls) from the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). We report statistically significant enrichment of association signal for the 'BMP signaling' and 'muscle cell differentiation' pathways, suggesting a possible role for these pathways onto the risk of colorectal cancer.

Altered microRNA profiles during early colon adenoma progression in a porcine model of familial adenomatous polyposis

MicroRNAs are dysregulated in various cancers including colorectal cancer, and are potential useful biomarkers of disease development. We used next generation sequencing to investigate miRNA expression profiles in low- and high-grade intraepithelial dysplastic polyps from pigs carrying a mutation in the adenomatous polyposis coli tumour suppressor (APC¹³¹¹ , orthologous to human APC¹³⁰⁹ ) that model an inherited predisposition to colorectal cancer, familial adenomatous polyposis. We identified several miRNAs and their isomiRs significantly (P < 0.05) differentially expressed between low and high-grade intraepithelial dysplastic polyps. Of these, ssc-let-7e, ssc-miR-98, ssc-miR-146a-5p, ssc-miR-146b, ssc-miR-183 and ssc-miR-196a were expressed at higher level and ssc-miR-126-3p at lower level in high-grade intraepithelial dysplastic polyps. Functional miRNA target analysis revealed significant (P < 0.001) over-representation of cancer-related pathways, including 'microRNAs in cancer', 'proteoglycans in cancer', 'pathways in cancer' and 'colorectal cancer'. This is the first study to reveal miRNAs associated with premalignant transformation of colon polyps.

Bayesian and frequentist analysis of an Austrian genome-wide association study of colorectal cancer and advanced adenomas

Most genome-wide association studies (GWAS) were analyzed using single marker tests in combination with stringent correction procedures for multiple testing. Thus, a substantial proportion of associated single nucleotide polymorphisms (SNPs) remained undetected and may account for missing heritability in complex traits. Model selection procedures present a powerful alternative to identify associated SNPs in high-dimensional settings. In this GWAS including 1060 colorectal cancer cases, 689 cases of advanced colorectal adenomas and 4367 controls we pursued a dual approach to investigate genome-wide associations with disease risk applying both, single marker analysis and model selection based on the modified Bayesian information criterion, mBIC2, implemented in the software package MOSGWA. For different case-control comparisons, we report models including between 1-14 candidate SNPs. A genome-wide significant association of rs17659990 (P=5.43×10^-9, DOCK3, chromosome 3p21.2) with colorectal cancer risk was observed. Furthermore, 56 SNPs known to influence susceptibility to colorectal cancer and advanced adenoma were tested in a hypothesis-driven approach and several of them were found to be relevant in our Austrian cohort. After correction for multiple testing (α=8.9×10^-4), the most significant associations were observed for SNPs rs10505477 (P=6.08×10^-4) and rs6983267 (P=7.35×10^-4) of CASC8, rs3802842 (P=8.98×10^-5, COLCA1,2), and rs12953717 (P=4.64×10^-4, SMAD7). All previously unreported SNPs demand replication in additional samples. Reanalysis of existing GWAS datasets using model selection as tool to detect SNPs associated with a complex trait may present a promising resource to identify further genetic risk variants not only for colorectal cancer.

Hemidesmosome integrity protects the colon against colitis and colorectal cancer

OBJECTIVE

Epidemiological and clinical data indicate that patients suffering from IBD with long-standing colitis display a higher risk to develop colorectal high-grade dysplasia. Whereas carcinoma invasion and metastasis rely on basement membrane (BM) disruption, experimental evidence is lacking regarding the potential contribution of epithelial cell/BM anchorage on inflammation onset and subsequent neoplastic transformation of inflammatory lesions. Herein, we analyse the role of the α6β4 integrin receptor found in hemidesmosomes that attach intestinal epithelial cells (IECs) to the laminin-containing BM.

DESIGN

We developed new mouse models inducing IEC-specific ablation of α6 integrin either during development (α6^ΔIEC) or in adults (α6^ΔIEC-TAM).

RESULTS

Strikingly, all α6^ΔIEC mutant mice spontaneously developed long-standing colitis, which degenerated overtime into infiltrating adenocarcinoma. The sequence of events leading to disease onset entails hemidesmosome disruption, BM detachment, IL-18 overproduction by IECs, hyperplasia and enhanced intestinal permeability. Likewise, IEC-specific ablation of α6 integrin induced in adult mice (α6^ΔIEC-TAM) resulted in fully penetrant colitis and tumour progression. Whereas broad-spectrum antibiotic treatment lowered tissue pathology and IL-1β secretion from infiltrating myeloid cells, it failed to reduce Th1 and Th17 response. Interestingly, while the initial intestinal inflammation occurred independently of the adaptive immune system, tumourigenesis required B and T lymphocyte activation.

CONCLUSIONS

We provide for the first time evidence that loss of IECs/BM interactions triggered by hemidesmosome disruption initiates the development of inflammatory lesions that progress into high-grade dysplasia and carcinoma. Colorectal neoplasia in our mouse models resemble that seen in patients with IBD, making them highly attractive for discovering more efficient therapies.

"Stealth dissemination" of macrophage-tumor cell fusions cultured from blood of patients with pancreatic ductal adenocarcinoma

Here we describe isolation and characterization of macrophage-tumor cell fusions (MTFs) from the blood of pancreatic ductal adenocarcinoma (PDAC) patients. The MTFs were generally aneuploidy, and immunophenotypic characterizations showed that the MTFs express markers characteristic of PDAC and stem cells, as well as M2-polarized macrophages. Single cell RNASeq analyses showed that the MTFs express many transcripts implicated in cancer progression, LINE1 retrotransposons, and very high levels of several long non-coding transcripts involved in metastasis (such as MALAT1). When cultured MTFs were transplanted orthotopically into mouse pancreas, they grew as obvious well-differentiated islands of cells, but they also disseminated widely throughout multiple tissues in "stealth" fashion. They were found distributed throughout multiple organs at 4, 8, or 12 weeks after transplantation (including liver, spleen, lung), occurring as single cells or small groups of cells, without formation of obvious tumors or any apparent progression over the 4 to 12 week period. We suggest that MTFs form continually during PDAC development, and that they disseminate early in cancer progression, forming "niches" at distant sites for subsequent colonization by metastasis-initiating cells.

Evaluation of genetic variants in association with colorectal cancer risk and survival in Asians

Genome-wide association studies (GWAS) have identified over 40 genetic loci associated with colorectal cancer (CRC) risk. The association of single nucleotide polymorphisms (SNPs) at these loci with CRC risk and survival has not been adequately evaluated in East Asians. GWAS-identified CRC risk variants were used to construct weighted genetic risk scores (GRSs). We evaluated these GRSs in association with CRC risk in 3,303 CRC cases and 3,553 controls using logistic regression models. Associations with overall and CRC-specific survival were assessed in 731 CRC patients using Cox regression models. The association between the GRSs (overall and Asian-specific) and CRC risk was approximately twofold (highest vs. lowest quintile), and the shape of the dose-response was linear (ptrend  = 1.24 × 10-13 and 3.02 × 10-14 for overall GRS and Asian-specific GRS, respectively). The association of the GRS with CRC risk was stronger among those with a family history of CRC (pinteraction  = 0.007). Asian-specific GRS using previously reported survival SNPs increased risk for mortality and the shape of the dose-response was linear for CRC-specific and all-cause mortality (ptrend  = 0.01 and 0.006, respectively). Furthermore, the minor alleles of rs6983267 and rs1957636 were associated with worse CRC-specific and overall survival. We show that GRSs constructed using GWAS-identified common variants are strongly associated with CRC risk in Asians. We confirm previous findings for the possible association between some SNPs with survival, and provide evidence for two additional CRC risk variants that may be related to CRC survival.

Investigating MicroRNA and transcription factor co-regulatory networks in colorectal cancer

Background

Colorectal cancer (CRC) is one of the most common malignancies worldwide with poor prognosis. Studies have showed that abnormal microRNA (miRNA) expression can affect CRC pathogenesis and development through targeting critical genes in cellular system. However, it is unclear about which miRNAs play central roles in CRC’s pathogenesis and how they interact with transcription factors (TFs) to regulate the cancer-related genes.

Results

To address this issue, we systematically explored the major regulation motifs, namely feed-forward loops (FFLs), that consist of miRNAs, TFs and CRC-related genes through the construction of a miRNA-TF regulatory network in CRC. First, we compiled CRC-related miRNAs, CRC-related genes, and human TFs from multiple data sources. Second, we identified 13,123 3-node FFLs including 25 miRNA-FFLs, 13,005 TF-FFLs and 93 composite-FFLs, and merged the 3-node FFLs to construct a CRC-related regulatory network. The network consists of three types of regulatory subnetworks (SNWs): miRNA-SNW, TF-SNW, and composite-SNW. To enhance the accuracy of the network, the results were filtered by using The Cancer Genome Atlas (TCGA) expression data in CRC, whereby we generated a core regulatory network consisting of 58 significant FFLs. We then applied a hub identification strategy to the significant FFLs and found 5 significant components, including two miRNAs (hsa-miR-25 and hsa-miR-31), two genes (ADAMTSL3 and AXIN1) and one TF (BRCA1). The follow up prognosis analysis indicated all of the 5 significant components having good prediction of overall survival of CRC patients.

Conclusions

In summary, we generated a CRC-specific miRNA-TF regulatory network, which is helpful to understand the complex CRC regulatory mechanisms and guide clinical treatment. The discovered 5 regulators might have critical roles in CRC pathogenesis and warrant future investigation.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1796-4) contains supplementary material, which is available to authorized users.

On estimation of time-dependent attributable fraction from population-based case-control studies

Population attributable fraction (PAF) is widely used to quantify the disease burden associated with a modifiable exposure in a population. It has been extended to a time-varying measure that provides additional information on when and how the exposure's impact varies over time for cohort studies. However, there is no estimation procedure for PAF using data that are collected from population-based case-control studies, which, because of time and cost efficiency, are commonly used for studying genetic and environmental risk factors of disease incidences. In this article, we show that time-varying PAF is identifiable from a case-control study and develop a novel estimator of PAF. Our estimator combines odds ratio estimates from logistic regression models and density estimates of the risk factor distribution conditional on failure times in cases from a kernel smoother. The proposed estimator is shown to be consistent and asymptotically normal with asymptotic variance that can be estimated empirically from the data. Simulation studies demonstrate that the proposed estimator performs well in finite sample sizes. Finally, the method is illustrated by a population-based case-control study of colorectal cancer.

Quantifying the Genetic Correlation between Multiple Cancer Types

Background: Many cancers share specific genetic risk factors, including both rare high-penetrance mutations and common SNPs identified through genome-wide association studies (GWAS). However, little is known about the overall shared heritability across cancers. Quantifying the extent to which two distinct cancers share genetic origin will give insights to shared biological mechanisms underlying cancer and inform design for future genetic association studies.Methods: In this study, we estimated the pair-wise genetic correlation between six cancer types (breast, colorectal, lung, ovarian, pancreatic, and prostate) using cancer-specific GWAS summary statistics data based on 66,958 case and 70,665 control subjects of European ancestry. We also estimated genetic correlations between cancers and 14 noncancer diseases and traits.Results: After adjusting for 15 pair-wise genetic correlation tests between cancers, we found significant (P < 0.003) genetic correlations between pancreatic and colorectal cancer (rg = 0.55, P = 0.003), lung and colorectal cancer (rg = 0.31, P = 0.001). We also found suggestive genetic correlations between lung and breast cancer (rg = 0.27, P = 0.009), and colorectal and breast cancer (rg = 0.22, P = 0.01). In contrast, we found no evidence that prostate cancer shared an appreciable proportion of heritability with other cancers. After adjusting for 84 tests studying genetic correlations between cancer types and other traits (Bonferroni-corrected P value: 0.0006), only the genetic correlation between lung cancer and smoking remained significant (rg = 0.41, P = 1.03 × 10-6). We also observed nominally significant genetic correlations between body mass index and all cancers except ovarian cancer.Conclusions: Our results highlight novel genetic correlations and lend support to previous observational studies that have observed links between cancers and risk factors.Impact: This study demonstrates modest genetic correlations between cancers; in particular, breast, colorectal, and lung cancer share some degree of genetic basis. Cancer Epidemiol Biomarkers Prev; 26(9); 1427-35. ©2017 AACR.

Complex genetic control of lung tumorigenesis in resistant mice strains

The SM/J mouse strain is resistant to chemically‐induced lung tumorigenesis despite having a haplotype, in the pulmonary adenoma susceptibility locus (Pas1) locus, that confers tumor susceptibility in other strains. To clarify this inconsistent genotype‐phenotype correlation, we crossed SM/J mice with another resistant strain and conducted genome‐wide linkage analysis in the (C57BL/6J × SM/J)F2 progeny exposed to urethane to induce lung tumors. Overall, >80% of F2 mice of both sexes developed from 1 to 20 lung tumors. Genotyping of 372 F2 mice for 744 informative non‐redundant SNPs dispersed over all autosomal chromosomes revealed four quantitative trait loci (QTLs) affecting lung tumor multiplicity, on chromosomes 3 (near rs13477379), 15 (rs6285067), 17 (rs33373629) and 18 (rs3706601), all with logarithm of the odds (LOD) scores >5. Four QTLs modulated total lung tumor volume, on chromosome 3 (rs13477379), 10 (rs13480702), 15 (rs6285067) and 17 (rs3682923), all with LOD scores >4. No QTL modulating lung tumor multiplicity or total volume was detected in Pas1 on chromosome 6. The present study demonstrates that the SM/J strain carries, at the Pas1 locus, the resistance allele: a finding that will facilitate identification of the Pas1 causal element. More generally, it demonstrates that lung tumorigenesis is under complex polygenic control even in a pedigree with low susceptibility to this neoplasia, suggesting that the genetics of lung tumorigenesis is much more complex than evidenced by the pulmonary adenoma susceptibility and resistance loci that have, so far, been mapped in a small number of crosses between a few inbred strains.

Inclusion of a Genetic Risk Score into a Validated Risk Prediction Model for Colorectal Cancer in Japanese Men Improves Performance

We previously developed and validated a risk prediction model for colorectal cancer in Japanese men using modifiable risk factors. To further improve risk prediction, we evaluated the degree of improvement obtained by adding a genetic risk score (GRS) using genome-wide association study (GWAS)-identified risk variants to our validated model. We examined the association between 36 risk variants identified by GWAS and colorectal cancer risk using a weighted Cox proportional hazards model in a nested case-control study within the Japan Public Health Center-based Prospective Study. GRS was constructed using six variants associated with risk in this study of the 36 tested. We assessed three models: a nongenetic model that included the same variables used in our previously validated model; a genetic model that used GRS; and an inclusive model, which included both. The c-statistic, integrated discrimination improvement (IDI), and net reclassification improvement (NRI) were calculated by the 5-fold cross-validation method. We estimated 10-year absolute risks for developing colorectal cancer. A statistically significant association was observed between the weighted GRS and colorectal cancer risk. The mean c-statistic for the inclusive model (0.66) was slightly greater than that for the nongenetic model (0.60). Similarly, the mean IDI and NRI showed improvement when comparing the nongenetic and inclusive models. These models for colorectal cancer were well calibrated. The addition of GRS using GWAS-identified risk variants to our validated model for Japanese men improved the prediction of colorectal cancer risk. Cancer Prev Res; 10(9); 535-41. ©2017 AACR.

Optimal sequencing strategies for identifying disease-associated singletons

With the increasing focus of genetic association on the identification of trait-associated rare variants through sequencing, it is important to identify the most cost-effective sequencing strategies for these studies. Deep sequencing will accurately detect and genotype the most rare variants per individual, but may limit sample size. Low pass sequencing will miss some variants in each individual but has been shown to provide a cost-effective alternative for studies of common variants. Here, we investigate the impact of sequencing depth on studies of rare variants, focusing on singletons—the variants that are sampled in a single individual and are hardest to detect at low sequencing depths. We first estimate the sensitivity to detect singleton variants in both simulated data and in down-sampled deep genome and exome sequence data. We then explore the power of association studies comparing burden of singleton variants in cases and controls under a variety of conditions. We show that the power to detect singletons increases with coverage, typically plateauing for coverage > ~25x. Next, we show that, when total sequencing capacity is fixed, the power of association studies focused on singletons is typically maximized for coverage of 15-20x, independent of relative risk, disease prevalence, singleton burden, and case-control ratio. Our results suggest sequencing depth of 15-20x as an appropriate compromise of singleton detection power and sample size for studies of rare variants in complex disease.

Genetic studies of rare variants can help us understand the biology of human disease. With modern techniques and sufficient effort, it is possible to very accurately resolve any human genome, identifying most of its unique features. When funding is limited, applying these techniques to study human disease often involves a trade-off between examining more samples, at reduced accuracy per sample, or fewer samples, each at greater accuracy. We evaluate these trade-offs for studies of very rare variants, using both simulation and real data. We propose cost effective strategies for increasing our understanding of human disease.

Genome-Wide Association Studies of Cancer in Diverse Populations

Genome-wide association studies (GWAS) of cancer have identified more than 700 risk loci, of which approximately 80% were first discovered in European ancestry populations, approximately 15% in East Asians, 3% in multiethnic scans, and less than 1% in African and Latin American populations. These percentages closely mirror the distribution of samples included in the discovery phase of cancer GWAS to date (84% European, 11% East Asian, 4% African, and 1% Latin American ancestry). GWAS in non-European ancestry populations have provided insight into ancestry-specific variation in cancer and have pointed to regions of susceptibility that are of particular importance in certain populations. Uncovering and characterizing cancer risk loci in diverse populations is critical for understanding underlying biological mechanisms and developing future genetic risk prediction models in non-European ancestry populations. New GWAS and continued collaborations will be required to eliminate population inequalities in the number of studies, sample sizes, and variant content on GWAS arrays, and to better align genetic research in cancer to the global distribution of race/ethnicity Cancer Epidemiol Biomarkers Prev; 27(4); 405-17. ©2018 AACRSee all articles in this CEBP Focus section, "Genome-Wide Association Studies in Cancer."

Genetic Variants in Epigenetic Pathways and Risks of Multiple Cancers in the GAME-ON Consortium

Background: Epigenetic disturbances are crucial in cancer initiation, potentially with pleiotropic effects, and may be influenced by the genetic background.Methods: In a subsets (ASSET) meta-analytic approach, we investigated associations of genetic variants related to epigenetic mechanisms with risks of breast, lung, colorectal, ovarian and prostate carcinomas using 51,724 cases and 52,001 controls. False discovery rate-corrected P values (q values < 0.05) were considered statistically significant.Results: Among 162,887 imputed or genotyped variants in 555 candidate genes, SNPs in eight genes were associated with risk of more than one cancer type. For example, variants in BABAM1 were confirmed as a susceptibility locus for squamous cell lung, overall breast, estrogen receptor (ER)-negative breast, and overall prostate, and overall serous ovarian cancer; the most significant variant was rs4808076 [OR = 1.14; 95% confidence interval (CI) = 1.10-1.19; q = 6.87 × 10-5]. DPF1 rs12611084 was inversely associated with ER-negative breast, endometrioid ovarian, and overall and aggressive prostate cancer risk (OR = 0.93; 95% CI = 0.91-0.96; q = 0.005). Variants in L3MBTL3 were associated with colorectal, overall breast, ER-negative breast, clear cell ovarian, and overall and aggressive prostate cancer risk (e.g., rs9388766: OR = 1.06; 95% CI = 1.03-1.08; q = 0.02). Variants in TET2 were significantly associated with overall breast, overall prostate, overall ovarian, and endometrioid ovarian cancer risk, with rs62331150 showing bidirectional effects. Analyses of subpathways did not reveal gene subsets that contributed disproportionately to susceptibility.Conclusions: Functional and correlative studies are now needed to elucidate the potential links between germline genotype, epigenetic function, and cancer etiology.Impact: This approach provides novel insight into possible pleiotropic effects of genes involved in epigenetic processes. Cancer Epidemiol Biomarkers Prev; 26(6); 816-25. ©2017 AACR.

In silico pathway analysis and tissue specific cis-eQTL for colorectal cancer GWAS risk variants

Background

Genome-wide association studies have identified 55 genetic variants associated with colorectal cancer risk to date. However, potential causal genes and pathways regulated by these risk variants remain to be characterized. Therefore, we performed gene ontology enrichment and pathway analyses to determine if there was an enrichment of genes in proximity to the colorectal cancer risk variants that could further elucidate the probable causal genes and pathways involved in colorectal cancer biology.

Results

For the 65 unique genes that either contained, or were immediately neighboring up- and downstream, of these variants there was a significant enrichment for the KEGG pathway, Pathways in Cancer (p-value = 2.67 × 10⁻⁵) and an enrichment for multiple biological processes (FDR < 0.05), such as cell junction organization, tissue morphogenesis, regulation of SMAD protein phosphorylation, and odontogenesis identified through Gene Ontology analysis. To identify potential causal genes, we conducted a cis-expression quantitative trait loci (cis-eQTL) analysis using gene expression and genotype data from the Genotype-Tissue Expression (GTEx) Project portal in normal sigmoid (n = 124) and transverse (n = 169) colon tissue. In addition, we also did a cis-eQTL analysis on colorectal tumor tissue (n = 147) from The Cancer Genome Atlas (TCGA). We identified two risk alleles that were significant cis-eQTLs for FADS2 (rs1535) and COLCA1 and 2 (rs3802842) genes in the normal transverse colon tissue and two risk alleles that were significant cis-eQTLs for the CABLES2 (rs2427308) and LIPG (rs7229639) genes in the normal sigmoid colon tissue, but not tumor tissue.

Conclusions

Our data reaffirm the potential to identify an enrichment for biological processes and candidate causal genes based on expression profiles correlated with genetic risk alleles of colorectal cancer, however, the identification of these significant cis-eQTLs is context and tissue specific.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3750-2) contains supplementary material, which is available to authorized users.

MicroRNAs: Potential candidates for diagnosis and treatment of colorectal cancer

Colorectal cancer (CRC) is known as the third common cancer worldwide and an important public health problem in different populations. Several genetics and environmental risk factors are involved in the development and progression of CRC including chromosomal abnormalities, epigenetic alterations, and unhealthy lifestyle. Identification of risk factors and biomarkers could lead to a better understanding of molecular pathways involved in CRC pathogenesis. MicroRNAs (miRNAs) are important regulatory molecules which could affect a variety of cellular and molecular targets in CRC. A large number of studies have indicated deregulations of some known tissue-specific miRNAs, for example, miR-21, miR-9, miR-155, miR-17, miR-19, let-7, and miR-24 as well as circulating miRNAs, for example, miR-181b, miR-21, miR-183, let-7g, miR-17, and miR-126, in patients with CRC. In the current review, we focus on the findings of preclinical and clinical studies performed on tissue-specific and circulating miRNAs as diagnostic biomarkers and therapeutic targets for the detection of patients at various stages of CRC.

Precancer Atlas to Drive Precision Prevention Trials

Cancer development is a complex process driven by inherited and acquired molecular and cellular alterations. Prevention is the holy grail of cancer elimination, but making this a reality will take a fundamental rethinking and deep understanding of premalignant biology. In this Perspective, we propose a national concerted effort to create a Precancer Atlas (PCA), integrating multi-omics and immunity - basic tenets of the neoplastic process. The biology of neoplasia caused by germline mutations has led to paradigm-changing precision prevention efforts, including: tumor testing for mismatch repair (MMR) deficiency in Lynch syndrome establishing a new paradigm, combinatorial chemoprevention efficacy in familial adenomatous polyposis (FAP), signal of benefit from imaging-based early detection research in high-germline risk for pancreatic neoplasia, elucidating early ontogeny in BRCA1-mutation carriers leading to an international breast cancer prevention trial, and insights into the intricate germline-somatic-immunity interaction landscape. Emerging genetic and pharmacologic (metformin) disruption of mitochondrial (mt) respiration increased autophagy to prevent cancer in a Li-Fraumeni mouse model (biology reproduced in clinical pilot) and revealed profound influences of subtle changes in mt DNA background variation on obesity, aging, and cancer risk. The elaborate communication between the immune system and neoplasia includes an increasingly complex cellular microenvironment and dynamic interactions between host genetics, environmental factors, and microbes in shaping the immune response. Cancer vaccines are in early murine and clinical precancer studies, building on the recent successes of immunotherapy and HPV vaccine immune prevention. Molecular monitoring in Barrett's esophagus to avoid overdiagnosis/treatment highlights an important PCA theme. Next generation sequencing (NGS) discovered age-related clonal hematopoiesis of indeterminate potential (CHIP). Ultra-deep NGS reports over the past year have redefined the premalignant landscape remarkably identifying tiny clones in the blood of up to 95% of women in their 50s, suggesting that potentially premalignant clones are ubiquitous. Similar data from eyelid skin and peritoneal and uterine lavage fluid provide unprecedented opportunities to dissect the earliest phases of stem/progenitor clonal (and microenvironment) evolution/diversity with new single-cell and liquid biopsy technologies. Cancer mutational signatures reflect exogenous or endogenous processes imprinted over time in precursors. Accelerating the prevention of cancer will require a large-scale, longitudinal effort, leveraging diverse disciplines (from genetics, biochemistry, and immunology to mathematics, computational biology, and engineering), initiatives, technologies, and models in developing an integrated multi-omics and immunity PCA - an immense national resource to interrogate, target, and intercept events that drive oncogenesis. Cancer Res; 77(7); 1510-41. ©2017 AACR.