Temporal Transcriptome Profiling of Pinus densiflora Infected with Pine Wood Nematode Reveals Genetically Programmed Changes upon Pine Wilt Disease.

Pine, an evergreen conifer, is widely distributed worldwide. It is economically, scientifically, and ecologically important. However, pine wilt disease (PWD) induced by the pine wood nematode (PWN) adversely affects pine trees. Many studies have been conducted on the PWN and its beetle vectors to prevent the spread of PWD. However, studies providing a comprehensive understanding of the pine tree transcriptome in response to PWN infection are lacking. Here, we performed temporal profiling of the pine tree transcriptome using PWD-infected red pine trees, Pinus densiflora, inoculated with the PWN by RNA sequencing. Our analysis revealed that defense-responsive genes involved in cell wall modification, jasmonic acid signaling, and phenylpropanoid-related processes were significantly enriched 2 weeks after PWD infection. Furthermore, some WRKY-type and MYB-type transcription factors were upregulated 2 weeks after PWD infection, suggesting that these transcription factors might be responsible for the genome-wide reprogramming of defense-responsive genes in the early PWD stage. Our comprehensive transcriptome analysis will assist in developing PWD-resistant pine trees and identifying genes to diagnose PWD at the early stage of infection, during which large-scale phenotypic changes are absent in PWD-infected pine trees.

Fetal neurodevelopmental spatio-temporal dynamic transcriptional landscape of maternal insult-induce autism spectrum disorder risk.

Maternal insults during pregnancy induces an increased risk of autism spectrum disorders (ASD) in offspring, but the neuropathological changes in this process remains not to be established. To shed light on this, the transcriptome datasets of maternal blood samples with children later diagnosed with ASD and typical development, and tissue samples of multiple brain regions from ASD patients and human neurodevelopment were conducted to identify the non-chasm differentially expressed genes (DEGs) to generate the spatio-temporal dynamic change. Combined enrichment and interaction network analysis revealed that non-chasm DEGs with similar expression trajectories in the same brain regions, were involved in neural, immune and metabolic GO functions and KEGG pathways, respectively, suggesting that did not performed exactly the same functions. Interestingly, our results found that non-chasm DEGs in frontal cortex and temporal cortex were associated with COVID-19, suggesting that as an environmental risk factor COVID-19 affects an increased risk of ASD.

Exogenous Melatonin Improves Seed Germination of Wheat (Triticum aestivum L.) under Salt Stress.

Melatonin (MT) can effectively reduce oxidative damage induced by abiotic stresses such as salt in plants. However, the effects of MT on physiological responses and molecular regulation during wheat germination remains largely elusive. In this study, the response of wheat seeds to MT under salt stress during germination was investigated at physiological and transcriptome levels. Our results revealed that application of MT significantly reduced the negative influence of salt stress on wheat seed germination. The oxidative load was reduced by inducing high activities of antioxidant enzymes. In parallel, the content of gibberellin A3 (GA3) and jasmonic acid (JA) increased in MT-treated seedling. RNA-seq analysis demonstrated that MT alters oxidoreductase activity and phytohormone-dependent signal transduction pathways under salt stress. Weighted correlation network analysis (WGCNA) revealed that MT participates in enhanced energy metabolism and protected seeds via maintained cell morphology under salt stress during wheat seed germination. Our findings provide a conceptual basis of the MT-mediated regulatory mechanism in plant adaptation to salt stress, and identify the potential candidate genes for salt-tolerant wheat molecular breeding.

Full-Length Transcriptomic Sequencing and Temporal Transcriptome Expression Profiling Analyses Offer Insights into Terpenoid Biosynthesis in Artemisia argyi.

Artemisiae argyi Folium is a traditional herbal medicine used for moxibustion heat therapy in China. The volatile oils in A.argyi leaves are closely related to its medicinal value. Records suggest that the levels of these terpenoids components within the leaves vary as a function of harvest time, with June being the optimal time for A. argyi harvesting, owing to the high levels of active ingredients during this month. However, the molecular mechanisms governing terpenoid biosynthesis and the time-dependent changes in this activity remain unclear. In this study, GC-MS analysis revealed that volatile oil levels varied across four different harvest months (April, May, June, and July) in A. argyi leaves, and the primarily terpenoids components (including both monoterpenes and sesquiterpenes) reached peak levels in early June. Through single-molecule real-time (SMRT) sequencing, corrected by Illumina RNA-sequencing (RNA-Seq), 44 full-length transcripts potentially involved in terpenoid biosynthesis were identified in this study. Differentially expressed genes (DEGs) exhibiting time-dependent expression patterns were divided into 12 coexpression clusters. Integrated chemical and transcriptomic analyses revealed distinct time-specific transcriptomic patterns associated with terpenoid biosynthesis. Subsequent hierarchical clustering and correlation analyses ultimately identified six transcripts that were closely linked to the production of these two types of terpenoid within A. argyi leaves, revealing that the structural diversity of terpenoid is related to the generation of the diverse terpene skeletons by prenyltransferase (TPS) family of enzymes. These findings can guide further studies of the molecular mechanisms underlying the quality of A. argyi leaves, aiding in the selection of optimal timing for harvests of A. argyi.

A spatial-temporal understanding of gene regulatory networks and NtARF-mediated regulation of potassium accumulation in tobacco.

MAIN CONCLUSION: After tobacco topping, changes in the auxin content could affect K+ uptake by inhibiting the activity of K+ uptake-related genes through the NtARF genes, thus causing changes in K+ content. Tobacco (Nicotiana tabacum) is a valuable industrial and commercial crop, and the leaf is its primary product. Topping (removing apical buds) is a common agronomic practice that significantly improves the yield of tobacco leaves. Potassium (K+) plays an important physiological role in tobacco growth and leaf traits, including combustibility, aroma, and safety in cigarette products, and its levels are significantly decreased after topping. Here, to present global spatial-temporal gene expression profiles and gene regulatory networks of the core elements of K+ uptake, leaves and roots from topped and untopped plants at short- and long-term time points after topping were sampled for transcriptome analysis. We found that the wounding response was initiated in leaves in the early stages after topping. Then, in the long term, processes related to metabolism and transcription regulation, as well as ion binding and transport, were altered. The expression profiles showed that core elements of K+ uptake and xylem loading were drastically suppressed in roots after topping. Finally, transient expression experiments confirmed that changes in the auxin content could affect K+ uptake by inhibiting the activity of K+ uptake-related genes through the tobacco auxin response factor (NtARF) genes, thus causing changes in the K+ content. These results suggest that some ARFs could be selected as targets to enhance the expressions of K+ uptake transporters, leading to increment of K+ contents and improvement of leaf quality in tobacco breeding.

Endoplasmic reticulum stress pathway mediates the early heat stress response of developing rice seeds.

A transient heat stress occurring during early seed development in rice (Oryza sativa) reduces seed size by altering endosperm development. However, the relationship between the timing of the stress and specific developmental stage on heat sensitivity is not well-understood. To address this, we imposed a series of non-overlapping heat stress treatments and found that young seeds are most sensitive during the first two days after flowering. Temporal transcriptome analysis of developing, heat stressed (35°C) seeds during this window shows that Inositol-requiring enzyme 1 (IRE1)-mediated endoplasmic reticulum (ER) stress response and jasmonic acid (JA) pathways are the early (1-3 h) drivers of heat stress response. We propose that increased JA levels under heat stress may precede ER stress response as JA application promotes the spliced form of OsbZIP50, an ER response marker gene linked to IRE1-specific pathway. This study presents temporal and mechanistic insights into the role of JA and ER stress signalling during early heat stress response of rice seeds that impact both grain size and quality. Modulating the heat sensitivity of the early sensing pathways and downstream endosperm development genes can enhance rice resilience to transient heat stress events.

Hybrid sequencing reveals insight into heat sensing and signaling of bread wheat.

Wheat (Triticum aestivum L.), a globally important crop, is challenged by increasing temperatures (heat stress, HS). However its polyploid nature, the incompleteness of its genome sequences and annotation, the lack of comprehensive HS-responsive transcriptomes and the unexplored heat sensing and signaling of wheat hinder our full understanding of its adaptations to HS. The recently released genome sequences of wheat, as well as emerging single-molecular sequencing technologies, provide an opportunity to thoroughly investigate the molecular mechanisms of the wheat response to HS. We generated a high-resolution spatio-temporal transcriptome map of wheat flag leaves and filling grain under HS at 0 min, 5 min, 10 min, 30 min, 1 h and 4 h by combining full-length single-molecular sequencing and Illumina short reads sequencing. This hybrid sequencing newly discovered 4947 loci and 70 285 transcripts, generating the comprehensive and dynamic list of HS-responsive full-length transcripts and complementing the recently released wheat reference genome. Large-scale analysis revealed a global landscape of heat adaptations, uncovering unexpected rapid heat sensing and signaling, significant changes of more than half of HS-responsive genes within 30 min, heat shock factor-dependent and -independent heat signaling, and metabolic alterations in early HS-responses. Integrated analysis also demonstrated the differential responses and partitioned functions between organs and subgenomes, and suggested a differential pattern of transcriptional and alternative splicing regulation in the HS response. This study provided comprehensive data for dissecting molecular mechanisms of early HS responses in wheat and highlighted the genomic plasticity and evolutionary divergence of polyploidy wheat.

Temporal transcriptome profiling of developing seeds reveals a concerted gene regulation in relation to oil accumulation in Pongamia (Millettia pinnata).

BACKGROUND: Pongamia (Millettia pinnata syn. Pongamia pinnata), an oilseed legume species, is emerging as potential feedstock for sustainable biodiesel production. Breeding Pongamia for favorable traits in commercial application will rely on a comprehensive understanding of molecular mechanism regulating oil accumulation during its seed development. To date, only limited genomic or transcript sequences are available for Pongamia, while a temporal transcriptome profiling of developing seeds is still lacking in this species. RESULTS: In this work, we conducted a time-series analysis of morphological and physiological characters, oil contents and compositions, as well as global gene expression profiles in developing Pongamia seeds. Firstly, three major developmental phases were characterized based on the combined evidences from embryonic shape, seed weight, seed moisture content, and seed color. Then, the gene expression levels at these three phases were quantified by RNA-Seq analyses with three biological replicates from each phase. Nearly 94% of unigenes were expressed at all three phases, whereas only less than 2% of unigenes were exclusively expressed at one of these phases. A total of 8881 differentially expressed genes (DEGs) were identified between phases. Furthermore, the qRT-PCR analyses for 10 DEGs involved in lipid metabolism demonstrated a good reliability of our RNA-Seq data in temporal gene expression profiling. We observed a dramatic increase in seed oil content from the embryogenesis phase to the early seed-filling phase, followed by a steady and moderate increase towards the maximum at the desiccation phase. We proposed that a highly active expression of most genes related to fatty acid (FA) and triacylglycerol (TAG) biosynthesis at the embryogenesis phase might trigger both the substantial oil accumulation and the membrane lipid synthesis for rapid cell proliferation at this phase, while a concerted reactivation of TAG synthesis-related genes at the desiccation phase might further promote storage lipid synthesis to achieve the maximum content of seed oils. CONCLUSIONS: This study not only built a bridge between gene expression profiles and oil accumulation in developing seeds, but also laid a foundation for future attempts on genetic engineering of Pongamia varieties to acquire higher oil yield or improved oil properties for biofuel applications.

Identification of novel anti-ZIKV drugs from viral-infection temporal gene expression profiles.

Zika virus (ZIKV) infections are typically asymptomatic but cause severe neurological complications (e.g. Guillain-Barré syndrome in adults, and microcephaly in newborns). There are currently no specific therapy or vaccine options available to prevent ZIKV infections. Temporal gene expression profiles of ZIKV-infected human brain microvascular endothelial cells (HBMECs) were used in this study to identify genes essential for viral replication. These genes were then used to identify novel anti-ZIKV agents and validated in publicly available data and functional wet-lab experiments. Here, we found that ZIKV effectively evaded activation of immune response-related genes and completely reprogrammed cellular transcriptional architectures. Knockdown of genes, which gradually upregulated during viral infection but showed distinct expression patterns between ZIKV- and mock infection, discovered novel proviral and antiviral factors. One-third of the 74 drugs found through signature-based drug repositioning and cross-reference with the Drug Gene Interaction Database (DGIdb) were known anti-ZIKV agents. In cellular assays, two promising antiviral candidates (Luminespib/NVP-AUY922, L-161982) were found to reduce viral replication without causing cell toxicity. Overall, our time-series transcriptome-based methods offer a novel and feasible strategy for antiviral drug discovery. Our strategies, which combine conventional and data-driven analysis, can be extended for other pathogens causing pandemics in the future.

Temporal transcriptome analysis reveals several key pathways involve in cadmium stress response in Nicotiana tabacum L.

Tobacco has a strong cadmium (Cd) enrichment capacity, meaning that it can absorb large quantities from the environment, but too much Cd will cause damage to the plant. It is not yet clear how the plant can dynamically respond to Cd stress. Here, we performed a temporal transcriptome analysis of tobacco roots under Cd treatment from 0 to 48 h. The number of differentially expressed genes (DEGs) was found to change significantly at 3 h of Cd treatment, which we used to define the early and middle stages of the Cd stress response. The gene ontology (GO) term analysis indicates that genes related to photosynthesis and fatty acid synthesis were enriched during the early phases of the stress response, and in the middle phase biological process related to metal ion transport, DNA damage repair, and metabolism were enriched. It was also found that plants use precursor mRNA (pre-mRNA) processes to first resist Cd stress, and with the increasing of Cd treatment time, the overlapped genes number of DEGs and DAS increased, suggesting the transcriptional levels and post-transcriptional level might influence each other. This study allowed us to better understand how plants dynamically respond to cadmium stress at the transcriptional and post-transcriptional levels and provided a reference for the screening of Cd-tolerant genes in the future.

Temporal transcriptome highlights the involvement of cytokine/JAK/STAT3 signaling pathway in the osteoinduction of BMSCs.

BACKGROUND: Mesenchymal stem cells (MSCs)-based therapy offers an effective strategy for bone regeneration to solve the clinical orthopedic problems. However, the transcriptional regulation of multiple transitional stages of continuous osteogenesis from MSCs has not been fully characterized. METHODS: Bone marrow mesenchymal stem cells (BMSCs) stimulated with osteogenic induction media were utilized to construct the in vitro osteogenic differentiation model. BMSCs were harvested after induction for 0, 7, 14 and 21 days, respectively, to perform the mRNA-sequencing (mRNA-Seq). The transcription factor networks and common molecules during the osteogenesis were revealed by using the temporal transcriptome. Further verification was performed by the quantitative real-time polymerase chain reaction (qRT-PCR), immunofluorescence and Western blotting. RESULTS: It showed that BMSCs could differentiate into osteogenic, and crucial regulator in the MAPK signaling pathway, the PPAR signaling pathway, the Toll-like receptor signaling and the Cytokine/JAK/STAT signaling pathway. PPI protein interaction analysis also suggested that three cytokines are involved in osteogenic differentiation as core genes, including leukemia inhibitory factor (LIF), interleukin-6 (IL6) and colony-stimulating factor 3 (CSF3). The osteogenic process was negatively affected by the inhibition of JAK/STAT3 signaling pathway. CONCLUSIONS: This work might provide new insights in the crucial features of the transcriptional regulation during the osteogenesis, as well as offer important clues about the activity and regulation of the relatively long-activated Cytokine/JAK/STAT3 signaling pathway in osteoinduction of BMSCs.

Dynamic Regulatory Event Mining by iDREM in Large-Scale Multi-omics Datasets During Biotic and Abiotic Stress in Plants.

The system-wide complexity of genome regulation encoding the organism phenotypic diversity is well understood. However, a major challenge persists about the appropriate method to describe the systematic dynamic genome regulation event utilizing enormous multi-omics datasets. Here, we describe Interactive Dynamic Regulatory Events Miner (iDREM) which reconstructs gene-regulatory networks from temporal transcriptome, proteome, and epigenome datasets during stress to envisage "master" regulators by simulating cascades of temporal transcription-regulatory and interactome events. The iDREM is a Java-based software that integrates static and time-series transcriptomics and proteomics datasets, transcription factor (TF)-target interactions, microRNA (miRNA)-target interaction, and protein-protein interactions to reconstruct temporal regulatory network and identify significant regulators in an unsupervised manner. The hidden Markov model detects specialized manipulated pathways as well as genes to recognize statistically significant regulators (TFs/miRNAs) that diverge in temporal activity. This method can be translated to any biotic or abiotic stress in plants and animals to predict the master regulators from condition-specific multi-omics datasets including host-pathogen interactions for comprehensive understanding of manipulated biological pathways.

NOCICEPTRA: Gene and microRNA Signatures and Their Trajectories Characterizing Human iPSC-Derived Nociceptor Maturation.

Nociceptors are primary afferent neurons serving the reception of acute pain but also the transit into maladaptive pain disorders. Since native human nociceptors are hardly available for mechanistic functional research, and rodent models do not necessarily mirror human pathologies in all aspects, human induced pluripotent stem cell-derived nociceptors (iDN) offer superior advantages as a human model system. Unbiased mRNA::microRNA co-sequencing, immunofluorescence staining, and qPCR validations, reveal expression trajectories as well as miRNA target spaces throughout the transition of pluripotent cells into iDNs. mRNA and miRNA candidates emerge as regulatory hubs for neurite outgrowth, synapse development, and ion channel expression. The exploratory data analysis tool NOCICEPTRA is provided as a containerized platform to retrieve experimentally determined expression trajectories, and to query custom gene sets for pathway and disease enrichments. Querying NOCICEPTRA for marker genes of cortical neurogenesis reveals distinct similarities and differences for cortical and peripheral neurons. The platform provides a public domain neuroresource to exploit the entire data sets and explore miRNA and mRNA as hubs regulating human nociceptor differentiation and function.