Sequencing of first-strand cDNA library reveals full-length transcriptomes

Biomolecules Interacting with Long Noncoding RNAs

This review explores the complex interactions between long noncoding RNAs (lncRNAs) and other biomolecules, highlighting their pivotal roles in gene regulation and cellular function. LncRNAs, defined as RNA transcripts exceeding 200 nucleotides without encoding proteins, are involved in diverse biological processes, from embryogenesis to pathogenesis. They interact with DNA through mechanisms like triplex structure formation, influencing chromatin organization and gene expression. LncRNAs also modulate RNA-mediated processes, including mRNA stability, translational control, and splicing regulation. Their versatility stems from their forming of complex structures that enable interactions with various biomolecules. This review synthesizes current knowledge on lncRNA functions, discusses emerging roles in development and disease, and evaluates potential applications in diagnostics and therapeutics. By examining lncRNA interactions, it provides insights into the intricate regulatory networks governing cellular processes, underscoring the importance of lncRNAs in molecular biology. Unlike the majority of previous reviews that primarily focused on individual aspects of lncRNA biology, this comprehensive review uniquely integrates structural, functional, and mechanistic perspectives on lncRNA interactions across diverse biomolecules. Additionally, this review critically evaluates cutting-edge methodologies for studying lncRNA interactions, bridges fundamental molecular mechanisms with potential clinical applications, and highlights their potential.

Coordinated neuron-specific splicing events restrict nucleosome engagement of the LSD1 histone demethylase complex

Chromatin regulatory proteins are expressed broadly and assumed to exert the same intrinsic function across cell types. Here, we report that 14 chromatin regulators undergo evolutionary-conserved neuron-specific splicing events involving microexons. Among them are two components of a histone demethylase complex: LSD1 H3K4 demethylase and the H3K4me0-reader PHF21A. We found that neuronal LSD1 splicing reduces the enzymes' affinity to the nucleosome. Meanwhile, neuronal PHF21A splicing significantly attenuates histone H3 binding and further ablates the DNA-binding function exerted by an AT-hook motif. Furthermore, in vitro reconstitution of the canonical and neuronal PHF21A-LSD1 complexes, combined with in vivo methylation mapping, identified the neuronal complex as a hypomorphic H3K4 demethylating machinery. The neuronal PHF21A, albeit with its weaker nucleosome binding, is necessary for normal gene expression and the H3K4 landscape in the developing brain. Thus, ubiquitously expressed chromatin regulatory complexes can exert neuron-specific functions via alternative splicing of their subunits.

Elucidating the regulatory role of long non-coding RNAs in drought stress response during seed germination in leaf mustard

Leaf mustard (Brassica juncea L. Czern & Coss), an important vegetable crop, experiences pronounced adversity due to seasonal drought stress, particularly at the seed germination stage. Although there is partial comprehension of drought-responsive genes, the role of long non-coding RNAs (lncRNAs) in adjusting mustard's drought stress response is largely unexplored. In this study, we showed that the drought-tolerant cultivar 'Weiliang' manifested a markedly lower base water potential (-1.073 MPa vs -0.437 MPa) and higher germination percentage (41.2% vs 0%) than the drought-susceptible cultivar 'Shuidong' under drought conditions. High throughput RNA sequencing techniques revealed a significant repertoire of lncRNAs from both cultivars during germination under drought stress, resulting in the identification of 2,087 differentially expressed lncRNAs (DELs) and their correspondingly linked 12,433 target genes. It was noted that 84 genes targeted by DEL exhibited enrichment in the photosynthesis pathway. Gene network construction showed that MSTRG.150397, a regulatory lncRNA, was inferred to potentially modulate key photosynthetic genes (Psb27, PetC, PetH, and PsbW), whilst MSTRG.107159 was indicated as an inhibitory regulator of six drought-responsive PIP genes. Further, weighted gene co-expression network analysis (WGCNA) corroborated the involvement of light intensity and stress response genes targeted by the identified DELs. The precision and regulatory impact of lncRNA were verified through qPCR. This study extends our knowledge of the regulatory mechanisms governing drought stress responses in mustard, which will help strategies to augment drought tolerance in this crop.

Type I/type III IFN and related factors regulate JEV infection and BBB endothelial integrity

BACKGROUND

Japanese encephalitis virus (JEV) remains a predominant cause of Japanese encephalitis (JE) globally. Its infection is usually accompanied by disrupted blood‒brain barrier (BBB) integrity and central nervous system (CNS) inflammation in a poorly understood pathogenesis. Productive JEV infection in brain microvascular endothelial cells (BMECs) is considered the initial event of the virus in penetrating the BBB. Type I/III IFN and related factors have been described as negative regulators in CNS inflammation, whereas their role in JE remains ambiguous.

METHODS

RNA-sequencing profiling (RNA-seq), real-time quantitative PCR, enzyme-linked immunosorbent assay, and Western blotting analysis were performed to analyze the gene and protein expression changes between mock- and JEV-infected hBMECs. Bioinformatic tools were used to cluster altered signaling pathway members during JEV infection. The shRNA-mediated immune factor-knockdown hBMECs and the in vitro transwell BBB model were utilized to explore the interrelation between immune factors, as well as between immune factors and BBB endothelial integrity.

RESULTS

RNA-Seq data of JEV-infected hBMECs identified 417, 1256, and 2748 differentially expressed genes (DEGs) at 12, 36, and 72 h post-infection (hpi), respectively. The altered genes clustered into distinct pathways in gene ontology (GO) terms and KEGG pathway enrichment analysis, including host antiviral immune defense and endothelial cell leakage. Further investigation revealed that pattern-recognition receptors (PRRs, including TLR3, RIG-I, and MDA5) sensed JEV and initiated IRF/IFN signaling. IFNs triggered the expression of interferon-induced proteins with tetratricopeptide repeats (IFITs) via the JAK/STAT pathway. Distinct PRRs exert different functions in barrier homeostasis, while treatment with IFN (IFN-β and IFN-λ1) in hBMECs stabilizes the endothelial barrier by alleviating exogenous destruction. Despite the complex interrelationship, IFITs are considered nonessential in the IFN-mediated maintenance of hBMEC barrier integrity.

CONCLUSIONS

This research provided the first comprehensive description of the molecular mechanisms of host‒pathogen interplay in hBMECs responding to JEV invasion, in which type I/III IFN and related factors strongly correlated with regulating the hBMEC barrier and restricting JEV infection. This might help with developing an attractive therapeutic strategy in JE.

Transcriptome profile analysis of Indian mustard (Brassica juncea L.) during seed germination reveals the drought stress-induced genes associated with energy, hormone, and phenylpropanoid pathways

Indian mustard (Brassica juncea L. Czern and Coss) is an important oil and vegetable crop frequently affected by seasonal drought stress during seed germination, which retards plant growth and causes yield loss considerably. However, the gene networks regulating responses to drought stress in leafy Indian mustard remain elusive. Here, we elucidated the underlying gene networks and pathways of drought response in leafy Indian mustard using next-generation transcriptomic techniques. Phenotypic analysis showed that the drought-tolerant leafy Indian mustard cv. 'WeiLiang' (WL) had a higher germination rate, antioxidant capacity, and better growth performance than the drought-sensitive cv. 'ShuiDong' (SD). Transcriptome analysis identified differentially expressed genes (DEGs) in both cultivars under drought stress during four germination time points (i.e., 0, 12, 24, and 36 h); most of which were classified as drought-responsive, seed germination, and dormancy-related genes. In the Kyoto Encyclopedia of Genes and Genome (KEGG) analyses, three main pathways (i.e., starch and sucrose metabolism, phenylpropanoid biosynthesis, and plant hormone signal transduction) were unveiled involved in response to drought stress during seed germination. Furthermore, Weighted Gene Co-expression Network Analysis (WGCNA) identified several hub genes (novel.12726, novel.1856, BjuB027900, BjuA003402, BjuA021578, BjuA005565, BjuB006596, novel.12977, and BjuA033308) associated with seed germination and drought stress in leafy Indian mustard. Taken together, these findings deepen our understanding of the gene networks for drought responses during seed germination in leafy Indian mustard and provide potential target genes for the genetic improvement of drought tolerance in this crop.

Tn5 tagments and transposes oligos to single-stranded DNA for strand-specific RNA sequencing

Tn5 transposon tagments double-stranded DNA and RNA/DNA hybrids to generate nucleic acids that are ready to be amplified for high-throughput sequencing. The nucleic acid substrates for the Tn5 transposon must be explored to increase the applications of Tn5. Here, we found that the Tn5 transposon can transpose oligos into the 5' end of single-stranded DNA longer than 140 nucleotides. Based on this property of Tn5, we developed a tagmentation-based and ligation-enabled single-stranded DNA sequencing method called TABLE-seq. Through a series of reaction temperature, time, and enzyme concentration tests, we applied TABLE-seq to strand-specific RNA sequencing, starting with as little as 30 pg of total RNA. Moreover, compared with traditional dUTP-based strand-specific RNA sequencing, this method detects more genes, has a higher strand specificity, and shows more evenly distributed reads across genes. Together, our results provide insights into the properties of Tn5 transposons and expand the applications of Tn5 in cutting-edge sequencing techniques.

Alcohol-specific transcriptional dynamics of memory reconsolidation and relapse

Relapse, a critical issue in alcohol addiction, can be attenuated by disruption of alcohol-associated memories. Memories are thought to temporarily destabilize upon retrieval during the reconsolidation process. Here, we provide evidence for unique transcriptional dynamics underpinning alcohol memory reconsolidation. Using a mouse place-conditioning procedure, we show that alcohol-memory retrieval increases the mRNA expression of immediate-early genes in the dorsal hippocampus and medial prefrontal cortex, and that alcohol seeking is abolished by post-retrieval non-specific inhibition of gene transcription, or by downregulating ARC expression using antisense-oligodeoxynucleotides. However, since retrieval of memories for a natural reward (sucrose) also increased the same immediate-early gene expression, we explored for alcohol-specific transcriptional changes using RNA-sequencing. We revealed a unique transcriptional fingerprint activated by alcohol memories, as the expression of this set of plasticity-related genes was not altered by sucrose-memory retrieval. Our results suggest that alcohol memories may activate two parallel transcription programs: one is involved in memory reconsolidation in general, and another is specifically activated during alcohol-memory processing.

Epitranscriptomic mapping of RNA modifications at single-nucleotide resolution using rhodamine sequencing (Rho-seq)

The recent development of epitranscriptomics revealed a new fundamental layer of gene expression, but the mapping of most RNA modifications remains technically challenging. Here, we describe our protocol for Rho-Seq, which enables the mapping of dihydrouridine RNA modification at single-nucleotide resolution. Rho-Seq relies on specific rhodamine-labeling of a subset of modified nucleotides that hinders reverse transcription. Although Rho-Seq was initially applied to the detection of dihydrouridine, we show here that it is applicable to other modifications including 7-methylguanosine or 4-thiouridine. For complete details on the use and execution of this protocol, please refer to Finet et al. (2022).

PTT-quant: a new method for direct identification and absolute quantification of premature transcription termination events, following the example of bacterial riboswitches

Regulation of gene expression by premature termination of transcription has been well described in all domains of life, including metazoans, yeast, plants, and bacteria. Although methods for identification of such regulatory events by sequencing are available, the focused biochemical studies of the mechanism are hampered by lack of highly sensitive and accurate experimental methods. Here, we propose a new method for absolute quantification of premature transcription termination events, PTT-quant. It is based on highly sensitive two-step digital droplet PCR protocol, coupled with normalized cDNA synthesis attained by site-specific pre-cleavage of investigated transcripts with RNase H. As a consequence, our method enables the reliable and sensitive quantification of both, prematurely terminated and full-length transcripts. By application of our method to investigation of transcriptional riboswitches in Bacillus subtilis, we were able to precisely measure the dynamics of S-adenosylmethionine (SAM) riboswitch induction, which turned to be ~ 23% higher in comparison the results obtained without cDNA synthesis normalization.

Key points

• A novel method for quantification of premature transcription termination events was established.

• PTT-quant measures absolute concentration of full-length and terminated transcripts.

• RNase H and the digital droplet PCR technique is used in PTT-quant.

Sequences to Differences in Gene Expression: Analysis of RNA-Seq Data

RNA-Seq is now a routinely employed assay to measure gene expression. As the technique matured over the last decade, so have dedicated analytic tools. In this chapter, we first describe the mainstream as well as the most up-to-date protocols and their implications on downstream analysis. We then detail the steps entailing RNA-Seq analysis in three main stages: (i) preprocessing and data preparation, (ii) upstream processing, and (iii) high-level analyses. We review the most recent and relevant tools as one workflow following a stepwise order. The chapter further encompasses in-depth features of these tools. Details of the required code are made available throughout the chapter, as well as of the underlying statistics. We illustrate these steps with analysis of publicly available RNA-Seq data.

An ultra low-input method for global RNA structure probing uncovers Regnase-1-mediated regulation in macrophages

To enable diverse functions and precise regulation, an RNA sequence often folds into complex yet distinct structures in different cellular states. Probing RNA in its native environment is essential to uncovering RNA structures of biological contexts. However, current methods generally require large amounts of input RNA and are challenging for physiologically relevant use. Here, we report smartSHAPE, a new RNA structure probing method that requires very low amounts of RNA input due to the largely reduced artefact of probing signals and increased efficiency of library construction. Using smartSHAPE, we showcased the profiling of the RNA structure landscape of mouse intestinal macrophages upon inflammation, and provided evidence that RNA conformational changes regulate immune responses. These results demonstrate that smartSHAPE can greatly expand the scope of RNA structure-based investigations in practical biological systems, and also provide a research paradigm for the study of post-transcriptional regulation.

Screening a Spliced Leader-Based Symbiodinium microadriaticum cDNA Library Using the Yeast-Two Hybrid System Reveals a Hemerythrin-Like Protein as a Putative SmicRACK1 Ligand

The dinoflagellate Symbiodiniaceae family plays a central role in the health of the coral reef ecosystem via the symbiosis that establishes with its inhabiting cnidarians and supports the host metabolism. In the last few decades, coral reefs have been threatened by pollution and rising temperatures which have led to coral loss. These events have raised interest in studying Symbiodiniaceae and their hosts; however, progress in understanding their metabolism, signal transduction pathways, and physiology in general, has been slow because dinoflagellates present peculiar characteristics. We took advantage of one of these peculiarities; namely, the post-transcriptional addition of a Dino Spliced Leader (Dino-SL) to the 5' end of the nuclear mRNAs, and used it to generate cDNA libraries from Symbiodinium microadriaticum. We compared sequences from two Yeast-Two Hybrid System cDNA Libraries, one based on the Dino-SL sequence, and the other based on the SMART technology (Switching Mechanism at 5' end of RNA Transcript) which exploits the template switching function of the reverse transcriptase. Upon comparison of the performance of both libraries, we obtained a significantly higher yield, number and length of sequences, number of transcripts, and better 5' representation from the Dino-SL based library than from the SMART library. In addition, we confirmed that the cDNAs from the Dino-SL library were adequately expressed in the yeast cells used for the Yeast-Two Hybrid System which resulted in successful screening for putative SmicRACK1 ligands, which yielded a putative hemerythrin-like protein.

KDM1A maintains genome-wide homeostasis of transcriptional enhancers

Transcriptional enhancers enable exquisite spatiotemporal control of gene expression in metazoans. Enrichment of monomethylation of histone H3 lysine 4 (H3K4me1) is a major chromatin signature of transcriptional enhancers. Lysine (K)-specific demethylase 1A (KDM1A, also known as LSD1), an H3K4me2/me1 demethylase, inactivates stem-cell enhancers during the differentiation of mouse embryonic stem cells (mESCs). However, its role in undifferentiated mESCs remains obscure. Here, we show that KDM1A actively maintains the optimal enhancer status in both undifferentiated and lineage-committed cells. KDM1A occupies a majority of enhancers in undifferentiated mESCs. KDM1A levels at enhancers exhibit clear positive correlations with its substrate H3K4me2, H3K27ac, and transcription at enhancers. In Kdm1a-deficient mESCs, a large fraction of these enhancers gains additional H3K4 methylation, which is accompanied by increases in H3K27 acetylation and increased expression of both enhancer RNAs (eRNAs) and target genes. In postmitotic neurons, loss of KDM1A leads to premature activation of neuronal activity-dependent enhancers and genes. Taken together, these results suggest that KDM1A is a versatile regulator of enhancers and acts as a rheostat to maintain optimal enhancer activity by counterbalancing H3K4 methylation at enhancers.

SL-quant: a fast and flexible pipeline to quantify spliced leader trans-splicing events from RNA-seq data

Background

The spliceosomal transfer of a short spliced leader (SL) RNA to an independent pre-mRNA molecule is called SL trans-splicing and is widespread in the nematode Caenorhabditis elegans. While RNA-sequencing (RNA-seq) data contain information on such events, properly documented methods to extract them are lacking.

Findings

To address this, we developed SL-quant, a fast and flexible pipeline that adapts to paired-end and single-end RNA-seq data and accurately quantifies SL trans-splicing events. It is designed to work downstream of read mapping and uses the reads left unmapped as primary input. Briefly, the SL sequences are identified with high specificity and are trimmed from the input reads, which are then remapped on the reference genome and quantified at the nucleotide position level (SL trans-splice sites) or at the gene level.

Conclusions

SL-quant completes within 10 minutes on a basic desktop computer for typical C. elegans RNA-seq datasets and can be applied to other species as well. Validating the method, the SL trans-splice sites identified display the expected consensus sequence, and the results of the gene-level quantification are predictive of the gene position within operons. We also compared SL-quant to a recently published SL-containing read identification strategy that was found to be more sensitive but less specific than SL-quant. Both methods are implemented as a bash script available under the MIT license [1]. Full instructions for its installation, usage, and adaptation to other organisms are provided.

Capturing the 'ome': the expanding molecular toolbox for RNA and DNA library construction

All sequencing experiments and most functional genomics screens rely on the generation of libraries to comprehensively capture pools of targeted sequences. In the past decade especially, driven by the progress in the field of massively parallel sequencing, numerous studies have comprehensively assessed the impact of particular manipulations on library complexity and quality, and characterized the activities and specificities of several key enzymes used in library construction. Fortunately, careful protocol design and reagent choice can substantially mitigate many of these biases, and enable reliable representation of sequences in libraries. This review aims to guide the reader through the vast expanse of literature on the subject to promote informed library generation, independent of the application.

Altered Gene-Regulatory Function of KDM5C by a Novel Mutation Associated With Autism and Intellectual Disability

Intellectual disability (ID) affects up to 2% of the population world-wide and often coincides with other neurological conditions such as autism spectrum disorders. Mutations in KDM5C cause Mental Retardation, X-linked, Syndromic, Claes-Jensen type (MRXSCJ, OMIM #300534) and are one of the most common causes of X-linked ID. KDM5C encodes a histone demethylase for di- and tri-methylated histone H3 lysine 4 (H3K4me2/3), which are enriched in transcriptionally engaged promoter regions. KDM5C regulates gene transcription; however, it remains unknown whether removal of H3K4me is fully responsible for KDM5C-mediated gene regulation. Most mutations functionally tested to date result in reduced enzymatic activity of KDM5C, indicating loss of demethylase function as the primary mechanism underlying MRXSCJ. Here, we report a novel KDM5C mutation, R1115H, identified in an individual displaying MRXSCJ-like symptoms. The carrier mother's cells exhibited a highly skewed X-inactivation pattern. The KDM5C-R1115H substitution does not have an impact on enzymatic activity nor protein stability. However, when overexpressed in post-mitotic neurons, KDM5C-R1115H failed to fully suppress expression of target genes, while the mutant also affected expression of a distinct set of genes compared to KDM5C-wildtype. These results suggest that KDM5C may have non-enzymatic roles in gene regulation, and alteration of these roles contributes to MRXSCJ in this patient.

Transcriptome Analysis Revealed Impaired cAMP Responsiveness in PHF21A-Deficient Human Cells

Potocki-Shaffer Syndrome is a rare neurodevelopmental syndrome associated with microdeletion of a region of Chromosome 11p11.2. Genetic evidence has implicated haploinsufficiency of PHF21A, a gene that encodes a histone-binding protein, as the likely cause of intellectual disability and craniofacial abnormalities in Potocki-Shaffer Syndrome. However, the molecular consequences of reduced PHF21A expression remain elusive. In this study, we analyzed by RNA-Sequencing (RNA-Seq) two patient-derived cell lines with heterozygous loss of PHF21A compared to unaffected individuals and identified 1,885 genes that were commonly misregulated. The patient cells displayed down-regulation of key pathways relevant to learning and memory, including Cyclic Adenosine Monophosphate (cAMP)-signaling pathway genes. We found that PHF21A is required for full induction of a luciferase reporter carrying cAMP-responsive elements (CRE) following stimulation by the cAMP analog, forskolin. Finally, PHF21A-deficient patient-derived cells exhibited a delayed induction of immediate early genes following forskolin stimulation. These results suggest that an impaired response to cAMP signaling might be involved in the pathology of PHF21A deficiency. This article is part of a Special Issue entitled: [SI: Molecules & Cognition].

A Mouse Model of X-linked Intellectual Disability Associated with Impaired Removal of Histone Methylation

Mutations in a number of chromatin modifiers are associated with human neurological disorders. KDM5C, a histone H3 lysine 4 di- and tri-methyl (H3K4me2/3)-specific demethylase, is frequently mutated in X-linked intellectual disability (XLID) patients. Here, we report that disruption of the mouse Kdm5c gene recapitulates adaptive and cognitive abnormalities observed in XLID, including impaired social behavior, memory deficits, and aggression. Kdm5c-knockout brains exhibit abnormal dendritic arborization, spine anomalies, and altered transcriptomes. In neurons, Kdm5c is recruited to promoters that harbor CpG islands decorated with high levels of H3K4me3, where it fine-tunes H3K4me3 levels. Kdm5c predominantly represses these genes, which include members of key pathways that regulate the development and function of neuronal circuitries. In summary, our mouse behavioral data strongly suggest that KDM5C mutations are causal to XLID. Furthermore, our findings suggest that loss of KDM5C function may impact gene expression in multiple regulatory pathways relevant to the clinical phenotypes.

Systematic identification of factors for provirus silencing in embryonic stem cells

Embryonic stem cells (ESCs) repress the expression of exogenous proviruses and endogenous retroviruses (ERVs). Here, we systematically dissected the cellular factors involved in provirus repression in embryonic carcinomas (ECs) and ESCs by a genome-wide siRNA screen. Histone chaperones (Chaf1a/b), sumoylation factors (Sumo2/Ube2i/Sae1/Uba2/Senp6), and chromatin modifiers (Trim28/Eset/Atf7ip) are key determinants that establish provirus silencing. RNA-seq analysis uncovered the roles of Chaf1a/b and sumoylation modifiers in the repression of ERVs. ChIP-seq analysis demonstrates direct recruitment of Chaf1a and Sumo2 to ERVs. Chaf1a reinforces transcriptional repression via its interaction with members of the NuRD complex (Kdm1a, Hdac1/2) and Eset, while Sumo2 orchestrates the provirus repressive function of the canonical Zfp809/Trim28/Eset machinery by sumoylation of Trim28. Our study reports a genome-wide atlas of functional nodes that mediate proviral silencing in ESCs and illuminates the comprehensive, interconnected, and multi-layered genetic and epigenetic mechanisms by which ESCs repress retroviruses within the genome.