SYSMut: decoding the functional significance of rare somatic mutations in cancer.

Current tailored-therapy efforts in cancer are largely focused on a small number of highly recurrently mutated driver genes but therapeutic targeting of these oncogenes remains challenging. However, the vast number of genes mutated infrequently across cancers has received less attention, in part, due to a lack of understanding of their biological significance. We present SYSMut, an extendable systems biology platform that can robustly infer the biologic consequences of somatic mutations by integrating routine multiomics profiles in primary tumors. We establish SYSMut's improved performance vis-à-vis state-of-the-art driver gene identification methodologies by recapitulating the functional impact of known driver genes, while additionally identifying novel functionally impactful mutated genes across 29 cancers. Subsequent application of SYSMut on low-frequency gene mutations in head and neck squamous cell (HNSC) cancers, followed by molecular and pharmacogenetic validation, revealed the lipidogenic network as a novel therapeutic vulnerability in aggressive HNSC cancers. SYSMut is thus a robust scalable framework that enables the discovery of new targetable avenues in cancer.

The Ephrin B2 Receptor Tyrosine Kinase Is a Regulator of Proto-oncogene MYC and Molecular Programs Central to Barrett's Neoplasia.

BACKGROUND & AIMS: Mechanisms contributing to the onset and progression of Barrett's (BE)-associated esophageal adenocarcinoma (EAC) remain elusive. Here, we interrogated the major signaling pathways deregulated early in the development of Barrett's neoplasia. METHODS: Whole-transcriptome RNA sequencing analysis was performed in primary BE, EAC, normal esophageal squamous, and gastric biopsy tissues (n = 89). Select pathway components were confirmed by quantitative polymerase chain reaction in an independent cohort of premalignant and malignant biopsy tissues (n = 885). Functional impact of selected pathway was interrogated using transcriptomic, proteomic, and pharmacogenetic analyses in mammalian esophageal organotypic and patient-derived BE/EAC cell line models, in vitro and/or in vivo. RESULTS: The vast majority of primary BE/EAC tissues and cell line models showed hyperactivation of EphB2 signaling. Transcriptomic/proteomic analyses identified EphB2 as an endogenous binding partner of MYC binding protein 2, and an upstream regulator of c-MYC. Knockdown of EphB2 significantly impeded the viability/proliferation of EAC and BE cells in vitro/in vivo. Activation of EphB2 in normal esophageal squamous 3-dimensional organotypes disrupted epithelial maturation and promoted columnar differentiation programs, notably including MYC. EphB2 and MYC showed selective induction in esophageal submucosal glands with acinar ductal metaplasia, and in a porcine model of BE-like esophageal submucosal gland spheroids. Clinically approved inhibitors of MEK, a protein kinase that regulates MYC, effectively suppressed EAC tumor growth in vivo. CONCLUSIONS: The EphB2 signaling is frequently hyperactivated across the BE-EAC continuum. EphB2 is an upstream regulator of MYC, and activation of EphB2-MYC axis likely precedes BE development. Targeting EphB2/MYC could be a promising therapeutic strategy for this often refractory and aggressive cancer.

Patients With Esophageal Adenocarcinoma With Prior Gastroesophageal Reflux Disease Symptoms Are Similar to Those Without Gastroesophageal Reflux Disease: A Cross-Sectional Study.

INTRODUCTION: A substantial proportion of patients with esophageal adenocarcinoma (EAC) do not report gastroesophageal reflux disease (GERD) symptoms. This study aimed to compare the risk factor profiles and cancer stage at presentation of patients with EAC with and without prior GERD. METHODS: In this retrospective cross-sectional study, patients with EAC were divided into 2 cohorts: (i) EAC with prior GERD: patients who reported typical GERD symptoms (heartburn or regurgitation) ≥1 year before cancer diagnosis and (ii) EAC without prior GERD: patients who did not report prior GERD symptoms or reported symptoms within 1 year of their cancer diagnosis. Baseline demographics, risk factors, and cancer stage at presentation were compared between the 2 cohorts. In addition, the distribution of patients based on numbers of BE/EAC-associated risk factors (1, 2, 3, 4, and 5 or more) was examined in the symptomatic and asymptomatic cohorts. RESULTS: Over 13 years, 388 patients with EAC with prior GERD and 245 patients with EAC without prior GERD were recruited. Both groups had similar baseline demographics and risk factors, but patients with EAC with prior GERD were more likely to have a history of BE. Asymptomatic patients had more advanced disease. Patients with 3 or more BE/EAC-related risk factors formed the largest proportion of patients in both the symptomatic and asymptomatic cohorts. DISCUSSION: Patients with EAC with and without prior GERD symptoms are phenotypically similar, suggesting that BE screening efforts to prevent or detect early EAC should not be restricted to just those with GERD.

Systems Biology Analyses Show Hyperactivation of Transforming Growth Factor-ß and JNK Signaling Pathways in Esophageal Cancer.

BACKGROUND & AIMS: Esophageal adenocarcinoma (EAC) is resistant to standard chemoradiation treatments, and few targeted therapies are available. We used large-scale tissue profiling and pharmacogenetic analyses to identify deregulated signaling pathways in EAC tissues that might be targeted to slow tumor growth or progression. METHODS: We collected 397 biopsy specimens from patients with EAC and nonmalignant Barrett's esophagus (BE), with or without dysplasia. We performed RNA-sequencing analyses and used systems biology approaches to identify pathways that are differentially activated in EAC vs nonmalignant dysplastic tissues; pathway activities were confirmed with immunohistochemistry and quantitative real-time polymerase chain reaction analyses of signaling components in patient tissue samples. Human EAC (FLO-1 and EsoAd1), dysplastic BE (CP-B, CP-C, CP-D), and nondysplastic BE (CP-A) cells were incubated with pharmacologic inhibitors or transfected with small interfering RNAs. We measured effects on proliferation, colony formation, migration, and/or growth of xenograft tumors in nude mice. RESULTS: Comparisons of EAC vs nondysplastic BE tissues showed hyperactivation of transforming growth factor-ß (TGFB) and/or Jun N-terminal kinase (JNK) signaling pathways in more than 80% of EAC samples. Immunohistochemical analyses showed increased nuclear localization of phosphorylated JUN and SMAD proteins in EAC tumor tissues compared with nonmalignant tissues. Genes regulated by the TGFB and JNK pathway were overexpressed specifically in EAC and dysplastic BE. Pharmacologic inhibition or knockdown of TGFB or JNK signaling components in EAC cells (FLO-1 or EsoAd1) significantly reduced cell proliferation, colony formation, cell migration, and/or growth of xenograft tumors in mice in a SMAD4-independent manner. Inhibition of the TGFB pathway in BE cell lines reduced the proliferation of dysplastic, but not nondysplastic, cells. CONCLUSIONS: In a transcriptome analysis of EAC and nondysplastic BE tissues, we found the TGFB and JNK signaling pathways to be hyperactivated in EACs and the genes regulated by these pathways to be overexpressed in EAC and dysplastic BE. Inhibiting these pathways in EAC cells reduces their proliferation, migration, and formation of xenograft tumors. Strategies to block the TGFB and JNK signaling pathways might be developed for treatment of EAC.

Evidence for GALNT12 as a moderate penetrance gene for colorectal cancer.

Characterizing moderate penetrance susceptibility genes is an emerging frontier in colorectal cancer (CRC) research. GALNT12 is a strong candidate CRC-susceptibility gene given previous linkage and association studies, and inactivating somatic and germline alleles in CRC patients. Previously, we found rare segregating germline GALNT12 variants in a clinic-based cohort (N = 118) with predisposition for CRC. Here, we screened a new population-based cohort of incident CRC cases (N = 479) for rare (MAF ≤1%) deleterious germline GALNT12 variants. GALNT12 screening revealed eight rare variants. Two variants were previously described (p.Asp303Asn, p.Arg297Trp), and additionally, we found six other rare variants: five missense (p.His101Gln, p.Ile142Thr, p.Glu239Gln, p.Thr286Met, p.Val290Phe) and one putative splice-altering variant (c.732-8 G>T). Sequencing of population-matched controls (N = 400) revealed higher burden of these variants in CRC cases compared with healthy controls (P = 0.0381). We then functionally characterized the impact of substitutions on GALNT12 enzyme activity using in vitro-derived peptide substrates. Three of the newly identified GALNT12 missense variants (p.His101Gln, p.Ile142Thr, p.Val290Phe) demonstrated a marked loss (>2-fold reduction) of enzymatic activity compared with wild-type (P ≤ 0.05), whereas p.Glu239Gln exhibited a ∼2-fold reduction in activity (P = 0.077). These findings provide strong, independent evidence for the association of GALNT12 defects with CRC-susceptibility; underscoring implications for glycosylation pathway defects in CRC.

Novel recurrently mutated genes in African American colon cancers.

We used whole-exome and targeted sequencing to characterize somatic mutations in 103 colorectal cancers (CRC) from African Americans, identifying 20 new genes as significantly mutated in CRC. Resequencing 129 Caucasian derived CRCs confirmed a 15-gene set as a preferential target for mutations in African American CRCs. Two predominant genes, ephrin type A receptor 6 (EPHA6) and folliculin (FLCN), with mutations exclusive to African American CRCs, are by genetic and biological criteria highly likely African American CRC driver genes. These previously unsuspected differences in the mutational landscapes of CRCs arising among individuals of different ethnicities have potential to impact on broader disparities in cancer behaviors.

Identification and functional characterization of paxillin as a target of protein tyrosine phosphatase receptor T.

Protein tyrosine phosphatase receptor-type T (PTPRT) is the most frequently mutated tyrosine phosphatase in human cancers. However, the cell signaling pathways regulated by PTPRT largely remain to be elucidated. Here, we show that paxillin is a direct substrate of PTPRT and that PTPRT specifically regulates paxillin phosphorylation at tyrosine residue 88 (Y88) in colorectal cancer (CRC) cells. We engineered CRC cells homozygous for a paxillin Y88F knock-in mutant and found that these cells exhibit significantly reduced cell migration and impaired anchorage-independent growth, fail to form xenograft tumors in nude mice, and have decreased phosphorylation of p130CAS, SHP2, and AKT. PTPRT knockout mice that we generated exhibit increased levels of colonic paxillin phosphorylation at residue Y88 and are highly susceptible to carcinogen azoxymethane-induced colon tumor, providing critical in vivo evidence that PTPRT normally functions as a tumor suppressor. Moreover, similarly increased paxillin pY88 is also found as a common feature of human colon cancers. These studies reveal an important signaling pathway that plays a critical role in colorectal tumorigenesis.

The microRNAs, MiR-31 and MiR-375, as candidate markers in Barrett's esophageal carcinogenesis.

There is a critical need to identify molecular markers that can reliably aid in stratifying esophageal adenocarcinoma (EAC) risk in patients with Barrett's esophagus. MicroRNAs (miRNA/miR) are one such class of biomolecules. In the present cross-sectional study, we characterized miRNA alterations in progressive stages of neoplastic development, i.e., metaplasia-dysplasia-adenocarcinoma, with an aim to identify candidate miRNAs potentially associated with progression. Using next generation sequencing (NGS) as an agnostic discovery platform, followed by quantitative real-time PCR (qPCR) validation in a total of 20 EACs, we identified 26 miRNAs that are highly and frequently deregulated in EACs (≥ 4-fold in >50% of cases) when compared to paired normal esophageal squamous (nSQ) tissue. We then assessed the 26 EAC-derived miRNAs in laser microdissected biopsy pairs of Barrett's metaplasia (BM)/nSQ (n = 15), and high-grade dysplasia (HGD)/nSQ (n = 14) by qPCR, to map the timing of deregulation during progression from BM to HGD and to EAC. We found that 23 of the 26 candidate miRNAs were deregulated at the earliest step, BM, and therefore noninformative as molecular markers of progression. Two miRNAs, miR-31 and -31*, however, showed frequent downregulation only in HGD and EAC cases suggesting association with transition from BM to HGD. A third miRNA, miR-375, showed marked downregulation exclusively in EACs and in none of the BM or HGD lesions, suggesting its association with progression to invasive carcinoma. Taken together, we propose miR-31 and -375 as novel candidate microRNAs specifically associated with early- and late-stage malignant progression, respectively, in Barrett's esophagus.

Global mutational profiling of formalin-fixed human colon cancers from a pathology archive.

The advent of Next-Generation sequencing technologies, which significantly increases the throughput and reduces the cost of large-scale sequencing efforts, provides an unprecedented opportunity for discovery of novel gene mutations in human cancers. However, it remains a challenge to apply Next-Generation technologies to DNA extracted from formalin-fixed paraffin-embedded cancer specimens. We describe here the successful development of a custom DNA capture method using Next-Generation for detection of 140 driver genes in five formalin-fixed paraffin-embedded human colon cancer samples using an improved extraction process to produce high-quality DNA. Isolated DNA was enriched for targeted exons and sequenced using the Illumina Next-Generation platform. An analytical pipeline using 3 software platforms to define single-nucleotide variants was used to evaluate the data output. Approximately 250 × average coverage was obtained with >96% of target bases having at least 30 sequence reads. Results were then compared with previously performed high-throughput Sanger sequencing. Using an algorithm of needing a positive call from all three callers to give a positive result, 98% of the verified Sanger sequencing somatic driver gene mutations were identified by our method with a specificity of 90%. In all, 13 insertions and deletions identified by Next-Generation were confirmed by Sanger sequencing. We also applied this technology to two components of a biphasic colon cancer, which had strikingly differing histology. Remarkably, no new driver gene mutation accumulation was identified in the more undifferentiated component. Applying this method to profiling of formalin-fixed paraffin-embedded colon cancer tissue samples yields equivalent sensitivity and specificity for mutation detection as Sanger sequencing of matched cell lines derived from these cancers. This method directly enables high-throughput comprehensive mutational profiling of colon cancer samples, and is easily extendable to enable targeted sequencing from formalin-fixed paraffin-embedded material for other tumor types.

Inactivating germ-line and somatic mutations in polypeptide N-acetylgalactosaminyltransferase 12 in human colon cancers.

Aberrant glycosylation is a pathological alteration that is widespread in colon cancer, and usually accompanies the onset and progression of the disease. To date, the molecular mechanisms underlying aberrant glycosylation remain largely unknown. In this study, we identify somatic and germ-line mutations in the gene encoding for polypeptide N-acetylgalactosaminyltransferase 12 (GALNT12) in individuals with colon cancer. Biochemical analyses demonstrate that each of the 8 GALNT12 mutations identified inactivates the normal function of the GALNT enzyme in initiating mucin type O-linked protein glycosylation. Two of these inactivating GALNT12 mutations were identified as acquired somatic mutations in a set of 30 microsatellite stable colon tumors. Relative to background gene mutation rates, finding these somatic GALNT12 mutations was statistically significant at P < 0.001. Six additional inactivating GALNT12 mutations were detected as germ-line changes carried by patients with colon cancer; however, no inactivating variants were detected among cancer-free controls (P = 0.005). Notably, in 3 of the 6 individuals harboring inactivating germ-line GALNT12 mutations, both a colon cancer and a second independent epithelial cancer had developed. These findings suggest that genetic defects in the O-glycosylation pathway in part underlie aberrant glycosylation in colon cancers, and they contribute to the development of a subset of these malignancies.

Allele-specific expression of TGFBR1 in colon cancer patients.

The genetic component of colorectal cancer (CRC) predisposition has been only partially explained. We recently suggested that a subtle decrease in the expression of one allele of the TGFBR1 gene was a heritable quantitative trait predisposing to CRC. Here, we refined the measurements of allele-specific expression (ASE) of TGFBR1 in a population-based series of CRC patients and controls. Five single-nucleotide polymorphisms (SNPs) in the 3'-untranslated region of the gene were genotyped and used for ASE determination by pyrosequencing. After eliminating non-informative samples and samples with RNA of insufficient quality 109 cases and 125 controls were studied. Allelic ratios ranged between 0.74 and 1.69 without evidence of bimodality or cutoff points for 'ASE' versus 'non-ASE'. Treating ASE as a continuous variable, cases had non-significantly different values than controls (P = 0.081 when comparing means by permutation test). However, cases had significantly higher ASE values when comparing medians by permutation test (P = 0.0027) and when using Wilcoxon test (P = 0.0094). We conclude that with the present-day technology, ASE differences between individuals and between cases and controls are too subtle to be used to assess CRC risk. More advanced technology is expected to resolve this issue as well as the low informativity caused by the limited heterozygosity of transcribed SNPs.

The noncoding RNA, miR-126, suppresses the growth of neoplastic cells by targeting phosphatidylinositol 3-kinase signaling and is frequently lost in colon cancers.

MicroRNAs (miRNA/miR) are a class of small noncoding RNAs implicated in the pathogenesis of various malignancies. In the current study, using micro(RNA) arrays, we found a ubiquitous loss of miR-126 expression in colon cancer lines when compared to normal human colon epithelia. Reconstitution of miR-126 in colon cancer cells resulted in a significant growth reduction as evidenced in clonogenic assays. A search for miR-126 gene targets revealed p85beta, a regulatory subunit involved in stabilizing and propagating the phosphatidylinositol 3-kinase (PI3K) signal, as one of the potential substrates. Restoration of miR-126 in cancer cells induced a > or =3-fold reduction in p85beta protein levels, with no concomitant change in p85alpha, a gene that is functionally related to p85beta but not a supposed target of miR-126. Additionally, using reporter constructs, we show that the p85beta-3' untranslated region is directly targeted by miR-126. Furthermore, this miR-126 mediated reduction of p85beta was accompanied by a substantial reduction in phosphorylated AKT levels in the cancer cells, suggesting an impairment in PI3K signaling. Finally, in a panel of matched normal colon and primary colon tumors, each of the tumors demonstrated miR-126 down-regulation together with an increase in the p85beta protein level. Taken together, we propose that miR-126 regulates PI3K signaling partly by targeting p85beta, and that the loss of miR-126 may provide a selective growth advantage during colon carcinogenesis.

Identifying DNA methylation biomarkers for non-endoscopic detection of Barrett's esophagus.

We report a biomarker-based non-endoscopic method for detecting Barrett's esophagus (BE) based on detecting methylated DNAs retrieved via a swallowable balloon-based esophageal sampling device. BE is the precursor of, and a major recognized risk factor for, developing esophageal adenocarcinoma. Endoscopy, the current standard for BE detection, is not cost-effective for population screening. We performed genome-wide screening to ascertain regions targeted for recurrent aberrant cytosine methylation in BE, identifying high-frequency methylation within the CCNA1 locus. We tested CCNA1 DNA methylation as a BE biomarker in cytology brushings of the distal esophagus from 173 individuals with or without BE. CCNA1 DNA methylation demonstrated an area under the curve of 0.95 for discriminating BE-related metaplasia and neoplasia cases versus normal individuals, performing identically to methylation of VIM DNA, an established BE biomarker. When combined, the resulting two biomarker panel was 95% sensitive and 91% specific. These results were replicated in an independent validation cohort of 149 individuals who were assayed using the same cutoff values for test positivity established in the training population. To progress toward non-endoscopic esophageal screening, we engineered a well-tolerated, swallowable, encapsulated balloon device able to selectively sample the distal esophagus within 5 min. In balloon samples from 86 individuals, tests of CCNA1 plus VIM DNA methylation detected BE metaplasia with 90.3% sensitivity and 91.7% specificity. Combining the balloon sampling device with molecular assays of CCNA1 plus VIM DNA methylation enables an efficient, well-tolerated, sensitive, and specific method of screening at-risk populations for BE.

Strain-specific homeostatic responses during early stages of Azoxymethane-induced colon tumorigenesis in mice.

Identifying molecular changes that predict the risk for developing colon cancer is critical for designing effective prevention strategies. In the present study, we determined early-stage molecular alterations within the colonic epithelium of A/J and AKR/J mice that are sensitive and resistant to Azoxymethane (AOM)-initiated tumor development, respectively. Six week-old male mice were injected intraperitoneally with AOM (10 mg/kg body weight) once a week for six weeks. One week after the last injection, distal colons from both strains were analyzed for cell proliferation using a proliferating cell nuclear antigen (PCNA) assay. Unlike AKR/J, a significant increase (2.5-fold, p<0.05) in the number of PCNA-positive cells within the upper third of the crypt compartment was observed in the A/J colons. This proliferative response was associated with a sizeable increase in the levels of c-myc mRNA, quantified by RNase protection assay. cDNA sequencing, protein expression and localization of beta-catenin, an upstream activator of c-myc, however, showed no aberrant changes within AOM-exposed A/J colons. Interestingly, TdT-mediated dUTP nick-end labeling assay revealed a significant increase (4-fold) in the number of apoptotic colonocytes in A/J mice following AOM treatment. Consistent with this finding, a modest increase in the expression of pro-apoptotic Bak was limited to the sensitive A/J colons. In summary, the current study suggests that a significant alteration in the rate of cell turnover in the normal appearing colonic mucosa, as observed in susceptible A/J mice, may be one of the earliest events predisposing the colon to neoplastic growth.

Genomic regions associated with susceptibility to Barrett's esophagus and esophageal adenocarcinoma in African Americans: The cross BETRNet admixture study.

BACKGROUND: Barrett's esophagus (BE) and esophageal adenocarcinoma (EAC) are far more prevalent in European Americans than in African Americans. Hypothesizing that this racial disparity in prevalence might represent a genetic susceptibility, we used an admixture mapping approach to interrogate disease association with genomic differences between European and African ancestry. METHODS: Formalin fixed paraffin embedded samples were identified from 54 African Americans with BE or EAC through review of surgical pathology databases at participating Barrett's Esophagus Translational Research Network (BETRNet) institutions. DNA was extracted from normal tissue, and genotyped on the Illumina OmniQuad SNP chip. Case-only admixture mapping analysis was performed on the data from both all 54 cases and also on a subset of 28 cases with high genotyping quality. Haplotype phases were inferred with Beagle 3.3.2, and local African and European ancestries were inferred with SABER plus. Disease association was tested by estimating and testing excess European ancestry and contrasting it to excess African ancestry. RESULTS: Both datasets, the 54 cases and the 28 cases, identified two admixture regions. An association of excess European ancestry on chromosome 11p reached a 5% genome-wide significance threshold, corresponding to -log10(P) = 4.28. A second peak on chromosome 8q reached -log10(P) = 2.73. The converse analysis examining excess African ancestry found no genetic regions with significant excess African ancestry associated with BE and EAC. On average, the regions on chromosomes 8q and 11p showed excess European ancestry of 15% and 20%, respectively. CONCLUSIONS: Chromosomal regions on 11p15 and 8q22-24 are associated with excess European ancestry in African Americans with BE and EAC. Because GWAS have not reported any variants in these two regions, low frequency and/or rare disease associated variants that confer susceptibility to developing BE and EAC may be driving the observed European ancestry association evidence.

Role of the alternating reading frame (P19)-p53 pathway in an in vivo murine colon tumor model.

Considering the importance of the oncogene checkpoint function of the alternating reading frame(ARF)-p53 pathway, studies were undertaken to evaluate the status of this pathway in azoxymethane (AOM)-induced mouse colon tumors. A PCR-based analysis of ARF and p53 cDNAs in normal colon tissues and AOM-induced colon tumors failed to detect mutations in either of these two critical tumor suppressor genes. In addition, laser capture microdissection of tumors followed by PCR-based sequencing of exons 5-7 of genomic p53 showed that even the most pleomorphic cancer cells were p53 normal. A marked increase in ARF mRNA and protein levels was observed in colon tumors, indicating activation of the ARF-p53 pathway in these tumors. High levels of ARF protein stabilized p53 protein in the tumors, but the p53 protein showed little biochemical activity. Compared with a mouse colonocyte cell line that expresses high levels of wild-type p53 (YAMC), the p53 protein in tumors had no detectable DNA binding activity nor did it activate p21 expression. In fact, p21 levels were lower in tumor tissue relative to normal mucosa, even though p53 levels were approximately 30-fold higher in tumors relative to control. Within the A/J tumors, we also used a cDNA microarray approach to screen a panel of genes that are transcriptionally up- or down-regulated by functional p53. The expression patterns of these p53-regulated genes were consistent with a lack of functional p53. This work demonstrates that the ARF-p53 oncogene checkpoint can be overcome without p53 mutations and that the mechanism used to overcome this checkpoint involves the suppression of p53 transcriptional activating activity. The AOM colon cancer model may be well suited for studying tumor promotion events that precede p53 disruption.