The DNA methylation landscape of advanced prostate cancer.

Although DNA methylation is a key regulator of gene expression, the comprehensive methylation landscape of metastatic cancer has never been defined. Through whole-genome bisulfite sequencing paired with deep whole-genome and transcriptome sequencing of 100 castration-resistant prostate metastases, we discovered alterations affecting driver genes that were detectable only with integrated whole-genome approaches. Notably, we observed that 22% of tumors exhibited a novel epigenomic subtype associated with hypermethylation and somatic mutations in TET2, DNMT3B, IDH1 and BRAF. We also identified intergenic regions where methylation is associated with RNA expression of the oncogenic driver genes AR, MYC and ERG. Finally, we showed that differential methylation during progression preferentially occurs at somatic mutational hotspots and putative regulatory regions. This study is a large integrated study of whole-genome, whole-methylome and whole-transcriptome sequencing in metastatic cancer that provides a comprehensive overview of the important regulatory role of methylation in metastatic castration-resistant prostate cancer.

Genomic Hallmarks and Structural Variation in Metastatic Prostate Cancer.

While mutations affecting protein-coding regions have been examined across many cancers, structural variants at the genome-wide level are still poorly defined. Through integrative deep whole-genome and -transcriptome analysis of 101 castration-resistant prostate cancer metastases (109X tumor/38X normal coverage), we identified structural variants altering critical regulators of tumorigenesis and progression not detectable by exome approaches. Notably, we observed amplification of an intergenic enhancer region 624 kb upstream of the androgen receptor (AR) in 81% of patients, correlating with increased AR expression. Tandem duplication hotspots also occur near MYC, in lncRNAs associated with post-translational MYC regulation. Classes of structural variations were linked to distinct DNA repair deficiencies, suggesting their etiology, including associations of CDK12 mutation with tandem duplications, TP53 inactivation with inverted rearrangements and chromothripsis, and BRCA2 inactivation with deletions. Together, these observations provide a comprehensive view of how structural variations affect critical regulators in metastatic prostate cancer.

Isaac: ultra-fast whole-genome secondary analysis on Illumina sequencing platforms.

SUMMARY: An ultrafast DNA sequence aligner (Isaac Genome Alignment Software) that takes advantage of high-memory hardware (>48 GB) and variant caller (Isaac Variant Caller) have been developed. We demonstrate that our combined pipeline (Isaac) is four to five times faster than BWA + GATK on equivalent hardware, with comparable accuracy as measured by trio conflict rates and sensitivity. We further show that Isaac is effective in the detection of disease-causing variants and can easily/economically be run on commodity hardware. AVAILABILITY: Isaac has an open source license and can be obtained at https://github.com/sequencing.

BAP1 loss defines a new class of renal cell carcinoma.

The molecular pathogenesis of renal cell carcinoma (RCC) is poorly understood. Whole-genome and exome sequencing followed by innovative tumorgraft analyses (to accurately determine mutant allele ratios) identified several putative two-hit tumor suppressor genes, including BAP1. The BAP1 protein, a nuclear deubiquitinase, is inactivated in 15% of clear cell RCCs. BAP1 cofractionates with and binds to HCF-1 in tumorgrafts. Mutations disrupting the HCF-1 binding motif impair BAP1-mediated suppression of cell proliferation but not deubiquitination of monoubiquitinated histone 2A lysine 119 (H2AK119ub1). BAP1 loss sensitizes RCC cells in vitro to genotoxic stress. Notably, mutations in BAP1 and PBRM1 anticorrelate in tumors (P = 3 × 10(-5)), [corrected] and combined loss of BAP1 and PBRM1 in a few RCCs was associated with rhabdoid features (q = 0.0007). BAP1 and PBRM1 regulate seemingly different gene expression programs, and BAP1 loss was associated with high tumor grade (q = 0.0005). Our results establish the foundation for an integrated pathological and molecular genetic classification of RCC, paving the way for subtype-specific treatments exploiting genetic vulnerabilities.

Genome-wide expression dynamics during mouse embryonic development reveal similarities to Drosophila development.

Gene transcription mediates many vital aspects of mammalian embryonic development. A comprehensive characterization and analysis of the dynamics of gene transcription in the embryo is therefore likely to provide significant insights into the basic mechanisms of this process. We used microarrays to map transcription in the mouse embryo in the important period from embryonic day 8 (e8.0) to postnatal day 1 (p1) during which the bulk of the differentiation and development of organ systems takes place. Analysis of these expression profiles revealed distinct patterns of gene expression which correlate with the differentiation of organs including the nervous system, liver, skin, lungs, and digestive system, among others. Statistical analysis of the data based on Gene Ontology (GO) group annotation showed that specific temporal sequence patterns in gene class utilization across development are very similar to patterns seen during the embryonic development of Drosophila, suggesting conservation of the temporal progression of these processes across 550 million years of evolution. The temporal profiles of gene expression and activation of processes revealed here provide intriguing insights into the mechanisms of mammalian development, embryogenesis, and organogenesis, as well as into the evolution of developmental processes.

Identification of neural genes using Xenopus DNA microarrays.

To isolate novel genes regulating neural induction, we used a DNA microarray approach. As neural induction is thought to occur by means of the inhibition of bone morphogenetic protein (BMP) signaling, BMP signaling was inhibited in ectodermal cells by overexpression of a dominant-negative receptor. RNAs were isolated from control animal cap explants and from dominant-negative BMP receptor expressing animal caps and subjected to a microarray experiment using newly generated high-density Xenopus DNA microarray chips representing over 17,000 unigenes. We have identified 77 genes that are induced in animal caps after inhibition of BMP signaling, and all of these genes were subjected to whole-mount in situ hybridization analysis. Thirty-two genes showed specific expression in neural tissues. Of the 32, 14 genes have never been linked to neural induction. Two genes that are highly induced by BMP inhibition are inhibitors of Wnt signaling, suggesting that a key step in neural induction is to produce Wnt antagonists to promote anterior neural plate development. Our current analysis also proves that a microarray approach is useful in identifying novel candidate factors involved in neural induction and patterning.

Transcriptional profiling of the heart reveals chamber-specific gene expression patterns.

Cardiac chamber-specific gene expression is critical for the normal development and function of the heart. To investigate the genetic basis of cardiac anatomical specialization, we have undertaken a nearly genome-wide transcriptional profiling of the four heart chambers and the interventricular septum. Rigorous statistical analysis has allowed the identification of known and novel members of gene families that are felt to be important in cardiac development and function, including LIM proteins, homeobox proteins, wnt and T-box pathway proteins, as well as structural proteins like actins and myosins. In addition, these studies have allowed the identification of thousands of additional differentially expressed genes, for which there is little structural or functional information. Clustering of genes with known and unknown functions provides insights into signaling pathways that are essential for development and maintenance of chamber-specific features. To facilitate future research in this area, a searchable internet database has been constructed that allows study of the chamber-specific expression of any gene represented on this comprehensive microarray. It is anticipated that further study of genes identified through this effort will provide insights into the specialization of heart chamber tissues, and their specific roles in cardiac development, aging, and disease.