Advances in vertebrate genomics have uncovered thousands of loci encoding long noncoding RNAs (lncRNAs). While progress has been made in elucidating the regulatory functions of lncRNAs, little is known about their origins and evolution. Here we explore the contribution of transposable elements (TEs) to the makeup and regulation of lncRNAs in human, mouse, and zebrafish. Surprisingly, TEs occur in more than two thirds of mature lncRNA transcripts and account for a substantial portion of total lncRNA sequence (∼30% in human), whereas they seldom occur in protein-coding transcripts. While TEs contribute less to lncRNA exons than expected, several TE families are strongly enriched in lncRNAs. There is also substantial interspecific variation in the coverage and types of TEs embedded in lncRNAs, partially reflecting differences in the TE landscapes of the genomes surveyed. In human, TE sequences in lncRNAs evolve under greater evolutionary constraint than their non–TE sequences, than their intronic TEs, or than random DNA. Consistent with functional constraint, we found that TEs contribute signals essential for the biogenesis of many lncRNAs, including ∼30,000 unique sites for transcription initiation, splicing, or polyadenylation in human. In addition, we identified ∼35,000 TEs marked as open chromatin located within 10 kb upstream of lncRNA genes. The density of these marks in one cell type correlate with elevated expression of the downstream lncRNA in the same cell type, suggesting that these TEs contribute to cis-regulation. These global trends are recapitulated in several lncRNAs with established functions. Finally a subset of TEs embedded in lncRNAs are subject to RNA editing and predicted to form secondary structures likely important for function. In conclusion, TEs are nearly ubiquitous in lncRNAs and have played an important role in the lineage-specific diversification of vertebrate lncRNA repertoires.

Human endogenous retrovirus subfamily H (HERVH) is a class of transposable elements expressed preferentially in human embryonic stem cells (hESCs). Here, we report that the long terminal repeats of HERVH function as enhancers and that HERVH is a nuclear long noncoding RNA required to maintain hESC identity. Furthermore, HERVH is associated with OCT4, coactivators and Mediator subunits. Together, these results uncover a new role of species-specific transposable elements in hESCs.

Single-cell technologies have enabled the characterization of the tumour microenvironment at unprecedented depth and have revealed vast cellular diversity among tumour cells and their niche. Anti-tumour immunity relies on cell–cell relationships within the tumour microenvironment1,2, yet many single-cell studies lack spatial context and rely on dissociated tissues3. Here we applied imaging mass cytometry to characterize the immunological landscape of 139 high-grade glioma and 46 brain metastasis tumours from patients. Single-cell analysis of more than 1.1 million cells across 389 high-dimensional histopathology images enabled the spatial resolution of immune lineages and activation states, revealing differences in immune landscapes between primary tumours and brain metastases from diverse solid cancers. These analyses revealed cellular neighbourhoods associated with survival in patients with glioblastoma, which we leveraged to identify a unique population of myeloperoxidase (MPO)-positive macrophages associated with long-term survival. Our findings provide insight into the biology of primary and metastatic brain tumours, reinforcing the value of integrating spatial resolution to single-cell datasets to dissect the microenvironmental contexture of cancer.

BackgroundA significant portion of expressed non-coding RNAs in human cells is derived from transposable elements (TEs). Moreover, it has been shown that various long non-coding RNAs (lncRNAs), which come from the human endogenous retrovirus subfamily H (HERVH), are not only expressed but required for pluripotency in human embryonic stem cells (hESCs).ResultsTo identify additional TE-derived functional non-coding transcripts, we generated RNA-seq data from induced pluripotent stem cells (iPSCs) of four primate species (human, chimpanzee, gorilla, and rhesus) and searched for transcripts whose expression was conserved. We observed that about 30% of TE instances expressed in human iPSCs had orthologous TE instances that were also expressed in chimpanzee and gorilla. Notably, our analysis revealed a number of repeat families with highly conserved expression profiles including HERVH but also MER53, which is known to be the source of a placental-specific family of microRNAs (miRNAs). We also identified a number of repeat families from all classes of TEs, including MLT1-type and Tigger families, that contributed a significant amount of sequence to primate lncRNAs whose expression was conserved.ConclusionsTogether, these results describe TE families and TE-derived lncRNAs whose conserved expression patterns can be used to identify what are likely functional TE-derived non-coding transcripts in primate iPSCs.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-017-3568-y) contains supplementary material, which is available to authorized users.

Purpose: Next-generation sequencing studies and CRISPR-Cas9 screens have established mutations in the IFNγ-JAK-STAT pathway as an immune checkpoint inhibitor (ICI) resistance mechanism in a subset of patients with melanoma. We hypothesized ICI resistance mutations in the IFNγ pathway would simultaneously render melanomas susceptible to oncolytic virus (OV) therapy. Experimental Design: Cytotoxicity experiments were performed with a number of OVs on a matched melanoma cell line pair generated from a baseline biopsy and a progressing lesion with complete JAK2 loss from a patient that relapsed on anti-PD-1 therapy, in melanoma lines following JAK1/2 RNA interference (RNAi) and pharmacologic inhibition and in Jak2 knockout (KO) B16-F10 mouse melanomas. Furthermore, we estimated the frequency of genetic alterations in the IFNγ-JAK-STAT pathway in human melanomas. Results: The melanoma line from an anti-PD-1 progressing lesion was 7- and 22-fold more sensitive to the modified OVs, herpes simplex virus 1 (HSV1-dICP0) and vesicular stomatitis virus (VSV-Δ51), respectively, compared with the line from the baseline biopsy. RNAi, JAK1/2 inhibitor studies, and in vivo studies of Jak2 KOs B16-F10 melanomas revealed a significant increase in VSV-Δ51 sensitivity with JAK/STAT pathway inhibition. Our analysis of The Cancer Genome Atlas data estimated that approximately 11% of ICI-naïve cutaneous melanomas have alterations in IFNγ pathway genes that may confer OV susceptibility. Conclusions: We provide mechanistic support for the use of OVs as a precision medicine strategy for both salvage therapy in ICI-resistant and first-line treatment in melanomas with IFNγ-JAK-STAT pathway mutations. Our study also supports JAK inhibitor–OV combination therapy for treatment-naïve melanomas without IFN signaling defects. See related commentary by Kaufman, p. 3278