Tumour genomic and microenvironmental heterogeneity for integrated prediction of 5-year biochemical recurrence of prostate cancer: a retrospective cohort study.

BACKGROUND Clinical prognostic groupings for localised prostate cancers are imprecise, with 30-50% of patients recurring after image-guided radiotherapy or radical prostatectomy. We aimed to test combined genomic and microenvironmental indices in prostate cancer to improve risk stratification and complement clinical prognostic factors. METHODS We used DNA-based indices alone or in combination with intra-prostatic hypoxia measurements to develop four prognostic indices in 126 low-risk to intermediate-risk patients (Toronto cohort) who will receive image-guided radiotherapy. We validated these indices in two independent cohorts of 154 (Memorial Sloan Kettering Cancer Center cohort [MSKCC] cohort) and 117 (Cambridge cohort) radical prostatectomy specimens from low-risk to high-risk patients. We applied unsupervised and supervised machine learning techniques to the copy-number profiles of 126 pre-image-guided radiotherapy diagnostic biopsies to develop prognostic signatures. Our primary endpoint was the development of a set of prognostic measures capable of stratifying patients for risk of biochemical relapse 5 years after primary treatment. FINDINGS Biochemical relapse was associated with indices of tumour hypoxia, genomic instability, and genomic subtypes based on multivariate analyses. We identified four genomic subtypes for prostate cancer, which had different 5-year biochemical relapse-free survival. Genomic instability is prognostic for relapse in both image-guided radiotherapy (multivariate analysis hazard ratio [HR] 4·5 [95% CI 2·1-9·8]; p=0·00013; area under the receiver operator curve [AUC] 0·70 [95% CI 0·65-0·76]) and radical prostatectomy (4·0 [1·6-9·7]; p=0·0024; AUC 0·57 [0·52-0·61]) patients with prostate cancer, and its effect is magnified by intratumoral hypoxia (3·8 [1·2-12]; p=0·019; AUC 0·67 [0·61-0·73]). A novel 100-loci DNA signature accurately classified treatment outcome in the MSKCC low-risk to intermediate-risk cohort (multivariate analysis HR 6·1 [95% CI 2·0-19]; p=0·0015; AUC 0·74 [95% CI 0·65-0·83]). In the independent MSKCC and Cambridge cohorts, this signature identified low-risk to high-risk patients who were most likely to fail treatment within 18 months (combined cohorts multivariate analysis HR 2·9 [95% CI 1·4-6·0]; p=0·0039; AUC 0·68 [95% CI 0·63-0·73]), and was better at predicting biochemical relapse than 23 previously published RNA signatures. INTERPRETATION This is the first study of cancer outcome to integrate DNA-based and microenvironment-based failure indices to predict patient outcome. Patients exhibiting these aggressive features after biopsy should be entered into treatment intensification trials. FUNDING Movember Foundation, Prostate Cancer Canada, Ontario Institute for Cancer Research, Canadian Institute for Health Research, NIHR Cambridge Biomedical Research Centre, The University of Cambridge, Cancer Research UK, Cambridge Cancer Charity, Prostate Cancer UK, Hutchison Whampoa Limited, Terry Fox Research Institute, Princess Margaret Cancer Centre Foundation, PMH-Radiation Medicine Program Academic Enrichment Fund, Motorcycle Ride for Dad (Durham), Canadian Cancer Society.

Epidermal growth factor receptor (EGFR) is an independent adverse prognostic factor in esophageal adenocarcinoma patients treated with cisplatin-based neoadjuvant chemotherapy.

Neoadjuvant platin-based therapy is accepted as a standard therapy for advanced esophageal adenocarcinoma (EAC). Patients who respond have a better survival prognosis, but still a significant number of responder patients die from tumor recurrence. Molecular markers for prognosis in neoadjuvantly treated EAC patients have not been identified yet. We investigated the epidermal growth factor receptor (EGFR) in prognosis and chemotherapy resistance in these patients. Two EAC patient cohorts, either treated by neoadjuvant cisplatin-based chemotherapy followed by surgery (n=86) or by surgical resection (n=46) were analyzed for EGFR protein expression and gene copy number. Data were correlated with clinical and histopathological response, disease-free and overall survival. In case of EGFR overexpression, the prognosis for neoadjuvant chemotherapy responders was poor as in non-responders. Responders had a significantly better disease-free survival than non-responders only if EGFR expression level (p=0.0152) or copy number (p=0.0050) was low. Comparing neoadjuvantly treated patients and primary resection patients, tumors of non-responder patients more frequently exhibited EGFR overexpression, providing evidence that EGFR is a factor for indicating chemotherapy resistance. EGFR overexpression and gene copy number are independent adverse prognostic factors for neoadjuvant chemotherapy-treated EAC patients, particularly for responders. Furthermore, EGFR overexpression is involved in resistance to cisplatin-based neoadjuvant chemotherapy.

16p11.2 600 kb Duplications confer risk for typical and atypical Rolandic epilepsy

Rolandic epilepsy (RE) is the most common idiopathic focal childhood epilepsy. Its molecular basis is largely unknown and a complex genetic etiology is assumed in the majority of affected individuals. The present study tested whether six large recurrent copy number variants at 1q21, 15q11.2, 15q13.3, 16p11.2, 16p13.11 and 22q11.2 previously associated with neurodevelopmental disorders also increase risk of RE. Our association analyses revealed a significant excess of the 600 kb genomic duplication at the 16p11.2 locus (chr16: 29.5-30.1 Mb) in 393 unrelated patients with typical (n = 339) and atypical (ARE; n = 54) RE compared with the prevalence in 65,046 European population controls (5/393 cases versus 32/65,046 controls; Fisher's exact test P = 2.83 × 10(-6), odds ratio = 26.2, 95% confidence interval: 7.9-68.2). In contrast, the 16p11.2 duplication was not detected in 1738 European epilepsy patients with either temporal lobe epilepsy (n = 330) and genetic generalized epilepsies (n = 1408), suggesting a selective enrichment of the 16p11.2 duplication in idiopathic focal childhood epilepsies (Fisher's exact test P = 2.1 × 10(-4)). In a subsequent screen among children carrying the 16p11.2 600 kb rearrangement we identified three patients with RE-spectrum epilepsies in 117 duplication carriers (2.6%) but none in 202 carriers of the reciprocal deletion. Our results suggest that the 16p11.2 duplication represents a significant genetic risk factor for typical and atypical RE.

CD41 T cell recovery during suppression of HIV replication: an international comparison of the immunological efficacy of antiretroviral therapy in North America, Asia and Africa.

BACKGROUND Even among HIV-infected patients who fully suppress plasma HIV RNA replication on antiretroviral therapy, genetic (e.g. CCL3L1 copy number), viral (e.g. tropism) and environmental (e.g. chronic exposure to microbial antigens) factors influence CD4 recovery. These factors differ markedly around the world and therefore the expected CD4 recovery during HIV RNA suppression may differ globally. METHODS We evaluated HIV-infected adults from North America, West Africa, East Africa, Southern Africa and Asia starting non-nucleoside reverse transcriptase inhibitorbased regimens containing efavirenz or nevirapine, who achieved at least one HIV RNA level <500/ml in the first year of therapy and observed CD4 changes during HIV RNA suppression. We used a piecewise linear regression to estimate the influence of region of residence on CD4 recovery, adjusting for socio-demographic and clinical characteristics. We observed 28 217 patients from 105 cohorts over 37 825 person-years. RESULTS After adjustment, patients from East Africa showed diminished CD4 recovery as compared with other regions. Three years after antiretroviral therapy initiation, the mean CD4 count for a prototypical patient with a pre-therapy CD4 count of 150/ml was 529/ml [95% confidence interval (CI): 517–541] in North America, 494/ml (95% CI: 429–559) in West Africa, 515/ml (95% CI: 508–522) in Southern Africa, 503/ml (95% CI: 478–528) in Asia and 437/ml (95% CI: 425–449) in East Africa. CONCLUSIONS CD4 recovery during HIV RNA suppression is diminished in East Africa as compared with other regions of the world, and observed differences are large enough to potentially influence clinical outcomes. Epidemiological analyses on a global scale can identify macroscopic effects unobservable at the clinical, national or individual regional level.

Differential spatial and structural organization of the X chromosome underlies dosage compensation in C. elegans.

The adjustment of X-linked gene expression to the X chromosome copy number (dosage compensation [DC]) has been widely studied as a model of chromosome-wide gene regulation. In Caenorhabditis elegans, DC is achieved by twofold down-regulation of gene expression from both Xs in hermaphrodites. We show that in males, the single X chromosome interacts with nuclear pore proteins, while in hermaphrodites, the DC complex (DCC) impairs this interaction and alters X localization. Our results put forward a structural model of DC in which X-specific sequences locate the X chromosome in transcriptionally active domains in males, while the DCC prevents this in hermaphrodites.

Detection of one VH antibody sequence in both healthy donors and urticaria patients.

We have previously isolated anti-FcepsilonRIalpha autoantibodies from phage libraries of healthy donors and urticaria patients. Strikingly, the same antibody, LTMalpha15, was isolated from both libraries. Sequence analysis revealed a germline configuration of the LTMalpha15 variable heavy (V(H)) chain with a slightly mutated variable light (V(L)) chain supporting its classification as a natural autoantibody. Distribution analysis of anti-FcepsilonRIalpha autoantibodies by functional or serological tests delivered conflicting data. For this reason we have developed a new real-time PCR to analyse the distribution of LTMalpha15V(H) in healthy donors and urticaria patients. Our new bioinformatic program permitted the design of a minor groove binder (MGB) TaqMan probe that specifically detected the LTMalpha15V(H). We were able to demonstrate a broad range of rearranged V(H) gene copy number without any correlation to the state of health. Monitoring LTMalpha15V(H) gene copy number in a single donor over a period of 70 days revealed a time-related fluctuation of circulating B cells carrying LTMalpha15V(H). We propose that our real-time PCR may serve as a model for the quantification of natural antibody sequences at a monoclonal level.

Melanotic tumors of the nervous system are characterized by distinct mutational, chromosomal and epigenomic profiles.

Melanotic tumors of the nervous system show overlapping histological characteristics but differ substantially in their biological behavior. In order to achieve a better delineation of such tumors, we performed an in-depth molecular characterization. Eighteen melanocytomas, 12 melanomas, and 14 melanotic and 14 conventional schwannomas (control group) were investigated for methylome patterns (450k array), gene mutations associated with melanotic tumors and copy number variants (CNVs). The methylome fingerprints assigned tumors to entity-specific groups. Methylation groups also showed a substantial overlap with histology-based diagnosis suggesting that they represent true biological entities. On the molecular level, melanotic schwannomas were characterized by a complex karyotype with recurrent monosomy of chromosome 22q and variable whole chromosomal gains and recurrent losses commonly involving chromosomes 1, 17p and 21. Melanocytomas carried GNAQ/11 mutations and presented with CNV involving chromosomes 3 and 6. Melanomas were frequently mutated in the TERT promoter, harbored additional oncogene mutations and showed recurrent chromosomal losses involving chromosomes 9, 10 and 6q, as well as gains of 22q. Together, melanotic nervous system tumors have several distinct mutational and chromosomal alterations and can reliably be distinguished by methylome profiling.

Electron cryomicroscopy of E. coli reveals filament bundles involved in plasmid DNA segregation.

Bipolar elongation of filaments of the bacterial actin homolog ParM drives movement of newly replicated plasmid DNA to opposite poles of a bacterial cell. We used a combination of vitreous sectioning and electron cryotomography to study this DNA partitioning system directly in native, frozen cells. The diffraction patterns from overexpressed ParM bundles in electron cryotomographic reconstructions were used to unambiguously identify ParM filaments in Escherichia coli cells. Using a low-copy number plasmid encoding components required for partitioning, we observed small bundles of three to five intracellular ParM filaments that were situated close to the edge of the nucleoid. We propose that this may indicate the capture of plasmid DNA within the periphery of this loosely defined, chromosome-containing region.

Spatial genomic heterogeneity within localized, multifocal prostate cancer

Herein we provide a detailed molecular analysis of the spatial heterogeneity of clinically localized, multifocal prostate cancer to delineate new oncogenes or tumor suppressors. We initially determined the copy number aberration (CNA) profiles of 74 patients with index tumors of Gleason score 7. Of these, 5 patients were subjected to whole-genome sequencing using DNA quantities achievable in diagnostic biopsies, with detailed spatial sampling of 23 distinct tumor regions to assess intraprostatic heterogeneity in focal genomics. Multifocal tumors are highly heterogeneous for single-nucleotide variants (SNVs), CNAs and genomic rearrangements. We identified and validated a new recurrent amplification of MYCL, which is associated with TP53 deletion and unique profiles of DNA damage and transcriptional dysregulation. Moreover, we demonstrate divergent tumor evolution in multifocal cancer and, in some cases, tumors of independent clonal origin. These data represent the first systematic relation of intraprostatic genomic heterogeneity to predicted clinical outcome and inform the development of novel biomarkers that reflect individual prognosis.

A certified plasmid reference material for the standardisation of BCR-ABL1 mRNA quantification by real-time quantitative PCR.

Serial quantification of BCR-ABL1 mRNA is an important therapeutic indicator in chronic myeloid leukaemia, but there is a substantial variation in results reported by different laboratories. To improve comparability, an internationally accepted plasmid certified reference material (CRM) was developed according to ISO Guide 34:2009. Fragments of BCR-ABL1 (e14a2 mRNA fusion), BCR and GUSB transcripts were amplified and cloned into pUC18 to yield plasmid pIRMM0099. Six different linearised plasmid solutions were produced with the following copy number concentrations, assigned by digital PCR, and expanded uncertainties: 1.08±0.13 × 10(6), 1.08±0.11 × 10(5), 1.03±0.10 × 10(4), 1.02±0.09 × 10(3), 1.04±0.10 × 10(2) and 10.0±1.5 copies/μl. The certification of the material for the number of specific DNA fragments per plasmid, copy number concentration of the plasmid solutions and the assessment of inter-unit heterogeneity and stability were performed according to ISO Guide 35:2006. Two suitability studies performed by 63 BCR-ABL1 testing laboratories demonstrated that this set of 6 plasmid CRMs can help to standardise a number of measured transcripts of e14a2 BCR-ABL1 and three control genes (ABL1, BCR and GUSB). The set of six plasmid CRMs is distributed worldwide by the Institute for Reference Materials and Measurements (Belgium) and its authorised distributors (https://ec.europa.eu/jrc/en/reference-materials/catalogue/; CRM code ERM-AD623a-f).

EXPRESSION OF THE HUMAN PLACENTAL LACTOGEN GENES

The placenta is the site of synthesis of various peptide and steroid hormones related to pregnancy. Human placental lactogen (hPL) is the predominant peptide hormone secreted by term placenta and its synthesis is tissue-specific and coupled to placenta development. The objective of this work was to study the structure and expression of the hPL.^ Poly(A('+))RNA from human term placenta was translated in a mouse-derived cell-free system. A major band corresponding to pre-hPL and a minor band comigrating with mature hPL, represent (TURN)15% of the total radioactively labeled proteins. Analysis of the poly(A('+))RNA showed a prominent band at approximately 860 nucleotides. A corresponding band was observed in Northern blots of total RNA, hybridized with {('32)P}-labeled recombinant plasmid containing a portion of hPL cDNA. Similar analyses of nuclear RNA showed at least four additional bands at 990, 1200, 1460 and 1760 nucleotides, respectively, which are likely precursors of hPL mRNA. Poly(A('+))RNA was used to construct a cDNA library, of which approximately 5% of the clones were found to hybridize to hPL DNA sequences. Heteroduplexes constructed between a clone containing a 815 bp hPL cDNA insert and a hPL genomic DNA clone revealed four small intervening sequences which can account for the lengths observed in hnRNA molecules.^ Recombinant plasmid HCS-pBR322 containing a 550 bp insert of a cDNA transcript of human placental lactogen (hPL) mRNA was ('3)H-labeled an hybridized in situ to human chromosome preparations. These experiments allowed assignment of the hPL and growth hormone (hGH) genes, which have over 90% nucleotide homology in their coding sequences, to band q22-24 of chromosome 17. A gene copy number experiment showed that both genes are present in (TURN)3 copies per haploid genome.^ Experiments were designed to determine if all members of the hPL gene cluster, consisting of four non-allelic genes, are transcribed in term placenta. Advantage was taken of differences in restriction endonuclease sites in the coding portions of the different hPL genes, to distinguish the putative cDNAs of the transcriptionally active genes. Two genes were found to be represented in the cDNA library and their cDNA transcripts were isolated and characterized. Three independent methods showed that their corresponding mRNAs are about equally represented in the hPL mRNA population. The two cDNAs code for prehPL proteins which differ at a single amino acid position. However the secreted hPLs have identical amino acid sequences. A tetramer insertion duplication was found in a palindrome area of the 3' untranslated region of one of the hPL mRNAs. ^

Molecular alterations in nucleotide excision repair genes in malignant glioma and their relevance to malignant progression and drug resistance

Recently, it has become apparent that DNA repair mechanisms are involved in the malignant progression and resistance to therapy of gliomas. Many investigators have shown that increased levels of O6-methyl guanine DNA alkyltransferase, a DNA monoalkyl adduct repair enzyme, are correlated with resistance of malignant glioma cell lines to nitrosourea-based chemotherapy. Three important DNA excision repair genes ERCC1 (excision repair cross complementation group 1), ERCC2 (excision repair cross complementation group 2), and ERCC6 (excision repair cross complementation group 6) have been studied in human tumors. Gene copy number variation of ERCC1 and ERCC2 has been observed in primary glioma tissues. A number of reports describing a relationship between ERCC1 gene alterations and resistance to anti-cancer drugs have been also described. The levels of ERCC1 gene expression, however, have not been correlated with drug resistance in gliomas. The expression of ERCC6 gene transcribes has been shown to vary with tissue types and to be highest in the brain. There have been no comprehensive studies so far, however, of ERCC6 gene expression and molecular alterations in malignant glioma. This project examined the ERCC1 expression levels and correlated them with cisplatin resistance in malignant glioma cell lines. We also examined the molecular alterations of ERCC6 gene in primary glioma tissues and cells and analyzed whether these alterations are related to tumor progression and chemotherapy resistance. Our results indicate the presence of mutations and/or deletions in exons II and V of the ERCC6 gene, and these alterations are more frequent in exon II. Furthermore, the mutations and/or deletions in exon II were shown to be associated with increased malignant grade of gliomas. The results on the Levels of ERCC1 gene transcripts showed that expression levels correlate with cisplatin resistance. The increase in ERCC1 mRNA induced by cisplatin could be down-regulated by cyclosporin A and herbimycin A. The results of this study are likely to provide useful information for clinical treatment of human gliomas. ^

High-resolution microarray analysis of chromosome 20q in human colon cancer metastasis model systems

Amplification of human chromosome 20q DNA is the most frequently occurring chromosomal abnormality detected in sporadic colorectal carcinomas and shows significant correlation with liver metastases. Through comprehensive high-resolution microarray comparative genomic hybridization and microarray gene expression profiling, we have characterized chromosome 20q amplicon genes associated with human colorectal cancer metastasis in two in vitro metastasis model systems. The results revealed increasing complexity of the 20q genomic profile from the primary tumor-derived cell lines to the lymph node and liver metastasis derived cell lines. Expression analysis of chromosome 20q revealed a subset of over expressed genes residing within the regions of genomic copy number gain in all the tumor cell lines, suggesting these are Chromosome 20q copy number responsive genes. Bases on their preferential expression levels in the model system cell lines and known biological function, four of the over expressed genes mapping to the common intervals of genomic copy gain were considered the most promising candidate colorectal metastasis-associated genes. Validation of genomic copy number and expression array data was carried out on these genes, with one gene, DNMT3B, standing out as expressed at a relatively higher levels in the metastasis-derived cell lines compared with their primary-derived counterparts in both the models systems analyzed. The data provide evidence for the role of chromosome 20q genes with low copy gain and elevated expression in the clonal evolution of metastatic cells and suggests that such genes may serve as early biomarkers of metastatic potential. The data also support the utility of the combined microarray comparative genomic hybridization and expression array analysis for identifying copy number responsive genes in areas of low DNA copy gain in cancer cells. ^

INTEGRATIVE BIOMARKER IDENTIFICATION AND CLASSIFICATION USING HIGH THROUGHPUT ASSAYS

It is well accepted that tumorigenesis is a multi-step procedure involving aberrant functioning of genes regulating cell proliferation, differentiation, apoptosis, genome stability, angiogenesis and motility. To obtain a full understanding of tumorigenesis, it is necessary to collect information on all aspects of cell activity. Recent advances in high throughput technologies allow biologists to generate massive amounts of data, more than might have been imagined decades ago. These advances have made it possible to launch comprehensive projects such as (TCGA) and (ICGC) which systematically characterize the molecular fingerprints of cancer cells using gene expression, methylation, copy number, microRNA and SNP microarrays as well as next generation sequencing assays interrogating somatic mutation, insertion, deletion, translocation and structural rearrangements. Given the massive amount of data, a major challenge is to integrate information from multiple sources and formulate testable hypotheses. This thesis focuses on developing methodologies for integrative analyses of genomic assays profiled on the same set of samples. We have developed several novel methods for integrative biomarker identification and cancer classification. We introduce a regression-based approach to identify biomarkers predictive to therapy response or survival by integrating multiple assays including gene expression, methylation and copy number data through penalized regression. To identify key cancer-specific genes accounting for multiple mechanisms of regulation, we have developed the integIRTy software that provides robust and reliable inferences about gene alteration by automatically adjusting for sample heterogeneity as well as technical artifacts using Item Response Theory. To cope with the increasing need for accurate cancer diagnosis and individualized therapy, we have developed a robust and powerful algorithm called SIBER to systematically identify bimodally expressed genes using next generation RNAseq data. We have shown that prediction models built from these bimodal genes have the same accuracy as models built from all genes. Further, prediction models with dichotomized gene expression measurements based on their bimodal shapes still perform well. The effectiveness of outcome prediction using discretized signals paves the road for more accurate and interpretable cancer classification by integrating signals from multiple sources.

Mechanisms of HER2/{\it neu\/} proto-oncogene overexpression in breast cancer cells

Overexpression and amplification of HER2/neu have been documented in many primary tumors, most notably in breast. Not only do approximately 30% of breast cancer patients carry tumors that overexpress the gene, but those that do generally have shorter overall and disease-free survival times than patients with tumors expressing low levels of HER2/neu. Thus, overexpression of HER2/neu plays an important role in the pathogenesis of breast cancer. We have examined the mechanisms that result in HER2/neu overexpression in breast cancer by using, as a model system, established breast cancer cell lines that express much higher levels of HER2/neu mRNA than normal breast tissue while maintaining a near normal HER2/neu gene copy number. Nuclear run-on experiments indicate that the breast cancer cell lines MDA-MB453, BT483, and BT474 have an increased HER2/neu gene transcription rate. By using HER2/neu promoter-CAT constructs, we have found that the enhanced HER2/neu transcription rate in MDA-MB453 cells is due to activation of the gene in trans, while the enhanced transcription rate in BT483 cells is due to activation of the gene in either trans or cis. In BT474 cells, transcriptional upregulation is primarily due to gene amplification. Since the levels of increased transcription are not as high as the levels of HER2/neu mRNA in any of these three lines, post-transcriptional deregulation that increases HER2/neu expression must also be functioning in these cells. The half-life of HER2/neu mRNA was measured and found to be equivalent in these lines as in a control. Thus, the post-transcriptional deregulation is not increased stability of the HER2/neu transcript.^ Much work has been performed in characterizing the altered trans-acting factor involved in increased HER2/neu transcription in MDA-MB453 cells. Using promoter deletion constructs linked to a reporter gene, the region responsive to this factor was localized in the rat neu promoter. When human HER2/neu promoter constructs were used, the homologous sequence in the human promoter was identified. Furthermore, a number of protein/DNA complexes are detected when these promoter regions are used in gel mobility shift assays. UV-crosslinking experiments indicate DNA-binding proteins of roughly 110 kDa, 70 kDa, and 35 kDa are capable of interacting with the human promoter element. ^