Landscape of somatic single nucleotide variants and indels in colorectal cancer and impact on survival.

Colorectal cancer (CRC) is a biologically heterogeneous disease. To characterize its mutational profile, we conduct targeted sequencing of 205 genes for 2,105 CRC cases with survival data. Our data shows several findings in addition to enhancing the existing knowledge of CRC. We identify PRKCI, SPZ1, MUTYH, MAP2K4, FETUB, and TGFBR2 as additional genes significantly mutated in CRC. We find that among hypermutated tumors, an increased mutation burden is associated with improved CRC-specific survival (HR = 0.42, 95% CI: 0.21-0.82). Mutations in TP53 are associated with poorer CRC-specific survival, which is most pronounced in cases carrying TP53 mutations with predicted 0% transcriptional activity (HR = 1.53, 95% CI: 1.21-1.94). Furthermore, we observe differences in mutational frequency of several genes and pathways by tumor location, stage, and sex. Overall, this large study provides deep insights into somatic mutations in CRC, and their potential relationships with survival and tumor features.

Correction to: Genetic variant predictors of gene expression provide new insight into risk of colorectal cancer.

Every author has erroneously been assigned to the affiliation "62". The affiliation 62 belongs to the author Graham Casey.

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.

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.

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.

Transcriptional profiling of androgen receptor (AR) mutants suggests instructive and permissive roles of AR signaling in germ cell development.

The androgen receptor (AR) is a transcription factor that plays a critical role in male sexual development, spermatogenesis, and maintenance of hormonal homeostasis. Despite the extensive knowledge of the phenotypic consequences of mutations in Ar, very little is known about the transcriptional targets of AR within the testis. To identify potential targets of androgen signaling in the testis, we have analyzed the transcriptional profile of adult testes from Ar hypomorphs alone or in combination with Sertoli cell-specific Ar ablation. Using Affymetrix MOE430A mouse genome arrays we interrogated more than 22,000 transcripts. We found the expression level of 62 transcripts in the Ar mutants differed by greater than 2-fold compared with wild type. We also found that more transcripts were up-regulated than down-regulated, highlighting AR's role as a transcriptional repressor in the testis. Twelve transcripts were uniquely affected, and 16 transcripts were more severely affected in Sertoli cell-specific Ar ablation compared with hypomorphic Ar mutants. Using a comparative genomic approach, we analyzed the 6 kb around the transcriptional start sites of affected transcripts for conserved AREs (androgen response elements). We identified at least one conserved ARE in 65% of the genes misregulated in our microarray analysis where clear mouse-human orthologs were available. We used a reporter assay in cell culture to functionally verify the AREs for the kallikrein 27 gene. This suggests that the majority of the misregulated transcripts have a high probability of being direct AR targets. The transcripts affected by these Ar mutations encode a diverse array of proteins whose molecular functions support the contention that AR supports spermatogenesis in both a permissive and instructive fashion.

Disruption of murine Tenr results in teratospermia and male infertility.

Gametes rely heavily on posttranscriptional control for their differentiation. Translational control, alternative splicing, and alternative processing of the 3' end of mRNAs are all common during spermatogenesis. Tenr, which encodes a highly conserved 72-kDa protein, is expressed solely in germ cells of the testis from the mid-pachytene stage until the elongating spermatid stage. TENR contains a double-stranded RNA binding domain, is localized to the nucleus, and is phylogenetically related to a family of adenosine deaminases involved in RNA editing. We show here that targeted mutation of the Tenr gene causes male sterility. Tenr mutant males have a reduced sperm count, and Tenr-/- sperm show a decrease in motility and an increase in malformed heads. Despite their sterility, some epididymal sperm from Tenr mutants have normal morphology. The ability of Tenr mutant sperm to fertilize zona pellucida-free oocytes and to bind to, but not fertilize, zona pellucida-intact oocytes suggests that the normal-appearing sperm are not able to penetrate the zona pellucida. These data suggest that TENR plays an essential function in spermatid morphogenesis.