DNA G-Quadruplex Is a Transcriptional Control Device That Regulates Memory.

The conformational state of DNA fine-tunes the transcriptional rate and abundance of RNA. Here, we report that G-quadruplex DNA (G4-DNA) accumulates in neurons, in an experience-dependent manner, and that this is required for the transient silencing and activation of genes that are critically involved in learning and memory in male C57/BL6 mice. In addition, site-specific resolution of G4-DNA by dCas9-mediated deposition of the helicase DHX36 impairs fear extinction memory. Dynamic DNA structure states therefore represent a key molecular mechanism underlying memory consolidation.One-Sentence Summary: G4-DNA is a molecular switch that enables the temporal regulation of the gene expression underlying the formation of fear extinction memory.

Divergent evolution of developmental timing in the neocortex revealed by marsupial and eutherian transcriptomes.

Only mammals evolved a neocortex, which integrates sensory-motor and cognitive functions. Significant diversifications in the cellular composition and connectivity of the neocortex occurred between the two main therian groups: marsupials and eutherians. However, the developmental mechanisms underlying these diversifications are largely unknown. Here, we compared the neocortical transcriptomes of Sminthopsis crassicaudata, a mouse-sized marsupial, with those of eutherian mice at two developmentally equivalent time points corresponding to deeper and upper layer neuron generation. Enrichment analyses revealed more mature gene networks in marsupials at the early stage, which reverted at the later stage, suggesting a more precocious but protracted neuronal maturation program relative to birth timing of cortical layers. We ranked genes expressed in different species and identified important differences in gene expression rankings between species. For example, genes known to be enriched in upper-layer cortical projection neuron subtypes, such as Cux1, Lhx2 and Satb2, likely relate to corpus callosum emergence in eutherians. These results show molecular heterochronies of neocortical development in Theria, and highlight changes in gene expression and cell type composition that may underlie neocortical evolution and diversification. This article has an associated 'The people behind the papers' interview.

Synapse-Enriched m6A-Modified Malat1 Interacts with the Novel m6A Reader, DPYSL2, and Is Required for Fear-Extinction Memory.

The RNA modification N6-methyladenosine (m6A) regulates the interaction between RNA and various RNA binding proteins within the nucleus and other subcellular compartments and has recently been shown to be involved in experience-dependent plasticity, learning, and memory. Using m6A RNA-sequencing, we have discovered a distinct population of learning-related m6A- modified RNAs at the synapse, which includes the long noncoding RNA metastasis-associated lung adenocarcinoma transcript 1 (Malat1). RNA immunoprecipitation and mass spectrometry revealed 12 new synapse-specific learning-induced m6A readers in the mPFC of male C57/BL6 mice, with m6A-modified Malat1 binding to a subset of these, including CYFIP2 and DPYSL2. In addition, a cell type- and synapse-specific, and state-dependent, reduction of m6A on Malat1 impairs fear-extinction memory; an effect that likely occurs through a disruption in the interaction between Malat1 and DPYSL2 and an associated decrease in dendritic spine formation. These findings highlight the critical role of m6A in regulating the functional state of RNA during the consolidation of fear-extinction memory, and expand the repertoire of experience-dependent m6A readers in the synaptic compartment.SIGNIFICANCE STATEMENT We have discovered that learning-induced m6A-modified RNA (including the long noncoding RNA, Malat1) accumulates in the synaptic compartment. We have identified several new m6A readers that are associated with fear extinction learning and demonstrate a causal relationship between m6A-modified Malat1 and the formation of fear-extinction memory. These findings highlight the role of m6A in regulating the functional state of an RNA during memory formation and expand the repertoire of experience-dependent m6A readers in the synaptic compartment.

Localised Cdr1as activity is required for fear extinction memory.

Circular RNAs (circRNAs) comprise a novel class of regulatory RNAs that are abundant in the brain, particularly within synapses. They are highly stable, dynamically regulated, and display a range of functions, including serving as decoys for microRNAs and proteins and, in some cases, circRNAs also undergo translation. Early work in animal models revealed an association between circRNAs and neurodegenerative and neuropsychiatric disorders; however, little is known about the link between circRNA function and memory. To address this, we examined circRNA in synaptosomes derived from the medial prefrontal cortex of fear extinction-trained male C57BL/6J mice and found 12,837 circRNAs that were enriched at the synapse, including cerebellar degeneration-related protein 1 antisense RNA (Cdr1as). Targeted knockdown of Cdr1as in the neural processes of the infralimbic cortex led to impaired fear extinction memory. These findings highlight the involvement of localised circRNA activity at the synapse in memory formation.

Splicing factor proline and glutamine rich intron retention, reduced expression and aggregate formation are pathological features of amyotrophic lateral sclerosis.

AIM: Splicing factor proline and glutamine rich (SFPQ) is an RNA-DNA binding protein that is dysregulated in Alzheimer's disease and frontotemporal dementia. Dysregulation of SFPQ, specifically increased intron retention and nuclear depletion, has been linked to several genetic subtypes of amyotrophic lateral sclerosis (ALS), suggesting that SFPQ pathology may be a common feature of this heterogeneous disease. Our study aimed to investigate this hypothesis by providing the first comprehensive assessment of SFPQ pathology in large ALS case-control cohorts. METHODS: We examined SFPQ at the RNA, protein and DNA levels. SFPQ RNA expression and intron retention were examined using RNA-sequencing and quantitative PCR. SFPQ protein expression was assessed by immunoblotting and immunofluorescent staining. At the DNA level, SFPQ was examined for genetic variation novel to ALS patients. RESULTS: At the RNA level, retention of SFPQ intron nine was significantly increased in ALS patients' motor cortex. In addition, SFPQ RNA expression was significantly reduced in the central nervous system, but not blood, of patients. At the protein level, neither nuclear depletion nor reduced expression of SFPQ was found to be a consistent feature of spinal motor neurons. However, SFPQ-positive ubiquitinated protein aggregates were observed in patients' spinal motor neurons. At the DNA level, our genetic screen identified two novel and two rare SFPQ sequence variants not previously reported in the literature. CONCLUSIONS: Our findings confirm dysregulation of SFPQ as a pathological feature of the central nervous system of ALS patients and indicate that investigation of the functional consequences of this pathology will provide insight into ALS biology.

Adolescent chronic intermittent toluene inhalation dynamically regulates the transcriptome and neuronal methylome within the rat medial prefrontal cortex.

Inhalants containing the volatile solvent toluene are misused to induce euphoria or intoxication. Inhalant abuse is most common during adolescence and can result in cognitive impairments during an important maturational period. Despite evidence suggesting that epigenetic modifications may underpin the cognitive effects of inhalants, no studies to date have thoroughly investigated toluene-induced regulation of the transcriptome or discrete epigenetic modifications within the brain. To address this, we investigated effects of adolescent chronic intermittent toluene (CIT) inhalation on gene expression and DNA methylation profiles within the rat medial prefrontal cortex (mPFC), which undergoes maturation throughout adolescence and has been implicated in toluene-induced cognitive deficits. Employing both RNA-seq and genome-wide Methyl CpG Binding Domain (MBD) Ultra-seq analysis, we demonstrate that adolescent CIT inhalation (10 000 ppm for 1 h/day, 3 days/week for 4 weeks) induces both transient and persistent changes to the transcriptome and DNA methylome within the rat mPFC for at least 2 weeks following toluene exposure. We demonstrate for the first time that adolescent CIT exposure results in dynamic regulation of the mPFC transcriptome likely relating to acute inflammatory responses and persistent deficits in synaptic plasticity. These adaptations may contribute to the cognitive deficits associated with chronic toluene exposure and provide novel molecular targets for preventing long-term neurophysiological abnormalities following chronic toluene inhalation.

Experience-Dependent Accumulation of N6-Methyladenosine in the Prefrontal Cortex Is Associated with Memory Processes in Mice.

UNLABELLED: The RNA modification N(6)-methyladenosine (m(6)A) influences mRNA stability and cell-type-specific developmental programming, and is highly abundant in the adult brain. However, it has not been determined whether m(6)A is dynamically regulated by experience. Based on transcriptome-wide profiling of m(6)A, we report that the level of m(6)A increases in the medial prefrontal cortex (mPFC) of mice in response to behavioral experience. The modulation was enriched near the stop codon of mRNAs, including genes related to neuronal plasticity. In primary cortical neurons, in vitro, modulation of m(6)A by the RNA demethylase FTO influenced the degradation profiles of a subset of transcripts with modulated sites. In vivo, the expression of Fto and the m(6)A methyltransferase, Mettl3 correlated with the observed increase in m(6)A levels post-training. Furthermore, targeted knockdown of FTO in the mPFC led to enhanced consolidation of cued fear memory. Thus, together with its role in early development, the dynamic regulation of m(6)A in the adult brain serves as an important epitranscriptomic mechanism associated with behavioral adaptation. SIGNIFICANCE STATEMENT: N(6)-methyladenosine (m(6)A) is the most prevalent internal modification on RNA, however, its cellular dynamics in vivo remains elusive. Here we provide the first demonstration of m(6)A upregulation in the mouse medial prefrontal cortex (mPFC) following behavioral training. Knocking down the m(6)A demethylase FTO in the mPFC, which increases total m(6)A level, results in enhanced consolidation of fear memory. Our findings suggest that m(6)A is regulated in an activity-dependent manner in the adult brain, and may function to fine-tune mRNA turnover during memory-related processes.

Neocortical Tet3-mediated accumulation of 5-hydroxymethylcytosine promotes rapid behavioral adaptation.

5-hydroxymethylcytosine (5-hmC) is a novel DNA modification that is highly enriched in the adult brain and dynamically regulated by neural activity. 5-hmC accumulates across the lifespan; however, the functional relevance of this change in 5-hmC and whether it is necessary for behavioral adaptation have not been fully elucidated. Moreover, although the ten-eleven translocation (Tet) family of enzymes is known to be essential for converting methylated DNA to 5-hmC, the role of individual Tet proteins in the adult cortex remains unclear. Using 5-hmC capture together with high-throughput DNA sequencing on individual mice, we show that fear extinction, an important form of reversal learning, leads to a dramatic genome-wide redistribution of 5-hmC within the infralimbic prefrontal cortex. Moreover, extinction learning-induced Tet3-mediated accumulation of 5-hmC is associated with the establishment of epigenetic states that promote gene expression and rapid behavioral adaptation.

Purification of neural precursor cells reveals the presence of distinct, stimulus-specific subpopulations of quiescent precursors in the adult mouse hippocampus.

The activity of neural precursor cells in the adult hippocampus is regulated by various stimuli; however, whether these stimuli regulate the same or different precursor populations remains unknown. Here, we developed a novel cell-sorting protocol that allows the purification to homogeneity of neurosphere-forming neural precursors from the adult mouse hippocampus and examined the responsiveness of individual precursors to various stimuli using a clonal assay. We show that within the Hes5-GFP(+)/Nestin-GFP(+)/EGFR(+) cell population, which comprises the majority of neurosphere-forming precursors, there are two distinct subpopulations of quiescent precursor cells, one directly activated by high-KCl depolarization, and the other activated by norepinephrine (NE). We then demonstrate that these two populations are differentially distributed along the septotemporal axis of the hippocampus, and show that the NE-responsive precursors are selectively regulated by GABA, whereas the KCl-responsive precursors are selectively modulated by corticosterone. Finally, based on RNAseq analysis by deep sequencing, we show that the progeny generated by activating NE-responsive versus KCl-responsive quiescent precursors are molecularly different. These results demonstrate that the adult hippocampus contains phenotypically similar but stimulus-specific populations of quiescent precursors, which may give rise to neural progeny with different functional capacity.

Global transcriptome and gene regulation network for secondary metabolite biosynthesis of tea plant (Camellia sinensis).

BACKGROUND: Major secondary metabolites, including flavonoids, caffeine, and theanine, are important components of tea products and are closely related to the taste, flavor, and health benefits of tea. Secondary metabolite biosynthesis in Camellia sinensis is differentially regulated in different tissues during growth and development. Until now, little was known about the expression patterns of genes involved in secondary metabolic pathways or their regulatory mechanisms. This study aimed to generate expression profiles for C. sinensis tissues and to build a gene regulation model of the secondary metabolic pathways. RESULTS: RNA sequencing was performed on 13 different tissue samples from various organs and developmental stages of tea plants, including buds and leaves of different ages, stems, flowers, seeds, and roots. A total of 43.7 Gbp of raw sequencing data were generated, from which 347,827 unigenes were assembled and annotated. There were 46,693, 8446, 3814, 10,206, and 4948 unigenes specifically expressed in the buds and leaves, stems, flowers, seeds, and roots, respectively. In total, 1719 unigenes were identified as being involved in the secondary metabolic pathways in C. sinensis, and the expression patterns of the genes involved in flavonoid, caffeine, and theanine biosynthesis were characterized, revealing the dynamic nature of their regulation during plant growth and development. The possible transcription factor regulation network for the biosynthesis of flavonoid, caffeine, and theanine was built, encompassing 339 transcription factors from 35 families, namely bHLH, MYB, and NAC, among others. Remarkably, not only did the data reveal the possible critical check points in the flavonoid, caffeine, and theanine biosynthesis pathways, but also implicated the key transcription factors and related mechanisms in the regulation of secondary metabolite biosynthesis. CONCLUSIONS: Our study generated gene expression profiles for different tissues at different developmental stages in tea plants. The gene network responsible for the regulation of the secondary metabolic pathways was analyzed. Our work elucidated the possible cross talk in gene regulation between the secondary metabolite biosynthetic pathways in C. sinensis. The results increase our understanding of how secondary metabolic pathways are regulated during plant development and growth cycles, and help pave the way for genetic selection and engineering for germplasm improvement.

Transcriptomic analyses reveal the adaptive features and biological differences of guts from two invasive whitefly species.

BACKGROUND: The gut of phloem feeding insects is critical for nutrition uptake and xenobiotics degradation. However, partly due to its tiny size, genomic information for the gut of phloem feeding insects is limited. RESULTS: In this study, the gut transcriptomes of two species of invasive whiteflies in the Bemisia tabaci complex, Middle East Asia Minor 1 (MEAM1) and Mediterranean (MED), were analyzed using the Illumina sequencing. A total of 12,879 MEAM1 transcripts and 11,246 MED transcripts were annotated with a significant Blastx hit. In addition, 7,000 and 5,771 gut specific genes were respectively identified for MEAM1 and MED. Functional analyses on these gut specific genes demonstrated the important roles of gut in metabolism of insecticides and secondary plant chemicals. To reveal the molecular difference between guts of MEAM1 and MED, a comparison between gut transcriptomes of the two species was conducted and 3,910 pairs of orthologous genes were identified. Based on the ratio of nonsynonymous and synonymous substitutions, 15 genes were found evolving under positive selection. Many of those genes are predicted to be involved in metabolism and insecticide resistance. Furthermore, many genes related to detoxification were expressed at an elevated level in the gut of MED compared to MEAM1, which might be responsible for the MED's higher resistance to insecticides and environmental stresses. CONCLUSION: The sequencing of MED and MEAM1 gut transcriptomes and extensive comparisons of MEAM1 and MED gut transcripts provide substantial sequence information for revealing the role of gut in whiteflies.

Alternative splicing and trans-splicing events revealed by analysis of the Bombyx mori transcriptome.

Alternative splicing and trans-splicing events have not been systematically studied in the silkworm Bombyx mori. Here, the silkworm transcriptome was analyzed by RNA-seq. We identified 320 novel genes, modified 1140 gene models, and found thousands of alternative splicing and 58 trans-splicing events. Studies of three SR proteins show that both their alternative splicing patterns and mRNA products are conserved from insect to human, and one isoform of Srsf6 with a retained intron is expressed sex-specifically in silkworm gonads. Trans-splicing of mod(mdg4) in silkworm was experimentally confirmed. We identified integrations from a common 5'-gene with 46 newly identified alternative 3'-exons that are located on both DNA strands over a 500-kb region. Other trans-splicing events in B. mori were predicted by bioinformatic analysis, in which 12 events were confirmed by RT-PCR, six events were further validated by chimeric SNPs, and two events were confirmed by allele-specific RT-PCR in F(1) hybrids from distinct silkworm lines of JS and L10, indicating that trans-splicing is more widespread in insects than previously thought. Analysis of the B. mori transcriptome by RNA-seq provides valuable information of regulatory alternative splicing events. The conservation of splicing events across species and newly identified trans-splicing events suggest that B. mori is a good model for future studies.

Wnt dose escalation during the exit from pluripotency identifies tranilast as a regulator of cardiac mesoderm.

Wnt signaling is a critical determinant of cell lineage development. This study used Wnt dose-dependent induction programs to gain insights into molecular regulation of stem cell differentiation. We performed single-cell RNA sequencing of hiPSCs responding to a dose escalation protocol with Wnt agonist CHIR-99021 during the exit from pluripotency to identify cell types and genetic activity driven by Wnt stimulation. Results of activated gene sets and cell types were used to build a multiple regression model that predicts the efficiency of cardiomyocyte differentiation. Cross-referencing Wnt-associated gene expression profiles to the Connectivity Map database, we identified the small-molecule drug, tranilast. We found that tranilast synergistically activates Wnt signaling to promote cardiac lineage differentiation, which we validate by in vitro analysis of hiPSC differentiation and in vivo analysis of developing quail embryos. Our study provides an integrated workflow that links experimental datasets, prediction models, and small-molecule databases to identify drug-like compounds that control cell differentiation.

The characteristics and expression profiles of the mitochondrial genome for the Mediterranean species of the Bemisia tabaci complex.

BACKGROUND: The whiteflies under the name Bemisia tabaci (Gennadius) (Aleyrodidae: Hemiptera) are species complex of at least 31 cryptic species some of which are globally invasive agricultural pests. Previously, the mitochondrial genome (mitogenome) of the indigenous New World B. tabaci species was sequenced and major differences of gene order from the postulated whitefly ancestral gene order were found. However, the sequence and gene order of mitogenomes in other B. tabaci species are unknown. In addition, the sequence divergences and gene expression profiles of mitogenomes in the B. tabaci species complex remain completely unexplored. RESULTS: In this study, we obtained the complete mitogenome (15,632 bp) of the invasive Mediterranean (MED), which has been identified as the type species of the B. tabaci complex. It encodes 37 genes, including 13 protein-coding genes (PCGs), 2 ribosomal RNAs and 22 transfer RNAs (tRNA). Comparative analyses of the mitogenomes from MED and New World (previously published) species reveal that there are no gene arrangements. Based on the Illumina sequencing data, the gene expression profile of the MED mitogenome was analyzed. We found that a number of genes were polyadenylated and the partial stop codons in cox1, cox2 and nd5 are completed via polyadenylation that changed T to the TAA stop codon. In addition, combining the transcriptome with the sequence alignment data, the possible termination site of some PCGs were defined. Our analyses also revealed that atp6 and atp8, nd4 and nd4l, nd6 and cytb were found on the same cistronic transcripts, whereas the other mature mitochondrial transcripts were monocistronic. Furthermore, RT-PCR analyses of the mitochondrial PCGs expression in different developmental stages revealed that the expression level of individual mitochondrial genes varied in each developmental stage of nymph, pupa and adult. Interestingly, mRNA levels showed significant differences among genes located in the same transcription unit suggesting that mitochondrial mRNA abundance is heavily modulated by post-transcriptional regulation. CONCLUSIONS: This work provides novel insights into the mitogenome evolution of B. tabaci species and demonstrates that utilizing RNA-seq data to obtain the mitogenome and analyze mitochondrial gene expression characteristics is practical.

Global transcriptome profiles of Camellia sinensis during cold acclimation.

BACKGROUND: Tea is the most popular non-alcoholic health beverage in the world. The tea plant (Camellia sinensis (L.) O. Kuntze) needs to undergo a cold acclimation process to enhance its freezing tolerance in winter. Changes that occur at the molecular level in response to low temperatures are poorly understood in tea plants. To elucidate the molecular mechanisms of cold acclimation, we employed RNA-Seq and digital gene expression (DGE) technologies to the study of genome-wide expression profiles during cold acclimation in tea plants. RESULTS: Using the Illumina sequencing platform, we obtained approximately 57.35 million RNA-Seq reads. These reads were assembled into 216,831 transcripts, with an average length of 356 bp and an N50 of 529 bp. In total, 1,770 differentially expressed transcripts were identified, of which 1,168 were up-regulated and 602 down-regulated. These include a group of cold sensor or signal transduction genes, cold-responsive transcription factor genes, plasma membrane stabilization related genes, osmosensing-responsive genes, and detoxification enzyme genes. DGE and quantitative RT-PCR analysis further confirmed the results from RNA-Seq analysis. Pathway analysis indicated that the "carbohydrate metabolism pathway" and the "calcium signaling pathway" might play a vital role in tea plants' responses to cold stress. CONCLUSIONS: Our study presents a global survey of transcriptome profiles of tea plants in response to low, non-freezing temperatures and yields insights into the molecular mechanisms of tea plants during the cold acclimation process. It could also serve as a valuable resource for relevant research on cold-tolerance and help to explore the cold-related genes in improving the understanding of low-temperature tolerance and plant-environment interactions.

Experimentally induced active and quiet sleep engage non-overlapping transcriptional programs in Drosophila.

Sleep in mammals can be broadly classified into two different physiological categories: rapid eye movement (REM) sleep and slow-wave sleep (SWS), and accordingly REM and SWS are thought to achieve a different set of functions. The fruit fly Drosophila melanogaster is increasingly being used as a model to understand sleep functions, although it remains unclear if the fly brain also engages in different kinds of sleep as well. Here, we compare two commonly used approaches for studying sleep experimentally in Drosophila: optogenetic activation of sleep-promoting neurons and provision of a sleep-promoting drug, gaboxadol. We find that these different sleep-induction methods have similar effects on increasing sleep duration, but divergent effects on brain activity. Transcriptomic analysis reveals that drug-induced deep sleep ('quiet' sleep) mostly downregulates metabolism genes, whereas optogenetic 'active' sleep upregulates a wide range of genes relevant to normal waking functions. This suggests that optogenetics and pharmacological induction of sleep in Drosophila promote different features of sleep, which engage different sets of genes to achieve their respective functions.

Experimentally induced active and quiet sleep engage non-overlapping transcriptional programs in Drosophila.

Sleep in mammals can be broadly classified into two different physiological categories: rapid eye movement (REM) sleep and slow wave sleep (SWS), and accordingly REM and SWS are thought to achieve a different set of functions. The fruit fly Drosophila melanogaster is increasingly being used as a model to understand sleep functions, although it remains unclear if the fly brain also engages in different kinds of sleep as well. Here, we compare two commonly used approaches for studying sleep experimentally in Drosophila: optogenetic activation of sleep-promoting neurons and provision of a sleep-promoting drug, Gaboxadol. We find that these different sleep-induction methods have similar effects on increasing sleep duration, but divergent effects on brain activity. Transcriptomic analysis reveals that drug-induced deep sleep ('quiet' sleep) mostly downregulates metabolism genes, whereas optogenetic 'active' sleep upregulates a wide range of genes relevant to normal waking functions. This suggests that optogenetics and pharmacological induction of sleep in Drosophila promote different features of sleep, which engage different sets of genes to achieve their respective functions.

Fear extinction is regulated by the activity of long noncoding RNAs at the synapse.

Long noncoding RNAs (lncRNAs) represent a multidimensional class of regulatory molecules that are involved in many aspects of brain function. Emerging evidence indicates that lncRNAs are localized to the synapse; however, a direct role for their activity in this subcellular compartment in memory formation has yet to be demonstrated. Using lncRNA capture-seq, we identified a specific set of lncRNAs that accumulate in the synaptic compartment within the infralimbic prefrontal cortex of adult male C57/Bl6 mice. Among these was a splice variant related to the stress-associated lncRNA, Gas5. RNA immunoprecipitation followed by mass spectrometry and single-molecule imaging revealed that this Gas5 isoform, in association with the RNA binding proteins G3BP2 and CAPRIN1, regulates the activity-dependent trafficking and clustering of RNA granules. In addition, we found that cell-type-specific, activity-dependent, and synapse-specific knockdown of the Gas5 variant led to impaired fear extinction memory. These findings identify a new mechanism of fear extinction that involves the dynamic interaction between local lncRNA activity and RNA condensates in the synaptic compartment.

Analysis of a native whitefly transcriptome and its sequence divergence with two invasive whitefly species.

BACKGROUND: Genomic divergence between invasive and native species may provide insight into the molecular basis underlying specific characteristics that drive the invasion and displacement of closely related species. In this study, we sequenced the transcriptome of an indigenous species, Asia II 3, of the Bemisia tabaci complex and compared its genetic divergence with the transcriptomes of two invasive whiteflies species, Middle East Asia Minor 1 (MEAM1) and Mediterranean (MED), respectively. RESULTS: More than 16 million reads of 74 base pairs in length were obtained for the Asia II 3 species using the Illumina sequencing platform. These reads were assembled into 52,535 distinct sequences (mean size: 466 bp) and 16,596 sequences were annotated with an E-value above 10-5. Protein family comparisons revealed obvious diversification among the transcriptomes of these species suggesting species-specific adaptations during whitefly evolution. On the contrary, substantial conservation of the whitefly transcriptomes was also evident, despite their differences. The overall divergence of coding sequences between the orthologous gene pairs of Asia II 3 and MEAM1 is 1.73%, which is comparable to the average divergence of Asia II 3 and MED transcriptomes (1.84%) and much higher than that of MEAM1 and MED (0.83%). This is consistent with the previous phylogenetic analyses and crossing experiments suggesting these are distinct species. We also identified hundreds of highly diverged genes and compiled sequence identify data into gene functional groups and found the most divergent gene classes are Cytochrome P450, Glutathione metabolism and Oxidative phosphorylation. These results strongly suggest that the divergence of genes related to metabolism might be the driving force of the MEAM1 and Asia II 3 differentiation. We also analyzed single nucleotide polymorphisms within the orthologous gene pairs of indigenous and invasive whiteflies which are helpful for the investigation of association between allelic and phenotypes. CONCLUSIONS: Our data present the most comprehensive sequences for the indigenous whitefly species Asia II 3. The extensive comparisons of Asia II 3, MEAM1 and MED transcriptomes will serve as an invaluable resource for revealing the genetic basis of whitefly invasion and the molecular mechanisms underlying their biological differences.

ADRAM is an experience-dependent long noncoding RNA that drives fear extinction through a direct interaction with the chaperone protein 14-3-3.

Here, we used RNA capture-seq to identify a large population of lncRNAs that are expressed in the infralimbic prefrontal cortex of adult male mice in response to fear-related learning. Combining these data with cell-type-specific ATAC-seq on neurons that had been selectively activated by fear extinction learning, we find inducible 434 lncRNAs that are derived from enhancer regions in the vicinity of protein-coding genes. In particular, we discover an experience-induced lncRNA we call ADRAM (activity-dependent lncRNA associated with memory) that acts as both a scaffold and a combinatorial guide to recruit the brain-enriched chaperone protein 14-3-3 to the promoter of the memory-associated immediate-early gene Nr4a2 and is required fear extinction memory. This study expands the lexicon of experience-dependent lncRNA activity in the brain and highlights enhancer-derived RNAs (eRNAs) as key players in the epigenomic regulation of gene expression associated with the formation of fear extinction memory.