Search

Purpose: Chordoma is a rare bone tumor with a high recurrence rate and limited treatment options. The aim of this study was to identify molecular subtypes of chordoma that may improve clinical management. Experimental Design: We conducted RNA sequencing in 48 tumors from patients with Chinese skull-base chordoma and identified two major molecular subtypes. We then replicated the classification using a NanoString panel in 48 patients with chordoma from North America. Results: Tumors in one subtype were more likely to have somatic mutations and reduced expression in chromatin remodeling genes, such as PBRM1 and SETD2, whereas the other subtype was characterized by the upregulation of genes in epithelial–mesenchymal transition and Sonic Hedgehog pathways. IHC staining of top differentially expressed genes between the two subtypes in 312 patients with Chinese chordoma with long-term follow-up data showed that the expression of some markers such as PTCH1 was significantly associated with survival outcomes. Conclusions: Our findings may improve the understanding of subtype-specific tumorigenesis of chordoma and inform clinical prognostication and targeted options.

PDF file - 132K, Affymetrix Probesets Differentially Expressed in Stage I Lung Adenocarcinoma Cases Compared to Controls in Current Smokers and All Subjects Combined, Italy, 2002-2005.

Affordable early screening in subjects with high risk of lung cancer has great potential to improve survival from this deadly disease. We measured gene expression from lung tissue and peripheral whole blood (PWB) from adenocarcinoma cases and controls to identify dysregulated lung cancer genes that could be tested in blood to improve identification of at-risk patients in the future. Genome-wide mRNA expression analysis was conducted in 153 subjects (73 adenocarcinoma cases, 80 controls) from the Environment And Genetics in Lung cancer Etiology study using PWB and paired snap-frozen tumor and noninvolved lung tissue samples. Analyses were conducted using unpaired t tests, linear mixed effects, and ANOVA models. The area under the receiver operating characteristic curve (AUC) was computed to assess the predictive accuracy of the identified biomarkers. We identified 50 dysregulated genes in stage I adenocarcinoma versus control PWB samples (false discovery rate ≤0.1, fold change ≥1.5 or ≤0.66). Among them, eight (TGFBR3, RUNX3, TRGC2, TRGV9, TARP, ACP1, VCAN, and TSTA3) differentiated paired tumor versus noninvolved lung tissue samples in stage I cases, suggesting a similar pattern of lung cancer–related changes in PWB and lung tissue. These results were confirmed in two independent gene expression analyses in a blood-based case–control study (n = 212) and a tumor–nontumor paired tissue study (n = 54). The eight genes discriminated patients with lung cancer from healthy controls with high accuracy (AUC = 0.81, 95% CI = 0.74–0.87). Our finding suggests the use of gene expression from PWB for the identification of early detection markers of lung cancer in the future. Cancer Prev Res; 4(10); 1599–608. ©2011 AACR.

PDF file - 222K, Affymetrix Quality Assessment for the Gene Expression Data from Peripheral Whole Blood in Cases and Controls, Italy, 2002-2005.

Figure S1. A. PBRM1 and SETD2 expression between CC1 and CC2; B. PBRM1 and SETD2 expression between NCC1 and NCC2.

Figure S2. A. Contributions of small indel (ID) signatures (A and E) between CC1 and CC2; B. Number of mutations in small indel (ID) signatures (A and E) between CC1 and CC2.

Figure S5. SNAI1, SNAI3, and TWIST expression in CC1 and CC2 tumors.

Supplementary Legend

Figure S3. Gene expression profiling based on the NanoString panel in 23 skull-base (top) and 25 non-skull-base (bottom) chordoma patients from North America.

Table S1. Genes included in the targeted DNA panel sequencing. Table S2. Differences in clinical characteristics by molecular subtype. Table S3. Distributions of mutation status in PBRM1 and SETD2 by molecular subtype in RNASeq and Nanostring datasets. Table S4. Differentially expressed genes among molecular subtypes. Table S5. GSEA results for up- and down-regulated gene sets comparing CC2 to CC1.

Purpose:Chordoma is a rare bone tumor with a high recurrence rate and limited treatment options. The aim of this study was to identify molecular subtypes of chordoma that may improve clinical management.Experimental Design:We conducted RNA sequencing in 48 tumors from patients with Chinese skull-base chordoma and identified two major molecular subtypes. We then replicated the classification using a NanoString panel in 48 patients with chordoma from North America.Results:Tumors in one subtype were more likely to have somatic mutations and reduced expression in chromatin remodeling genes, such as PBRM1 and SETD2, whereas the other subtype was characterized by the upregulation of genes in epithelial–mesenchymal transition and Sonic Hedgehog pathways. IHC staining of top differentially expressed genes between the two subtypes in 312 patients with Chinese chordoma with long-term follow-up data showed that the expression of some markers such as PTCH1 was significantly associated with survival outcomes.Conclusions:Our findings may improve the understanding of subtype-specific tumorigenesis of chordoma and inform clinical prognostication and targeted options.

Figure S4. Expression of selected markers in relation to overall survival (left) and recurrence-free survival (right) in Chinese chordoma patients with long-term follow-up data. A. CLDN11; B. GLI1; C. EPHA3; D. LEF1.

Primers used for validating somatic SNVs

Purpose: Epstein–Barr virus (EBV) infection and a family history of nasopharyngeal carcinoma (NPC) are associated with NPC risk. We examined the risk associated with EBV markers and their clinical utility to identify NPC susceptibles within high-risk NPC families.Experimental Design: We evaluated antibody titers against viral capsid antigen (VCA) IgA, EBV nuclear antigen-1 (EBNA1) IgA, and DNase among unaffected relatives of NPC cases from 358 multiplex families in Taiwan. Incident NPC cases were identified via linkage to the National Cancer Registry. Clinical examinations of 924 individuals were also done to identify occult, asymptomatic NPC. Baseline EBV serology was used to estimate NPC risk using rate ratios with 95% CI. Associated sensitivity/specificity and receiver operating characteristic (ROC) curves were calculated.Results: A total of 2,444 unaffected individuals with 15,519 person-years (6.5 years median follow-up) yielded 14 incident NPC cases (nearly 11 times the general population rate). The absolute rate of NPC among anti-EBV EBNA1 IgA seropositives using a standard positivity cutoff versus an optimized cutoff point defined by ROC analyses was 265/100,000 person-years with a 4.7-fold increased risk of NPC (95% CI: 1.4–16) and 166/100,000 person-years with a 6.6-fold increase (95% CI: 1.5–61), respectively. Sensitivity and specificity using the optimized positivity cutoff points were 85.7% and 51.2%, respectively. It is estimated that active evaluation of 49% of individuals from high-risk NPC families seropositive for this marker could lead to earlier detection of up to 86% of NPC cases. Risks associated with the other three EBV markers were weaker.Conclusions: Future efforts are needed to identify susceptibility markers among high-risk NPC families that maximize both sensitivity and specificity. Clin Cancer Res; 17(7); 1906–14. ©2011 AACR.

Primers used for SV validation experiments

Somatic mutations validated by Sanger sequencing

Background: Germline genetic variants in PLCE1 (10q23) have demonstrated consistent associations with risk of esophageal squamous cell carcinoma (ESCC) and gastric cancer among Chinese. We evaluated PLCE1 mRNA and protein expression in paired tumor-normal tissues, and their relationship with survival.Methods:PLCE1 mRNA was profiled using three probes in the Affymetrix GeneChip U133 for paired tumor-normal tissues of ESCC (n = 132), gastric cardia adenocarcinoma (GCA, n = 62), and gastric noncardia adenocarcinoma (GNCA, n = 72). We used immunohistochemistry to detect PLCE1 protein on slides from tissue microarrays in paired tumor-normal tissues of ESCC (n = 303), and tumors of GCA (n = 298) and GNCA (n = 124).Results: Compared with normal tissues, PLCE1 mRNA expression was significantly reduced in ESCC tumors (P = 0.03, probe_205112_at), as well as in GCA and GNCA tumors (P < 0.0001, each probe). Protein expression was nonsignificantly reduced in ESCC tumors (P = 0.51). Increased tumor-normal mRNA fold change (probe_205112_at) was associated with longer survival in ESCC (9.6 months for highest vs. lowest quartile; Ptrend = 0.02). Increased mRNA tumor-normal fold change (probe_205111_at) was associated with longer survival for GCA (10.7 months for highest quartile; Ptrend = 0.04), but not for GNCA cases (P = 0.72). Similar to mRNA, elevated tumor-normal fold change for protein in ESCC was also associated with improved survival (8.1 months for highest quartile; Ptrend = 0.04).Conclusions: Dysregulated PLCE1 mRNA expression was observed for both ESCC (one probe only) and GCA tumors, and the altered PLCE1 expression seems to be associated with cancer prognosis.Impact: A potential role for PLCE1 in the early detection and/or therapy of ESCC and GCA warrants further investigation. Cancer Epidemiol Biomarkers Prev; 23(8); 1579–88. ©2014 AACR.

Purpose: The molecular drivers that determine histology in lung cancer are largely unknown. We investigated whether microRNA (miR) expression profiles can differentiate histologic subtypes and predict survival for non–small cell lung cancer.Experimental Design: We analyzed miR expression in 165 adenocarcinoma and 125 squamous cell carcinoma (SQ) tissue samples from the Environment And Genetics in Lung cancer Etiology (EAGLE) study using a custom oligo array with 440 human mature antisense miRs. We compared miR expression profiles using t tests and F tests and accounted for multiple testing using global permutation tests. We assessed the association of miR expression with tobacco smoking using Spearman correlation coefficients and linear regression models, and with clinical outcome using log-rank tests, Cox proportional hazards, and survival risk prediction models, accounting for demographic and tumor characteristics.Results: MiR expression profiles strongly differed between adenocarcinoma and SQ (Pglobal < 0.0001), particularly in the early stages, and included miRs located on chromosome loci most often altered in lung cancer (e.g., 3p21-22). Most miRs, including all members of the let-7 family, were downregulated in SQ. Major findings were confirmed by quantitative real time-polymerase chain reaction (qRT-PCR) in EAGLE samples and in an independent set of lung cancer cases. In SQ, the low expression of miRs that are downregulated in the histology comparison was associated with 1.2- to 3.6-fold increased mortality risk. A five-miR signature significantly predicted survival for SQ.Conclusions: We identified a miR expression profile that strongly differentiated adenocarcinoma from SQ and had prognostic implications. These findings may lead to histology-based therapeutic approaches. Clin Cancer Res; 16(2); 430–41

PDF - 194K, Supplementary Figure S1. Immunohistochemistry for PLCE1 protein expression with areas of annotation selected for image quantification (20x.)Supplementary Table S1. PLCE1 mRNA expression in GCA, GNCA and ESCC tumor tissues and tumor-normal expression fold change by characteristics of the subjects. Supplementary Table S2. PLCE1 protein expression in GCA, GNCA and ESCC tumor tissues and tumor-normal expression fold change in ESCC by characteristics of the participants. Supplementary Table S3. The probe-specific association between PLCE1 mRNA expression in tumors and tumor-normal fold change and mortality of ESCC, GCA, and GNCAs. Supplementary Table S4. Spearman correlation coefficients (P values) between PLCE1 mRNA expression in normal tissues of ESCC, GCA, and GNCA and selected PLCE1 SNPs.

Summary of the WGS mapping data

Summary of demographic and clinical data for 15 sequenced cancer cases.

Gastric cancer and esophageal cancer are the second and sixth leading causes of cancer-related death worldwide. Multiple genomic alterations underlying gastric cancer and esophageal squamous cell carcinoma (ESCC) have been identified, but the full spectrum of genomic structural variations and mutations have yet to be uncovered. Here, we report the results of whole-genome sequencing of 30 samples comprising tumor and blood from 15 patients, four of whom presented with ESCC, seven with gastric cardia adenocarcinoma (GCA), and four with gastric noncardia adenocarcinoma. Analyses revealed that an A>C mutation was common in GCA, and in addition to the preferential nucleotide sequence of A located 5 prime to the mutation as noted in previous studies, we found enrichment of T in the 5 prime base. The A>C mutations in GCA suggested that oxidation of guanine may be a potential mechanism underlying cancer mutagenesis. Furthermore, we identified genes with mutations in gastric cancer and ESCC, including well-known cancer genes, TP53, JAK3, BRCA2, FGF2, FBXW7, MSH3, PTCH, NF1, ERBB2, and CHEK2, and potentially novel cancer-associated genes, KISS1R, AMH, MNX1, WNK2, and PRKRIR. Finally, we identified recurrent chromosome alterations in at least 30% of tumors in genes, including MACROD2, FHIT, and PARK2 that were often intragenic deletions. These structural alterations were validated using the The Cancer Genome Atlas dataset. Our studies provide new insights into understanding the genomic landscape, genome instability, and mutation profile underlying gastric cancer and ESCC development. Cancer Res; 76(7); 1714–23. ©2016 AACR.

Summary of somatic SVs

It has legends for the supplementary tables.

Summary of agreement between SNP5 and CGI data

IVA disease analysis report

Supplementary Data from MicroRNA Expression Differentiates Histology and Predicts Survival of Lung Cancer

Supplementary Tables S1-S2.

Supplementary Tables 1A-1C from Genomic Characterization of Esophageal Squamous Cell Carcinoma from a High-Risk Population in China

Genomic instability plays an important role in most human cancers. To characterize genomic instability in esophageal squamous cell carcinoma (ESCC), we examined loss of heterozygosity (LOH), copy number (CN) loss, CN gain, and gene expression using the Affymetrix GeneChip Human Mapping 500K (n = 30 cases) and Human U133A (n = 17 cases) arrays in ESCC cases from a high-risk region of China. We found that genomic instability measures varied widely among cases and separated them into two groups: a high-frequency instability group (two-thirds of all cases with one or more instability category of ≥10%) and a low-frequency instability group (one-third of cases with instability of

Esophageal squamous cell carcinoma (ESCC) shows remarkable variation in incidence that is not fully explained by known lifestyle and environmental risk factors. It has been speculated that an unknown exogenous exposure(s) could be responsible. Here we combine the fields of mutational signature analysis with cancer epidemiology to study 552 ESCC genomes from eight countries with varying incidence rates. Mutational profiles were similar across all countries studied. Associations between specific mutational signatures and ESCC risk factors were identified for tobacco, alcohol, opium and germline variants, with modest impacts on mutation burden. We find no evidence of a mutational signature indicative of an exogenous exposure capable of explaining differences in ESCC incidence. Apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like (APOBEC)-associated mutational signatures single-base substitution (SBS)2 and SBS13 were present in 88% and 91% of cases, respectively, and accounted for 25% of the mutation burden on average, indicating that APOBEC activation is a crucial step in ESCC tumor development.

Genome-wide association studies (GWASs) have identified a melanoma-associated locus on chromosome band 7p21.1 with rs117132860 as the lead SNP and a secondary independent signal marked by rs73069846. rs117132860 is also associated with tanning ability and cutaneous squamous cell carcinoma (cSCC). Because ultraviolet radiation (UVR) is a key environmental exposure for all three traits, we investigated the mechanisms by which this locus contributes to melanoma risk, focusing on cellular response to UVR. Fine-mapping of melanoma GWASs identified four independent sets of candidate causal variants. A GWAS region-focused Capture-C study of primary melanocytes identified physical interactions between two causal sets and the promoter of the aryl hydrocarbon receptor (AHR). Subsequent chromatin state annotation, eQTL, and luciferase assays identified rs117132860 as a functional variant and reinforced AHR as a likely causal gene. Because AHR plays critical roles in cellular response to dioxin and UVR, we explored links between this SNP and AHR expression after both 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and ultraviolet B (UVB) exposure. Allele-specific AHR binding to rs117132860-G was enhanced following both, consistent with predicted weakened AHR binding to the risk/poor-tanning rs117132860-A allele, and allele-preferential AHR expression driven from the protective rs117132860-G allele was observed following UVB exposure. Small deletions surrounding rs117132860 introduced via CRISPR abrogates AHR binding, reduces melanocyte cell growth, and prolongs growth arrest following UVB exposure. These data suggest AHR is a melanoma susceptibility gene at the 7p21.1 risk locus and rs117132860 is a functional variant within a UVB-responsive element, leading to allelic AHR expression and altering melanocyte growth phenotypes upon exposure.

Pancreatic ductal adenocarcinoma (PDAC) is a component of familial melanoma due to germline pathogenic variants (GPVs) inWe analyzed 189 cancer predisposition genes using parametric rare-variant association (RVA) tests and nonparametric permutation tests to identify gene-level associations in PDAC for patients with (In RVA tests,These results suggest that variants in other genes likely play a role in PDAC in all patients and that PDAC in

Background Few studies have evaluated the relationship between CDKN2A germline pathogenic variants (GPV), transcript (p16/p14ARF) alteration, and cancer risk. Methods Standardized incidence ratios (SIRs) comparing cancer risk with the general population were calculated for 385 CDKN2A GPV carriers from 2 large cohorts (259 United States and 126 Swedish individuals) using Poisson regression; statistical significance was defined as P less than .002 (Bonferroni correction). Cumulative incidence is reported for melanoma and nonmelanoma cancer. Results Incidence was increased for melanoma (SIR = 159.8, 95% confidence interval [CI] = 132.1 to 193.2), pancreatic cancer (SIR = 24.1, 95% CI = 14.7 to 39.4), head and neck squamous cell carcinoma (SIR = 16.2, 95% CI = 9.5 to 27.6), and lung cancer (SIR = 5.6, 95% CI = 3.4 to 9.1) in GPV carriers. Similar associations were observed with p16 alteration. Combined p16 and p14ARF alteration was associated with increased incidence of esophageal cancer (SIR = 16.7, 95% CI = 5.7 to 48.9) and malignant peripheral nerve sheath tumor (SIR = 113.0, 95% CI = 16.4 to 780.9), although cancer events were limited (n Conclusion These findings highlight the impact of p16 and p14ARF alteration on cancer risk. Smoking was an important risk factor for smoking-related cancers in our study.

Background: Human Cub and Sushi Multiple Domains 1 (CSMD1) is a novel candidate tumor-suppressor gene that codes for multiple domains, including complement regulatory and adhesion proteins, and has recently been shown to have alterations in multiple cancers. We investigated CSMD1 in esophageal squamous cell carcinoma (ESCC) by performing an integrated analysis on somatic copy number alterations (CNAs), including copy-number gain or loss, allelic imbalance (AI), loss of heterozygosity (LOH), and the expressions of mRNA and its target miRNAs on specimens from the same patients with ESCC. Results: (i) Two-thirds of ESCC patients had all three types of alterations studied—somatic DNA alterations in 70%, and abnormal expressions of CSMD1 RNA in 69% and in target miRNAs in 66%; patterns among these alterations were complex. (ii) In total, 97% of 888 CSMD1 SNPs studied showed somatic DNA alterations, with most located near exons 4–11, 24–25, 39–40, 55–56, and 69–70. (iii) In total, 68% of SNPs with a CNA were correlated with expression of CSMD1. (iv) A total of 33 correlations between non-coding SNPs and expression of CSMD1 target miRs were found. Conclusions: Our results indicate that the CSMD1 gene may play a role in ESCC through complex patterns of DNA alterations and RNA and miRNA expressions. Alterations in some somatic SNPs in non-coding regions of CSMD1 appear to influence expression of this gene and its target miRNAs.

While several high-penetrance melanoma risk genes are known, variation in these genes fail to explain melanoma susceptibility in a large proportion of high-risk families. As part of a melanoma family sequencing study, including 435 families from Mediterranean populations, we identified a novel NRAS variant (c.170A>C, p.D57A) in a melanoma-prone family. This variant is absent in exomes in gnomAD, ESP, UKBiobank, and the 1000 Genomes Project, as well as in 11 273 Mediterranean individuals and 109 melanoma-prone families from the US and Australia. This variant occurs in the GTP-binding pocket of NRAS. Differently from other RAS activating alterations, NRAS D57A expression is unable to activate MAPK-pathway both constitutively and after stimulation but enhances EGF-induced PI3K-pathway signaling in serum starved conditions in vitro. Consistent with in vitro data demonstrating that NRAS D57A does not enrich GTP binding, molecular modeling suggests that the D57A substitution would be expected to impair Mg2+ binding and decrease nucleotide-binding and GTPase activity of NRAS. While we cannot firmly establish NRAS c.170A>C (p.D57A) as a melanoma susceptibility variant, further investigation of NRAS as a familial melanoma gene is warranted.

Chordoma is a rare bone tumor with genetic risk factors largely unknown. We conducted a whole-exome sequencing (WES) analysis of germline DNA from 19 familial chordoma cases in five pedigrees and 137 sporadic chordoma patients and identified 17 rare germline variants in PALB2 and BRCA2, whose products play essential roles in homologous recombination (HR) and tumor suppression. One PALB2 variant showed disease cosegregation in a family with four affected people or obligate gene carrier. Chordoma cases had a significantly increased burden of rare variants in both genes when compared to population-based controls. Four of the six PALB2 variants identified from chordoma patients modestly affected HR function and three of the 11 BRCA2 variants caused loss of function in experimental assays. These results, together with previous reports of abnormal morphology and Brachyury expression of the notochord in Palb2 knockout mouse embryos and genomic signatures associated with HR defect and HR gene mutations in advanced chordomas, suggest that germline mutations in PALB2 and BRCA2 may increase chordoma susceptibility. Our data shed light on the etiology of chordoma and support the previous finding that PARP-1 inhibitors may be a potential therapy for some chordoma patients.