CpG Island Methylation in Colorectal Cancer: Past, Present and Future

The concept of a CpG island methylator phenotype, or CIMP, quickly became the focus of several colorectal cancer studies describing its clinical and pathological features after its introduction in 1999 by Toyota and colleagues. Further characterization of CIMP in tumors lead to widespread acceptance of the concept, as expressed by Shen and Issa in their 2005 editorial, “CIMP, at last.” Since that time, extensive research efforts have brought great insights into the epidemiology and prognosis of CIMP+ tumors and other epigenetic mechanisms underlying tumorigenesis. With the advances in technology and subsequent cataloging of the human methylome in cancer and normal tissue, new directions in research to understand CIMP and its role in complex biological systems yield hope for future epigenetically based diagnostics and treatments

Power in Pairs: Assessing the Statistical Value of Paired Samples in Tests for Differential Expression

Background: When genomics researchers design a high-throughput study to test for differential expression, some biological systems and research questions provide opportunities to use paired samples from subjects, and researchers can plan for a certain proportion of subjects to have paired samples. We consider the effect of this paired samples proportion on the statistical power of the study, using characteristics of both count (RNA-Seq) and continuous (microarray) expression data from a colorectal cancer study. Results: We demonstrate that a higher proportion of subjects with paired samples yields higher statistical power, for various total numbers of samples, and for various strengths of subject-level confounding factors. In the design scenarios considered, the statistical power in a fully-paired design is substantially (and in many cases several times) greater than in an unpaired design. Conclusions: For the many biological systems and research questions where paired samples are feasible and relevant, substantial statistical power gains can be achieved at the study design stage when genomics researchers plan on using paired samples from the largest possible proportion of subjects. Any cost savings in a study design with unpaired samples are likely accompanied by underpowered and possibly biased results

Changing population characteristics, effect-measure modification, and cancer risk factor identification

Epidemiologic studies have identified a number of lifestyle factors, e.g. diet, obesity, and use of certain medications, which affect risk of colon cancer. However, the magnitude and significance of risk factor-disease associations differ among studies. We propose that population trends of changing prevalence of risk factors explains some of the variability between studies when factors that change prevalence also modify the effect of other risk factors. We used data collected from population-based control who were selected as study participants for two time periods, 1991–1994 and 1997–2000, along with data from the literature, to examine changes in the population prevalence of aspirin and non-steroidal anti-inflammatory medication (NSAID) use, obesity, and hormone replacement therapy (HRT) over time. Data from a population-based colon cancer case-control study were used to estimate effect-measurement modification among these factors. Sizeable changes in aspirin use, HRT use, and the proportion of the population that is obese were observed between the 1980s and 2000. Use of NSAIDs interacted with BMI and HRT; HRT use interacted with body mass index (BMI). We estimate that as the prevalence of NSAIDs use changed from 10% to almost 50%, the colon cancer relative risk associated with BMI >30 would change from 1.3 to 1.9 because of the modifying effect of NSAIDs. Similarly, the relative risk estimated for BMI would increase as the prevalence of use of HRT among post-menopausal women increased. In conclusion, as population characteristics change over time, these changes may have an influence on relative risk estimates for colon cancer for other exposures because of effect-measure modification. The impact of population changes on comparability between epidemiologic studies can be kept to a minimum if investigators assess exposure-disease associations within strata of other exposures, and present results in a manner that allows comparisons across studies. Effect-measure modification is an important component of data analysis that should be evaluated to obtain a complete understanding of disease etiology

Genetic Variation in Selenoprotein Genes, Lifestyle, and Risk of Colon and Rectal Cancer

BACKGROUND: Associations between selenium and cancer have directed attention to role of selenoproteins in the carcinogenic process. METHODS: We used data from two population-based case-control studies of colon (n = 1555 cases, 1956 controls) and rectal (n = 754 cases, 959 controls) cancer. We evaluated the association between genetic variation in TXNRD1, TXNRD2, TXNRD3, C11orf31 (SelH), SelW, SelN1, SelS, SepX, and SeP15 with colorectal cancer risk. RESULTS: After adjustment for multiple comparisons, several associations were observed. Two SNPs in TXNRD3 were associated with rectal cancer (rs11718498 dominant OR 1.42 95% CI 1.16,1.74 pACT 0.0036 and rs9637365 recessive 0.70 95% CI 0.55,0.90 pACT 0.0208). Four SNPs in SepN1 were associated with rectal cancer (rs11247735 recessive OR 1.30 95% CI 1.04,1.63 pACT 0.0410; rs2072749 GGvsAA OR 0.53 95% CI 0.36,0.80 pACT 0.0159; rs4659382 recessive OR 0.58 95% CI 0.39,0.86 pACT 0.0247; rs718391 dominant OR 0.76 95% CI 0.62,0.94 pACT 0.0300). Interaction between these genes and exposures that could influence these genes showed numerous significant associations after adjustment for multiple comparisons. Two SNPs in TXNRD1 and four SNPs in TXNRD2 interacted with aspirin/NSAID to influence colon cancer; one SNP in TXNRD1, two SNPs in TXNRD2, and one SNP in TXNRD3 interacted with aspirin/NSAIDs to influence rectal cancer. Five SNPs in TXNRD2 and one in SelS, SeP15, and SelW1 interacted with estrogen to modify colon cancer risk; one SNP in SelW1 interacted with estrogen to alter rectal cancer risk. Several SNPs in this candidate pathway influenced survival after diagnosis with colon cancer (SeP15 and SepX1 increased HRR) and rectal cancer (SepX1 increased HRR). CONCLUSIONS: Findings support an association between selenoprotein genes and colon and rectal cancer development and survival after diagnosis. Given the interactions observed, it is likely that the impact of cancer susceptibility from genotype is modified by lifestyle

Exploring multilocus associations of inflammation genes and colorectal cancer risk using hapConstructor

<p>Abstract</p> <p>Background</p> <p>In candidate-gene association studies of single nucleotide polymorphisms (SNPs), multilocus analyses are frequently of high dimensionality when considering haplotypes or haplotype pairs (diplotypes) and differing modes of expression. Often, while candidate genes are selected based on their biological involvement in a given pathway, little is known about the functionality of SNPs to guide association studies. Investigators face the challenge of exploring multiple SNP models to elucidate which variants, independently or in combination, might be associated with a disease of interest. A data mining module, hapConstructor (freely-available in Genie software) performs systematic construction and association testing of multilocus genotype data in a Monte Carlo framework. Our objective was to assess its utility to guide statistical analyses of haplotypes within a candidate region (or combined genotypes across candidate genes) beyond that offered by a standard logistic regression approach.</p> <p>Methods</p> <p>We applied the hapConstructor method to a multilocus investigation of candidate genes involved in pro-inflammatory cytokine IL6 production, <it>IKBKB</it>, <it>IL6</it>, and <it>NFKB1 </it>(16 SNPs total) hypothesized to operate together to alter colorectal cancer risk. Data come from two U.S. multicenter studies, one of colon cancer (1,556 cases and 1,956 matched controls) and one of rectal cancer (754 cases and 959 matched controls).</p> <p>Results</p> <p>HapConstrcutor enabled us to identify important associations that were further analyzed in logistic regression models to simultaneously adjust for confounders. The most significant finding (nominal <it>P </it>= 0.0004; false discovery rate <it>q </it>= 0.037) was a combined genotype association across <it>IKBKB </it>SNP rs5029748 (1 or 2 variant alleles), <it>IL6 </it>rs1800797 (1 or 2 variant alleles), and <it>NFKB1 </it>rs4648110 (2 variant alleles) which conferred an ~80% decreased risk of colon cancer.</p> <p>Conclusions</p> <p>Strengths of hapConstructor were: systematic identification of multiple loci within and across genes important in CRC risk; false discovery rate assessment; and efficient guidance of subsequent logistic regression analyses.</p

IκBKβ and NFκB1 , NSAID use and risk of colorectal cancer in the Colon Cancer Family Registry

The NFκB-signaling pathway regulates cell proliferation and inflammation. Activation of the pathway is implicated in the etiology of colorectal cancer (CRC). NSAIDs may reduce CRC risk partially through a nuclear factor-kappa B (NFκB)-dependent pathway. In this study, we investigated associations between 34 NFκB1 and 8 IκBKβ tagSNPs and CRC risk and examined interactions with non-steroidal anti-inflammatory drug (NSAID) use. Using conditional logistic regression, we investigated these associations among 1584 incident CRC cases and 2516 sibling controls from the Colon Cancer Family Registry. Three IκBKβ SNPs were associated with a statistically significant lower colorectal or colon cancer risk: rs9694958 (A>G intron 5) (colorectal: ORhzv = 0.26(0.07–0.99), Ptrend = 0.048, Padj = 0.25), rs10958713 (A>C intron 19) (colon: ORhzv = 0.62(0.42–0.92), Ptrend = 0.005, Padj = 0.03) and rs5029748 (C>A intron 2) (colon: ORhet = 0.72(0.56–0.91), Ptrend = 0.01, Padj = 0.08). We replicated trends associated with NFκB1 and IκBKβ variants identified in a previous study (rs4648110 (T>A intron 22), rs13117745 (G>A intron 5) and rs3747811 (T>A intron 1)). IκBKβ’s rs6474387 (C>T intron 20) and rs11986055 (A>C intron 2) showed substantially lower colon cancer risk among current NSAID users (Pinteraction = 0.01 and Pinteraction = 0.045, respectively), whereas NFκB1’s rs230490 (G>A 5ʹ (outside UTR)) and rs997476 (C>A 3ʹ (outside UTR)) showed higher CRC risk among current NSAID users (Pinteraction = 0.01 and Pinteraction = 0.03, respectively). These findings suggest that variants in NFκB1 and IκBKβ are associated with CRC risk and NSAIDs may function partially through an NFκB-dependent pathway. The SNPs identified here should be considered for future functional studies and may be useful in designing a pharmacogenetic approach to preventive NSAID use

Postmenopausal Hormone Therapy and Colorectal Cancer Risk by Molecularly Defined Subtypes and Tumor Location

Background: Postmenopausal hormone therapy (HT) is associated with a decreased colorectal cancer (CRC) risk. As CRC is a heterogeneous disease, we evaluated whether the association of HT and CRC differs across etiologically relevant, molecularly defined tumor subtypes and tumor location. Methods: We pooled data on tumor subtypes (microsatellite instability status, CpG island methylator phenotype status, BRAF and KRAS mutations, pathway: adenoma-carcinoma, alternate, serrated), tumor location (proximal colon, distal colon, rectum), and HT use among 8220 postmenopausal women (3898 CRC cases and 4322 controls) from 8 observational studies. We used multinomial logistic regression to estimate odds ratios (OR) and 95% confidence intervals (CIs) for the association of ever vs never HT use with each tumor subtype compared with controls. Models were adjusted for study, age, body mass index, smoking status, and CRC family history. All statistical tests were 2-sided. Results: Among postmenopausal women, ever HT use was associated with a 38% reduction in overall CRC risk (OR =0.62, 95% CI = 0.56 to 0.69). This association was similar according to microsatellite instability, CpG island methylator phenotype and BRAF or KRAS status. However, the association was attenuated for tumors arising through the serrated pathway (OR = 0.81, 95% CI = 0.66 to 1.01) compared with the adenoma-carcinoma pathway (OR = 0.63, 95% CI = 0.55 to 0.73; P het =.04) and alternate pathway (OR = 0.61, 95% CI = 0.51 to 0.72). Additionally, proximal colon tumors had a weaker association (OR = 0.71, 95% CI = 0.62 to 0.80) compared with rectal (OR = 0.54, 95% CI = 0.46 to 0.63) and distal colon (OR = 0.57, 95% CI = 0.49 to 0.66; P het =.01) tumors. Conclusions: We observed a strong inverse association between HT use and overall CRC risk, which may predominantly reflect a benefit of HT use for tumors arising through the adenoma-carcinoma and alternate pathways as well as distal colon and rectal tumors

Association between germline variants and somatic mutations in colorectal cancer

Colorectal cancer (CRC) is a heterogeneous disease with evidence of distinct tumor types that develop through different somatically altered pathways. To better understand the impact of the host genome on somatically mutated genes and pathways, we assessed associations of germline variations with somatic events via two complementary approaches. We first analyzed the association between individual germline genetic variants and the presence of non-silent somatic mutations in genes in 1375 CRC cases with genome-wide SNPs data and a tumor sequencing panel targeting 205 genes. In the second analysis, we tested if germline variants located within previously identified regions of somatic allelic imbalance were associated with overall CRC risk using summary statistics from a recent large scale GWAS (n similar or equal to 125 k CRC cases and controls). The first analysis revealed that a variant (rs78963230) located within a CNA region associated with TLR3 was also associated with a non-silent mutation within gene FBXW7. In the secondary analysis, the variant rs2302274 located in CDX1/PDGFRB frequently gained/lost in colorectal tumors was associated with overall CRC risk (OR = 0.96, p = 7.50e-7). In summary, we demonstrate that an integrative analysis of somatic and germline variation can lead to new insights about CRC

Association Between Smoking and Molecular Subtypes of Colorectal Cancer

Background: Smoking is associated with colorectal cancer (CRC) risk. Previous studies suggested this association may be restricted to certain molecular subtypes of CRC, but large-scale comprehensive analysis is lacking. Methods: A total of 9789 CRC cases and 11 231 controls of European ancestry from 11 observational studies were included. We harmonized smoking variables across studies and derived sex study-specific quartiles of pack-years of smoking for analysis. Four somatic colorectal tumor markers were assessed individually and in combination, including BRAF mutation, KRAS mutation, CpG island methylator phenotype (CIMP), and microsatellite instability (MSI) status. A multinomial logistic regression analysis was used to assess the association between smoking and risk of CRC subtypes by molecular characteristics, adjusting for age, sex, and study. All statistical tests were 2-sided and adjusted for Bonferroni correction. Results: Heavier smoking was associated with higher risk of CRC overall and stratified by individual markers (P-trend <.001). The associations differed statistically significantly between all molecular subtypes, which was the most statistically significant for CIMP and BRAF. Compared with never-smokers, smokers in the fourth quartile of pack-years had a 90% higher risk of CIMP-positive CRC (odds ratio = 1.90, 95% confidence interval = 1.60 to 2.26) but only 35% higher risk for CIMP-negative CRC (odds ratio = 1.35, 95% confidence interval = 1.22 to 1.49; P-difference = 2.1 x 10(-6)). The association was also stronger in tumors that were CIMP positive, MSI high, or KRAS wild type when combined (P-difference <.001). Conclusion: Smoking was associated with differential risk of CRC subtypes defined by molecular characteristics. Heavier smokers had particularly higher risk of CRC subtypes that were CIMP positive and MSI high in combination, suggesting that smoking may be involved in the development of colorectal tumors via the serrated pathway