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.

Abstract PR03: Genetic mechanisms of immune evasion in colorectal cancer

Background: Immune checkpoint blockade has shown activity in microsatellite-instability high (MSI-H) colorectal cancers (CRCs). However, despite a very high mutational and neoantigen load among virtually all MSI-H tumors, the response rate is around 40-50%. In addition, for the majority of CRCs, which are microsatellite stable (MSS), immune checkpoint inhibitors have so far proven ineffective. Thus, to better understand the genetic drivers of immune evasion in CRC, we integrated next generation sequencing data from over 1200 tumors with transcriptional and immunohistochemical measures of immune infiltration.

Methods: We molecularly characterized 1,211 colorectal cancers, including 592 tumors from The Cancer Genome Atlas with whole exome sequencing (WES) and whole transcriptome (RNAseq) data and 619 cancers from two prospective cohort studies with WES and immunohistochemical (IHC) annotations. To identify driver gene alterations and their selection pressures specific to MSI-high tumors, we developed a statistical method to identify significantly mutated microsatellite tracts and we further developed a method to identify copy neutral loss of heterozygosity (CN-LOH) events. We used an established immune-gene transcriptional signature as well as IHC stains against specific subsets of immune-infiltrating cells to identify genetic events associated with immune evasion.

Results: We demonstrated that WNT-signaling and immune-related genes are significantly mutated in colorectal cancer. In MSI-H CRCs, we found biallelic antigen-presentation machinery (APM) mutations in the HLA Class I genes, B2M and TAP2, in addition to recurrent mutations in NLRC5 and RFX5, which downregulated HLA Class I expression. In all CRCs, we showed that WNT-signaling activity and APC-biallelic mutations are inversely associated with both transcriptional and IHC measures of T-cell infiltration. Specifically, nuclear CTNNB1 expression was inversely correlated with overall tumor-infiltrating lymphocytes (p = 0.027), CD8+ subset of T-cells (p = 0.0019) and CD45RO+ subset of T-cells (0.0080). Meanwhile, colorectal tumors with biallelic disruptive mutations in APC had a significantly decreased activated T-cell transcriptional signature (p = 4e-12) relative to samples with no disruptive mutations in APC. We further showed that in MSS tumors, AXIN2 (a key WNT-signaling effector) super-enhancer hypomethylation, independent of the APC mutation status, was further associated with decreased T-cell activity.

Conclusions: In this largest CRC genomic analysis to date, we identify genetic events that are associated with immune evasion in this disease. Specifically, we find evidence of immuno-editing in MSI-H tumors through disruptive mutations in APM. In the MSS and MSI-H subtype of CRCs, we use transcriptional and immunohistochemical orthogonal analyses to demonstrate exclusion of an effective immune infiltrate in CRC through an active WNT-signaling pathway. Our results shed light on the underlying molecular mechanisms of immune exclusion in CRC and have direct implications for novel combination immunotherapy trials for patients with this disease.

This abstract is also being presented as Poster B23.

Citation Format: Marios Giannakis, Catherine Grasso, Daniel Wells, Tsuyoshi Hamada, Xinmeng Jasmine Mu, Michael Quist, Jonathan Nowak, Reiko Nishihara, Charles M. Connolly, Sachet Shukla, William M. Grady, David Wheeler, Catherine J. Wu, Jesse Zaretsky, Levi Garraway, Thomas Hudson, Charles Fuchs, Antoni Ribas, Riki Peters, Shuji Ogino. Genetic mechanisms of immune evasion in colorectal cancer [abstract]. In: Proceedings of the AACR Special Conference on Tumor Immunology and Immunotherapy; 2017 Oct 1-4; Boston, MA. Philadelphia (PA): AACR; Cancer Immunol Res 2018;6(9 Suppl):Abstract nr PR03.

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.

RNA polymerase III interferes with Ty3 integration

Ty3, a gypsylike retrotransposon of budding yeast, integrates at the transcription initiation site of genes transcribed by RNA polymerase III (pol III). It was previously shown that integration in vitro requires intact promoter elements and the pol III transcription factors TFIIIB and TFIIIC. In order to test the effect of pol III on integration, increasing amounts of a pol III-containing fraction were added to Ty3 in vitro integration reactions. The pol III-containing fraction was inhibitory to integration. These results are consistent with a model where the Ty3 integration complex and pol III recognize similar features of the stable transcription complex and compete with each other for access to the transcription initiation site.

Requirement of RNA polymerase III transcription factors for in vitro position-specific integration of a retroviruslike element

The yeast retroviruslike element Ty3 inserts at the transcription initiation sites of genes transcribed by RNA polymerase III (Pol III). An in vitro integration assay was developed with the use of Ty3 viruslike particles and a modified SUP2 tyrosine transfer RNA (tRNA(Tyr)) gene target. Integration was position-specific and required Ty3 integrase, Pol III transcription factor (TF) IIIB-, TFIIIC-, and Pol III-containing fractions showed that TFIIIB and TFIIIC, together, were sufficient for position-specific Ty3 integration, but not for transcription. This report demonstrates that in vitro integration of a retroelement can be targeted by cellular proteins.

Selective decontamination with nystatin for control of a Candida outbreak in a neonatal intensive care unit

Selective decontamination of the digestive tract (SDD) with oral nystatin was evaluated as a measure to control an outbreak of Candida infection in a neonatal intensive care unit (NICU). Seventy-six out of 106 neonates who carried Candida spp. received the main study manoeuvre (the application of oral nystatin in the throat and stomach) during the 12-month open trial. One third of the neonates weighed < 1500 g whilst about half were being ventilated. The mean stay was 33.2 d (SD +/- 46.9). Two cases with candidaemia within a fortnight were associated with a yeast carriage rate in the NICU of about 50%; more than 80% of the isolates were Candida parapsilosis. During the implementation period there were four new neonates with fungaemia caused by C. parapsilosis. Once the carriage rate dropped below 5% (P < 0.001), no new cases of systemic infection with the outbreak strain were recognized in the following 8 months. It took 3.5 months to control the outbreak. The observation that all other clinical diagnostic samples were free from Candida suggests that translocation from throat or gut into the systemic circulation occurred. SDD with oral nystatin was effective in reducing the yeast carriage index (mean index 1.93, before SDD; 0.45, after SDD; P < 0.001). A significant reduction of carriage, both in rates and indices, is thought to have contributed to the control of this candida outbreak.

Sequence of pedal peptide: a novel neuropeptide from the central nervous system of Aplysia

We report the identification of a novel neuropeptide from Aplysia nervous tissue. The peptide was termed Pedal peptide (Pep) because it was predominantly synthesized in the pedal ganglia. Pep was purified and sequenced from pooled extracts of pedal ganglia. The following sequence was proposed: Pro-Leu-Asp-Ser-Val-Tyr-Gly-Thr-His-Gly-Met-Ser-Gly-Phe-Ala. Enzymatic hydrolysis procedures indicated that Pep had a free carboxyl terminal. A peptide with the proposed sequence was synthesized and compared with the native peptide. Chromatographic properties of the 2 peptides under 3 different conditioned were compared and found to be identical. Electrophysiological responses to the 2 peptides were compared on an identified neuron in the abdominal ganglia and found to be qualitatively and quantitatively very similar. Both peptides produced net inward currents that were associated with a decrease in membrane conductance. The results from these 2 procedures confirmed that the proposed Pep sequence was correct. Quantitative measurements of the incorporation of 35S-methionine into Pep suggest that cell bodies that synthesize Pep were present predominantly in the pedal ganglia but should also be found in other central ganglia as well. Pep-like immunoreactive neurons are found predominantly in the pedal ganglia and less frequently in the other ganglia (Pearson and Lloyd, 1989). Quantitatively, Pep constitutes one of the predominant peptides in the nervous system of Aplysia. Pep does not appear to be a member of any other previously identified invertebrate or vertebrate peptide family.