Molecular dissection of CRC primary tumors and their matched liver metastases reveals critical role of immune microenvironment, EMT and angiogenesis in cancer metastasis.

Metastasis is the primary cause of cancer mortality. The primary tumors of colorectal cancer (CRC) often metastasize to the liver. In this study, we have collected 122 samples from 45 CRC patients. Among them, 32 patients have primary tumors, adjacent normal tissues, and matched liver metastases. Thirteen patients have primary tumors without distant metastasis and matched normal tissues. Characterization of these samples was conducted by whole-exome and RNA sequencing and SNP6.0 analysis. Our results revealed no significant difference in genetic alterations including common oncogenic mutations, whole genome mutations and copy number variations between primary and metastatic tumors. We then assembled gene co-expression networks and identified metastasis-correlated gene networks of immune-suppression, epithelial-mesenchymal transition (EMT) and angiogenesis as the key events and potentially synergistic drivers associated with CRC metastasis. Further independent cohort validation using published datasets has verified that these specific gene networks are up regulated throughout the tumor progression. The gene networks of EMT, angiogenesis, immune-suppression and T cell exhaustion are closely correlated with the poor patient outcome and intrinsic anti-PD-1 resistance. These results offer insights of combinational strategy for the treatment of metastatic CRC.

Genomic characterization of intrinsic and acquired resistance to cetuximab in colorectal cancer patients.

Anti-EGFR antibodies are effective in therapies for late-stage colorectal cancer (CRC); however, many tumours are unresponsive or develop resistance. We performed genomic analysis of intrinsic and acquired resistance to anti-EGFR therapy in prospectively collected tumour samples from 25 CRC patients receiving cetuximab (an EGFR inhibitor). Of 25 CRC patients, 13 displayed intrinsic resistance to cetuximab; 12 were intrinsically sensitive. We obtained six re-biopsy samples at acquired resistance from the intrinsically sensitive patients. NCOA4-RET and LMNA-NTRK1 fusions and NRG1 and GNAS amplifications were found in intrinsic-resistant patients. In cetuximab-sensitive patients, we found KRAS K117N and A146T mutations in addition to BRAF V600E, AKT1 E17K, PIK3CA E542K, and FGFR1 or ERBB2 amplifications. The comparison between baseline and acquired-resistant tumours revealed an extreme shift in variant allele frequency of somatic variants, suggesting that cetuximab exposure dramatically selected for rare resistant subclones that were initially undetectable. There was also an increase in epithelial-to-mesenchymal transition at acquired resistance, with a reduction in the immune infiltrate. Furthermore, characterization of an acquired-resistant, patient-derived cell line showed that PI3K/mTOR inhibition could rescue cetuximab resistance. Thus, we uncovered novel genomic alterations that elucidate the mechanisms of sensitivity and resistance to anti-EGFR therapy in metastatic CRC patients.

Mouse PDX Trial Suggests Synergy of Concurrent Inhibition of RAF and EGFR in Colorectal Cancer with BRAF or KRAS Mutations.

Purpose: To evaluate the antitumor efficacy of cetuximab in combination with LSN3074753, an analog of LY3009120 and pan-RAF inhibitor in 79 colorectal cancer patient-derived xenograft (PDX) models.Experimental Design: Seventy-nine well-characterized colorectal cancer PDX models were employed to conduct a single mouse per treatment group (n = 1) trial.Results: Consistent with clinical results, cetuximab was efficacious in wild-type KRAS and BRAF PDX models, with an overall response rate of 6.3% and disease control rate (DCR) of 20.3%. LSN3074753 was active in a small subset of PDX models that harbored KRAS or BRAF mutations. However, the combination treatment displayed the enhanced antitumor activity with DCR of 35.4%. Statistical analysis revealed that BRAF and KRAS mutations were the best predictors of the combinatorial activity and were significantly associated with synergistic effect with a P value of 0.01 compared with cetuximab alone. In 12 models with BRAF mutations, the combination therapy resulted in a DCR of 41.7%, whereas either monotherapy had a DCR of 8.3%. Among 44 KRAS mutation models, cetuximab or LSN3074753 monotherapy resulted in a DCR of 13.6% or 11.4%, respectively, and the combination therapy increased DCR to 34.1%. Molecular analysis suggests that EGFR activation is a potential feedback and resistant mechanism of pan-RAF inhibition.Conclusions: MAPK and EGFR pathway activations are two major molecular hallmarks of colorectal cancer. This mouse PDX trial recapitulated clinical results of cetuximab. Concurrent EGFR and RAF inhibition demonstrated synergistic antitumor activity for colorectal cancer PDX models with a KRAS or BRAF mutation. Clin Cancer Res; 23(18); 5547-60. ©2017 AACR.

Comprehensive Characterization of Oncogenic Drivers in Asian Lung Adenocarcinoma.

INTRODUCTION: The incidence rate of lung adenocarcinoma (LUAD), the predominant histological subtype of lung cancer, is elevated in Asians, particularly in female nonsmokers. The mutation patterns in LUAD in Asians might be distinct from those in LUAD in whites. METHODS: We profiled 271 resected LUAD tumors (mainly stage I) to characterize the genomic landscape of LUAD in Asians with a focus on female nonsmokers. RESULTS: Mutations in EGFR, KRAS, erb-b2 receptor tyrosine kinase 2 gene (ERBB2), and BRAF; gene fusions involving anaplastic lymphoma receptor tyrosine kinase gene (ALK), ROS1, and ret proto-oncogene (RET); and Met Proto-Oncogene Tyrosine Kinase (MET) exon 14 skipping were the major drivers in LUAD in Asians, exhibiting mutually exclusive and differing prevalence from those reported in studies of LUAD in non-Asians. In addition, we identified a novel mutational signature of XNX (the mutated base N in the middle flanked by two identical bases at the 5' and 3' positions) that was overrepresented in LUAD tumors in nonsmokers and negatively correlated with the overall mutational frequency. CONCLUSIONS: In this cohort, approximately 85% of individuals have known driver mutations (EGFR 59.4%, KRAS 7.4%, ALK 7.4%, ERBB2 2.6%, ROS1 2.2%, RET 2.2%, MET 1.8%, BRAF 1.1%, and NRAS 0.4%). Seventy percent of smokers and 90% of nonsmokers had defined oncogenic drivers matching the U.S. Food and Drug Administration-approved targeted therapies.

Molecular analysis of gastric cancer identifies subtypes associated with distinct clinical outcomes.

Gastric cancer, a leading cause of cancer-related deaths, is a heterogeneous disease. We aim to establish clinically relevant molecular subtypes that would encompass this heterogeneity and provide useful clinical information. We use gene expression data to describe four molecular subtypes linked to distinct patterns of molecular alterations, disease progression and prognosis. The mesenchymal-like type includes diffuse-subtype tumors with the worst prognosis, the tendency to occur at an earlier age and the highest recurrence frequency (63%) of the four subtypes. Microsatellite-unstable tumors are hyper-mutated intestinal-subtype tumors occurring in the antrum; these have the best overall prognosis and the lowest frequency of recurrence (22%) of the four subtypes. The tumor protein 53 (TP53)-active and TP53-inactive types include patients with intermediate prognosis and recurrence rates (with respect to the other two subtypes), with the TP53-active group showing better prognosis. We describe key molecular alterations in each of the four subtypes using targeted sequencing and genome-wide copy number microarrays. We validate these subtypes in independent cohorts in order to provide a consistent and unified framework for further clinical and preclinical translational research.

Genomic landscape and genetic heterogeneity in gastric adenocarcinoma revealed by whole-genome sequencing.

Gastric cancer (GC) is the second most common cause of cancer-related deaths. It is known to be a heterogeneous disease with several molecular and histological subtypes. Here we perform whole-genome sequencing of 49 GCs with diffuse (N=31) and intestinal (N=18) histological subtypes and identify three mutational signatures, impacting TpT, CpG and TpCp[A/T] nucleotides. The diffuse-type GCs show significantly lower clonality and smaller numbers of somatic and structural variants compared with intestinal subtype. We further divide the diffuse subtype into one with infrequent genetic changes/low clonality and another with relatively higher clonality and mutations impacting TpT dinucleotide. Notably, we discover frequent and exclusive mutations in Ephrins and SLIT/ROBO signalling pathway genes. Overall, this study delivers new insights into the mutational heterogeneity underlying distinct histologic subtypes of GC that could have important implications for future research in the diagnosis and treatment of GC.

Whole-genome sequencing of asian lung cancers: second-hand smoke unlikely to be responsible for higher incidence of lung cancer among Asian never-smokers.

Asian nonsmoking populations have a higher incidence of lung cancer compared with their European counterparts. There is a long-standing hypothesis that the increase of lung cancer in Asian never-smokers is due to environmental factors such as second-hand smoke. We analyzed whole-genome sequencing of 30 Asian lung cancers. Unsupervised clustering of mutational signatures separated the patients into two categories of either all the never-smokers or all the smokers or ex-smokers. In addition, nearly one third of the ex-smokers and smokers classified with the never-smoker-like cluster. The somatic variant profiles of Asian lung cancers were similar to that of European origin with G.C>T.A being predominant in smokers. We found EGFR and TP53 to be the most frequently mutated genes with mutations in 50% and 27% of individuals, respectively. Among the 16 never-smokers, 69% had an EGFR mutation compared with 29% of 14 smokers/ex-smokers. Asian never-smokers had lung cancer signatures distinct from the smoker signature and their mutation profiles were similar to European never-smokers. The profiles of Asian and European smokers are also similar. Taken together, these results suggested that the same mutational mechanisms underlie the etiology for both ethnic groups. Thus, the high incidence of lung cancer in Asian never-smokers seems unlikely to be due to second-hand smoke or other carcinogens that cause oxidative DNA damage, implying that routine EGFR testing is warranted in the Asian population regardless of smoking status.

Locations and patterns of meiotic recombination in two-generation pedigrees.

BACKGROUND: Meiotic crossovers are the major mechanism by which haplotypes are shuffled to generate genetic diversity. Previously available methods for the genome-wide, high-resolution identification of meiotic crossover sites are limited by the laborious nature of the assay (as in sperm typing). METHODS: Several methods have been introduced to identify crossovers using high density single nucleotide polymorphism (SNP) array technologies, although programs are not widely available to implement such analyses. RESULTS: Here we present a two-generation "reverse pedigree analysis" method (analyzing the genotypes of two children relative to each parent) and a web-accessible tool to determine and visualize inheritance differences among siblings and crossover locations on each parental gamete. This approach is complementary to existing methods and uses informative markers which provide high resolution for locating meiotic crossover sites. We introduce a segmentation algorithm to identify crossover sites, and used a synthetic data set to determine that the segmentation algorithm specificity was 92% and sensitivity was 89%. The use of reverse pedigrees allows the inference of crossover locations on the X chromosome in a maternal gamete through analysis of two sons and their father. We further analyzed genotypes from eight multiplex autism families, observing a 1.462 maternal to paternal recombination ratio and no significant differences between affected and unaffected children. Meiotic recombination results from pediSNP can also be used to identify haplotypes that are shared by probands within a pedigree, as we demonstrated with a multiplex autism family. CONCLUSION: Using "reverse pedigrees" and defining unique sets of genotype markers within pedigree data, we introduce a method that identifies inherited allelic differences and meiotic crossovers. We implemented the method in the pediSNP software program, and we applied it to several data sets. This approach uses data from two generations to identify crossover sites, facilitating studies of recombination in disease. pediSNP is available online at http://pevsnerlab.kennedykrieger.org/pediSNP.

Visualization of uniparental inheritance, Mendelian inconsistencies, deletions, and parent of origin effects in single nucleotide polymorphism trio data with SNPtrio.

A variety of alterations occur in chromosomal DNA, many of which can be detected using high density single nucleotide polymorphism (SNP) microarrays. These include deletions and duplications (assessed by observing changes in copy number) and regions of homozygosity. The analysis of SNP data from trios can provide an additional category of information about the nature and origin of inheritance patterns, including uniparental disomy (UPD), loss of transmitted allele (LTA), and nonparental relationship. The main purpose of SNPtrio is to locate regions of uniparental inheritance (UPI) and Mendelian inconsistency (MI), identify the type (paternal vs. maternal, iso- vs. hetero-), and assess the associated statistical probability of occurrence by chance. SNPtrio's schema permits the identification of hemizygous or homozygous deletions as well as UPD. We validated the performance of SNPtrio on three platforms (Affymetrix 10 K and 100 K arrays and Illumina 550 K arrays) using SNP data obtained from DNA samples of patients known to have UPD and diagnosed with Prader-Willi syndrome, Angelman syndrome, Beckwith-Wiedemann syndrome, pseudohypoparathyroidism, and a complex chromosome 2 abnormality. We further validated SNPtrio using DNA from patients previously shown to have microdeletions that were verified by fluorescence in situ hybridization (FISH). SNPtrio successfully identified previously known UPD and deletion regions, and generated associated probability values. SNPtrio analysis of trisomy 21 (Down syndrome) cases and their parents permitted identification of the parent of origin of the extra chromosomal copy. SNPtrio is freely accessible at http://pevsnerlab.kennedykrieger.org/SNPtrio.htm (Last accessed: 20 June 2007).

SNPchip: R classes and methods for SNP array data.

UNLABELLED: High-density single nucleotide polymorphism microarrays (SNP chips) provide information on a subject's genome, such as copy number and genotype (heterozygosity/homozygosity) at a SNP. While fluorescence in situ hybridization and karyotyping reveal many abnormalities, SNP chips provide a higher resolution map of the human genome that can be used to detect, e.g., aneuploidies, microdeletions, microduplications and loss of heterozygosity (LOH). As a variety of diseases are linked to such chromosomal abnormalities, SNP chips promise new insights for these diseases by aiding in the discovery of such regions, and may suggest targets for intervention. The R package SNPchip contains classes and methods useful for storing, visualizing and analyzing high density SNP data. Originally developed from the SNPscan web-tool, SNPchip utilizes S4 classes and extends other open source R tools available at Bioconductor. This has numerous advantages, including the ability to build statistical models for SNP-level data that operate on instances of the class, and to communicate with other R packages that add additional functionality. AVAILABILITY: The package is available from the Bioconductor web page at www.bioconductor.org. SUPPLEMENTARY INFORMATION: The supplementary material as described in this article (case studies, installation guidelines and R code) is available from http://biostat.jhsph.edu/~iruczins/publications/sm/

Analysis and visualization of chromosomal abnormalities in SNP data with SNPscan.

BACKGROUND: A variety of diseases are caused by chromosomal abnormalities such as aneuploidies (having an abnormal number of chromosomes), microdeletions, microduplications, and uniparental disomy. High density single nucleotide polymorphism (SNP) microarrays provide information on chromosomal copy number changes, as well as genotype (heterozygosity and homozygosity). SNP array studies generate multiple types of data for each SNP site, some with more than 100,000 SNPs represented on each array. The identification of different classes of anomalies within SNP data has been challenging. RESULTS: We have developed SNPscan, a web-accessible tool to analyze and visualize high density SNP data. It enables researchers (1) to visually and quantitatively assess the quality of user-generated SNP data relative to a benchmark data set derived from a control population, (2) to display SNP intensity and allelic call data in order to detect chromosomal copy number anomalies (duplications and deletions), (3) to display uniparental isodisomy based on loss of heterozygosity (LOH) across genomic regions, (4) to compare paired samples (e.g. tumor and normal), and (5) to generate a file type for viewing SNP data in the University of California, Santa Cruz (UCSC) Human Genome Browser. SNPscan accepts data exported from Affymetrix Copy Number Analysis Tool as its input. We validated SNPscan using data generated from patients with known deletions, duplications, and uniparental disomy. We also inspected previously generated SNP data from 90 apparently normal individuals from the Centre d'Etude du Polymorphisme Humain (CEPH) collection, and identified three cases of uniparental isodisomy, four females having an apparently mosaic X chromosome, two mislabelled SNP data sets, and one microdeletion on chromosome 2 with mosaicism from an apparently normal female. These previously unrecognized abnormalities were all detected using SNPscan. The microdeletion was independently confirmed by fluorescence in situ hybridization, and a region of homozygosity in a UPD case was confirmed by sequencing of genomic DNA. CONCLUSION: SNPscan is useful to identify chromosomal abnormalities based on SNP intensity (such as chromosomal copy number changes) and heterozygosity data (including regions of LOH and some cases of UPD). The program and source code are available at the SNPscan website http://pevsnerlab.kennedykrieger.org/snpscan.htm.

Primary and secondary transcriptional effects in the developing human Down syndrome brain and heart.

BACKGROUND: Down syndrome, caused by trisomic chromosome 21, is the leading genetic cause of mental retardation. Recent studies demonstrated that dosage-dependent increases in chromosome 21 gene expression occur in trisomy 21. However, it is unclear whether the entire transcriptome is disrupted, or whether there is a more restricted increase in the expression of those genes assigned to chromosome 21. Also, the statistical significance of differentially expressed genes in human Down syndrome tissues has not been reported. RESULTS: We measured levels of transcripts in human fetal cerebellum and heart tissues using DNA microarrays and demonstrated a dosage-dependent increase in transcription across different tissue/cell types as a result of trisomy 21. Moreover, by having a larger sample size, combining the data from four different tissue and cell types, and using an ANOVA approach, we identified individual genes with significantly altered expression in trisomy 21, some of which showed this dysregulation in a tissue-specific manner. We validated our microarray data by over 5,600 quantitative real-time PCRs on 28 genes assigned to chromosome 21 and other chromosomes. Gene expression values from chromosome 21, but not from other chromosomes, accurately classified trisomy 21 from euploid samples. Our data also indicated functional groups that might be perturbed in trisomy 21. CONCLUSIONS: In Down syndrome, there is a primary transcriptional effect of disruption of chromosome 21 gene expression, without a pervasive secondary effect on the remaining transcriptome. The identification of dysregulated genes and pathways suggests molecular changes that may underlie the Down syndrome phenotypes.