LEDGF/p75 interacts with mRNA splicing factors and targets HIV-1 integration to highly spliced genes.

The host chromatin-binding factor LEDGF/p75 interacts with HIV-1 integrase and directs integration to active transcription units. To understand how LEDGF/p75 recognizes transcription units, we sequenced 1 million HIV-1 integration sites isolated from cultured HEK293T cells. Analysis of integration sites showed that cancer genes were preferentially targeted, raising concerns about using lentivirus vectors for gene therapy. Additional analysis led to the discovery that introns and alternative splicing contributed significantly to integration site selection. These correlations were independent of transcription levels, size of transcription units, and length of the introns. Multivariate analysis with five parameters previously found to predict integration sites showed that intron density is the strongest predictor of integration density in transcription units. Analysis of previously published HIV-1 integration site data showed that integration density in transcription units in mouse embryonic fibroblasts also correlated strongly with intron number, and this correlation was absent in cells lacking LEDGF. Affinity purification showed that LEDGF/p75 is associated with a number of splicing factors, and RNA sequencing (RNA-seq) analysis of HEK293T cells lacking LEDGF/p75 or the LEDGF/p75 integrase-binding domain (IBD) showed that LEDGF/p75 contributes to splicing patterns in half of the transcription units that have alternative isoforms. Thus, LEDGF/p75 interacts with splicing factors, contributes to exon choice, and directs HIV-1 integration to transcription units that are highly spliced.

Transposable element insertions in fission yeast drive adaptation to environmental stress.

Cells are regularly exposed to a range of naturally occurring stress that can restrict growth or cause lethality. In response, cells activate expression networks with hundreds of genes that together increase resistance to common environmental insults. However, stress response networks can be insufficient to ensure survival, which raises the question of whether cells possess genetic programs that can promote adaptation to novel forms of stress. We found transposable element (TE) mobility in Schizosaccharomyces pombe was greatly increased when cells were exposed to unusual forms of stress such as heavy metals, caffeine, and the plasticizer phthalate. By subjecting TE-tagged cells to CoCl2, we found the TE integration provided the major path to resistance. Groups of insertions that provided resistance were linked to TOR regulation and metal response genes. We extended our study of adaptation by analyzing TE positions in 57 genetically distinct wild strains. The genomic positions of 1048 polymorphic LTRs were strongly associated with a range of stress response genes, indicating TE integration promotes adaptation in natural conditions. These data provide strong support for the idea, first proposed by Barbara McClintock, that TEs provide a system to modify the genome in response to stress.

Resource: A multi-species multi-timepoint transcriptome database and webpage for the pineal gland and retina.

The website and database https://snengs.nichd.nih.gov provides RNA sequencing data from multi-species analysis of the pineal glands from zebrafish (Danio rerio), chicken (White Leghorn), rat (Rattus novegicus), mouse (Mus musculus), rhesus macaque (Macaca mulatta), and human (Homo sapiens); in most cases, retinal data are also included along with results of the analysis of a mixture of RNA from tissues. Studies cover day and night conditions; in addition, a time series over multiple hours, a developmental time series and pharmacological experiments on rats are included. The data have been uniformly re-processed using the latest methods and assemblies to allow for comparisons between experiments and to reduce processing differences. The website presents search functionality, graphical representations, Excel tables, and track hubs of all data for detailed visualization in the UCSC Genome Browser. As more data are collected from investigators and improved genomes become available in the future, the website will be updated. This database is in the public domain and elements can be reproduced by citing the URL and this report. This effort makes the results of 21st century transcriptome profiling widely available in a user-friendly format that is expected to broadly influence pineal research.

Host factors that promote retrotransposon integration are similar in distantly related eukaryotes.

Retroviruses and Long Terminal Repeat (LTR)-retrotransposons have distinct patterns of integration sites. The oncogenic potential of retrovirus-based vectors used in gene therapy is dependent on the selection of integration sites associated with promoters. The LTR-retrotransposon Tf1 of Schizosaccharomyces pombe is studied as a model for oncogenic retroviruses because it integrates into the promoters of stress response genes. Although integrases (INs) encoded by retroviruses and LTR-retrotransposons are responsible for catalyzing the insertion of cDNA into the host genome, it is thought that distinct host factors are required for the efficiency and specificity of integration. We tested this hypothesis with a genome-wide screen of host factors that promote Tf1 integration. By combining an assay for transposition with a genetic assay that measures cDNA recombination we could identify factors that contribute differentially to integration. We utilized this assay to test a collection of 3,004 S. pombe strains with single gene deletions. Using these screens and immunoblot measures of Tf1 proteins, we identified a total of 61 genes that promote integration. The candidate integration factors participate in a range of processes including nuclear transport, transcription, mRNA processing, vesicle transport, chromatin structure and DNA repair. Two candidates, Rhp18 and the NineTeen complex were tested in two-hybrid assays and were found to interact with Tf1 IN. Surprisingly, a number of pathways we identified were found previously to promote integration of the LTR-retrotransposons Ty1 and Ty3 in Saccharomyces cerevisiae, indicating the contribution of host factors to integration are common in distantly related organisms. The DNA repair factors are of particular interest because they may identify the pathways that repair the single stranded gaps flanking the sites of strand transfer following integration of LTR retroelements.

Serial number tagging reveals a prominent sequence preference of retrotransposon integration.

Transposable elements (TE) have both negative and positive impact on the biology of their host. As a result, a balance is struck between the host and the TE that relies on directing integration to specific genome territories. The extraordinary capacity of DNA sequencing can create ultra dense maps of integration that are being used to study the mechanisms that position integration. Unfortunately, the great increase in the numbers of insertion sites detected comes with the cost of not knowing which positions are rare targets and which sustain high numbers of insertions. To address this problem we developed the serial number system, a TE tagging method that measures the frequency of integration at single nucleotide positions. We sequenced 1 million insertions of retrotransposon Tf1 in the genome of Schizosaccharomyces pombe and obtained the first profile of integration with frequencies for each individual position. Integration levels at individual nucleotides varied over two orders of magnitude and revealed that sequence recognition plays a key role in positioning integration. The serial number system is a general method that can be applied to determine precise integration maps for retroviruses and gene therapy vectors.

The role of LEDGF in transcription is exploited by HIV-1 to position integration.

HIV-1 integration occurs across actively transcribed genes due to the interaction of integrase with host chromatin factor LEDGF. Although LEDGF was originally isolated as a co-activator that stimulates promoter activity in purified systems, this role is inconsistent with LEDGF-mediated integration across gene bodies and with data indicating LEDGF is a histone chaperone that promotes transcriptional elongation. We found LEDGF is enriched in pronounced peaks that match the enrichments of H3K4me3 and RNA Pol II at transcription start sites (TSSs) of active promoters. Our genome-wide chromatin mapping revealed that MLL1 had a dominant role in recruiting LEDGF to promoters and the presence of LEDGF recruits RNA Pol II. Enrichment of LEDGF at TSSs correlates strongly with levels of integration across the transcribed sequences, indicating that LEDGF at TSSs contributed to integration across gene bodies. Although the N-terminal Pro-Trp-Trp-Pro (PWWP) domain of LEDGF interacts with nucleosomes containing H3K36me3, a modification thought to recruit LEDGF to chromatin, we found H3K36me3 does not contribute to gene specificity of integration. These data support a dual role model of LEDGF where it is tethered to promoters by MLL1 and recruits RNA Pol II. Subsequently, LEDGF travels across genes to effect HIV-1 integration. Our data also provides a mechanistic context for the contribution made by LEDGF to MLL1-based infant acute leukemia and acute myeloid leukemia in adults.

tirant, a newly discovered active endogenous retrovirus in Drosophila simulans.

Endogenous retroviruses have the ability to become permanently integrated into the genomes of their host, and they are generally transmitted vertically from parent to progeny. With the exception of gypsy, few endogenous retroviruses have been identified in insects. In this study, we describe the tirant endogenous retrovirus in a subset of Drosophila simulans natural populations. By focusing on the envelope gene, we show that the entire retroviral cycle (transcription, translation, and retrotransposition) can be completed for tirant within one population of this species.

Extensive diversity in RNA termination and regulation revealed by transcriptome mapping for the Lyme pathogen B. burgdorferi.

Transcription termination is an essential and dynamic process that can tune gene expression in response to diverse molecular signals. Yet, the genomic positions, molecular mechanisms, and regulatory consequences of termination have only been studied thoroughly in model bacteria. We employed complementary RNA-seq approaches to map RNA ends for the transcriptome of the spirochete Borrelia burgdorferi - the etiological agent of Lyme disease. By systematically mapping B. burgdorferi RNA ends at single nucleotide resolution, we delineated complex gene arrangements and operons and mapped untranslated regions (UTRs) and small RNAs (sRNAs). We experimentally tested modes of B. burgdorferi transcription termination and compared our findings to observations in E. coli , P. aeruginosa , and B. subtilis . We discovered 63% of B. burgdorferi RNA 3' ends map upstream or internal to open reading frames (ORFs), suggesting novel mechanisms of regulation. Northern analysis confirmed the presence of stable 5' derived RNAs from mRNAs encoding gene products involved in the unique infectious cycle of B. burgdorferi . We suggest these RNAs resulted from premature termination and regulatory events, including forms of cis- acting regulation. For example, we documented that the polyamine spermidine globally influences the generation of truncated mRNAs. In one case, we showed that high spermidine concentrations increased levels of RNA fragments derived from an mRNA encoding a spermidine import system, with a concomitant decrease in levels of the full- length mRNA. Collectively, our findings revealed new insight into transcription termination and uncovered an abundance of potential RNA regulators.

Extensive diversity in RNA termination and regulation revealed by transcriptome mapping for the Lyme pathogen Borrelia burgdorferi.

Transcription termination is an essential and dynamic process that can tune gene expression in response to diverse molecular signals. Yet, the genomic positions, molecular mechanisms, and regulatory consequences of termination have only been studied thoroughly in model bacteria. Here, we use several RNA-seq approaches to map RNA ends for the transcriptome of the spirochete Borrelia burgdorferi - the etiological agent of Lyme disease. We identify complex gene arrangements and operons, untranslated regions and small RNAs. We predict intrinsic terminators and experimentally test examples of Rho-dependent transcription termination. Remarkably, 63% of RNA 3' ends map upstream of or internal to open reading frames (ORFs), including genes involved in the unique infectious cycle of B. burgdorferi. We suggest these RNAs result from premature termination, processing and regulatory events such as cis-acting regulation. Furthermore, the polyamine spermidine globally influences the generation of truncated mRNAs. Collectively, our findings provide insights into transcription termination and uncover an abundance of potential RNA regulators in B. burgdorferi.

Integrated spatiotemporal transcriptomic resolution of embryonic palate osteogenesis.

The differentiation of osteoblasts and the subsequent formation of bone marks an important terminal phase in palate formation that leads to the separation of the oral and nasal cavities. While the developmental events that precede palatal osteogenesis are well explored, major gaps remain in our understanding of the molecular mechanisms that lead to the bony union of fusing palatal shelves. Herein, the timeline of osteogenic transcriptional programming is unveiled in the embryonic palate by way of integrated bulk, single-cell, and spatially resolved RNA-seq analyses. We define spatially restricted expression patterns of key marker genes, both regulatory and structural, that are differentially expressed during palatal fusion, including the identification of several novel genes ( Deup1, Dynlrb2, Lrrc23 ) spatially restricted in expression to the palate, providing a relevant framework for future studies that identify new candidate genes for cleft palate anomalies in humans as well as the timing of mammalian embryonic palatal osteogenesis.

Multimodal spatiotemporal transcriptomic resolution of embryonic palate osteogenesis.

The terminal differentiation of osteoblasts and subsequent formation of bone marks an important phase in palate development that leads to the separation of the oral and nasal cavities. While the morphogenetic events preceding palatal osteogenesis are well explored, major gaps remain in our understanding of the molecular mechanisms driving the formation of this bony union of the fusing palate. Through bulk, single-nucleus, and spatially resolved RNA-sequencing analyses of the developing secondary palate, we identify a shift in transcriptional programming between embryonic days 14.5 and 15.5 pinpointing the onset of osteogenesis. We define spatially restricted expression patterns of key osteogenic marker genes that are differentially expressed between these developmental timepoints. Finally, we identify genes in the palate highly expressed by palate nasal epithelial cells, also enriched within palatal osteogenic mesenchymal cells. This investigation provides a relevant framework to advance palate-specific diagnostic and therapeutic biomarker discovery.

Identification of an integrase-independent pathway of retrotransposition.

Retroviruses and long terminal repeat retrotransposons rely on integrase (IN) to insert their complementary DNA (cDNA) into the genome of host cells. Nevertheless, in the absence of IN, retroelements can retain notable levels of insertion activity. We have characterized the IN-independent pathway of Tf1 and found that insertion sites had homology to the primers of reverse transcription, which form single-stranded DNAs at the termini of the cDNA. In the absence of IN activity, a similar bias was observed with HIV-1. Our studies showed that the Tf1 insertions result from single-strand annealing, a noncanonical form of homologous recombination mediated by Rad52. By expanding our analysis of insertions to include repeat sequences, we found most formed tandem elements by inserting at preexisting copies of a related transposable element. Unexpectedly, we found that wild-type Tf1 uses the IN-independent pathway as an alternative mode of insertion.

A protocol for transposon insertion sequencing in Schizosaccharomyces pombe to identify factors that maintain heterochromatin.

Transposon insertion sequencing (TIS) is a highly effective method used with bacteria to identify genes important for growth in any condition of interest. Previously, we adapted this method to identify essential genes of the yeast Schizosaccharomyces pombe. Here, we describe modifications used to identify genes necessary for the formation of centromeric heterochromatin. For complete details on the use and execution of this protocol, please refer to Lee et al. (2020).

Regulatory roles of Escherichia coli 5' UTR and ORF-internal RNAs detected by 3' end mapping.

Many bacterial genes are regulated by RNA elements in their 5´ untranslated regions (UTRs). However, the full complement of these elements is not known even in the model bacterium Escherichia coli. Using complementary RNA-sequencing approaches, we detected large numbers of 3´ ends in 5´ UTRs and open reading frames (ORFs), suggesting extensive regulation by premature transcription termination. We documented regulation for multiple transcripts, including spermidine induction involving Rho and translation of an upstream ORF for an mRNA encoding a spermidine efflux pump. In addition to discovering novel sites of regulation, we detected short, stable RNA fragments derived from 5´ UTRs and sequences internal to ORFs. Characterization of three of these transcripts, including an RNA internal to an essential cell division gene, revealed that they have independent functions as sRNA sponges. Thus, these data uncover an abundance of cis- and trans-acting RNA regulators in bacterial 5´ UTRs and internal to ORFs.

In most organisms, specific segments of a cell's genetic information are copied to form single-stranded molecules of various sizes and purposes. Each of these RNA molecules, as they are known, is constructed as a chain that starts at the 5´ end and terminates at the 3´ end. Certain RNAs carry the information present in a gene, which provides the instructions that a cell needs to build proteins. Some, however, are 'non-coding' and instead act to fine-tune the activity of other RNAs. These regulatory RNAs can be separate from the RNAs they control, or they can be embedded in the very sequences they regulate; new evidence also shows that certain regulatory RNAs can act in both ways. Many regulatory RNAs are yet to be catalogued, even in simple, well-studied species such as the bacterium Escherichia coli. Here, Adams et al. aimed to better characterize the regulatory RNAs present in E. coli by mapping out the 3´ ends of every RNA molecule in the bacterium. This revealed many new regulatory RNAs and offered insights into where these sequences are located. For instance, the results show that several of these RNAs were embedded within RNA produced from larger genes. Some were nested in coding RNAs, and were parts of a longer RNA sequence that is adjacent to the protein coding segment. Others, however, were present within the instructions that code for a protein. The work by Adams et al. reveals that regulatory RNAs can be located in unexpected places, and provides a method for identifying them. This can be applied to other types of bacteria, in particular in species with few known RNA regulators.

NMDARs Drive the Expression of Neuropsychiatric Disorder Risk Genes Within GABAergic Interneuron Subtypes in the Juvenile Brain.

Medial ganglionic eminence (MGE)-derived parvalbumin (PV)+, somatostatin (SST)+and Neurogliaform (NGFC)-type cortical and hippocampal interneurons, have distinct molecular, anatomical, and physiological properties. However, the molecular mechanisms regulating their maturation remain poorly understood. Here, via single-cell transcriptomics, we show that the obligate NMDA-type glutamate receptor (NMDAR) subunit gene Grin1 mediates transcriptional regulation of gene expression in specific subtypes of MGE-derived interneurons, leading to altered subtype abundances. Notably, MGE-specific early developmental Grin1 loss results in a broad downregulation of diverse transcriptional, synaptogenic and membrane excitability regulatory programs in the juvenile brain. These widespread gene expression abnormalities mirror aberrations that are typically associated with neurodevelopmental disorders. Our study hence provides a road map for the systematic examination of NMDAR signaling in interneuron subtypes, revealing potential MGE-specific genetic targets that could instruct future therapies of psychiatric disorders.

Dense Transposon Integration Reveals Essential Cleavage and Polyadenylation Factors Promote Heterochromatin Formation.

Heterochromatin functions as a scaffold for factors responsible for gene silencing and chromosome segregation. Heterochromatin can be assembled by multiple pathways, including RNAi and RNA surveillance. We identified factors that form heterochromatin using dense profiles of transposable element integration in Schizosaccharomyces pombe. The candidates include a large number of essential proteins such as four canonical mRNA cleavage and polyadenylation factors. We find that Iss1, a subunit of the poly(A) polymerase module, plays a role in forming heterochromatin in centromere repeats that is independent of RNAi. Genome-wide maps reveal that Iss1 accumulates at genes regulated by RNA surveillance. Iss1 interacts with RNA surveillance factors Mmi1 and Rrp6, and importantly, Iss1 contributes to RNA elimination that forms heterochromatin at meiosis genes. Our profile of transposable element integration supports the model that a network of mRNA cleavage and polyadenylation factors coordinates RNA surveillance, including the mechanism that forms heterochromatin at meiotic genes.

Nucleoporin TPR is an integral component of the TREX-2 mRNA export pathway.

Nuclear pore complexes (NPCs) are important for cellular functions beyond nucleocytoplasmic trafficking, including genome organization and gene expression. This multi-faceted nature and the slow turnover of NPC components complicates investigations of how individual nucleoporins act in these diverse processes. To address this question, we apply an Auxin-Induced Degron (AID) system to distinguish roles of basket nucleoporins NUP153, NUP50 and TPR. Acute depletion of TPR causes rapid and pronounced changes in transcriptomic profiles. These changes are dissimilar to shifts observed after loss of NUP153 or NUP50, but closely related to changes caused by depletion of mRNA export receptor NXF1 or the GANP subunit of the TRanscription-EXport-2 (TREX-2) mRNA export complex. Moreover, TPR depletion disrupts association of TREX-2 subunits (GANP, PCID2, ENY2) to NPCs and results in abnormal RNA transcription and export. Our findings demonstrate a unique and pivotal role of TPR in gene expression through TREX-2- and/or NXF1-dependent mRNA turnover.

The Long Terminal Repeat Retrotransposons Tf1 and Tf2 of Schizosaccharomyces pombe.

The long terminal repeat (LTR) retrotransposons Tf1 and Tf2 of Schizosaccharomyces pombe are active mobile elements of the Ty3/gypsy family. The mobilization of these retrotransposons depends on particle formation, reverse transcription and integration, processes typical of other LTR retrotransposons. However, Tf1 and Tf2 are distinct from other LTR elements in that they assemble virus-like particles from a single primary translation product, initiate reverse transcription with an unusual self-priming mechanism, and, in the case of Tf1, integrate with a pattern that favors specific promoters of RNA pol II-transcribed genes. To avoid the chromosome instability and genome damage that results from increased copy number, S. pombe applies a variety of defense mechanisms that restrict Tf1 and Tf2 activity. The mRNA of the Tf elements is eliminated by an exosome-based pathway when cells are in favorable conditions whereas nutrient deprivation triggers an RNA interference-dependent pathway that results in the heterochromatization of the elements. Interestingly, Tf1 integrates into the promoters of stress-induced genes and these insertions are capable of increasing the expression of adjacent genes. These properties of Tf1 transposition raise the possibility that Tf1 benefits cells with specific insertions by providing resistance to environmental stress.

Single-Nucleotide-Specific Targeting of the Tf1 Retrotransposon Promoted by the DNA-Binding Protein Sap1 of Schizosaccharomyces pombe.

Transposable elements (TEs) constitute a substantial fraction of the eukaryotic genome and, as a result, have a complex relationship with their host that is both adversarial and dependent. To minimize damage to cellular genes, TEs possess mechanisms that target integration to sequences of low importance. However, the retrotransposon Tf1 of Schizosaccharomyces pombe integrates with a surprising bias for promoter sequences of stress-response genes. The clustering of integration in specific promoters suggests that Tf1 possesses a targeting mechanism that is important for evolutionary adaptation to changes in environment. We report here that Sap1, an essential DNA-binding protein, plays an important role in Tf1 integration. A mutation in Sap1 resulted in a 10-fold drop in Tf1 transposition, and measures of transposon intermediates support the argument that the defect occurred in the process of integration. Published ChIP-Seq data on Sap1 binding combined with high-density maps of Tf1 integration that measure independent insertions at single-nucleotide positions show that 73.4% of all integration occurs at genomic sequences bound by Sap1. This represents high selectivity because Sap1 binds just 6.8% of the genome. A genome-wide analysis of promoter sequences revealed that Sap1 binding and amounts of integration correlate strongly. More important, an alignment of the DNA-binding motif of Sap1 revealed integration clustered on both sides of the motif and showed high levels specifically at positions +19 and -9. These data indicate that Sap1 contributes to the efficiency and position of Tf1 integration.