Person:

Szallasi, Zoltan

Gene expression profiles of clinical cohorts can be used to identify genes that are correlated with a clinical variable of interest such as patient outcome or response to a particular drug. However, expression measurements are susceptible to technical bias caused by variation in extraneous factors such as RNA quality and array hybridization conditions. If such technical bias is correlated with the clinical variable of interest, the likelihood of identifying false positive genes is increased. Here we describe a method to visualize an expression matrix as a projection of all genes onto a plane defined by a clinical variable and a technical nuisance variable. The resulting plot indicates the extent to which each gene is correlated with the clinical variable or the technical variable. We demonstrate this method by applying it to three clinical trial microarray data sets, one of which identified genes that may have been driven by a confounding technical variable. This approach can be used as a quality control step to identify data sets that are likely to yield false positive results.

Purpose Quantifying chromosomal instability (CIN) has both prognostic and predictive clinical utility in breast cancer. In order to establish a robust and clinically applicable gene expression-based measure of CIN, we assessed the ability of four qPCR quantified genes selected from the 70-gene Chromosomal Instability (CIN70) expression signature to stratify outcome in patients with grade 2 breast cancer. Methods: AURKA, FOXM1, TOP2A and TPX2 (CIN4), were selected from the CIN70 signature due to their high level of correlation with histological grade and mean CIN70 signature expression in silico. We assessed the ability of CIN4 to stratify outcome in an independent cohort of patients diagnosed between 1999 and 2002. 185 formalin-fixed, paraffin-embedded (FFPE) samples were included in the qPCR measurement of CIN4 expression. In parallel, ploidy status of tumors was assessed by flow cytometry. We investigated whether the categorical CIN4 score derived from the CIN4 signature was correlated with recurrence-free survival (RFS) and ploidy status in this cohort. Results: We observed a significant association of tumor proliferation, defined by Ki67 and mitotic index (MI), with both CIN4 expression and aneuploidy. The CIN4 score stratified grade 2 carcinomas into good and poor prognostic cohorts (mean RFS: 83.8±4.9 and 69.4±8.2 months, respectively, p = 0.016) and its predictive power was confirmed by multivariate analysis outperforming MI and Ki67 expression. Conclusions: The first clinically applicable qPCR derived measure of tumor aneuploidy from FFPE tissue, stratifies grade 2 tumors into good and poor prognosis groups.

Background: Interpretation of gene expression microarrays requires a mapping from probe set to gene. On many Affymetrix gene expression microarrays, a given gene may be detected by multiple probe sets, which may deliver inconsistent or even contradictory measurements. Therefore, obtaining an unambiguous expression estimate of a pre-specified gene can be a nontrivial but essential task. Results: We developed scoring methods to assess each probe set for specificity, splice isoform coverage, and robustness against transcript degradation. We used these scores to select a single representative probe set for each gene, thus creating a simple one-to-one mapping between gene and probe set. To test this method, we evaluated concordance between protein measurements and gene expression values, and between sets of genes whose expression is known to be correlated. For both test cases, we identified genes that were nominally detected by multiple probe sets, and we found that the probe set chosen by our method showed stronger concordance. Conclusions: This method provides a simple, unambiguous mapping to allow assessment of the expression levels of specific genes of interest.

Breast cancer in younger patients often presents with adverse histopathological features, including increased frequency of estrogen receptor negative and lymph node positive disease status. Chromosomal instability (CIN) is increasingly recognised as an important prognostic variable in solid tumours. In a breast cancer meta-analysis of 2423 patients we examine the relationship between clinicopathological parameters and two distinct chromosomal instability gene expression signatures in order to address whether younger age at diagnosis is associated with increased tumour genome instability. We find that CIN, assessed by the two independently derived CIN expression signatures, is significantly associated with increased tumour size, ER negative or HER2 positive disease, higher tumour grade and younger age at diagnosis in ER negative breast cancer. These data support the hypothesis that chromosomal instability may be a defining feature of breast cancer biology and clinical outcome.

Our bloodstream is considered to be an environment well separated from the outside world and the digestive tract. According to the standard paradigm large macromolecules consumed with food cannot pass directly to the circulatory system. During digestion proteins and DNA are thought to be degraded into small constituents, amino acids and nucleic acids, respectively, and then absorbed by a complex active process and distributed to various parts of the body through the circulation system. Here, based on the analysis of over 1000 human samples from four independent studies, we report evidence that meal-derived DNA fragments which are large enough to carry complete genes can avoid degradation and through an unknown mechanism enter the human circulation system. In one of the blood samples the relative concentration of plant DNA is higher than the human DNA. The plant DNA concentration shows a surprisingly precise log-normal distribution in the plasma samples while non-plasma (cord blood) control sample was found to be free of plant DNA.

Thrombocytosis portends adverse prognostic significance in many types of cancers including ovarian and lung carcinoma. In this study, we determined the prevalence and prognostic significance of thrombocytosis (defined as platelet count in excess of 400 × 10 3/μl) in patients with colorectal cancer. We performed a retrospective analysis of 310 consecutive patients diagnosed at our Institution between 2004 and 2013. The patients (48.7% male and 51.3% female) had a mean age of 69.9 years (+/- 12.7 years) at diagnosis. Thrombocytosis was found in a total of 25 patients, with a higher incidence in those with stage III and IV disease (14.4% of patients). Although the mean platelet count increased with the depth of tumor invasion (pT), its values remained within normal limits in the whole patient cohort. No patient with stage I cancer (n=57) had elevated platelet count at diagnosis. By contrast, five of the 78 patients (6.4%) with stage II cancer showed thrombocytosis, and four of these patients showed early recurrence and/or metastatic disease, resulting in shortened survival (they died within one year after surgery). The incidence of thrombocytosis increased to 12.2% and 20.6%, respectively, in patients with stage III and IV disease. The overall survival rate of patients with thrombocytosis was lower than those without thrombocytosis in the stage II and III disease groups, but this difference disappeared in patients with stage IV cancer who did poorly regardless of their platelet count. We concluded that thrombocytosis at diagnosis indicates adverse clinical outcome in colorectal cancer patients with stage II or III disease. This observation is especially intriguing in stage II patients because the clinical management of these patients is controversial. If our data are confirmed in larger studies, stage II colon cancer patients with thrombocytosis may be considered for adjuvant chemotherapy.

The chicken DT40 cell line is a widely used model system in the study of multiple cellular processes due to the efficiency of homologous gene targeting. The cell line was derived from a bursal lymphoma induced by avian leukosis virus infection. In this study we characterized the genome of the cell line using whole genome shotgun sequencing and single nucleotide polymorphism array hybridization. The results indicate that wild-type DT40 has a relatively normal karyotype, except for whole chromosome copy number gains, and no karyotype variability within stocks. In a comparison to two domestic chicken genomes and the Gallus gallus reference genome, we found no unique mutational processes shaping the DT40 genome except for a mild increase in insertion and deletion events, particularly deletions at tandem repeats. We mapped coding sequence mutations that are unique to the DT40 genome; mutations inactivating the PIK3R1 and ATRX genes likely contributed to the oncogenic transformation. In addition to a known avian leukosis virus integration in the MYC gene, we detected further integration sites that are likely to de-regulate gene expression. The new findings support the hypothesis that DT40 is a typical transformed cell line with a relatively intact genome; therefore, it is well-suited to the role of a model system for DNA repair and related processes. The sequence data generated by this study, including a searchable de novo genome assembly and annotated lists of mutated genes, will support future research using this cell line.

Antigen-specific immune responses against peptides derived from missense gene mutations have been identified in multiple cancers. The application of personalized peptide vaccines based on the tumor mutation repertoire of each cancer patient is a near-term clinical reality. These peptides can be identified for pre-validation by leveraging the results of massive gene sequencing efforts in cancer. In this study, we utilized NetMHC 3.2 to predict nanomolar peptide binding affinity to 57 human HLA-A and B alleles. All peptides were derived from 5,685 missense mutations in 312 genes annotated as functionally relevant in the Cancer Genome Project. Of the 26,672,189 potential 8–11 mer peptide-HLA pairs evaluated, 0.4% (127,800) display binding affinities < 50 nM, predicting high affinity interactions. These peptides can be segregated into two groups based on the binding affinity to HLA proteins relative to germline-encoded sequences: peptides for which both the mutant and wild-type forms are high affinity binders, and peptides for which only the mutant form is a high affinity binder. Current evidence directs the attention to mutations that increase HLA binding affinity, as compared with cognate wild-type peptide sequences, as these potentially are more relevant for vaccine development from a clinical perspective. Our analysis generated a database including all predicted HLA binding peptides and the corresponding change in binding affinity as a result of point mutations. Our study constitutes a broad foundation for the development of personalized peptide vaccines that hone-in on functionally relevant targets in multiple cancers in individuals with diverse HLA haplotypes.

Background: Increased number of single nucleotide substitutions is seen in breast and ovarian cancer genomes carrying disease-associated mutations in BRCA1 or BRCA2. The significance of these genome-wide mutations is unknown. We hypothesize genome-wide mutation burden mirrors deficiencies in DNA repair and is associated with treatment outcome in ovarian cancer. Methods and Results: The total number of synonymous and non-synonymous exome mutations (Nmut), and the presence of germline or somatic mutation in BRCA1 or BRCA2 (mBRCA) were extracted from whole-exome sequences of high-grade serous ovarian cancers from The Cancer Genome Atlas (TCGA). Cox regression and Kaplan-Meier methods were used to correlate Nmut with chemotherapy response and outcome. Higher Nmut correlated with a better response to chemotherapy after surgery. In patients with mBRCA-associated cancer, low Nmut was associated with shorter progression-free survival (PFS) and overall survival (OS), independent of other prognostic factors in multivariate analysis. Patients with mBRCA-associated cancers and a high Nmut had remarkably favorable PFS and OS. The association with survival was similar in cancers with either BRCA1 or BRCA2 mutations. In cancers with wild-type BRCA, tumor Nmut was associated with treatment response in patients with no residual disease after surgery. Conclusions: Tumor Nmut was associated with treatment response and with both PFS and OS in patients with high-grade serous ovarian cancer carrying BRCA1 or BRCA2 mutations. In the TCGA cohort, low Nmut predicted resistance to chemotherapy, and for shorter PFS and OS, while high Nmut forecasts a remarkably favorable outcome in mBRCA-associated ovarian cancer. Our observations suggest that the total mutation burden coupled with BRCA1 or BRCA2 mutations in ovarian cancer is a genomic marker of prognosis and predictor of treatment response. This marker may reflect the degree of deficiency in BRCA-mediated pathways, or the extent of compensation for the deficiency by alternative mechanisms.

An unbiased genome-scale screen for unmutated genes that drive cancer growth when overexpressed identified MECP2 as a novel oncogene. MECP2 resides in a region of the Xchromosome that is significantly amplified across 18% of cancers, and many cancer cell lines have amplified, overexpressed MECP2 and are dependent on MECP2 expression for growth. MECP2 copy number gain and RAS family member alterations are mutually exclusive in several cancer types. The MECP2 splicing isoforms activate the major growth factor pathways targeted by activated RAS, the MAPK and PI3K pathways. MECP2 rescued the growth of a KRASG12Caddicted cell line after KRAS down-regulation, and activated KRAS rescues the growth of an MECP2-addicted cell line after MECP2 downregulation. MECP2 binding to the epigenetic modification 5-hydroxymethylcytosine is required for efficient transformation. These observations suggest that MECP2 is a commonly amplified oncogene with an unusual epigenetic mode of action.