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.

Genomic heterogeneity of small ruminant lentiviruses detected by PCR

In order to detect a large spectrum of small ruminant lentiviruses, primers for PCR were chosen in conserved parts of the LTR and GAG genes of Icelandic Visna virus 1514 and of the POL gene of caprine arthritis-encephalitis virus. This set of primers was tested in six different caprine arthritis-encephalitis virus (CAEV)- and Maedi-Visna virus isolates of Dutch, American and Swiss origin. The LTR primers allowed the detection of the corresponding fragments of all isolates. The GAG primers allowed amplification of the corresponding fragments of all but the Swiss Maedi-Visna virus strain OLV. Using the POL primers, one Maedi-Visna- and two caprine arthritis-encephalitis virus strains were detected after one round of amplification. Sequencing of the GAG and POL amplification products and comparison to Icelandic Visna virus and CAEV strain CO revealed total heterogeneity of 38% for the GAG- and 28% for the POL fragment. The virus strains studied fall into two groups which are more closely related to one another than to Icelandic Visna virus.

A genomic perspective on a new bacterial genus and species from the Alcaligenaceae family, Basilea psittacipulmonis

BACKGROUND A novel Gram-negative, non-haemolytic, non-motile, rod-shaped bacterium was discovered in the lungs of a dead parakeet (Melopsittacus undulatus) that was kept in captivity in a petshop in Basel, Switzerland. The organism is described with a chemotaxonomic profile and the nearly complete genome sequence obtained through the assembly of short sequence reads. RESULTS Genome sequence analysis and characterization of respiratory quinones, fatty acids, polar lipids, and biochemical phenotype is presented here. Comparison of gene sequences revealed that the most similar species is Pelistega europaea, with BLAST identities of only 93% to the 16S rDNA gene, 76% identity to the rpoB gene, and a similar GC content (~43%) as the organism isolated from the parakeet, DSM 24701 (40%). The closest full genome sequences are those of Bordetella spp. and Taylorella spp. High-throughput sequencing reads from the Illumina-Solexa platform were assembled with the Edena de novo assembler to form 195 contigs comprising the ~2 Mb genome. Genome annotation with RAST, construction of phylogenetic trees with the 16S rDNA (rrs) gene sequence and the rpoB gene, and phylogenetic placement using other highly conserved marker genes with ML Tree all suggest that the bacterial species belongs to the Alcaligenaceae family. Analysis of samples from cages with healthy parakeets suggested that the newly discovered bacterial species is not widespread in parakeet living quarters. CONCLUSIONS Classification of this organism in the current taxonomy system requires the formation of a new genus and species. We designate the new genus Basilea and the new species psittacipulmonis. The type strain of Basilea psittacipulmonis is DSM 24701 (= CIP 110308 T, 16S rDNA gene sequence Genbank accession number JX412111 and GI 406042063).

An arranged marriage for precision medicine: hypoxia and genomic assays in localized prostate cancer radiotherapy

Prostate cancer (CaP) is the most commonly diagnosed malignancy in males in the Western world with one in six males diagnosed in their lifetime. Current clinical prognostication groupings use pathologic Gleason score, pre-treatment prostatic-specific antigen and Union for International Cancer Control-TNM staging to place patients with localized CaP into low-, intermediate- and high-risk categories. These categories represent an increasing risk of biochemical failure and CaP-specific mortality rates, they also reflect the need for increasing treatment intensity and justification for increased side effects. In this article, we point out that 30-50% of patients will still fail image-guided radiotherapy or surgery despite the judicious use of clinical risk categories owing to interpatient heterogeneity in treatment response. To improve treatment individualization, better predictors of prognosis and radiotherapy treatment response are needed to triage patients to bespoke and intensified CaP treatment protocols. These should include the use of pre-treatment genomic tests based on DNA or RNA indices and/or assays that reflect cancer metabolism, such as hypoxia assays, to define patient-specific CaP progression and aggression. More importantly, it is argued that these novel prognostic assays could be even more useful if combined together to drive forward precision cancer medicine for localized CaP.

DNA hypomethylation and oxidative stress-mediated increase in genomic instability in equine sarcoid-derived fibroblasts

It is widely accepted that equine sarcoid disease, the most common skin associated neoplasm in equids, is induced by bovine papillomavirus (BPV-1). Although BPV-1 DNA has been found in almost all examined sarcoids so far, its detailed impact on the horse's host cell metabolism is largely unknown. We used equine fibroblast cell lines originating from sarcoid biopsies to study BPV-1-associated changes on DNA methylation status and oxidative stress parameters. Sarcoid-derived fibroblasts manifested increased proliferation in vitro, transcriptional rDNA activity (NORs expression) and DNA hypomethylation compared to control cells. Cells isolated from equine sarcoids suffered from oxidative stress: the expression of antioxidant enzymes was decreased and the superoxide production was increased. Moreover, increased ploidy, oxidative DNA damage and micronuclei formation was monitored in sarcoid cells. We postulate that both altered DNA methylation status and redox milieu may affect genomic stability in BPV-1-infected cells and in turn contribute to sarcoid pathology.

The genomic substrate for adaptive radiation in African cichlid fish

Cichlid fishes are famous for large, diverse and replicated adaptive radiations in the Great Lakes of East Africa. To understand the molecular mechanisms underlying cichlid phenotypic diversity, we sequenced the genomes and transcriptomes of five lineages of African cichlids: the Nile tilapia (Oreochromis niloticus), an ancestral lineage with low diversity; and four members of the East African lineage: Neolamprologus brichardi/pulcher (older radiation, Lake Tanganyika), Metriaclima zebra (recent radiation, Lake Malawi), Pundamilia nyererei (very recent radiation, Lake Victoria), and Astatotilapia burtoni (riverine species around Lake Tanganyika). We found an excess of gene duplications in the East African lineage compared to tilapia and other teleosts, an abundance of non-coding element divergence, accelerated coding sequence evolution, expression divergence associated with transposable element insertions, and regulation by novel microRNAs. In addition, we analysed sequence data from sixty individuals representing six closely related species from Lake Victoria, and show genome-wide diversifying selection on coding and regulatory variants, some of which were recruited from ancient polymorphisms. We conclude that a number of molecular mechanisms shaped East African cichlid genomes, and that amassing of standing variation during periods of relaxed purifying selection may have been important in facilitating subsequent evolutionary diversification.

Impact of range expansions on current human genomic diversity

The patterns of population genetic diversity depend to a large extent on past demographic history. Most human populations are known to have gone recently through a series of range expansions within and out of Africa, but these spatial expansions are rarely taken into account when interpreting observed genomic diversity, possibly because they are difficult to model. Here we review available evidence in favour of range expansions out of Africa, and we discuss several of their consequences on neutral and selected diversity, including some recent observations on an excess of rare neutral and selected variants in large samples. We further show that in spatially subdivided populations, the sampling strategy can severely impact the resulting genetic diversity and be confounded by past demography. We conclude that ignoring the spatial structure of human population can lead to some misinterpretations of extant genetic diversity.

Identification of a new genomic hot spot of evolutionary diversification of protein function.

Establishment of phylogenetic relationships remains a challenging task because it is based on computational analysis of genomic hot spots that display species-specific sequence variations. Here, we identify a species-specific thymine-to-guanine sequence variation in the Glrb gene which gives rise to species-specific splice donor sites in the Glrb genes of mouse and bushbaby. The resulting splice insert in the receptor for the inhibitory neurotransmitter glycine (GlyR) conveys synaptic receptor clustering and specific association with a particular synaptic plasticity-related splice variant of the postsynaptic scaffold protein gephyrin. This study identifies a new genomic hot spot which contributes to phylogenetic diversification of protein function and advances our understanding of phylogenetic relationships.

Genomic knock-out of rancRNAs in the archaeon H. volcanii

As translation is the final step in gene expression it is particularly important to understand the processes involved in translation regulation. It was shown in the last years that a class of RNA, the non-protein-coding RNAs (ncRNAs), is involved in regulation of gene expression via various mechanisms [1]. Herein included is the prominent example of gene silencing caused by micro RNAs (miRNAs) and small interfering RNAs (siRNAs). Almost all of these ncRNA discovered so far target the mRNA in order to modulate protein biosynthesis, this is rather unexpected considering the crucial role of the ribosome during gene expression. However, recent data from our laboratory showed that there is a new class of RNAs among the well-studied ncRNAs that target the ribosome itself [2,3]. These so called ribosome-associated ncRNAs (rancRNAs) have an impact on translation regulation, mainly by interfering / modulating the rate of protein biosynthesis. Recent studies show the presence of small regulatory RNAs (sRNAs) in archaea which are involved in many biological processes including stress response and metabolic regulation [4]. To date the biological function and the targets of these archaeal sRNAs are only described for a few examples. There are reports of sRNAs binding to the 5’ as well as to the 3’ of mRNAs [5,6]. In addition to these findings, a tRNA derived fragment (tRF) of Valine tRNA was found in a genomic screen of RNAs associated with the ribosome in H. volcanii in our laboratory [3]. This Valine tRF seems to be processed in a stress-dependent manner and showed in vitro binding to the ribosome and inhibited in vitro translation. These results showed that Valine tRF is capable to regulate translation in H. volcanii by targeting the ribosome. The main goal of this project is to identify and describe novel potential regulatory rancRNAs in H. volcanii with the focus on intergenic candidates. Northern blot analyses already revealed interactions with the ribosome and showed differential expression patterns in response to stress conditions. To investigate the biological relevance of some of the ribosome-associated ncRNA candidates, knock-out and phenotypic characterization studies are done. The genomic knock out of a hypothetical ORF (198nt), where one putative rancRNA candidate (46nt) named IG33 was detected in the library at the beginning of the ORF, showed interesting growth phenotype under specific stress conditions. Furthermore a strain with an introduced start to stop codon mutation in this hypothetical ORF still shows the same phenotype indicating that rather the missing protein than the missing sRNA causes this growth phenotype.

Profiling invasiveness in head and neck cancer: recent contributions of genomic and transcriptomic approaches.

High-throughput molecular profiling approaches have emerged as precious research tools in the field of head and neck translational oncology. Such approaches have identified and/or confirmed the role of several genes or pathways in the acquisition/maintenance of an invasive phenotype and the execution of cellular programs related to cell invasion. Recently published new-generation sequencing studies in head and neck squamous cell carcinoma (HNSCC) have unveiled prominent roles in carcinogenesis and cell invasion of mutations involving NOTCH1 and PI3K-patwhay components. Gene-expression profiling studies combined with systems biology approaches have allowed identifying and gaining further mechanistic understanding into pathways commonly enriched in invasive HNSCC. These pathways include antigen-presenting and leucocyte adhesion molecules, as well as genes involved in cell-extracellular matrix interactions. Here we review the major insights into invasiveness in head and neck cancer provided by high-throughput molecular profiling approaches.

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.

Privacy-preserving genomic testing in the clinic: a model using HIV treatment

PURPOSE The implementation of genomic-based medicine is hindered by unresolved questions regarding data privacy and delivery of interpreted results to health-care practitioners. We used DNA-based prediction of HIV-related outcomes as a model to explore critical issues in clinical genomics. METHODS We genotyped 4,149 markers in HIV-positive individuals. Variants allowed for prediction of 17 traits relevant to HIV medical care, inference of patient ancestry, and imputation of human leukocyte antigen (HLA) types. Genetic data were processed under a privacy-preserving framework using homomorphic encryption, and clinical reports describing potentially actionable results were delivered to health-care providers. RESULTS A total of 230 patients were included in the study. We demonstrated the feasibility of encrypting a large number of genetic markers, inferring patient ancestry, computing monogenic and polygenic trait risks, and reporting results under privacy-preserving conditions. The average execution time of a multimarker test on encrypted data was 865 ms on a standard computer. The proportion of tests returning potentially actionable genetic results ranged from 0 to 54%. CONCLUSIONS The model of implementation presented herein informs on strategies to deliver genomic test results for clinical care. Data encryption to ensure privacy helps to build patient trust, a key requirement on the road to genomic-based medicine.Genet Med advance online publication 14 January 2016Genetics in Medicine (2016); doi:10.1038/gim.2015.167.