Retrospective clinical trial experimentally validates glioblastoma genome-wide pattern of DNA copy-number alterations predictor of survival.

Modeling of genomic profiles from the Cancer Genome Atlas (TCGA) by using recently developed mathematical frameworks has associated a genome-wide pattern of DNA copy-number alterations with a shorter, roughly one-year, median survival time in glioblastoma (GBM) patients. Here, to experimentally test this relationship, we whole-genome sequenced DNA from tumor samples of patients. We show that the patients represent the U.S. adult GBM population in terms of most normal and disease phenotypes. Intratumor heterogeneity affects ≈ 11 % and profiling technology and reference human genome specifics affect <1% of the classifications of the tumors by the pattern, where experimental batch effects normally reduce the reproducibility, i.e., precision, of classifications based upon between one to a few hundred genomic loci by >30%. With a 2.25-year Kaplan-Meier median survival difference, a 3.5 univariate Cox hazard ratio, and a 0.78 concordance index, i.e., accuracy, the pattern predicts survival better than and independent of age at diagnosis, which has been the best indicator since 1950. The prognostic classification by the pattern may, therefore, help to manage GBM pseudoprogression. The diagnostic classification may help drugs progress to regulatory approval. The therapeutic predictions, of previously unrecognized targets that are correlated with survival, may lead to new drugs. Other methods missed this relationship in the roughly 3B-nucleotide genomes of the small, order of magnitude of 100, patient cohorts, e.g., from TCGA. Previous attempts to associate GBM genotypes with patient phenotypes were unsuccessful. This is a proof of principle that the frameworks are uniquely suitable for discovering clinically actionable genotype-phenotype relationships.

Differential gene expression patterns in vein regions susceptible versus resistant to neointimal hyperplasia.

Arteriovenous hemodialysis graft (AVG) stenosis results in thrombosis and AVG failure, but prevention of stenosis has been unsuccessful due in large part to our limited understanding of the molecular processes involved in neointimal hyperplasia (NH) formation. AVG stenosis develops chiefly as a consequence of highly localized NH formation in the vein-graft anastomosis region. Surprisingly, the vein region just downstream of the vein-graft anastomosis (herein termed proximal vein region) is relatively resistant to NH. We hypothesized that the gene expression profiles of the NH-prone and NH-resistant regions will be different from each other after graft placement, and analysis of their genomic profiles may yield potential therapeutic targets to prevent AVG stenosis. To test this, we evaluated the vein-graft anastomosis (NH-prone) and proximal vein (NH-resistant) regions in a porcine model of AVG stenosis with a porcine microarray. Gene expression changes in these two distinct vein regions, relative to the gene expression in unoperated control veins, were examined at early (5 days) and later (14 days) time points following graft placement. Global genomic changes were much greater in the NH-prone region than in the NH-resistant region at both time points. In the NH-prone region, genes related to regulation of cell proliferation and osteo-/chondrogenic vascular remodeling were most enriched among the significantly upregulated genes, and genes related to smooth muscle phenotype were significantly downregulated. These results provide insights into the spatial and temporal genomic modulation underlying NH formation in AVG and suggest potential therapeutic strategies to prevent and/or limit AVG stenosis.

Temporal Dysynchrony in brain connectivity gene expression following hypoxia.

BACKGROUND: Despite the fundamental biological importance and clinical relevance of characterizing the effects of chronic hypoxia exposure on central nervous system (CNS) development, the changes in gene expression from hypoxia are unknown. It is not known if there are unifying principles, properties, or logic in the response of the developing CNS to hypoxic exposure. Here, we use the small vertebrate zebrafish (Danio rerio) to study the effects of hypoxia on connectivity gene expression across development. We perform transcriptional profiling at high temporal resolution to systematically determine and then experimentally validate the response of CNS connectivity genes to hypoxia exposure. RESULTS: We characterized mRNA changes during development, comparing the effects of chronic hypoxia exposure at different time-points. We focused on changes in expression levels of a subset of 1270 genes selected for their roles in development of CNS connectivity, including axon pathfinding and synapse formation. We found that the majority of CNS connectivity genes were unaffected by hypoxia. However, for a small subset of genes hypoxia significantly affected their gene expression profiles. In particular, hypoxia appeared to affect both the timing and levels of expression, including altering expression of interacting gene pairs in a fashion that would potentially disrupt normal function. CONCLUSIONS: Overall, our study identifies the response of CNS connectivity genes to hypoxia exposure during development. While for most genes hypoxia did not significantly affect expression, for a subset of genes hypoxia changed both levels and timing of expression. Importantly, we identified that some genes with interacting proteins, for example receptor/ligand pairs, had dissimilar responses to hypoxia that would be expected to interfere with their function. The observed dysynchrony of gene expression could impair the development of normal CNS connectivity maps.

Spontaneous shaker rat mutant - a new model for X-linked tremor/ataxia.

The shaker rat is an X-linked recessive spontaneous model of progressive Purkinje cell (PC) degeneration exhibiting a shaking ataxia and wide stance. Generation of Wistar Furth (WF)/Brown Norwegian (BN) F1 hybrids and genetic mapping of F2 sib-sib offspring using polymorphic markers narrowed the candidate gene region to 26 Mbp denoted by the last recombinant genetic marker DXRat21 at 133 Mbp to qter (the end of the long arm). In the WF background, the shaker mutation has complete penetrance, results in a stereotypic phenotype and there is a narrow window for age of disease onset; by contrast, the F2 hybrid phenotype was more varied, with a later age of onset and likely non-penetrance of the mutation. By deep RNA-sequencing, five variants were found in the candidate region; four were novel without known annotation. One of the variants caused an arginine (R) to cysteine (C) change at codon 35 of the ATPase, Ca(2+) transporting, plasma membrane 3 (Atp2b3) gene encoding PMCA3 that has high expression in the cerebellum. The variant was well supported by hundreds of overlapping reads, and was found in 100% of all affected replicas and 0% of the wild-type (WT) replicas. The mutation segregated with disease in all affected animals and the amino acid change was found in an evolutionarily conserved region of PMCA3. Despite strong genetic evidence for pathogenicity, in vitro analyses of PMCA3(R35C) function did not show any differences to WT PMCA3. Because Atp2b3 mutation leads to congenital ataxia in humans, the identified Atp2b3 missense change in the shaker rat presents a good candidate for the shaker rat phenotype based on genetic criteria, but cannot yet be considered a definite pathogenic variant owing to lack of functional changes.

Zinc-dependent regulation of the Adh1 antisense transcript in fission yeast.

In yeast, Adh1 (alcohol dehydrogenase 1) is an abundant zinc-binding protein that is required for the conversion of acetaldehyde to ethanol. Through transcriptome profiling of the Schizosaccharomyces pombe genome, we identified a natural antisense transcript at the adh1 locus that is induced in response to zinc limitation. This antisense transcript (adh1AS) shows a reciprocal expression pattern to that of the adh1 mRNA partner. In this study, we show that increased expression of the adh1AS transcript in zinc-limited cells is necessary for the repression of adh1 gene expression and that the increased level of the adh1AS transcript in zinc-limited cells is a result of two mechanisms. At the transcriptional level, the adh1AS transcript is expressed at a high level in zinc-limited cells. In addition to this transcriptional control, adh1AS transcripts preferentially accumulate in zinc-limited cells when the adh1AS transcript is expressed from a constitutive promoter. This secondary mechanism requires the simultaneous expression of adh1. Our studies reveal how multiple mechanisms can synergistically control the ratio of sense to antisense transcripts and highlight a novel mechanism by which adh1 gene expression can be controlled by cellular zinc availability.

Chemical genetic screen reveals a role for desmosomal adhesion in mammary branching morphogenesis.

During the process of branching morphogenesis, the mammary gland undergoes distinct phases of remodeling to form an elaborate ductal network that ultimately produces and delivers milk to newborn animals. These developmental events rely on tight regulation of critical cellular pathways, many of which are probably disrupted during initiation and progression of breast cancer. Transgenic mouse and in vitro organoid models previously identified growth factor signaling as a key regulator of mammary branching, but the functional downstream targets of these pathways remain unclear. Here, we used purified primary mammary epithelial cells stimulated with fibroblast growth factor-2 (FGF2) to model mammary branching morphogenesis in vitro. We employed a forward chemical genetic approach to identify modulators of this process and describe a potent compound, 1023, that blocks FGF2-induced branching. In primary mammary epithelial cells, we used lentivirus-mediated knockdown of the aryl hydrocarbon receptor (AHR) to demonstrate that 1023 acts through AHR to block branching. Using 1023 as a tool, we identified desmosomal adhesion as a novel target of AHR signaling and show that desmosomes are critical for AHR agonists to block branching. Our findings support a functional role for desmosomes during mammary morphogenesis and also in blocking FGF-induced invasion.

Ccdc94 protects cells from ionizing radiation by inhibiting the expression of p53.

DNA double-strand breaks (DSBs) represent one of the most deleterious forms of DNA damage to a cell. In cancer therapy, induction of cell death by DNA DSBs by ionizing radiation (IR) and certain chemotherapies is thought to mediate the successful elimination of cancer cells. However, cancer cells often evolve to evade the cytotoxicity induced by DNA DSBs, thereby forming the basis for treatment resistance. As such, a better understanding of the DSB DNA damage response (DSB-DDR) pathway will facilitate the design of more effective strategies to overcome chemo- and radioresistance. To identify novel mechanisms that protect cells from the cytotoxic effects of DNA DSBs, we performed a forward genetic screen in zebrafish for recessive mutations that enhance the IR-induced apoptotic response. Here, we describe radiosensitizing mutation 7 (rs7), which causes a severe sensitivity of zebrafish embryonic neurons to IR-induced apoptosis and is required for the proper development of the central nervous system. The rs7 mutation disrupts the coding sequence of ccdc94, a highly conserved gene that has no previous links to the DSB-DDR pathway. We demonstrate that Ccdc94 is a functional member of the Prp19 complex and that genetic knockdown of core members of this complex causes increased sensitivity to IR-induced apoptosis. We further show that Ccdc94 and the Prp19 complex protect cells from IR-induced apoptosis by repressing the expression of p53 mRNA. In summary, we have identified a new gene regulating a dosage-sensitive response to DNA DSBs during embryonic development. Future studies in human cancer cells will determine whether pharmacological inactivation of CCDC94 reduces the threshold of the cancer cell apoptotic response.

Antibody detection of translocations in Ewing sarcoma.

The detection of chromosomal translocations has important implications in the diagnosis, prognosis and treatment of patients with cancer. Current approaches to translocation detection have significant shortcomings, including limited sensitivity and/or specificity, and difficulty in application to formalin-fixed paraffin-embedded (FFPE) clinical samples. We developed a new approach called antibody detection of translocations (ADOT) that avoids the shortcomings of current techniques. ADOT combines a transcriptional microarray-based approach with a novel antibody-based detection method. ADOT allows for the accurate and sensitive identification of translocations and provides exon-level information about the fusion transcript. ADOT can detect translocations in poor-quality unprocessed total ribonucleic acid (RNA). Furthermore, the technique is readily generalizable to detect any potential fusion transcript, including previously undescribed fusions. We demonstrate the feasibility of ADOT by examples in which both known and unknown Ewing sarcoma translocations are identified from cell lines, tumour xenografts and FFPE primary tumours. These results demonstrate that ADOT may be an effective approach for translocation analysis in clinical specimens with significant RNA degradation and may offer a novel diagnostic tool for translocation-based cancers.

Deep sequencing for the detection of virus-like sequences in the brains of patients with multiple sclerosis: detection of GBV-C in human brain.

Multiple sclerosis (MS) is a demyelinating disease of unknown origin that affects the central nervous system of an estimated 400,000 Americans. GBV-C or hepatitis G is a flavivirus that is found in the serum of 1-2% of blood donors. It was originally associated with hepatitis, but is now believed to be a relatively non-pathogenic lymphotropic virus. Fifty frozen specimens from the brains of deceased persons affected by MS were obtained along with 15 normal control brain specimens. RNA was extracted and ribosomal RNAs were depleted before sequencing on the Illumina GAII. These 36 bp reads were compared with a non-redundant database derived from the 600,000+ viral sequences in GenBank organized into 4080 taxa. An individual read successfully aligned to the viral database was considered to be a "hit". Normalized MS specimen hit rates for each viral taxon were compared to the distribution of hits in the normal controls. Seventeen MS and 11 control brain extracts were sequenced, yielding 4-10 million sequences ("reads") each. Over-representation of sequence from at least one of 12 viral taxa was observed in 7 of the 17 MS samples. Sequences resembling other viruses previously implicated in the pathogenesis of MS were not significantly enriched in any of the diseased brain specimens. Sequences from GB virus C (GBV-C), a flavivirus not previously isolated from brain, were enriched in one of the MS samples. GBV-C in this brain specimen was confirmed by specific amplification in this single MS brain specimen, but not in the 30 other MS brain samples available. The entire 9.4 kb sequence of this GBV-C isolate is reported here. This study shows the feasibility of deep sequencing for the detection of occult viral infections in the brains of deceased persons with MS. The first isolation of GBV-C from human brain is reported here.

Novel venom peptides from the cone snail Conus pulicarius discovered through next-generation sequencing of its venom duct transcriptome.

The venom peptides (i.e., conotoxins or conopeptides) that species in the genus Conus collectively produce are remarkably diverse, estimated to be around 50,000 to 140,000, but the pace of discovery and characterization of these peptides have been rather slow. To date, only a minor fraction have been identified and studied. However, the advent of next-generation DNA sequencing technologies has opened up opportunities for expediting the exploration of this diversity. The whole transcriptome of a venom duct from the vermivorous marine snail C. pulicarius was sequenced using the 454 sequencing platform. Analysis of the data set resulted in the identification of over eighty unique putative conopeptide sequences, the highest number discovered so far from a Conus venom duct transcriptome. More importantly, majority of the sequences were potentially novel, many with unexpected structural features, hinting at the vastness of the diversity of Conus venom peptides that remains to be explored. The sequences represented at least 14 major superfamilies/types (disulfide- and non-disulfide-rich), indicating the structural and functional diversity of conotoxins in the venom of C. pulicarius. In addition, the contryphans were surprisingly more diverse than what is currently known. Comparative analysis of the O-superfamily sequences also revealed insights into the complexity of the processes that drive the evolution and diversification of conotoxins.

Targeting the PyMT Oncogene to Diverse Mammary Cell Populations Enhances Tumor Heterogeneity and Generates Rare Breast Cancer Subtypes.

Human breast cancer is a heterogeneous disease composed of different histologies and molecular subtypes, many of which are not replicated in animal models. Here, we report a mouse model of breast cancer that generates unique tumor histologies including tubular, adenosquamous, and lipid-rich carcinomas. Utilizing a nononcogenic variant of polyoma middle T oncogene (PyMT) that requires a spontaneous base-pair deletion to transform cells, in conjunction with lentiviral transduction and orthotopic transplantation of primary mammary epithelial cells, this model sporadically induces oncogene expression in both the luminal and myoepithelial cell lineages of the normal mouse mammary epithelium. Microarray and hierarchical analyses using an intrinsic subtype gene set revealed that lentiviral PyMT generates both luminal and basal-like tumors. Cumulatively, these results show that low-level expression of PyMT in a broad range of cell types significantly increases tumor heterogeneity and establishes a mouse model of several rare human breast cancer subtypes.

Tumor grafts derived from women with breast cancer authentically reflect tumor pathology, growth, metastasis and disease outcomes.

Development and preclinical testing of new cancer therapies is limited by the scarcity of in vivo models that authentically reproduce tumor growth and metastatic progression. We report new models for breast tumor growth and metastasis in the form of transplantable tumors derived directly from individuals undergoing treatment for breast cancer. These tumor grafts illustrate the diversity of human breast cancer and maintain essential features of the original tumors, including metastasis to specific sites. Co-engraftment of primary human mesenchymal stem cells maintains phenotypic stability of the grafts and increases tumor growth by promoting angiogenesis. We also report that tumor engraftment is a prognostic indicator of disease outcome for women with newly diagnosed breast cancer; orthotopic breast tumor grafting is a step toward individualized models for tumor growth, metastasis and prognosis. This bank of tumor grafts also serves as a publicly available resource for new models in which to study the biology of breast cancer.

Polymorphisms in the NPY2R gene show significant associations with BMI that are additive to FTO, MC4R, and NPFFR2 gene effects.

Neuropeptide Y (NPY) is an appetite hormone that acts centrally to control feeding behavior. The 5' and exon 2 regions of NPY2R, one of five NPY receptor genes, have been weakly and inconsistently implicated with obesity. With the ATG start site of the gene at the beginning of exon 2, single-nucleotide polymorphisms (SNPs) across intron 1 may show stronger associations with obesity than expected. Two 5' SNPs, three intron 1 SNPs, and one synonymous exon 2 SNP were genotyped on 2,985 white Utah subjects. Previously associated FTO, NPY, NPY1R, MC4R, PPARGC1A, OR7D4, and four NPFFR2 SNPs were also genotyped and related to BMI. One NPY2R 5' SNP (rs12649641, P = 0.008), an exon 2 SNP (rs2880415, P = 0.009), and an intron 1 SNP (rs17376826, P = 7 × 10(-6)) were each significantly associated with BMI. All three SNPs, plus FTO (rs9939609, P = 1.5 × 10(-6)) and two NPFFR2 SNPs (rs4129733, P = 3.7 × 10(-13) and rs11940196, 4.2 × 10(-10)) remained significant in a multiple regression additive model. Diplotypes using the estimated haplotypes of NPY2R, NPFFR2, and MC4R were significantly associated with BMI (P = 1.0 × 10(-10), 3.2 × 10(-8), and 1.1 × 10(-4), respectively). Haplotypes of NPY2R, NPFFR2, and MC4R, plus the FTO SNP, explained 9.6% of the BMI variance. SNP effect sizes per allele for the four genes ranged from 0.8 to 3.5 kg/m(2). We conclude that haplotypes containing the rs17376826 SNP in intron 1 of NPY2R have strong associations with BMI, some NPFFR2 haplotypes are strongly protective against or increase risk of obesity, and both NPY2R and NPFFR2 play important roles in obesity predisposition independent of FTO and MC4R.

MicroRNAs and colon and rectal cancer: differential expression by tumor location and subtype.

MicroRNAs are thought to have an impact on cell proliferation, apoptosis, stress responses, maintenance of stem cell potency, and metabolism and are, therefore, important in the carcinogenic process. In this study, we examined 40 colon tumors, 30 rectal tumors, and 30 normal tissue samples (10 proximal colon, 10 distal colon, and 10 rectal paired with cancer cases) to examine miRNA expression profiles in colon and rectal tumors. MiRNA expression levels were adjusted for multiple comparisons; tumor tissue was compared with noncancerous tissue from the same site. A comparison of normal tissue showed 287 unique miRNAs that were significantly differentially expressed at the 1.5-fold level and 73 with over a two-fold difference in expression between colon and rectal tissue. Examination of miRNAs that were significantly differentially expressed at the 1.5-fold level by tumor phenotype showed 143 unique miRNAs differentially expression for microsatellite instability positive (MSI+) colon tumors; 129 unique miRNAs differentially expressed for CpG Island Methylator Phenotype positive (CIMP+) colon tumors; 135 miRNAs were differentially expressed for KRAS2-mutated colon tumors, and 139 miRNAs were differentially expressed for TP53-mutated colon tumors. Similar numbers of differentially expressed miRNAs were observed for rectal tumors, although the miRNAs differentially expressed differed. There were 129 unique miRNAs for CIMP+, 143 unique miRNAs for KRAS2-mutated, and 136 unique miRNAs for TP53-mutated rectal tumors. These results suggest the importance of miRNAs in colorectal cancer and the need for studies that can confirm these results and provide insight into the diet, lifestyle, and genetic factors that influence miRNA expression.

Next generation tools for genomic data generation, distribution, and visualization.

BACKGROUND: With the rapidly falling cost and availability of high throughput sequencing and microarray technologies, the bottleneck for effectively using genomic analysis in the laboratory and clinic is shifting to one of effectively managing, analyzing, and sharing genomic data. RESULTS: Here we present three open-source, platform independent, software tools for generating, analyzing, distributing, and visualizing genomic data. These include a next generation sequencing/microarray LIMS and analysis project center (GNomEx); an application for annotating and programmatically distributing genomic data using the community vetted DAS/2 data exchange protocol (GenoPub); and a standalone Java Swing application (GWrap) that makes cutting edge command line analysis tools available to those who prefer graphical user interfaces. Both GNomEx and GenoPub use the rich client Flex/Flash web browser interface to interact with Java classes and a relational database on a remote server. Both employ a public-private user-group security model enabling controlled distribution of patient and unpublished data alongside public resources. As such, they function as genomic data repositories that can be accessed manually or programmatically through DAS/2-enabled client applications such as the Integrated Genome Browser. CONCLUSIONS: These tools have gained wide use in our core facilities, research laboratories and clinics and are freely available for non-profit use. See http://sourceforge.net/projects/gnomex/, http://sourceforge.net/projects/genoviz/, and http://sourceforge.net/projects/useq.

Reactive astrocytosis, microgliosis and inflammation in rats with neonatal hydrocephalus.

The deleterious effects of hydrocephalus, a disorder that primarily affects children, include reactive astrocytosis, microgliosis and inflammatory responses; however, the roles that these mechanisms play in the pathophysiology of hydrocephalus are still not clear in terms of cytopathology and gene expression. Therefore we have examined neuroinflammation at both the cellular and the molecular levels in an experimental model of neonatal obstructive hydrocephalus. On post-natal day 1, rats received an intracisternal injection of kaolin to induce hydrocephalus; control animals received saline injections. Prior to sacrifice on post-natal day 22, animals underwent magnetic resonance imaging to quantify ventricular enlargement, and the parietal cortex was harvested for analysis. Immunohistochemistry and light microscopy were performed on 5 hydrocephalic and 5 control animals; another set of 5 hydrocephalic and 5 control animals underwent molecular testing with Western blots and a gene microarray. Scoring of immunoreactivity on a 4-point ranking scale for GFAP and Iba-1 demonstrated an increase in reactive astrocytes and reactive microglia respectively in the hydrocephalic animals compared to controls (2.90±0.11 vs. 0.28±0.26; 2.91±0.11 vs. 0.58±0.23, respectively). Western blots confirmed these results. Microarray analysis identified significant (1.5-fold) changes in 1729 of 33,951 genes, including 26 genes out of 185 genes (26/185) in the cytokine-cytokine receptor interaction pathway, antigen processing and presentation pathways (15/66), and the apoptosis pathway (10/69). Collectively, these results demonstrate alterations in normal physiology and an up-regulation of the inflammatory response. These findings lead to a better understanding of neonatal hydrocephalus and begin to form a baseline for future treatments that may reverse these effects.

Radiosensitization of cervical cancer cells via double-strand DNA break repair inhibition.

PURPOSE: LY294002, a phosphatidylinositol 3-kinase (PI3K) inhibitor, has been found to radiosensitize various human cancer cells. However, its potential to act as an effective therapeutic agent is diminished by its toxicity levels. The purposes of this study were to determine the mechanism by which LY294002 radiosensitizes. MATERIALS AND METHODS: Cell growth curves and clonogenic assays were performed with increasing LY294002 exposure times proximate to the radiation dose. Protein levels of downstream PI3K effectors were analyzed. Detection of phosphorylated histone H2AX (gammaH2AX) was used to identify DNA double-strand breaks at various time points post-radiation. RESULTS: LY294002 significantly radiosensitized HeLa cervical cancer cells when administered for just 12 h following radiation. Cell growth curves also decreased with brief LY294002 application. DNA double-strand breaks are typically repaired within 2-6 h following radiation. Interestingly, at 48, 72, and 96 h post-irradiation, gammaH2AX was still significantly elevated in cells radiated in combination with LY294002. Protein expressions of ATM and ATR downstream effectors showed no differences among the treated groups, however, DNA-PK activity was significantly inhibited by LY294002. CONCLUSIONS: These results lead us to conclude that the central mechanism by which LY294002 radiosensitizes is via DNA-PK inhibition which induces DNA double-strand break repair inhibition. We are currently investigating radiosensitization induced by DNA-PK-specific inhibition in efforts to find a less toxic, yet equally effective, chemotherapeutic agent than LY294002.

Dynamic transcriptome of Schizosaccharomyces pombe shown by RNA-DNA hybrid mapping.

We have determined the high-resolution strand-specific transcriptome of the fission yeast S. pombe under multiple growth conditions using a novel RNA-DNA hybridization mapping (HybMap) technique. HybMap uses an antibody against an RNA-DNA hybrid to detect RNA molecules hybridized to a high-density DNA oligonucleotide tiling microarray. HybMap showed exceptional dynamic range and reproducibility, and allowed us to identify strand-specific coding, noncoding and structural RNAs, as well as previously unknown RNAs conserved in distant yeast species. Notably, we found that virtually the entire euchromatic genome (including intergenics) is transcribed, with heterochromatin dampening intergenic transcription. We identified features including large numbers of condition-specific noncoding RNAs, extensive antisense transcription, new properties of antisense transcripts and induced divergent transcription. Furthermore, our HybMap data informed the efficiency and locations of RNA splicing genome-wide. Finally, we observed strand-specific transcription islands around tRNAs at heterochromatin boundaries inside centromeres. Here, we discuss these new features in terms of organism fitness and transcriptome evolution.