Full-length transcriptome sequences of ephemeral plant Arabidopsis pumila provides insight into gene expression dynamics during continuous salt stress

The ldp1 Mutation Affects the Expression of Auxin-Related Genes and Enhances SAM Size in Rice

Panicle type is one of the important factors affecting rice (Oryza sativa L.) yield, and the identification of regulatory genes in panicle development can provide significant insights into the molecular network involved. This study identified a large and dense panicle 1 (ldp1) mutant produced from the Wuyunjing 7 (WYJ7) genotype, which displayed significant relative increases in panicle length, number of primary and secondary branches, number of grains per panicle, grain width, and grain yield per plant. Scanning electron microscopy results showed that the shoot apical meristem (SAM) of ldp1 was relatively larger at the bract stage (BM), with a significantly increased number of primary (PBM) and secondary branch (SBM) meristematic centers, indicating that the ldp1 mutation affects early stages in SAM development Comparative RNA-Seq analysis of meristem tissues from WYJ7 and ldp1 at the BM, PBM, and SBM developmental stages indicated that the number of differentially expressed genes (DEGs) were highest (1407) during the BM stage. Weighted gene coexpression network analysis (WGCNA) revealed that genes in one module (turquoise) are associated with the ldp1 phenotype and highly expressed during the BM stage, suggesting their roles in the identity transition and branch differentiation stages of rice inflorescences. Hub genes involved in auxin synthesis and transport pathways, such as OsAUX1, OsAUX4, and OsSAUR25, were identified. Moreover, GO and KEGG analysis of the DEGs in the turquoise module and the 1407 DEGs in the BM stage revealed that a majority of genes involved in tryptophan metabolism and auxin signaling pathway were differentially expressed between WYJ and ldp1. The genetic analysis indicated that the ldp1 phenotype is controlled by a recessive monogene (LDP1), which was mapped to a region between 16.9 and 18.1 Mb on chromosome seven. This study suggests that the ldp1 mutation may affect the expression of key genes in auxin synthesis and signal transduction, enhance the size of SAM, and thus affect panicle development. This study provides insights into the molecular regulatory network underlying rice panicle morphogenesis and lays an important foundation for further understanding the function and molecular mechanism of LDP1 during panicle development.

Molecular evolution and interaction of 14-3-3 proteins with H+-ATPases in plant abiotic stresses

Environmental stresses severely affect plant growth and crop productivity. Regulated by 14-3-3 proteins (14-3-3s), H+-ATPases (AHAs) are important proton pumps that can induce diverse secondary transport via channels and co-transporters for the abiotic stress response of plants. Many studies demonstrated the roles of 14-3-3s and AHAs in coordinating the processes of plant growth, phytohormone signaling, and stress responses. However, the molecular evolution of 14-3-3s and AHAs has not been summarized in parallel with evolutionary insights across multiple plant species. Here, we comprehensively review the roles of 14-3-3s and AHAs in cell signaling to enhance plant responses to diverse environmental stresses. We analyzed the molecular evolution of key proteins and functional domains that are associated with 14-3-3s and AHAs in plant growth and hormone signaling. The results revealed evolution, duplication, contraction, and expansion of 14-3-3s and AHAs in green plants. We also discussed the stress-specific expression of those 14-3-3and AHA genes in a eudicotyledon (Arabidopsis thaliana), a monocotyledon (Hordeum vulgare), and a moss (Physcomitrium patens) under abiotic stresses. We propose that 14-3-3s and AHAs respond to abiotic stresses through many important targets and signaling components of phytohormones, which could be promising to improve plant tolerance to single or multiple environmental stresses.

CAP-RNAseq: an integrated pipeline for functional annotation and prioritization of co-expression clusters

Cluster analysis is one of the most widely used exploratory methods for visualization and grouping of gene expression patterns across multiple samples or treatment groups. Although several existing online tools can annotate clusters with functional terms, there is no all-in-one webserver to effectively prioritize genes/clusters using gene essentiality as well as congruency of mRNA-protein expression. Hence, we developed CAP-RNAseq that makes possible (1) upload and clustering of bulk RNA-seq data followed by identification, annotation and network visualization of all or selected clusters; and (2) prioritization using DepMap gene essentiality and/or dependency scores as well as the degree of correlation between mRNA and protein levels of genes within an expression cluster. In addition, CAP-RNAseq has an integrated primer design tool for the prioritized genes. Herein, we showed using comparisons with the existing tools and multiple case studies that CAP-RNAseq can uniquely aid in the discovery of co-expression clusters enriched with essential genes and prioritization of novel biomarker genes that exhibit high correlations between their mRNA and protein expression levels. CAP-RNAseq is applicable to RNA-seq data from different contexts including cancer and available at http://konulabapps.bilkent.edu.tr:3838/CAPRNAseq/ and the docker image is downloadable from https://hub.docker.com/r/konulab/caprnaseq.

Integrating Full-Length and Second-Generation Transcriptomics to Reveal Differentially Expressed Genes Associated with the Development of Corydalis yanhusuo Tuber

Corydalis yanhusuo is a medicinal herb in China that has been widely used to treat various kinds of pain. The tuber is the main organ of C. yanhusuo used for medicinal purposes, but changes in related genes during the development of the tuber have rarely been reported. To identify the differentially expressed genes during tuber development, C. yanhusuo full-length transcriptomic sequencing was performed using single-molecule real-time technology, and tubers at three development stages were selected for comparative transcriptome analysis. A total of 90,496 full-length non-chimeric transcripts were obtained, and 19,341 transcripts were annotated in at least one public database. A total of 9221 differentially expressed genes were identified during the swelling process of C. yanhusuo tuber. A Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis revealed that differentially expressed genes associated with a "starch and sucrose metabolism pathway", "phenylpropanoid biosynthesis pathway", "isoquinoline alkaloid biosynthesis pathway", "zeatin biosynthesis pathway", and "brassinosteroid biosynthesis pathway" were predominantly enriched. In addition, the genes involved in cell wall metabolism were potentially associated with tuber swelling. These processes regulated and were involved in C. yanhusuo tuber development. The results provide a foundation for further research on tuber formation in medicinal plants.

Comparative transcriptome analysis of gene responses of salt-tolerant and salt-sensitive rice cultivars to salt stress

Salt stress is one unfavorable factor of global climate change that adversely affects rice plant growth and yield. To identify novel salt-tolerant genes and new varieties of salt-tolerant rice, a better understanding of the molecular regulation mechanism of salt tolerance in rice is needed. In this study we used transcriptome analyses to examine changes in gene expression of salt-tolerant and salt-sensitive rice plants. The salt-tolerant cultivar HH11 and salt-sensitive cultivar IR29 were treated with 200 mM NaCl solution for 0 h, 6 h, 24 h and 48 h at the three leaf stage. Physiological parameters and transcriptome were measured and analyzed after each treatment. Activity of SOD and POD, as well as the MDA and protein content of the two rice cultivars generally increased with increasing time of exposure to NaCl. Meanwhile, the APX activity first increased, then decreased in both cultivars, with maximum values seen at 6 h for IR29 and at 24 h for HH11. The GR and GPX activity of HH11 were stronger than that of IR29 in response to salt stress. The H2O2 content first increased at 0-6 h, then decreased at 6-24 h, and then increased again at 24-48 h under salt stress. Compared with IR29, SOD, POD and APX activity of HH11 was more sluggish in response to salt stress, reaching the maximum at 24 h or 48 h. The MDA, H2O2 and proline content of HH11 was lower than that of IR29 under salt stress. Relative to untreated HH11 plants (0 h) and those exposed to salt for 6 h, 24 h, and 48 h (H0-H6, H0-H24 and H0-H48), 7462, 6363 and 6636, differentially expressed genes (DEGs), respectively, were identified. For IR29, the respective total DEGs were 7566, 6075 and 6136. GO and KEGG enrichment analysis showed that metabolic pathways related to antioxidative responses and osmotic balance played vital roles in salt stress tolerance. Sucrose and starch metabolism, in addition to flavonoid biosynthesis and glutathione metabolism, showed positive responses to salt stress. Expression of two SPS genes (LOC_Os01g69030 and LOC_Os08g20660) and two GST genes (LOC_Os06g12290 and LOC_Os10g38740) was up-regulated in both HH11 and IR29, whereas expression of LOC_Os09g12660, a glucose-1-phosphate adenylyltransferase gene, and two SS genes (LOC_Os04g17650 and LOC_Os04g24430) was up-regulated differential expression in HH11. The results showed that HH11 had more favorable adjustment in antioxidant and osmotic activity than IR29 upon exposure to salt stress, and highlighted candidate genes that could play roles in the function and regulation mechanism of salt tolerance in rice.

Finding Balance in Adversity: Nitrate Signaling as the Key to Plant Growth, Resilience, and Stress Response

To effectively adapt to changing environments, plants must maintain a delicate balance between growth and resistance or tolerance to various stresses. Nitrate, a significant inorganic nitrogen source in soils, not only acts as an essential nutrient but also functions as a critical signaling molecule that regulates multiple aspects of plant growth and development. In recent years, substantial advancements have been made in understanding nitrate sensing, calcium-dependent nitrate signal transmission, and nitrate-induced transcriptional cascades. Mounting evidence suggests that the primary response to nitrate is influenced by environmental conditions, while nitrate availability plays a pivotal role in stress tolerance responses. Therefore, this review aims to provide an overview of the transcriptional and post-transcriptional regulation of key components in the nitrate signaling pathway, namely, NRT1.1, NLP7, and CIPK23, under abiotic stresses. Additionally, we discuss the specificity of nitrate sensing and signaling as well as the involvement of epigenetic regulators. A comprehensive understanding of the integration between nitrate signaling transduction and abiotic stress responses is crucial for developing future crops with enhanced nitrogen-use efficiency and heightened resilience.

Integrative Analysis of Metabolomic and Transcriptomic Data Reveals the Mechanism of Color Formation in Corms of Pinellia ternata

Pinellia ternata (Thunb.) Breit. (P. ternata) is a very important plant that is commonly used in traditional Chinese medicine. Its corms can be used as medicine and function to alleviate cough, headache, and phlegm. The epidermis of P. ternata corms is often light yellow to yellow in color; however, within the range of P. ternata found in JingZhou City in Hubei Province, China, there is a form of P. ternata in which the epidermis of the corm is red. We found that the total flavonoid content of red P. ternata corms is significantly higher than that of yellow P. ternata corms. The objective of this study was to understand the molecular mechanisms behind the difference in epidermal color between the two forms of P. ternata. The results showed that a high content of anthocyanidin was responsible for the red epidermal color in P. ternata, and 15 metabolites, including cyanidin-3-O-rutinoside-5-O-glucoside, cyanidin-3-O-glucoside, and cyanidin-3-O-rutinoside, were screened as potential color markers in P. ternata through metabolomic analysis. Based on an analysis of the transcriptome, seven genes, including PtCHS1, PtCHS2, PtCHI1, PtDFR5, PtANS, PtUPD-GT2, and PtUPD-GT3, were found to have important effects on the biosynthesis of anthocyanins in the P. ternata corm epidermis. Furthermore, two transcription factors (TFs), bHLH1 and bHLH2, may have regulatory functions in the biosynthesis of anthocyanins in red P. ternata corms. Using an integrative analysis of the metabolomic and transcriptomic data, we identified five genes, PtCHI, PtDFR2, PtUPD-GT1, PtUPD-GT2, and PtUPD-GT3, that may play important roles in the presence of the red epidermis color in P. ternata corms.

Screening and validating of endogenous reference genes in Chlorella sp. TLD 6B under abiotic stress

Chlorella sp. TLD 6B, a microalgae growing in the Taklamakan Desert, Xinjiang of China, is a good model material for studying the physiological and environmental adaptation mechanisms of plants in their arid habitats, as its adaptation to the harsh desert environment has led to its strong resistance. However, when using real-time quantitative polymerase chain reaction (RT-qPCR) to analyze the gene expression of this algae under abiotic stress, it is essential to find the suitable endogenous reference genes so to obtain reliable results. This study assessed the expression stability of 9 endogenous reference genes of Chlorella sp. TLD 6B under four abiotic stresses (drought, salt, cold and heat). These genes were selected based on the analysis results calculated by the three algorithmic procedures of geNorm, NormFinder, and BestKeeper, which were ranked by refinder. Our research showed that 18S and GTP under drought stress, 18S and IDH under salt stress, CYP and 18S under cold stress, GTP and IDH under heat stress were the most stable endogenous reference genes. Moreover, UBC and 18S were the most suitable endogenous reference gene combinations for all samples. In contrast, GAPDH and α-TUB were the two least stable endogenous reference genes in all experimental samples. Additionally, the selected genes have been verified to be durable and reliable by detecting POD and PXG3 genes using above endogenous reference genes. The identification of reliable endogenous reference genes guarantees more accurate RT-qPCR quantification for Chlorella sp. TLD 6B, facilitating functional genomics studies of deserts Chlorella as well as the mining of resistance genes.

Transcriptome Level Reveals the Triterpenoid Saponin Biosynthesis Pathway of Bupleurum falcatum L

Bupleurum falcatum L. is frequently used in traditional herbal medicine in Asia. Saikosaponins (SSs) are the main bioactive ingredients of B. falcatum, but the biosynthetic pathway of SSs is unclear, and the biosynthesis of species-specific phytometabolites is little known. Here we resolved the transcriptome profiles of B. falcatum to identify candidate genes that might be involved in the biosynthesis of SSs. By isoform sequencing (Iso-Seq) analyses of the whole plant, a total of 26.98 Gb of nucleotides were obtained and 124,188 unigenes were identified, and 81,594 unigenes were successfully annotated. A total of 1033 unigenes of 20 families related to the mevalonate (MVA) pathway and methylerythritol phosphate (MEP) pathway of the SS biosynthetic pathway were identified. The WGCNA (weighted gene co-expression network analysis) of these unigenes revealed that only the co-expression module of MEmagenta, which contained 343 unigenes, was highly correlated with the biosynthesis of SSs. Comparing differentially expressed gene analysis and the WGCNA indicated that 130 out of 343 genes of the MEmagenta module exhibited differential expression levels, and genes with the most "hubness" within this module were predicted. Manipulation of these genes might improve the biosynthesis of SSs.

H2S Enhanced the Tolerance of Malus hupehensis to Alkaline Salt Stress through the Expression of Genes Related to Sulfur-Containing Compounds and the Cell Wall in Roots

Malus is an economically important plant that is widely cultivated worldwide, but it often encounters saline-alkali stress. The composition of saline-alkali land is a variety of salt and alkali mixed with the formation of alkaline salt. Hydrogen sulfide (H₂S) has been reported to have positive effects on plant responses to abiotic stresses. Our previous study showed that H₂S pretreatment alleviated the damage caused by alkaline salt stress to Malus hupehensis Rehd. var. pingyiensis Jiang (Pingyi Tiancha, PYTC) roots by regulating Na⁺/K⁺ homeostasis and oxidative stress. In this study, transcriptome analysis was used to investigate the overall mechanism through which H₂S alleviates alkaline salt stress in PYTC roots. Simultaneously, differentially expressed genes (DEGs) were explored. Transcriptional profiling of the Control-H₂S, Control-AS, Control-H₂S + AS, and AS-H₂S + AS comparison groups identified 1618, 18,652, 16,575, and 4314 DEGs, respectively. Further analysis revealed that H₂S could alleviate alkaline salt stress by increasing the energy maintenance capacity and cell wall integrity of M. hupehensis roots and by enhancing the capacity for reactive oxygen species (ROS) metabolism because more upregulated genes involved in ROS metabolism and sulfur-containing compounds were identified in M. hupehensis roots after H₂S pretreatment. qRT-PCR analysis of H₂S-induced and alkaline salt-response genes showed that these genes were consistent with the RNA-seq analysis results, which indicated that H₂S alleviation of alkaline salt stress involves the genes of the cell wall and sulfur-containing compounds in PYTC roots.

Phenotypic and transcriptomic responses of cultivated sunflower seedlings (Helianthus annuus L.) to four abiotic stresses

Plants encounter and respond to numerous abiotic stresses during their lifetimes. These stresses are often related and could therefore elicit related responses. There are, however, relatively few detailed comparisons between multiple different stresses at the molecular level. Here, we investigated the phenotypic and transcriptomic response of cultivated sunflower (Helianthus annuus L.) seedlings to three water-related stresses (i.e., dry-down, an osmotic challenge, and salt stress), as well as a generalized low-nutrient stress. All four stresses negatively impacted seedling growth, with the nutrient stress having a more divergent response from control as compared to the water-related stresses. Phenotypic responses were consistent with expectations for growth in low-resource environments, including increased (i.e., less negative) carbon fractionation values and leaf C:N ratios, as well as increased belowground biomass allocation. The number of differentially expressed genes (DEGs) under stress was greater in leaf tissue, but roots exhibited a higher proportion of DEGs unique to individual stresses. Overall, the three water-related stresses had a more similar transcriptomic response to each other vs. nutrient stress, though this pattern was more pronounced in root vs. leaf tissue. In contrast to our DEG analyses, co-expression network analysis revealed that there was little indication of a shared response between the four stresses in despite the majority of DEGs being shared between multiple stresses. Importantly, osmotic stress, which is often used to simulate drought stress in experimental settings, had little transcriptomic resemblance to true water limitation (i.e., dry-down) in our study, calling into question its utility as a means for simulating drought.

Comparative Analysis on the Evolution of Flowering Genes in Sugar Pathway in Brassicaceae

Sugar plays an important role in regulating the flowering of plants. However, studies of genes related to flowering regulation by the sugar pathway of Brassicaceae plants are scarce. In this study, we performed a comprehensive comparative genomics analysis of the flowering genes in the sugar pathway from seven members of the Brassicaceae, including: Arabidopsis thaliana, Arabidopsis lyrata, Astelia pumila, Camelina sativa, Brassica napus, Brassica oleracea, and Brassica rapa. We identified 105 flowering genes in the sugar pathway of these plants, and they were categorized into nine groups. Protein domain analysis demonstrated that the IDD8 showed striking structural variations in different Brassicaceae species. Selection pressure analysis revealed that sugar pathway genes related to flowering were subjected to strong purifying selection. Collinearity analysis showed that the identified flowering genes expanded to varying degrees, but SUS4 was absent from the genomes of Astelia pumila, Camelina sativa, Brassica napus, Brassica oleracea, and Brassica rapa. Tissue-specific expression of ApADG indicated functional differentiation. To sum up, genome-wide identification revealed the expansion, contraction, and diversity of flowering genes in the sugar pathway during Brassicaceae evolution. This study lays a foundation for further study on the evolutionary characteristics and potential biological functions of flowering genes in the sugar pathway of Brassicaceae.

Exploring epitranscriptomics for crop improvement and environmental stress tolerance

Climate change and stressful environmental conditions severely hamper crop growth, development and yield. Plants respond to environmental perturbations, through their plasticity provided by key-genes, governed at post-/transcriptional levels. Gene-regulation in plants is a multilevel process controlled by diverse cellular entities that includes transcription factors (TF), epigenetic regulators and non-coding RNAs beside others. There are successful studies confirming the role of epigenetic modifications (DNA-methylation/histone-modifications) in gene expression. Recent years have witnessed emergence of a highly specialized field the "Epitranscriptomics". Epitranscriptomics deals with investigating post-transcriptional RNA chemical-modifications present across the life forms that change structural, functional and biological characters of RNA. However, deeper insights on of epitranscriptomic modifications, with >140 types known so far, are to be understood fully. Researchers have identified epitranscriptome marks (writers, erasers and readers) and mapped the site-specific RNA modifications (m6A, m⁵C, 3' uridylation, etc.) responsible for fine-tuning gene expression in plants. Simultaneous advancement in sequencing platforms, upgraded bioinformatic tools and pipelines along with conventional labelled techniques have further given a statistical picture of these epitranscriptomic modifications leading to their potential applicability in crop improvement and developing climate-smart crops. We present herein the insights on epitranscriptomic machinery in plants and how epitranscriptome and epitranscriptomic modifications underlying plant growth, development and environmental stress responses/adaptations. Third-generation sequencing technology, advanced bioinformatics tools and databases being used in plant epitranscriptomics are also discussed. Emphasis is given on potential exploration of epitranscriptome engineering for crop-improvement and developing environmental stress tolerant plants covering current status, challenges and future directions.

Full-Length Transcriptomics Reveals Complex Molecular Mechanism of Salt Tolerance in Bromus inermis L

Bromus inermis L. (commonly known as smooth bromegrass) is a grass species with high nutritional value, great palatability, cold tolerance, and grazing resistance, which has been widely cultivated for pasture and sand fixation in northern and northwestern China. Salt stress is a main environmental factor limiting growth and production of smooth bromegrass. In this study, we performed PacBio Iso-Seq to construct the first full-length transcriptome database for smooth bromegrass under 300 mM NaCl treatment at different time points. Third-generation full-length transcriptome sequencing yielded 19.67 G polymerase read bases, which were assembled into 355,836 full-length transcripts with an average length of 2,542 bp. A total of 116,578 differentially expressed genes were obtained by comparing the results of third-generation sequencing and second-generation sequencing. GO and KEGG enrichment analyses revealed that multiple pathways were differently activated in leaves and roots. In particular, a number of genes participating in the molecular network of plant signal perception, signal transduction, transcription regulation, antioxidant defense, and ion regulation were affected by NaCl treatment. In particular, the CBL-CIPK, MAPK, ABA signaling network, and SOS core regulatory pathways of Ca²⁺ signal transduction were activated to regulate salt stress response. In addition, the expression patterns of 10 salt-responsive genes were validated by quantitative real-time PCR, which were consistent with those detected by RNA-Seq. Our results reveal the molecular regulation of smooth bromegrass in response to salt stress, which are important for further investigation of critical salt responsive genes and molecular breeding of salt-tolerant smooth bromegrass.

Full-length transcriptome sequencing reveals the molecular mechanism of potato seedlings responding to low-temperature

BACKGROUND

Potato (Solanum tuberosum L.) is one of the world's most important crops, the cultivated potato is frost-sensitive, and low-temperature severely influences potato production. However, the mechanism by which potato responds to low-temperature stress is unclear. In this research, we apply a combination of second-generation sequencing and third-generation sequencing technologies to sequence full-length transcriptomes in low-temperature-sensitive cultivars to identify the important genes and main pathways related to low-temperature resistance.

RESULTS

In this study, we obtained 41,016 high-quality transcripts, which included 15,189 putative new transcripts. Amongst them, we identified 11,665 open reading frames, 6085 simple sequence repeats out of the potato dataset. We used public available genomic contigs to analyze the gene features, simple sequence repeat, and alternative splicing event of 24,658 non-redundant transcript sequences, predicted the coding sequence and identified the alternative polyadenylation. We performed cluster analysis, GO, and KEGG functional analysis of 4518 genes that were differentially expressed between the different low-temperature treatments. We examined 36 transcription factor families and identified 542 transcription factors in the differentially expressed genes, and 64 transcription factors were found in the AP2 transcription factor family which was the most. We measured the malondialdehyde, soluble sugar, and proline contents and the expression genes changed associated with low temperature resistance in the low-temperature treated leaves. We also tentatively speculate that StLPIN10369.5 and StCDPK16 may play a central coordinating role in the response of potatoes to low temperature stress.

CONCLUSIONS

Overall, this study provided the first large-scale full-length transcriptome sequencing of potato and will facilitate structure-function genetic and comparative genomics studies of this important crop.

Long-read transcriptome sequencing provides insight into lignan biosynthesis during fruit development in Schisandra chinensis

BACKGROUND

Schisandra chinensis, an ancient member of the most basal angiosperm lineage which is known as the ANITA, is a fruit-bearing vine with the pharmacological effects of a multidrug system, such as antioxidant, anti-inflammatory, cardioprotective, neuroprotective, anti-osteoporosis effects. Its major bioactive compound is represented by lignans such as schisandrin. Molecular characterization of lignan biosynthesis in S. chinensis is of great importance for improving the production of this class of active compound. However, the biosynthetic mechanism of schisandrin remains largely unknown.

RESULTS

To understand the potential key catalytic steps and their regulation of schisandrin biosynthesis, we generated genome-wide transcriptome data from three different tissues of S. chinensis cultivar Cheongsoon, including leaf, root, and fruit, via long- and short-read sequencing technologies. A total of 132,856 assembled transcripts were generated with an average length of 1.9 kb and high assembly completeness. Overall, our data presented effective, accurate gene annotation in the prediction of functional pathways. In particular, the annotation revealed the abundance of transcripts related to phenylpropanoid biosynthesis. Remarkably, transcriptome profiling during fruit development of S. chinensis cultivar Cheongsoon revealed that the phenylpropanoid biosynthetic pathway, specific to coniferyl alcohol biosynthesis, showed a tendency to be upregulated at the postfruit development stage. Further the analysis also revealed that the pathway forms a transcriptional network with fruit ripening-related genes, especially the ABA signaling-related pathway. Finally, candidate unigenes homologous to isoeugenol synthase 1 (IGS1) and dirigent-like protein (DIR), which are subsequently activated by phenylpropanoid biosynthesis and thus catalyze key upstream steps in schisandrin biosynthesis, were identified. Their expression was increased at the postfruit development stage, suggesting that they may be involved in the regulation of schisandrin biosynthesis in S. chinensis.

CONCLUSIONS

Our results provide new insights into the production and accumulation of schisandrin in S. chinensis berries and will be utilized as a valuable transcriptomic resource for improving the schisandrin content.

Plant Dehydrins: Expression, Regulatory Networks, and Protective Roles in Plants Challenged by Abiotic Stress

Dehydrins, also known as Group II late embryogenesis abundant (LEA) proteins, are classic intrinsically disordered proteins, which have high hydrophilicity. A wide range of hostile environmental conditions including low temperature, drought, and high salinity stimulate dehydrin expression. Numerous studies have furnished evidence for the protective role played by dehydrins in plants exposed to abiotic stress. Furthermore, dehydrins play important roles in seed maturation and plant stress tolerance. Hence, dehydrins might also protect plasma membranes and proteins and stabilize DNA conformations. In the present review, we discuss the regulatory networks of dehydrin gene expression including the abscisic acid (ABA), mitogen-activated protein (MAP) kinase cascade, and Ca²⁺ signaling pathways. Crosstalk among these molecules and pathways may form a complex, diverse regulatory network, which may be implicated in regulating the same dehydrin.

Network response of two cherry tomato (Lycopersicon esculentum) cultivars to Cadmium stress as revealed by transcriptome analysis

Soil cadmium (Cd) contamination severely threatens human health. Therefore, screening and breeding low-Cd absorption cultivars of cherry tomato (Solanum lycopersicum L.) is essential to restrict human Cd intake. In this study, a hydroponic experiment was conducted to perform a comparative transcriptome analysis of the leaves of two cherry tomato cultivars with different Cd contents under different Cd stress (0, 10, and 50 μM), for the purpose of exploring the differences in the transcriptional responses to Cd stress between the two cultivars. Our results revealed that the Cd content in the leaves of HLZ (Hanluzhe; a low-Cd accumulation cultivar) was significantly lower than that in the leaves of LFC (Lvfeicui; a high-Cd accumulation cultivar). Transcriptome analysis showed that the different expression genes (DEGs) were mainly involved in plant hormone signal transduction, antioxidant enzymes, cell wall biosynthesis, and metal transportation. In the LFC leaves, DEGs in the IAA signal transduction and antioxidant enzymes exhibited higher transcription levels. However, the DEGs in the ETH signal transduction demonstrated a lower transcription level compared to that of HLZ. Over-expressed genes in the pectin biosynthesis and pectin methylesterase (PME) of the LFC leaves might result in the trapping of Cd by increased levels of low-methylated pectin around the cell wall. Furthermore, Cd transporter genes, such as HMA5, NRAMP6, CAX3, ABCC3, and PDR1, were up-regulated in the HLZ leaves, indicating that the HLZ cultivar comprised an active Cd transport capacity from apoplast to vacuolar. This may contribute to the low Cd concentration observed in the HLZ leaves. Overall, our study provides a molecular basis for tomato screening and breeding.

Transcriptome sequencing of wild soybean revealed gene expression dynamics under low nitrogen stress

Nitrogen is one of the essential elements for plant growth. Wild soybeans (Glycine soja) have strong abilities to survive in harsh and barren environments, and hence become ideal plant model for studying plant adaptability to low nitrogen (LN) conditions. In this study, we analyzed and compared the transcriptomes of wild soybean subjected to LN treatments. We totally identified 1095 (681 up and 414 down) and 5490 (2998 up and 2492 down) differentially expressed genes (DEGs) in the aerial parts (leaf and stem, LS) and roots, respectively. Gene ontology classification analysis revealed that the categories related to LN stress (including oxidation reduction, transcriptional regulation, membrane, and protein phosphorylation) were highly enriched among DEGs. In addition, a total of 784 transcription factor (TF) and 84 transporter protein (TP) genes were determined in LS DEGs, of which some TF genes (NAC1, NAC35, ZFP1, CIM1, and WRKY25) and TP genes like NRT2.5 (nitrate transporter) and ABCC12 (ABC transporter) were widely upregulated under LN stress. Nevertheless, a total of 3859 TF and 370 TP genes were identified in root DEGs, of which some TF genes (NAC6, NAC14, MYB29, MYB92, bZIP62, bZIP72, WRKY60, WRKY58) and TP genes like NRT2.4 and HAK5 (potassium transporter) were upregulated under LN stress. These findings suggest that the identified DEGs may play vital roles in plant responses to LN stress, providing important genetic resources for further functional dissection of plant molecular mechanisms to LN stress.

Third-generation sequencing and metabolome analysis reveal candidate genes and metabolites with altered levels in albino jackfruit seedlings

BACKGROUND

Most plants rely on photosynthesis; therefore, albinism in plants with leaves that are white instead of green causes slow growth, dwarfing, and even death. Although albinism has been characterized in annual model plants, little is known about albino trees. Jackfruit (Artocarpus heterophyllus) is an important tropical fruit tree species. To gain insight into the mechanisms underlying the differential growth and development between albino jackfruit mutants and green seedlings, we analyzed root, stem, and leaf tissues by combining PacBio single-molecule real-time (SMRT) sequencing, high-throughput RNA-sequencing (RNA-seq), and metabolomic analysis.

RESULTS

We identified 8,202 differentially expressed genes (DEGs), including 225 genes encoding transcription factors (TFs), from 82,572 full-length transcripts. We also identified 298 significantly changed metabolites (SCMs) in albino A. heterophyllus seedlings from a set of 692 metabolites in A. heterophyllus seedlings. Pathway analysis revealed that these DEGs were highly enriched in metabolic pathways such as 'photosynthesis', 'carbon fixation in photosynthetic organisms', 'glycolysis/gluconeogenesis', and 'TCA cycle'. Analysis of the metabolites revealed 76 SCMs associated with metabolic pathways in the albino mutants, including L-aspartic acid, citric acid, succinic acid, and fumaric acid. We selected 225 differentially expressed TF genes, 333 differentially expressed metabolic pathway genes, and 76 SCMs to construct two correlation networks. Analysis of the TF-DEG network suggested that basic helix-loop-helix (bHLH) and MYB-related TFs regulate the expression of genes involved in carbon fixation and energy metabolism to affect light responses or photomorphogenesis and normal growth. Further analysis of the DEG-SCM correlation network and the photosynthetic carbon fixation pathway suggested that NAD-ME2 (encoding a malic enzyme) and L-aspartic acid jointly inhibit carbon fixation in the albino mutants, resulting in reduced photosynthetic efficiency and inhibited plant growth.

CONCLUSIONS

Our preliminarily screening identified candidate genes and metabolites specifically affected in albino A. heterophyllus seedlings, laying the foundation for further study of the regulatory mechanism of carbon fixation during photosynthesis and energy metabolism. In addition, our findings elucidate the way genes and metabolites respond in albino trees.

Full-length transcriptome analysis provides new insights into the early bolting occurrence in medicinal Angelica sinensis

Angelica sinensis (Oliv.) Diels root part is an integral component of traditional Chinese medicine, widely prescribed to improve blood circulation and blood stasis. However, early bolting of A. sinensis compromises the quality of the roots and hence is a major limitation for yield of medicinal materials. To date, little information about the molecular mechanisms underlying bolting is available for this important medicinal plant. To identify genes putatively involved in early bolting, we have conducted the transcriptome analysis of the shoot tips of the early-bolting plants and non-bolting (normal) plants of A. sinensis, respectively, using a combination of third-generation sequencing and next-generation sequencing. A total of 43,438 non-redundant transcripts were collected and 475 unique differentially expressed genes (DEGs) were identified. Gene annotation and functional analyses revealed that DEGs were highly involved in plant hormone signaling and biosynthesis pathways, three main flowering pathways, pollen formation, and very-long-chain fatty acids biosynthesis pathways. The levels of endogenous hormones were also changed significantly in the early bolting stage of A. sinensis. This study provided new insights into the transcriptomic control of early bolting in A. sinensis, which could be further applied to enhance the yield of medicinally important raw materials.

Full-length SMRT transcriptome sequencing and microsatellite characterization in Paulownia catalpifolia

Paulownia catalpifolia is an important, fast-growing timber species known for its high density, color and texture. However, few transcriptomic and genetic studies have been conducted in P. catalpifolia. In this study, single-molecule real-time sequencing technology was applied to obtain the full-length transcriptome of P. catalpifolia leaves treated with varying degrees of drought stress. The sequencing data were then used to search for microsatellites, or simple sequence repeats (SSRs). A total of 28.83 Gb data were generated, 25,969 high-quality (HQ) transcripts with an average length of 1624 bp were acquired after removing the redundant reads, and 25,602 HQ transcripts (98.59%) were annotated using public databases. Among the HQ transcripts, 16,722 intact coding sequences, 149 long non-coding RNAs and 179 alternative splicing events were predicted, respectively. A total of 7367 SSR loci were distributed throughout 6293 HQ transcripts, of which 763 complex SSRs and 6604 complete SSRs. The SSR appearance frequency was 28.37%, and the average distribution distance was 5.59 kb. Among the 6604 complete SSR loci, 1-3 nucleotide repeats were dominant, occupying 97.85% of the total SSR loci, of which mono-, di- and tri-nucleotide repeats were 44.68%, 33.86% and 19.31%, respectively. We detected 112 repeat motifs, of which A/T (42.64%), AG/CT (12.22%), GA/TC (9.63%), GAA/TTC (1.57%) and CCA/TGG (1.54%) were most common in mono-, di- and tri-nucleotide repeats, respectively. The length of the repeat SSR motifs was 10-88 bp, and 4997 (75.67%) were ≤ 20 bp. This study provides a novel full-length transcriptome reference for P. catalpifolia and will facilitate the identification of germplasm resources and breeding of new drought-resistant P. catalpifolia varieties.

Comparative physiological and full-length transcriptome analyses reveal the molecular mechanism of melatonin-mediated salt tolerance in okra (Abelmoschus esculentus L.)

BACKGROUND

Melatonin, a multifunctional signal molecule, has been reported to play crucial roles in growth and development and stress responses in various plant species. Okra (Abelmoschus esculentus L.) is a food crop with extremely high values of nutrition and healthcare. Recent reports have revealed the protective role of melatonin in alleviating salt stress. However, little is known about its regulatory mechanisms in response to salt stress in okra.

RESULTS

In this study, we explored whether exogenous melatonin pretreatment could alleviate salt stress (300 mM NaCl) of okra plants. Results showed that exogenous application of melatonin (50 μM) significantly enhanced plant tolerance to salt stress, as demonstrated by the plant resistant phenotype, as well as by the higher levels of the net photosynthetic rate, chlorophyll fluorescence and chlorophyll content in comparison with nontreated salt-stressed plants. Additionally, melatonin pretreatment remarkably decreased the levels of lipid peroxidation and H₂O₂ content and scavenged O₂^•- in melatonin-pretreated plants, which may be attributed to the higher levels of enzyme activities including POD and GR. Moreover, a combination of third- (PacBio) and second-generation (Illumina) sequencing technologies was applied to sequence full-length transcriptomes of okra. A total of 121,360 unigenes was obtained, and the size of transcript lengths ranged from 500 to 6000 bp. Illumina RNA-seq analysis showed that: Comparing with control, 1776, 1063 and 1074 differential expression genes (DEGs) were identified from the three treatments (NaCl, MT50 and MT + NaCl, respectively). These genes were enriched in more than 10 GO terms and 34 KEGG pathways. Nitrogen metabolism, sulfur metabolism, and alanine, aspartate and glutamate metabolism were significantly enriched in all three treatments. Many transcription factors including MYB, WRKY, NAC etc., were also identified as DEGs.

CONCLUSIONS

Our preliminary results suggested that melatonin pretreatment enhanced salt tolerance of okra plants for the first time. These data provide the first set of full-length isoforms in okra and more comprehensive insights into the molecular mechanism of melatonin responses to salt stress.

De novo full length transcriptome analysis of Arachis glabrata provides insights into gene expression dynamics in response to biotic and abiotic stresses

The perennial ornamental peanut Arachis glabrata represents one of the most adaptable wild Arachis species. This study used PacBio combined with BGISEQ-500 RNA-seq technology to study the transcriptome and gene expression dynamics of A. glabrata. Of the total 109,747 unique transcripts obtained, >90,566 transcripts showed significant homology to known proteins and contained the complete coding sequence (CDS). RNA-seq revealed that 1229, 1039, 1671, 3923, 1521 and 1799 transcripts expressed specifically in the root, stem, leaf, flower, peg and pod, respectively. We also identified thousands of differentially expressed transcripts in response to drought, salt, cold and leaf spot disease. Furthermore, we identified 30 polyphenol oxidase encoding genes associated with the quality of forage, making A. glabrata suitable as a forage crop. Our findings presented the first transcriptome study of A. glabrata which will facilitate genetic and genomics studies and lays the groundwork for a deeper understanding of the A. glabrata genome.

Plastidial transporters KEA1 and KEA2 at the inner envelope membrane adjust stromal pH in the dark

Photosynthesis and carbon fixation depend critically on the regulation of pH in chloroplast compartments in the daylight and at night. While it is established that an alkaline stroma is required for carbon fixation, it is not known how alkaline stromal pH is formed, maintained or regulated. We tested whether two envelope transporters, AtKEA1 and AtKEA2, directly affected stromal pH in isolated Arabidopsis chloroplasts using the fluorescent probe 2',7'-bis-(2-carboxyethyl)-5-(and-6)-carboxyfluorescein (BCECF). External K⁺ -induced alkalinization of the stroma was observed in chloroplasts from wild-type (WT) plants but not from kea1kea2 mutants, suggesting that KEA1 and KEA2 mediate K⁺ uptake/H⁺ loss to modulate stromal pH. While light-stimulated alkalinization of the stroma was independent of KEA1 and KEA2, the rate of decay to neutral pH in the dark is delayed in kea1kea2 mutants. However, the dark-induced loss of a pH gradient across the thylakoid membrane was similar in WT and mutant chloroplasts. This indicates that proton influx from the cytosol mediated by envelope K⁺ /H⁺ antiporters contributes to adjustment of stromal pH upon light to dark transitions.

A Tissue-Specific Landscape of Alternative Polyadenylation, lncRNAs, TFs, and Gene Co-expression Networks in Liriodendron chinense

Liriodendron chinense is an economically and ecologically important deciduous tree species. Although the reference genome has been revealed, alternative polyadenylation (APA), transcription factors (TFs), long non-coding RNAs (lncRNAs), and co-expression networks of tissue-specific genes remain incompletely annotated. In this study, we used the bracts, petals, sepals, stamens, pistils, leaves, and shoot apex of L. chinense as materials for hybrid sequencing. On the one hand, we improved the annotation of the genome. We detected 13,139 novel genes, 7,527 lncRNAs, 1,791 TFs, and 6,721 genes with APA sites. On the other hand, we found that tissue-specific genes play a significant role in maintaining tissue characteristics. In total, 2,040 tissue-specific genes were identified, among which 9.2% of tissue-specific genes were affected by APA, and 1,809 tissue-specific genes were represented in seven specific co-expression modules. We also found that bract-specific hub genes were associated plant defense, leaf-specific hub genes were involved in energy metabolism. Moreover, we also found that a stamen-specific hub TF Lchi25777 may be involved in the determination of stamen identity, and a shoot-apex-specific hub TF Lchi05072 may participate in maintaining meristem characteristic. Our study provides a landscape of APA, lncRNAs, TFs, and tissue-specific gene co-expression networks in L. chinense that will improve genome annotation, strengthen our understanding of transcriptome complexity, and drive further research into the regulatory mechanisms of tissue-specific genes.

RNA-seq reveals the salt tolerance of Ipomoea pes-caprae, a wild relative of sweet potato

Wild relatives of crops are often rich in genetic resources and provide great possibilities for crop improvement. Ipomoea pes-caprae is one of the wild relatives of sweet potato and has high salt tolerance. Transcriptomes in the treatment and control groups at various times were sequenced to identify salt tolerance genes and salt response pathways. A total of 40,525 genes were obtained, of which 2478 and 3334 were differentially expressed in the roots and leaves of I. pes-caprae under salt stress, respectively. Identification of candidate genes revealed that the mitogen-activated protein kinase (MAPK) signaling pathway of plants and plant hormone signal transduction participates in the salt signal of I. pes-caprae under salt stress. Homology to ABI2 (HAB2) and Clade A protein phosphatases type 2C (HAI1), which encode two protein phosphatases 2C (PP2C) in the abscisic acid (ABA) signal pathway, were continuously up-regulated upon salt stress, indicating their key role in the salt signal transduction pathway of I. pes-caprae. The expression of EIN3-binding F-box protein 1 (EBF1) in the ethylene signaling pathway was also up-regulated, revealing that the salt tolerance of I. pes-caprae was related to the scavenging of reactive oxygen species (ROS). This study provides insights into the mechanism of salt-tolerant plants and the mining of salt-tolerant genes in sweet potato for the innovation of germplasm resources.

Differential Expression of Genes Involved in Saikosaponin Biosynthesis Between Bupleurum chinense DC. and Bupleurum scorzonerifolium Willd

Radix Bupleuri (roots of Bupleurum spp.) is an important medicinal herb. Triterpenoid saponins of saikosaponins generally constitute the main class of secondary metabolites of plants in the Bupleurum genus. However, the molecular regulatory mechanism underlying their biosynthesis remains elusive. In this study, we observed significantly different saikosaponin biosynthesis between Bupleurum chinense and Bupleurum scorzonerifolium at the seedling stage. The sequential and expression characterization of 232 genes in the triterpenoid saponin biosynthetic pathway, which includes the mevalonate (MVA) pathway and methylerythritol phosphate (MEP) pathway, between B. chinense and B. scorzonerifolium was also investigated. Sixty of these genes may be involved in saikosaponin biosynthesis. Manipulation of these genes, especially those of the β-AS, P450, and UGT families, may improve saikosaponin production.

The Role of N6-Methyladenosine Methylation in the Progression of Endometrial Cancer

Purpose: N6-methyladenosine (m6A) methylation was the most abundant internal modification on messenger RNAs in eukaryotes. This study intended to explore the role of m6A methylation in endometrial cancer (EC). Materials and Methods: The m6A-sequencing data "GSE93911" of human EC were downloaded from Gene Expression Omnibus database. Hisat2 software and MACS2 were used to perform the alignment of reads and m6A methylation peak calling, and the peaks were annotated using Chipseeker. Then, differential m6A methylation peaks between normal and tumor samples were analyzed, followed by the functional enrichment analysis of the differentially methylated genes in promoter and 3' untranslated region (UTR) using Clusterprofiler. Based on the 450K methylated chip data, gene expression and clinical data in The Cancer Genome Atlas, the differentially methylated genes were verified, followed by Cox univariate/multivariate regression analysis and survival analysis. Finally, a risk prognosis model was constructed. Results: The m6A peak number was decreased in EC. The distribution of m6A peaks was highly enriched near transcriptional start site, in promoter, UTR, intron and exon, followed by distal intergenic. A total of 581 differentially methylated genes (361 hyper- and 220 hypomethylated genes) were identified in promoter and UTR regions that were enriched in insulin resistance (IR) and extracellular matrix (ECM). A total of 181 genes with significant differential expressions and differential methylation site in EC were selected. Of which, 31 genes were correlated with survival, and an 11-gene risk prognosis model was identified, including GDF7, BNC2, SLC8A1, B4GALNT3, DHCR24, ESRP1, HOXB9, IGSF9, KIAA1324, MSnX1, and PHGDH. Conclusion: The m6A methylation regulated EC progression by targeting the genes related to IR and ECM. A 11-gene risk prognosis model was identified to predict survival of patients with EC.

Single-Molecule Real-Time Transcript Sequencing of Turnips Unveiling the Complexity of the Turnip Transcriptome

To generate the full-length transcriptome of Xinjiang green and purple turnips, Brassica rapa var. Rapa, using single-molecule real-time (SMRT) sequencing. The samples of two varieties of Brassica rapa var. Rapa at five developmental stages were collected and combined to perform SMRT sequencing. Meanwhile, next generation sequencing was performed to correct SMRT sequencing data. A series of analyses were performed to investigate the transcript structure. Finally, the obtained transcripts were mapped to the genome of Brassica rapa ssp. pekinesis Chiifu to identify potential novel transcripts. For green turnip (F01), a total of 19.54 Gb clean data were obtained from 8 cells. The number of reads of insert (ROI) and full-length non-chimeric (FLNC) reads were 510,137 and 267,666. In addition, 82,640 consensus isoforms were obtained in the isoform sequences clustering, of which 69,480 were high-quality, and 13,160 low-quality sequences were corrected using Illumina RNA seq data. For purple turnip (F02), there were 20.41 Gb clean data, 552,829 ROIs, and 274,915 FLNC sequences. A total of 93,775 consensus isoforms were obtained, of which 78,798 were high-quality, and the 14,977 low-quality sequences were corrected. Following the removal of redundant sequences, there were 46,516 and 49,429 non-redundant transcripts for F01 and F02, respectively; 7,774 and 9,385 alternative splicing events were predicted for F01 and F02; 63,890 simple sequence repeats, 59,460 complete coding sequences, and 535 long-non coding RNAs were predicted. Moreover, 5,194 and 5,369 novel transcripts were identified by mapping to Brassica rapa ssp. pekinesis Chiifu. The obtained transcriptome data may improve turnip genome annotation and facilitate further study of the Brassica rapa var. Rapa genome and transcriptome.

Characterization and Analysis of the Full-Length Transcriptomes of Multiple Organs in Pseudotaxus chienii (W.C.Cheng) W.C.Cheng

Pseudotaxus chienii, a rare tertiary relict species with economic and ecological value, is a representative of the monotypic genus Pseudotaxus that is endemic to China. P. chienii can adapt well to habitat isolation and ecological heterogeneity under a variety of climate and soil conditions, and is able to survive in harsh environments. However, little is known about the molecular and genetic resources of this long-lived conifer. Herein, we sequenced the transcriptomes of four organs of P. chienii using the PacBio Isoform Sequencing and Illumina RNA Sequencing platforms. Based on the PacBio Iso-Seq data, we obtained 44,896, 58,082, 50,485, and 67,638 full-length unigenes from the root, stem, leaf, and strobilus, respectively, with a mean length of 2692 bp, and a mean N50 length of 3010.75 bp. We then comprehensively annotated these unigenes. The number of organ-specific expressed unigenes ranged from 4393 in leaf to 9124 in strobilus, suggesting their special roles in physiological processes, organ development, and adaptability in the different four organs. A total of 16,562 differentially expressed genes (DEGs) were identified among the four organs and clustered into six subclusters. The gene families related to biotic/abiotic factors, including the TPS, CYP450, and HSP families, were characterized. The expression levels of most DEGs in the phenylpropanoid biosynthesis pathway and plant–pathogen interactions were higher in the root than in the three other organs, suggesting that root constitutes the main organ of defensive compound synthesis and accumulation and has a stronger ability to respond to stress. The sequences were analyzed to predict transcription factors, long non-coding RNAs, and alternative splicing events. The expression levels of most DEGs of C2H2, C3H, bHLH, and bZIP families in the root and stem were higher than those in the leaf and strobilus, indicating that these TFs may play a crucial role in the survival of the root and stem. These results comprise the first comprehensive gene expression profiles obtained for different organs of P. chienii. Our findings will facilitate further studies on the functional genomics, adaptive evolution, and phylogeny of P. chienii, and lay the foundation for the development of conservation strategies for this endangered conifer.

Transcriptome profiling of spike provides expression features of genes related to terpene biosynthesis in lavender

Lavender (Lavandula angustifolia) is an important economic plant because of the value of its essential oil (EO). The Yili Valley in Xinjiang has become the largest lavender planting base in China. However, there is a lack of research on the gene expression regulation of EO biosynthesis and metabolism in local varieties. Here, de novo transcriptome analysis of inflorescence of three development stages from initial flower bud to flowering stage 50% from two lavender cultivars with contrasting EO production revealed the dynamics of 100,177 differentially expressed transcripts (DETs) in various stages of spike development within and across the cultivars. The lavender transcriptome contained 77 DETs with annotations related to terpenoid biosynthesis. The expression profiles of the 27 genes involved in the methylerythritol phosphate (MEP) pathway, 22 genes in the mevalonate (MVA) pathway, 28 genes related to monoterpene and sesquiterpene biosynthesis during inflorescence development were comprehensively characterized, and possible links between the expression changes of genes and contents of EO constituents were explored. The upregulated genes were mainly concentrated in the MEP pathway, while most genes in the MVA pathway were downregulated during flower development, and cultivars with a higher EO content presented higher expression of genes in the MEP pathway, indicating that EOs were chiefly produced through the MEP pathway. Additionally, MYB transcription factors constituted the largest number of transcripts in all samples, suggesting their potential roles in regulating EO biosynthesis. The sequences and transcriptional patterns of the transcripts will be helpful for understanding the molecular basis of lavender terpene biosynthesis.

Identification of drought resistant miRNA in Macleaya cordata by high-throughput sequencing

Drought is one of the most serious factors affecting crop yields in the world. Macleaya cordata (Willd.) is a draught-tolerant medicinal plant that has been proposed as a pioneer crop to be cultivated in arid areas. However, the exact molecular mechanisms through which M. cordata responds to draught stress remain elusive. In recent years, microRNA (miRNAs) in plants have been associated with stress response. Based on these findings, the current study aimed to shed light on the potential regulatory roles of miRNAs in the draught tolerance of M. cordata by employing high-throughput RNA sequencing and degradation sequencing. Six M. cordata plants were randomly divided into two equal experiment groups, including one draught group and one control group. High-throughput sequencing of the M. cordata samples led to the identification of 895 miRNAs, of which 18 showed significantly different expression levels between the two groups. PsRobot analysis and degradation sequencing predicted the differential miRNAs to target 59 and 36 genes, respectively. Functional analysis showed that 38 of the predicted genes could be implicated in the modulation of stress response. Four miRNAs and eight target genes were selected for quantitative real-time polymerase chain reaction (qRT-PCR) validation. The expression trend of each miRNA analyzed by qRT-PCR was consistent with that determined by sequencing, and was negatively correlated with those of its target genes. The results of our current study supported the involvement of miRNAs in the draught tolerance of M. cordata and could pave the way for further investigation into the related regulatory mechanisms.

Full-length transcriptome sequencing reveals the low-temperature-tolerance mechanism of Medicago falcata roots

BACKGROUND

Low temperature is one of the main environmental factors that limits crop growth, development, and production. Medicago falcata is an important leguminous herb that is widely distributed worldwide. M. falcata is related to alfalfa but is more tolerant to low temperature than alfalfa. Understanding the low temperature tolerance mechanism of M. falcata is important for the genetic improvement of alfalfa.

RESULTS

In this study, we explored the transcriptomic changes in the roots of low-temperature-treated M. falcata plants by combining SMRT sequencing and NGS technologies. A total of 115,153 nonredundant sequences were obtained, and 8849 AS events, 73,149 SSRs, and 4189 lncRNAs were predicted. A total of 111,587 genes from SMRT sequencing were annotated, and 11,369 DEGs involved in plant hormone signal transduction, protein processing in endoplasmic reticulum, carbon metabolism, glycolysis/gluconeogenesis, starch and sucrose metabolism, and endocytosis pathways were identified. We characterized 1538 TF genes into 45 TF gene families, and the most abundant TF family was the WRKY family, followed by the ERF, MYB, bHLH and NAC families. A total of 134 genes, including 101 whose expression was upregulated and 33 whose expression was downregulated, were differentially coexpressed at all five temperature points. PB40804, PB75011, PB110405 and PB108808 were found to play crucial roles in the tolerance of M. falcata to low temperature. WGCNA revealed that the MEbrown module was significantly correlated with low-temperature stress in M. falcata. Electrolyte leakage was correlated with most genetic modules and verified that electrolyte leakage can be used as a direct stress marker in physiological assays to indicate cell membrane damage from low-temperature stress. The consistency between the qRT-PCR results and RNA-seq analyses confirmed the validity of the RNA-seq data and the analysis of the regulatory mechanism of low-temperature stress on the basis of the transcriptome.

CONCLUSIONS

The full-length transcripts generated in this study provide a full characterization of the transcriptome of M. falcata and may be useful for mining new low-temperature stress-related genes specific to M. falcata. These new findings could facilitate the understanding of the low-temperature-tolerance mechanism of M. falcata.

Full-length transcriptome analysis of shade-induced promotion of tuber production in Pinellia ternata

BACKGROUND

Pinellia ternata is native to China and has been used as a traditional herb due to its antiemetic, antitussive, analgesic, and anxiolytic effects. When exposed to strong light intensity and high temperature during the reproductive growth process, P. ternata withers in a phenomenon known as "sprout tumble", which largely limits tuber production. Shade was previously found to delay sprout tumble formation (STF); however, no information exists regarding this process at the molecular level. Hence, we determined the genes involved in tuber development and STF in P. ternata.

RESULTS

Compared to that with natural sun-light (control), shade significantly induced chlorophyll accumulation, increased chlorophyll fluorescence parameters including initial fluorescence, maximal fluorescence, and qP, and dramatically repressed chlorophyll a:b and NPQ. Catalase (CAT) activity was largely induced by shade, and tuber products were largely increased in this environment. Transcriptome profiles of P. ternata grown in natural sun-light and shaded environments were analyzed by a combination of next generation sequencing (NGS) and third generation single-molecule real-time (SMRT) sequencing. Corrections of SMRT long reads based on NGS short reads yielded 136,163 non-redundant transcripts, with an average N50 length of 2578 bp. In total, 6738 deferentially-expressed genes (DEGs) were obtained from the comparisons, specifically D5S vs D5CK, D20S vs D20CK, D20S vs D5S, and D20CK vs D5CK, of which, 6384 DEGs (94.8%) were generated from the D20S vs D20CK comparison. Gene annotation and functional analyses revealed that these genes were related to auxin signal transduction, polysaccharide and sugar metabolism, phenylpropanoid biosynthesis, and photosynthesis. Moreover, the expression of genes enriched in photosynthesis appeared to be significantly altered by shade. The expression patterns of 16 candidate genes were consistent with changes in their transcript abundance as identified by RNA-Seq, and these might contribute to STF and tuber production.

CONCLUSION

The full-length transcripts identified in this study have provided a more accurate depiction of P. ternata gene transcription. Further, we identified potential genes involved in STF and tuber growth. Such data could serve as a genetic resource and a foundation for further research on this important traditional herb.

Characterization and Activity Analyses of the FLOWERING LOCUS T Promoter in Gossypium Hirsutum

Flowering transition is a crucial development process in cotton (Gossypium hirsutum L.), and the flowering time is closely correlated with the timing of FLOWERING LOCUS T (FT) expression. However, the mechanism underlying the coordination of various cis-regulatory elements in the FT promoter of cotton has not been determined. In this study, a 5.9-kb promoter of FT was identified from cotton. A bioinformatics analysis showed that multiple insertion-deletion sites existed in the 5.9-kb promoter. Different expression levels of a reporter gene, and the induction by sequential deletions in GhFT promoter, demonstrated that 1.8-kb of the GhFT promoter was stronger than 4.2-, 4.8-, and 5.9-kb promoter fragments. The binding sites of the CONSTANS (CO) and NUCLEAR FACTOR Y transcription factors were located within the 1.0-kb sequence upstream of the FT transcription start site. A large number of repeat segments were identified in proximal promoter regions (-1.1 to -1.4 kb). A complementation analysis of deletion constructs between 1.0 and 1.8 kb of G. hirsutum, Gossypium arboretum, and Gossypium raimondii FT promoters revealed that the 1.0-kb fragment significantly rescued the late-flowering phenotype of the Arabidopsis FT loss-of-function mutant ft-10, whereas the 1.8-kb promoter only slightly rescued the late-flowering phenotype. Furthermore, the conserved CORE motif in the cotton FT promoter is an atypical TGTG(N2-3)ATG, but the number of arbitrary bases between TGTG and ATG is uncertain. Thus, the proximal FT promoter region might play an important role affecting the activity levels of FT promoters in cotton flowering.

Analysis of Chromosomal Numbers, Mitochondrial Genome, and Full-Length Transcriptome of Onychostoma brevibarba

Onychostoma brevibarba is a new discovered species which is distributed in Xiang Jiang River of the middle Chang Jiang basin in Hunan Province, South China. In this study, the ploidy levels of O. brevibarba were confirmed by counting chromosomal numbers and analyzing karyotype. The complete mitochondrial genome of O. brevibarba was determined and analyzed. Besides, we firstly performed the full-length transcriptome of O. brevibarba derived from 5 different tissues using the PacBio SMRT sequencing. The result shows that O. brevibarba was a diploid with 50 chromosomes [corrected]. The complete mitogenome of O. brevibarba was 16,602 bp in size and very similar (89.1-91.3%) to that of other Onychostoma species but was distinct from all congeners. The full-length transcriptome dataset of O. brevibarba comprised 120,239 unigenes. Among the unigenes, 91,542 were functionally annotated, whereas 26,794 were found to have two or more isoforms. This study could provide many new insights into cytology and molecular characteristics of O. brevibarba; it laid the foundation for further exploration of the genomic signatures of species of Onychostoma.

Identification of reliable reference genes for qRT-PCR in the ephemeral plant Arabidopsis pumila based on full-length transcriptome data

Arabidopsis pumila, an annual ephemeral plant, plays important roles in preventing wind and sand erosion, water and soil conservation, and microhabitat improvement in the North of Xinjiang, China. Studies of adaptive mechanisms in harsh desert environments at the genetic and genomic levels can be used to more effectively develop and protect this species. The quantitative real-time polymerase chain reaction (qRT-PCR) method is one of the essential means to achieve these goals, and the selection of an appropriate reference gene is the prerequisite for qRT-PCR. In this study, 10 candidate reference genes were identified from the full-length transcriptome data of A. pumila, and their expression stabilities under four abiotic stresses (drought, heat, cold and salt) and in seven different tissues (roots, hypocotyl, cotyledon, leaves, stems, flowers and siliques) were evaluated with four programmes geNorm, NormFinder, Bestkeeper and RefFinder. Although the most stable reference genes were variable under different treatments using different software, comprehensive ranking revealed that UEP and HAF1 under drought stress, UBQ9 and GAPDH under heat stress, UBC35 and GAPDH under cold stress, GAPDH and ACT1 under salt stress, and ACT1 and GAPDH in different tissues were the most stable reference genes. Moreover, GAPDH and UBQ9 were the most suitable reference gene combinations for all samples. The expression pattern of the K⁺ uptake permease gene KUP9 further validated that the selected reference genes were suitable for normalization of gene expression. The identification of reliable reference genes guarantees more accurate qRT-PCR quantification for A. pumila and facilitates functional genomics studies of ephemeral plants.

Analysis of Transcriptome and Epitranscriptome in Plants Using PacBio Iso-Seq and Nanopore-Based Direct RNA Sequencing

Nanopore sequencing from Oxford Nanopore Technologies (ONT) and Pacific BioSciences (PacBio) single-molecule real-time (SMRT) long-read isoform sequencing (Iso-Seq) are revolutionizing the way transcriptomes are analyzed. These methods offer many advantages over most widely used high-throughput short-read RNA sequencing (RNA-Seq) approaches and allow a comprehensive analysis of transcriptomes in identifying full-length splice isoforms and several other post-transcriptional events. In addition, direct RNA-Seq provides valuable information about RNA modifications, which are lost during the PCR amplification step in other methods. Here, we present a comprehensive summary of important applications of these technologies in plants, including identification of complex alternative splicing (AS), full-length splice variants, fusion transcripts, and alternative polyadenylation (APA) events. Furthermore, we discuss the impact of the newly developed nanopore direct RNA-Seq in advancing epitranscriptome research in plants. Additionally, we summarize computational tools for identifying and quantifying full-length isoforms and other co/post-transcriptional events and discussed some of the limitations with these methods. Sequencing of transcriptomes using these new single-molecule long-read methods will unravel many aspects of transcriptome complexity in unprecedented ways as compared to previous short-read sequencing approaches. Analysis of plant transcriptomes with these new powerful methods that require minimum sample processing is likely to become the norm and is expected to uncover novel co/post-transcriptional gene regulatory mechanisms that control biological outcomes during plant development and in response to various stresses.