Caveolin-1-driven membrane remodelling regulates hnRNPK-mediated exosomal microRNA sorting in cancer.

BACKGROUND: Caveolae proteins play diverse roles in cancer development and progression. In prostate cancer, non-caveolar caveolin-1 (CAV1) promotes metastasis, while CAVIN1 attenuates CAV1-induced metastasis. Here, we unveil a novel mechanism linking CAV1 to selective loading of exosomes with metastasis-promoting microRNAs. RESULTS: We identify hnRNPK as a CAV1-regulated microRNA binding protein. In the absence of CAVIN1, non-caveolar CAV1 drives localisation of hnRPNK to multi-vesicular bodies (MVBs), recruiting AsUGnA motif-containing miRNAs and causing their release within exosomes. This process is dependent on the lipid environment of membranes as shown by cholesterol depletion using methyl-ß-cyclodextrin or by treatment with n-3 polyunsaturated fatty acids. Consistent with a role in bone metastasis, knockdown of hnRNPK in prostate cancer PC3 cells abolished the ability of PC3 extracellular vesicles (EV) to induce osteoclastogenesis, and biofluid EV hnRNPK is elevated in metastatic prostate and colorectal cancer. CONCLUSIONS: Taken together, these results support a novel pan-cancer mechanism for CAV1-driven exosomal release of hnRNPK and associated miRNA in metastasis, which is modulated by the membrane lipid environment.

Genome-wide gene expression changes in postpartum depression point towards an altered immune landscape.

Maternal postpartum depression (PPD) is a significant public health concern due to the severe negative impact on maternal and child health and well-being. In this study, we aimed to identify genes associated with PPD. To do this, we investigated genome-wide gene expression profiles of pregnant women during their third trimester of pregnancy and tested the association of gene expression with perinatal depressive symptoms. A total of 137 women from a cohort from the University of North Carolina, USA were assessed. The main phenotypes analysed were Edinburgh Postnatal Depression Scale (EPDS) scores at 2 months postpartum and PPD (binary yes/no) based on an EPDS cutoff of 10. Illumina NextSeq500/550 transcriptomic sequencing from whole blood was analysed using the edgeR package. We identified 71 genes significantly associated with postpartum depression scores at 2 months, after correction for multiple testing at 5% FDR. These included several interesting candidates including TNFRSF17, previously reported to be significantly upregulated in women with PPD and MMP8, a matrix metalloproteinase gene, associated with depression in a genome-wide association study. Functional annotation of differentially expressed genes revealed an enrichment of immune response-related biological processes. Additional analysis of genes associated with changes in depressive symptoms from recruitment to 2 months postpartum identified 66 genes significant at an FDR of 5%. Of these genes, 33 genes were also associated with depressive symptoms at 2 months postpartum. Comparing the results with previous studies, we observed that 15.4% of genes associated with PPD in this study overlapped with 700 core maternal genes that showed significant gene expression changes across multiple brain regions (P = 7.9e-05) and 29-53% of the genes were also associated with estradiol changes in a pharmacological model of depression (P values range = 1.2e-4-2.1e-14). In conclusion, we identified novel genes and validated genes previously associated with oestrogen sensitivity in PPD. These results point towards the role of an altered immune transcriptomic landscape as a vulnerability factor for PPD.

The NK cell granule protein NKG7 regulates cytotoxic granule exocytosis and inflammation.

Immune-modulating therapies have revolutionized the treatment of chronic diseases, particularly cancer. However, their success is restricted and there is a need to identify new therapeutic targets. Here, we show that natural killer cell granule protein 7 (NKG7) is a regulator of lymphocyte granule exocytosis and downstream inflammation in a broad range of diseases. NKG7 expressed by CD4+ and CD8+ T cells played key roles in promoting inflammation during visceral leishmaniasis and malaria-two important parasitic diseases. Additionally, NKG7 expressed by natural killer cells was critical for controlling cancer initiation, growth and metastasis. NKG7 function in natural killer and CD8+ T cells was linked with their ability to regulate the translocation of CD107a to the cell surface and kill cellular targets, while NKG7 also had a major impact on CD4+ T cell activation following infection. Thus, we report a novel therapeutic target expressed on a range of immune cells with functions in different immune responses.

Cumulative influence of parity-related genomic changes in multiple sclerosis.

Pregnancy reduces the frequency of relapses in Multiple Sclerosis (MS) and parity also has a beneficial long term effect on disease outcome. We aimed to uncover the biological mechanisms underlying the beneficial long-term effects of parity in MS. Genome-wide gene expression revealed 574 genes associated with parity; 38.3% showed significant DNA methylation changes (enrichment p = 0.029). These genes overlapped with previous MS genes in humans and a rat MS model and were overrepresented within axon guidance (P = 1.6e-05), developmental biology (P = 0.0094) and cell-cell communication (P = 0.019) pathways. This gene regulation could provide a basis for a protective effect of parity on the long-term outcome of MS.

Long Noncoding RNAs CUPID1 and CUPID2 Mediate Breast Cancer Risk at 11q13 by Modulating the Response to DNA Damage.

Breast cancer risk is strongly associated with an intergenic region on 11q13. We have previously shown that the strongest risk-associated SNPs fall within a distal enhancer that regulates CCND1. Here, we report that, in addition to regulating CCND1, this enhancer regulates two estrogen-regulated long noncoding RNAs, CUPID1 and CUPID2. We provide evidence that the risk-associated SNPs are associated with reduced chromatin looping between the enhancer and the CUPID1 and CUPID2 bidirectional promoter. We further show that CUPID1 and CUPID2 are predominantly expressed in hormone-receptor-positive breast tumors and play a role in modulating pathway choice for the repair of double-strand breaks. These data reveal a mechanism for the involvement of this region in breast cancer.

Shutdown of HIV-1 Transcription in T Cells by Nullbasic, a Mutant Tat Protein.

UNLABELLED: Nullbasic is a derivative of the HIV-1 transactivator of transcription (Tat) protein that strongly inhibits HIV-1 replication in lymphocytes. Here we show that lentiviral vectors that constitutively express a Nullbasic-ZsGreen1 (NB-ZSG1) fusion protein by the eEF1α promoter led to robust long-term inhibition of HIV-1 replication in Jurkat cells. Although Jurkat-NB-ZSG1 cells were infected by HIV-1, no virus production could be detected and addition of phorbol ester 12-myristate 13-acetate (PMA) and JQ1 had no effect, while suberanilohydroxamic acid (SAHA) modestly stimulated virus production but at levels 300-fold lower than those seen in HIV-1-infected Jurkat-ZSG1 cells. Virus replication was not recovered by coculture of HIV-1-infected Jurkat-NB-ZSG1 cells with uninfected Jurkat cells. Latently infected Jurkat latent 6.3 and ACH2 cells treated with latency-reversing agents produced measurable viral capsid (CA), but little or none was made when they expressed NB-ZSG1. When Jurkat cells chronically infected with HIV-1 were transduced with lentiviral virus-like particles conveying NB-ZSG1, a >3-log reduction in CA production was observed. Addition of PMA increased virus CA production but at levels 500-fold lower than those seen in nontransduced Jurkat cells. Transcriptome sequencing analysis confirmed that HIV-1 mRNA was strongly inhibited by NB-ZSG1 but indicated that full-length viral mRNA was made. Analysis of HIV-1-infected Jurkat cells expressing NB-ZSG1 by chromatin immunoprecipitation assays indicated that recruitment of RNA polymerase II (RNAPII) and histone 3 lysine 9 acetylation were inhibited. The reduction of HIV-1 promoter-associated RNAPII and epigenetic changes in viral nucleosomes indicate that Nullbasic can inhibit HIV-1 replication by enforcing viral silencing in cells. IMPORTANCE: HIV-1 infection is effectively controlled by antiviral therapy that inhibits virus replication and reduces measurable viral loads in patients below detectable levels. However, therapy interruption leads to viral rebound due to latently infected cells that serve as a source of continued viral infection. Interest in strategies leading to a functional cure of HIV infection by permanent viral suppression, which may be achievable, is growing. Here we show that a mutant form of the HIV-1 Tat protein, referred to as Nullbasic, can inhibit HIV-1 transcription in infected Jurkat T cell to undetectable levels. Analysis shows that Nullbasic alters the epigenetic state of the HIV-1 long terminal repeat promoter, inhibiting its association with RNA polymerase II. This study indicates that key cellular proteins and pathways targeted here can silence HIV-1 transcription. Further elucidation could lead to functional-cure strategies by suppression of HIV transcription, which may be achievable by a pharmacological method.

Genomic analyses identify molecular subtypes of pancreatic cancer.

Integrated genomic analysis of 456 pancreatic ductal adenocarcinomas identified 32 recurrently mutated genes that aggregate into 10 pathways: KRAS, TGF-ß, WNT, NOTCH, ROBO/SLIT signalling, G1/S transition, SWI-SNF, chromatin modification, DNA repair and RNA processing. Expression analysis defined 4 subtypes: (1) squamous; (2) pancreatic progenitor; (3) immunogenic; and (4) aberrantly differentiated endocrine exocrine (ADEX) that correlate with histopathological characteristics. Squamous tumours are enriched for TP53 and KDM6A mutations, upregulation of the TP63∆N transcriptional network, hypermethylation of pancreatic endodermal cell-fate determining genes and have a poor prognosis. Pancreatic progenitor tumours preferentially express genes involved in early pancreatic development (FOXA2/3, PDX1 and MNX1). ADEX tumours displayed upregulation of genes that regulate networks involved in KRAS activation, exocrine (NR5A2 and RBPJL), and endocrine differentiation (NEUROD1 and NKX2-2). Immunogenic tumours contained upregulated immune networks including pathways involved in acquired immune suppression. These data infer differences in the molecular evolution of pancreatic cancer subtypes and identify opportunities for therapeutic development.

Integrated genomic and transcriptomic analysis of human brain metastases identifies alterations of potential clinical significance.

Treatment options for patients with brain metastases (BMs) have limited efficacy and the mortality rate is virtually 100%. Targeted therapy is critically under-utilized, and our understanding of mechanisms underpinning metastatic outgrowth in the brain is limited. To address these deficiencies, we investigated the genomic and transcriptomic landscapes of 36 BMs from breast, lung, melanoma and oesophageal cancers, using DNA copy-number analysis and exome- and RNA-sequencing. The key findings were as follows. (a) Identification of novel candidates with possible roles in BM development, including the significantly mutated genes DSC2, ST7, PIK3R1 and SMC5, and the DNA repair, ERBB-HER signalling, axon guidance and protein kinase-A signalling pathways. (b) Mutational signature analysis was applied to successfully identify the primary cancer type for two BMs with unknown origins. (c) Actionable genomic alterations were identified in 31/36 BMs (86%); in one case we retrospectively identified ERBB2 amplification representing apparent HER2 status conversion, then confirmed progressive enrichment for HER2-positivity across four consecutive metastatic deposits by IHC and SISH, resulting in the deployment of HER2-targeted therapy for the patient. (d) In the ERBB/HER pathway, ERBB2 expression correlated with ERBB3 (r(2) = 0.496; p < 0.0001) and HER3 and HER4 were frequently activated in an independent cohort of 167 archival BM from seven primary cancer types: 57.6% and 52.6% of cases were phospho-HER3(Y1222) or phospho-HER4(Y1162) membrane-positive, respectively. The HER3 ligands NRG1/2 were barely detectable by RNAseq, with NRG1 (8p12) genomic loss in 63.6% breast cancer-BMs, suggesting a microenvironmental source of ligand. In summary, this is the first study to characterize the genomic landscapes of BM. The data revealed novel candidates, potential clinical applications for genomic profiling of resectable BMs, and highlighted the possibility of therapeutically targeting HER3, which is broadly over-expressed and activated in BMs, independent of primary site and systemic therapy.

Whole genomes redefine the mutational landscape of pancreatic cancer.

Pancreatic cancer remains one of the most lethal of malignancies and a major health burden. We performed whole-genome sequencing and copy number variation (CNV) analysis of 100 pancreatic ductal adenocarcinomas (PDACs). Chromosomal rearrangements leading to gene disruption were prevalent, affecting genes known to be important in pancreatic cancer (TP53, SMAD4, CDKN2A, ARID1A and ROBO2) and new candidate drivers of pancreatic carcinogenesis (KDM6A and PREX2). Patterns of structural variation (variation in chromosomal structure) classified PDACs into 4 subtypes with potential clinical utility: the subtypes were termed stable, locally rearranged, scattered and unstable. A significant proportion harboured focal amplifications, many of which contained druggable oncogenes (ERBB2, MET, FGFR1, CDK6, PIK3R3 and PIK3CA), but at low individual patient prevalence. Genomic instability co-segregated with inactivation of DNA maintenance genes (BRCA1, BRCA2 or PALB2) and a mutational signature of DNA damage repair deficiency. Of 8 patients who received platinum therapy, 4 of 5 individuals with these measures of defective DNA maintenance responded.

A workflow to increase verification rate of chromosomal structural rearrangements using high-throughput next-generation sequencing.

Somatic rearrangements, which are commonly found in human cancer genomes, contribute to the progression and maintenance of cancers. Conventionally, the verification of somatic rearrangements comprises many manual steps and Sanger sequencing. This is labor intensive when verifying a large number of rearrangements in a large cohort. To increase the verification throughput, we devised a high-throughput workflow that utilizes benchtop next-generation sequencing and in-house bioinformatics tools to link the laboratory processes. In the proposed workflow, primers are automatically designed. PCR and an optional gel electrophoresis step to confirm the somatic nature of the rearrangements are performed. PCR products of somatic events are pooled for Ion Torrent PGM and/or Illumina MiSeq sequencing, the resulting sequence reads are assembled into consensus contigs by a consensus assembler, and an automated BLAT is used to resolve the breakpoints to base level. We compared sequences and breakpoints of verified somatic rearrangements between the conventional and high-throughput workflow. The results showed that next-generation sequencing methods are comparable to conventional Sanger sequencing. The identified breakpoints obtained from next-generation sequencing methods were highly accurate and reproducible. Furthermore, the proposed workflow allows hundreds of events to be processed in a shorter time frame compared with the conventional workflow.

Imperfect centered miRNA binding sites are common and can mediate repression of target mRNAs.

BACKGROUND: MicroRNAs (miRNAs) bind to mRNAs and target them for translational inhibition or transcriptional degradation. It is thought that most miRNA-mRNA interactions involve the seed region at the 5' end of the miRNA. The importance of seed sites is supported by experimental evidence, although there is growing interest in interactions mediated by the central region of the miRNA, termed centered sites. To investigate the prevalence of these interactions, we apply a biotin pull-down method to determine the direct targets of ten human miRNAs, including four isomiRs that share centered sites, but not seeds, with their canonical partner miRNAs. RESULTS: We confirm that miRNAs and their isomiRs can interact with hundreds of mRNAs, and that imperfect centered sites are common mediators of miRNA-mRNA interactions. We experimentally demonstrate that these sites can repress mRNA activity, typically through translational repression, and are enriched in regions of the transcriptome bound by AGO. Finally, we show that the identification of imperfect centered sites is unlikely to be an artifact of our protocol caused by the biotinylation of the miRNA. However, the fact that there was a slight bias against seed sites in our protocol may have inflated the apparent prevalence of centered site-mediated interactions. CONCLUSIONS: Our results suggest that centered site-mediated interactions are much more frequent than previously thought. This may explain the evolutionary conservation of the central region of miRNAs, and has significant implications for decoding miRNA-regulated genetic networks, and for predicting the functional effect of variants that do not alter protein sequence.

Somatic point mutation calling in low cellularity tumors.

Somatic mutation calling from next-generation sequencing data remains a challenge due to the difficulties of distinguishing true somatic events from artifacts arising from PCR, sequencing errors or mis-mapping. Tumor cellularity or purity, sub-clonality and copy number changes also confound the identification of true somatic events against a background of germline variants. We have developed a heuristic strategy and software (http://www.qcmg.org/bioinformatics/qsnp/) for somatic mutation calling in samples with low tumor content and we show the superior sensitivity and precision of our approach using a previously sequenced cell line, a series of tumor/normal admixtures, and 3,253 putative somatic SNVs verified on an orthogonal platform.

MicroRNA-182-5p targets a network of genes involved in DNA repair.

MicroRNAs are noncoding regulators of gene expression, which act by repressing protein translation and/or degrading mRNA. Many have been shown to drive tumorigenesis in cancer, but functional studies to understand their mode of action are typically limited to single-target genes. In this study, we use synthetic biotinylated miRNA to pull down endogenous targets of miR-182-5p. We identified more than 1000 genes as potential targets of miR-182-5p, most of which have a known function in pathways underlying tumor biology. Specifically, functional enrichment analysis identified components of both the DNA damage response pathway and cell cycle to be highly represented in this target cohort. Experimental validation confirmed that miR-182-5p-mediated disruption of the homologous recombination (HR) pathway is a consequence of its ability to target multiple components in that pathway. Although there is a strong enrichment for the cell cycle ontology, we do not see primary proliferative defects as a consequence of miR-182-5p overexpression. We highlight targets that could be responsible for miR-182-5p-mediated disruption of other biological processes attributed in the literature so far. Finally, we show that miR-182-5p is highly expressed in a panel of human breast cancer samples, highlighting its role as a potential oncomir in breast cancer.

Pancreatic cancer genomes reveal aberrations in axon guidance pathway genes.

Pancreatic cancer is a highly lethal malignancy with few effective therapies. We performed exome sequencing and copy number analysis to define genomic aberrations in a prospectively accrued clinical cohort (n = 142) of early (stage I and II) sporadic pancreatic ductal adenocarcinoma. Detailed analysis of 99 informative tumours identified substantial heterogeneity with 2,016 non-silent mutations and 1,628 copy-number variations. We define 16 significantly mutated genes, reaffirming known mutations (KRAS, TP53, CDKN2A, SMAD4, MLL3, TGFBR2, ARID1A and SF3B1), and uncover novel mutated genes including additional genes involved in chromatin modification (EPC1 and ARID2), DNA damage repair (ATM) and other mechanisms (ZIM2, MAP2K4, NALCN, SLC16A4 and MAGEA6). Integrative analysis with in vitro functional data and animal models provided supportive evidence for potential roles for these genetic aberrations in carcinogenesis. Pathway-based analysis of recurrently mutated genes recapitulated clustering in core signalling pathways in pancreatic ductal adenocarcinoma, and identified new mutated genes in each pathway. We also identified frequent and diverse somatic aberrations in genes described traditionally as embryonic regulators of axon guidance, particularly SLIT/ROBO signalling, which was also evident in murine Sleeping Beauty transposon-mediated somatic mutagenesis models of pancreatic cancer, providing further supportive evidence for the potential involvement of axon guidance genes in pancreatic carcinogenesis.

Deep-transcriptome and ribonome sequencing redefines the molecular networks of pluripotency and the extracellular space in human embryonic stem cells.

Recent RNA-sequencing studies have shown remarkable complexity in the mammalian transcriptome. The ultimate impact of this complexity on the predicted proteomic output is less well defined. We have undertaken strand-specific RNA sequencing of multiple cellular RNA fractions (>20 Gb) to uncover the transcriptional complexity of human embryonic stem cells (hESCs). We have shown that human embryonic stem (ES) cells display a high degree of transcriptional diversity, with more than half of active genes generating RNAs that differ from conventional gene models. We found evidence that more than 1000 genes express long 5' and/or extended 3'UTRs, which was confirmed by "virtual Northern" analysis. Exhaustive sequencing of the membrane-polysome and cytosolic/untranslated fractions of hESCs was used to identify RNAs encoding peptides destined for secretion and the extracellular space and to demonstrate preferential selection of transcription complexity for translation in vitro. The impact of this newly defined complexity on known gene-centric network models such as the Plurinet and the cell surface signaling machinery in human ES cells revealed a significant expansion of known transcript isoforms at play, many predicting possible alternative functions based on sequence alterations within key functional domains.

Refining transcriptional programs in kidney development by integration of deep RNA-sequencing and array-based spatial profiling.

BACKGROUND: The developing mouse kidney is currently the best-characterized model of organogenesis at a transcriptional level. Detailed spatial maps have been generated for gene expression profiling combined with systematic in situ screening. These studies, however, fall short of capturing the transcriptional complexity arising from each locus due to the limited scope of microarray-based technology, which is largely based on "gene-centric" models. RESULTS: To address this, the polyadenylated RNA and microRNA transcriptomes of the 15.5 dpc mouse kidney were profiled using strand-specific RNA-sequencing (RNA-Seq) to a depth sufficient to complement spatial maps from pre-existing microarray datasets. The transcriptional complexity of RNAs arising from mouse RefSeq loci was catalogued; including 3568 alternatively spliced transcripts and 532 uncharacterized alternate 3' UTRs. Antisense expressions for 60% of RefSeq genes was also detected including uncharacterized non-coding transcripts overlapping kidney progenitor markers, Six2 and Sall1, and were validated by section in situ hybridization. Analysis of genes known to be involved in kidney development, particularly during mesenchymal-to-epithelial transition, showed an enrichment of non-coding antisense transcripts extended along protein-coding RNAs. CONCLUSION: The resulting resource further refines the transcriptomic cartography of kidney organogenesis by integrating deep RNA sequencing data with locus-based information from previously published expression atlases. The added resolution of RNA-Seq has provided the basis for a transition from classical gene-centric models of kidney development towards more accurate and detailed "transcript-centric" representations, which highlights the extent of transcriptional complexity of genes that direct complex development events.

MicroRNAs and their isomiRs function cooperatively to target common biological pathways.

BACKGROUND: Variants of microRNAs (miRNAs), called isomiRs, are commonly reported in deep-sequencing studies; however, the functional significance of these variants remains controversial. Observational studies show that isomiR patterns are non-random, hinting that these molecules could be regulated and therefore functional, although no conclusive biological role has been demonstrated for these molecules. RESULTS: To assess the biological relevance of isomiRs, we have performed ultra-deep miRNA-seq on ten adult human tissues, and created an analysis pipeline called miRNA-MATE to align, annotate, and analyze miRNAs and their isomiRs. We find that isomiRs share sequence and expression characteristics with canonical miRNAs, and are generally strongly correlated with canonical miRNA expression. A large proportion of isomiRs potentially derive from AGO2 cleavage independent of Dicer. We isolated polyribosome-associated mRNA, captured the mRNA-bound miRNAs, and found that isomiRs and canonical miRNAs are equally associated with translational machinery. Finally, we transfected cells with biotinylated RNA duplexes encoding isomiRs or their canonical counterparts and directly assayed their mRNA targets. These studies allow us to experimentally determine genome-wide mRNA targets, and these experiments showed substantial overlap in functional mRNA networks suppressed by both canonical miRNAs and their isomiRs. CONCLUSIONS: Together, these results find isomiRs to be biologically relevant and functionally cooperative partners of canonical miRNAs that act coordinately to target pathways of functionally related genes. This work exposes the complexity of the miRNA-transcriptome, and helps explain a major miRNA paradox: how specific regulation of biological processes can occur when the specificity of miRNA targeting is mediated by only 6 to 11 nucleotides.

Gene-specific methylation: potential markers for colorectal cancer.

PURPOSE: Aberrant methylation of the promoter region is associated with silencing of many genes in neoplasia. CpG island methylation is an epigenetic mechanism for transcriptional silencing that occurs at various stages of colon tumorigenesis. In this study, we tested the promoter methylation and expression of seven genes from various pathways of DNA repair, apoptosis and inflammation, i.e., sFRP1, MLH1, RASSF1A, CDA, v-fgr, LYN-B, and TNFR10d. METHOD: The genes were analyzed by quantitative polymerase chain reaction for the level of gene expression. The promoter methylation status of the genes was studied by methylation-specific polymerase chain reaction. RESULT: The correlation of promoter methylation status with suppressed gene expression patterns suggested a potential role for the silencing these genes in colon cancer progression. CONCLUSION: Promoter methylations of the studied genes could be explored as promising biomarkers for new diagnostic, prognostic and therapeutic targets of colorectal cancer.

The regulated retrotransposon transcriptome of mammalian cells.

Although repetitive elements pervade mammalian genomes, their overall contribution to transcriptional activity is poorly defined. Here, as part of the FANTOM4 project, we report that 6-30% of cap-selected mouse and human RNA transcripts initiate within repetitive elements. Analysis of approximately 250,000 retrotransposon-derived transcription start sites shows that the associated transcripts are generally tissue specific, coincide with gene-dense regions and form pronounced clusters when aligned to full-length retrotransposon sequences. Retrotransposons located immediately 5' of protein-coding loci frequently function as alternative promoters and/or express noncoding RNAs. More than a quarter of RefSeqs possess a retrotransposon in their 3' UTR, with strong evidence for the reduced expression of these transcripts relative to retrotransposon-free transcripts. Finally, a genome-wide screen identifies 23,000 candidate regulatory regions derived from retrotransposons, in addition to more than 2,000 examples of bidirectional transcription. We conclude that retrotransposon transcription has a key influence upon the transcriptional output of the mammalian genome.