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

Phenotypic Physiological and Metabolomic Analyses Reveal Crucial Metabolic Pathways in Quinoa (Chenopodium quinoa Willd.) in Response to PEG-6000 Induced Drought Stress

Drought stress seriously threatens human food security, and enhancing crops' drought tolerance is an urgent problem to be solved in breeding. Quinoa is known for its high nutritional value and strong drought tolerance, but its molecular mechanism in response to drought stress is still unclear. In this study, we used drought-tolerant (D2) and drought-sensitive (ZK1) quinoa varieties, and PEG-6000 was used to simulate drought stress in quinoa seedlings. Phenotypic and physiological biochemical indicators were measured during the seedling stage, and LC-MS was used for a metabolite analysis of drought stress to explore the drought tolerance mechanism of quinoa under drought stress. With the intensification of drought stress, chlorophyll content gradually increased, and D2 reached its maximum at W4, an increase of 49.85% compared with W1. The total chlorophyll content, photosynthesis rate, and stomatal conductance of ZK1 were significantly lower than D2 under moderate and severe drought stress. Metabolomic results showed that a total of 1295 positive ion mode (pos) metabolites and 914 negative ion mode (neg) metabolites were identified. Of these, 12(R)-HETE, phosphatidylcholine, monogalactose diester (MGDG), and stachyose up-regulated expression under drought stress. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that unsaturated fatty acid biosynthesis and glycerophospholipid metabolism pathways were significantly enriched. In summary, our results elucidate that quinoa responds to drought stress by accumulating chlorophyll and sugars, activating unsaturated fatty acid metabolism, and protecting the photosynthetic system. These findings provide new insights for the breeding of drought-tolerant quinoa varieties and the study of drought tolerance mechanisms.

Full-Length Transcriptomes Reconstruction Reveals Intraspecific Diversity in Hairy Vetch (Vicia villosa Roth) and Smooth Vetch (V. villosa Roth var. glabrescens)

Hairy vetch (Vicia villosa Roth) and smooth vetch (V. villosa Roth var. glabrescens) are important cover crops and legume forage with great economic and ecological values. Due to the large and highly heterozygous genome, full-length transcriptome reconstruction is a cost-effective route to mining their genetic resources. In this study, a hybrid sequencing approach combining SMRT and NGS technologies was applied. The results showed that 28,747 and 40,600 high-quality non-redundant transcripts with an average length of 1808 bp and 1768 bp were generated from hairy vetch and smooth vetch, including 24,864 and 35,035 open reading frames (ORFs), respectively. More than 96% of transcripts were annotated to the public databases, and around 25% of isoforms underwent alternative splicing (AS) events. In addition, 987 and 1587 high-confidence lncRNAs were identified in two vetches. Interestingly, smooth vetch contains more specific transcripts and orthologous clusters than hairy vetch, revealing intraspecific transcript diversity. The phylogeny revealed that they were clustered together and closely related to the genus Pisum. Furthermore, the estimation of Ka/Ks ratios showed that purifying selection was the predominant force. A putative 3-dehydroquinate dehydratase/shikimate dehydrogenase (DHD/SDH) gene underwent strong positive selection and might regulate phenotypic differences between hairy vetch and smooth vetch. Overall, our study provides a vital characterization of two full-length transcriptomes in Vicia villosa, which will be valuable for their molecular research and breeding.

Metabolomics combined with transcriptomics and physiology reveals the regulatory responses of soybean plants to drought stress

Drought, a prevalent environmental stressor, has had significant consequences on soybean (Glycine max L.), notably impeding its growth and production. Therefore, it is crucial to gain insight into the regulatory responses of soybean plants exposed to drought stress during soybean flowering in the field. In this study, the cultivar 'Liaodou 15' was performed light drought (LD, 24.3% soil moisture content), moderate drought (MD, 20.6% soil moisture content) and severe drought (SD, 16.9% soil moisture content) treatments at flowering stages of soybean and then rehydrated (30% soil moisture content) until harvest. The yield-related indicators were measured and revealed that MD and SD treatments significantly reduced 6.3% and 10.8% of the 100-grain weight. Soybean plants subjected to three drought stresses showed that net photosynthetic rates were 20.8%, 51.5% and 71.8% lower in LD, MD and SD than that of CK. The WUE increased by 31.8%, 31.5% and 18.8% under three drought stress treatments compared to CK. In addition, proline content was 25.94%, 41.01% and 65.43% greater than that of CK under three drought stress treatments. The trend of the MDA content was consistent with that of the proline content. SOD activity was significantly increasing by 10.86%, 46.73% and 14.54% under three drought stress treatments. The activity of CAT in the SD treatment increased by 49.28%. All the indices recovered after rehydration. Furthermore, 54,78 and 51 different expressed metabolomics (DEMs) were identified in the LDCK/LD, MDCK/MD and SDCK/SD groups, respectively. There were 1,211, 1,265 and 1,288 different expressed genes (DEGs) were upregulated and 1,003, 1,819 and 1,747 DEGs were downregulated. Finally, combined transcriptomic and metabolomic analysis suggested that 437 DEGs and 24 DEMs of LDCK/LD group, 741 DEGs and 35 DEMs of MDCK/MD group, 633 DEGs and 23 DEMs of SDCK/SD group, were highly positively correlated in soybean plants under drought stress. Drought stress induced the expression of the PAO1, PAO4, PAO5 and P5CS genes to promote the accumulation of spermidine and proline. Our study elucidates the responses of drought-stressed soybean plants in the field and provides a genetic basis for the breeding of drought-tolerant soybean plants.

Preferentially expressed endosperm genes reveal unique activities in wheat endosperm during grain filling

BACKGROUND

Bread wheat (Triticum aestivum L.) endosperm contains starch and proteins, which determine the final yield, quality, and nutritional value of wheat grain. The preferentially expressed endosperm genes can precisely provide targets in the endosperm for improving wheat grain quality and nutrition using modern bioengineering technologies. However, the genes specifically expressed in developing endosperms remain largely unknown.

RESULTS

In this study, 315 preferentially expressed endosperm genes (PEEGs) in the spring wheat landrace, Chinese Spring, were screened using data obtained from an open bioinformatics database, which reveals a unique grain reserve deposition process and special signal transduction in a developing wheat endosperm. Furthermore, transcription and accumulation of storage proteins in the wheat cultivar, XC26 were evaluated. The results revealed that 315 PEEG plays a critical role in storage protein fragment deposition and is a potential candidate for modifying grain quality and nutrition.

CONCLUSION

These results provide new insights into endosperm development and candidate genes and promoters for improving wheat grain quality through genetic engineering and plant breeding techniques.

Differentially Expressed Genes Identification of Kohlrabi Seedlings (Brassica oleracea var. caulorapa L.) under Polyethylene Glycol Osmotic Stress and AP2/ERF Transcription Factor Family Analysis

Osmotic stress is a condition in which plants do not get enough water due to changes in environmental factors. Plant response to osmotic stress is a complex process involving the interaction of different stress-sensitive mechanisms. Differentially expressed genes and response mechanisms of kohlrabi have not been reported under osmotic stress. A total of 196,642 unigenes and 33,040 differentially expressed unigenes were identified in kohlrabi seedlings under polyethylene glycol osmotic stress. AP2/ERF, NAC and eight other transcription factor family members with a high degree of interaction with CAT and SOD antioxidant enzyme activity were identified. Subsequently, 151 AP2/ERF genes were identified and analyzed. Twelve conserved motifs were searched and all AP2/ERF genes were clustered into four groups. A total of 149 AP2/ERF genes were randomly distributed on the chromosome, and relative expression level analysis showed that BocAP2/ERF genes of kohlrabi have obvious specificity in different tissues. This study lays a foundation for explaining the osmotic stress resistance mechanism of kohlrabi and provides a theoretical basis for the functional analysis of BocAP2/ERF transcription factor family members.

Genome-wide expression analysis of novel heat-responsive microRNAs and their targets in contrasting wheat genotypes at reproductive stage under terminal heat stress

Introduction

Heat stress at terminal stage of wheat is critical and leads to huge yield losses worldwide. microRNAs (miRNAs) play significant regulatory roles in gene expression associated with abiotic and biotic stress at the post-transcriptional level.

Methods

In the present study, we carried out a comparative analysis of miRNAs and their targets in flag leaves as well as developing seeds of heat tolerant (RAJ3765) and heat susceptible (HUW510) wheat genotypes under heat stress and normal conditions using small RNA and degradome sequencing.

Results and discussion

A total of 84 conserved miRNAs belonging to 35 miRNA families and 93 novel miRNAs were identified in the 8 libraries. Tae-miR9672a-3p, tae-miR9774, tae-miR9669-5p, and tae-miR5048-5p showed the highest expression under heat stress. Tae-miR9775, tae-miR9662b-3p, tae-miR1120a, tae-miR5084, tae-miR1122a, tae-miR5085, tae-miR1118, tae-miR1130a, tae-miR9678-3p, tae-miR7757-5p, tae-miR9668-5p, tae-miR5050, tae-miR9652-5p, and tae-miR9679-5p were expressed only in the tolerant genotype, indicating their role in heat tolerance. Comparison between heat-treated and control groups revealed that 146 known and 57 novel miRNAs were differentially expressed in the various tissues. Eight degradome libraries sequence identified 457 targets of the differentially expressed miRNAs. Functional analysis of the targets indicated their involvement in photosynthesis, spliceosome, biosynthesis of nucleotide sugars and protein processing in the endoplasmic reticulum, arginine and proline metabolism and endocytosis.

Conclusion

This study increases the number of identified and novel miRNAs along with their roles involved in heat stress response in contrasting genotypes at two developing stages of wheat.

Polyploidization: A Biological Force That Enhances Stress Resistance

Organisms with three or more complete sets of chromosomes are designated as polyploids. Polyploidy serves as a crucial pathway in biological evolution and enriches species diversity, which is demonstrated to have significant advantages in coping with both biotic stressors (such as diseases and pests) and abiotic stressors (like extreme temperatures, drought, and salinity), particularly in the context of ongoing global climate deterioration, increased agrochemical use, and industrialization. Polyploid cultivars have been developed to achieve higher yields and improved product quality. Numerous studies have shown that polyploids exhibit substantial enhancements in cell size and structure, physiological and biochemical traits, gene expression, and epigenetic modifications compared to their diploid counterparts. However, some research also suggested that increased stress tolerance might not always be associated with polyploidy. Therefore, a more comprehensive and detailed investigation is essential to complete the underlying stress tolerance mechanisms of polyploids. Thus, this review summarizes the mechanism of polyploid formation, the polyploid biochemical tolerance mechanism of abiotic and biotic stressors, and molecular regulatory networks that confer polyploidy stress tolerance, which can shed light on the theoretical foundation for future research.

Deciphering the evolution and complexity of wheat germplasm from a genomic perspective

Bread wheat provides an essential fraction of the daily calorific intake for humanity. Due to its huge and complex genome, progress in studying on the wheat genome is substantially trailed behind those of the other two major crops, rice and maize, for at least a decade. With rapid advances in genome assembling and reduced cost of high-throughput sequencing, emerging de novo genome assemblies of wheat and whole-genome sequencing data are leading to a paradigm shift in wheat research. Here, we review recent progress in dissecting the complex genome and germplasm evolution of wheat since the release of the first high-quality wheat genome. New insights have been gained in the evolution of wheat germplasm during domestication and modern breeding progress, genomic variations at multiple scales contributing to the diversity of wheat germplasm, and complex transcriptional and epigenetic regulations of functional genes in polyploid wheat. Genomics databases and bioinformatics tools meeting the urgent needs of wheat genomics research are also summarized. The ever-increasing omics data, along with advanced tools and well-structured databases, are expected to accelerate deciphering the germplasm and gene resources in wheat for future breeding advances.

Emerging roles of long non-coding RNAs in regulating agriculturally important seed traits

Seeds have enormous economic importance as source of calories, nutrition, edible oil, and biofuels. Therefore, seed traits like seed size and shape, weight, micronutrient content, oil content, quality, post-harvest sprouting, etc., are some of the main targets in crop improvement. Designing the strategies for their improvement benefits heavily from understanding the regulatory aspects of seed development. Recent studies indicate that long non-coding RNAs (lncRNAs) are one of the important regulators of seed development. They played a significant role in crop domestication by influencing seed traits. LncRNAs are conventionally defined as non-coding RNAs greater than 200 bp in length but lacking protein coding potential. Here we highlight the emerging pieces of evidence of lncRNA-mediated regulation of seed development through diverse mechanisms, for instance, by acting as target mimics or precursors of regulatory small RNAs or through chromatin remodeling and post-transcriptional repression. We also enumerate the insights from high-throughput transcriptomic studies from developing seeds of cereal, oilseed, biofuel, and pulse crops. We highlight the lncRNA candidates and lncRNA-mediated regulatory networks regulating seed development and related agronomic traits. Further, we discuss the potential of lncRNAs for improvement of agriculturally important seed traits through marker-assisted breeding and/or transgenic approaches.

Integration of transcriptomic and metabolomic analysis unveils the response mechanism of sugar metabolism in Cyclocarya paliurus seedlings subjected to PEG-induced drought stress

Cyclocarya paliurus (Batal.) Iljinskaja is a multiple function tree species used for functional food and valued timber production. Carbohydrates, especially water-soluble carbohydrates, play an important role in osmotic protection, signal transduction and carbon storage. Under the circumstance of global climate change the abiotic stress would restrict the development of C. paliurus plantation, whereas there is few knowledge on the regulatory mechanisms of sugar metabolism under drought stress in C. paliurus. To investigate the drought response of C. paliurus at molecular level, we conducted an integrated analysis of transcriptomic and metabolomic of C. paliurus at three PEG-induced drought stress levels (0%: control; 15%: moderate drought; 25%: severe drought) in short term. Both moderate and severe drought treatments activated the chemical defense with lowering relative water content, and enhancing the contents of soluble protein, proline and malondialdehyde in the leaves. Meanwhile, alterations in the expression of differentially expressed genes and carbohydrate metabolism profiles were observed among the treatments. Weighted gene co-expression network analysis (WGCNA) showed 3 key modules, 8 structural genes (such as genes encoding beta-fructofuranosidase (INV), sucrose synthase (SUS), raffinose synthase (RS)) and 14 regulatory transcription factors were closely linked to sugar metabolism. Our results provided the foundation to understand the response mechanism of sugar metabolism in C. paliurus under drought stress, and would drive progress in breeding of drought-tolerant varieties and plantation development of the species.

PacBio single-molecule long-read sequencing provides new insights into the complexity of full-length transcripts in oriental river prawn, macrobrachium nipponense

BACKGROUND

Oriental river prawn (Macrobrachium nipponense) is one of the most dominant species in shrimp farming in China, which is a rich source of protein and contributes to a significant impact on the quality of human life. Thus, more complete and accurate annotation of gene models are important for the breeding research of oriental river prawn.

RESULTS

A full-length transcriptome of oriental river prawn muscle was obtained using the PacBio Sequel platform. Then, 37.99 Gb of subreads were sequenced, including 584,498 circular consensus sequences, among which 512,216 were full length non-chimeric sequences. After Illumina-based correction of long PacBio reads, 6,599 error-corrected isoforms were identified. Transcriptome structural analysis revealed 2,263 and 2,555 alternative splicing (AS) events and alternative polyadenylation (APA) sites, respectively. In total, 620 novel genes (NGs), 197 putative transcription factors (TFs), and 291 novel long non-coding RNAs (lncRNAs) were identified.

CONCLUSIONS

In summary, this study offers novel insights into the transcriptome complexity and diversity of this prawn species, and provides valuable information for understanding the genomic structure and improving the draft genome annotation of oriental river prawn.

A wheat integrative regulatory network from large-scale complementary functional datasets enables trait-associated gene discovery for crop improvement

Gene regulation is central to all aspects of organism growth, and understanding it using large-scale functional datasets can provide a whole view of biological processes controlling complex phenotypic traits in crops. However, the connection between massive functional datasets and trait-associated gene discovery for crop improvement is still lacking. In this study, we constructed a wheat integrative gene regulatory network (wGRN) by combining an updated genome annotation and diverse complementary functional datasets, including gene expression, sequence motif, transcription factor (TF) binding, chromatin accessibility, and evolutionarily conserved regulation. wGRN contains 7.2 million genome-wide interactions covering 5947 TFs and 127 439 target genes, which were further verified using known regulatory relationships, condition-specific expression, gene functional information, and experiments. We used wGRN to assign genome-wide genes to 3891 specific biological pathways and accurately prioritize candidate genes associated with complex phenotypic traits in genome-wide association studies. In addition, wGRN was used to enhance the interpretation of a spike temporal transcriptome dataset to construct high-resolution networks. We further unveiled novel regulators that enhance the power of spike phenotypic trait prediction using machine learning and contribute to the spike phenotypic differences among modern wheat accessions. Finally, we developed an interactive webserver, wGRN (http://wheat.cau.edu.cn/wGRN), for the community to explore gene regulation and discover trait-associated genes. Collectively, this community resource establishes the foundation for using large-scale functional datasets to guide trait-associated gene discovery for crop improvement.

Transcriptome Analysis of Heat Shock Factor C2a Over-Expressing Wheat Roots Reveals Ferroptosis-like Cell Death in Heat Stress Recovery

Wheat (Triticum aestivum L.) growing areas in many regions of the world are subject to heat waves which are predicted to increase in frequency because of climate change. The engineering of crop plants can be a useful strategy to mitigate heat stress-caused yield losses. Previously, we have shown that heat shock factor subclass C (TaHsfC2a-B)-overexpression significantly increased the survival of heat-stressed wheat seedlings. Although previous studies have shown that the overexpression of Hsf genes enhanced the survival of plants under heat stress, the molecular mechanisms are largely unknown. To understand the underlying molecular mechanisms involved in this response, a comparative analysis of the root transcriptomes of untransformed control and TaHsfC2a-overexpressing wheat lines by RNA-sequencing have been performed. The results of RNA-sequencing indicated that the roots of TaHsfC2a-overexpressing wheat seedlings showed lower transcripts of hydrogen peroxide-producing peroxidases, which corresponds to the reduced accumulation of hydrogen peroxide along the roots. In addition, suites of genes from iron transport and nicotianamine-related gene ontology categories showed lower transcript abundance in the roots of TaHsfC2a-overexpressing wheat roots than in the untransformed control line following heat stress, which are in accordance with the reduction in iron accumulation in the roots of transgenic plants under heat stress. Overall, these results suggested the existence of ferroptosis-like cell death under heat stress in wheat roots, and that TaHsfC2a is a key player in this mechanism. To date, this is the first evidence to show that a Hsf gene plays a key role in ferroptosis under heat stress in plants. In future, the role of Hsf genes could be further studied on ferroptosis in plants to identify root-based marker genes to screen for heat-tolerant genotypes.

Regulatory Mechanisms of Heat Stress Response and Thermomorphogenesis in Plants

As worldwide warming intensifies, the average temperature of the earth continues to increase. Temperature is a key factor for the growth and development of all organisms and governs the distribution and seasonal behavior of plants. High temperatures lead to various biochemical, physiological, and morphological changes in plants and threaten plant productivity. As sessile organisms, plants are subjected to various hostile environmental factors and forced to change their cellular state and morphological architecture to successfully deal with the damage they suffer. Therefore, plants have evolved multiple strategies to cope with an abnormal rise in temperature. There are two main mechanisms by which plants respond to elevated environmental temperatures. One is the heat stress response, which is activated under extremely high temperatures; the other is the thermomorphogenesis response, which is activated under moderately elevated temperatures, below the heat-stress range. In this review, we summarize recent progress in the study of these two important heat-responsive molecular regulatory pathways mediated, respectively, by the Heat Shock Transcription Factor (HSF)-Heat Shock Protein (HSP) pathway and PHYTOCHROME INTER-ACTING FACTOR 4 (PIF4) pathways in plants and elucidate the regulatory mechanisms of the genes involved in these pathways to provide comprehensive data for researchers studying the heat response. We also discuss future perspectives in this field.

Heat Stress Tolerance 2 confers basal heat stress tolerance in allohexaploid wheat (Triticum aestivum L.)

Heat stress substantially reduces the yield potential of wheat (Triticum aestivum L.), one of the most widely cultivated staple crops, and greatly threatens global food security in the context of global warming. However, few studies have explored the heat stress tolerance (HST)-related genetic resources in wheat. Here, we identified and fine-mapped a wheat HST locus, TaHST2, which is indispensable for HST in both the vegetative and reproductive stages of the wheat life cycle. The studied pair of near isogenic lines (NILs) exhibited diverse morphologies under heat stress, based on which we mapped TaHST2 to a 485 kb interval on chromosome arm 4DS. Under heat stress, TaHST2 confers a superior conversion rate from soluble sugars to starch in wheat grains, resulting in faster grain filling and a higher yield potential. A further exploration of genetic resources indicated that TaHST2 underwent strong artificial selection during wheat domestication, suggesting it is an essential locus for basal HST in wheat. Our findings provide deeper insights into the genetic basis of wheat HST and might be useful for global efforts to breed heat-stress-tolerant cultivars.

Evolution of isoform-level gene expression patterns across tissues during lotus species divergence

Both gene duplication and alternative splicing (AS) drive the functional diversity of gene products in plants, yet the relative contributions of the two key mechanisms to the evolution of gene function are largely unclear. Here, we studied AS in two closely related lotus plants, Nelumbo lutea and Nelumbo nucifera, and the outgroup Arabidopsis thaliana, for both single-copy and duplicated genes. We show that most splicing events evolved rapidly between orthologs and that the origin of lineage-specific splice variants or isoforms contributed to gene functional changes during species divergence within Nelumbo. Single-copy genes contain more isoforms, have more AS events conserved across species, and show more complex tissue-dependent expression patterns than their duplicated counterparts. This suggests that expression divergence through isoforms is a mechanism to extend the expression breadth of genes with low copy numbers. As compared to isoforms of local, small-scale duplicates, isoforms of whole-genome duplicates are less conserved and display a less conserved tissue bias, pointing towards their contribution to subfunctionalization. Through comparative analysis of isoform expression networks, we identified orthologous genes of which the expression of at least some of their isoforms displays a conserved tissue bias across species, indicating a strong selection pressure for maintaining a stable expression pattern of these isoforms. Overall, our study shows that both AS and gene duplication contributed to the diversity of gene function during the evolution of lotus.

Distribution, expression of hexaploid wheat Fes1s and functional characterization of two TaFes1As in Arabidopsis

Hexaploid wheat is a major food crop and is sensitive to heat stress. It is necessary to discover genes related to thermotolerance in wheat. Fes1s is a class of nucleotide exchange factor of heat shock protein 70s, proven to be participated in heat response in human, yeast, and Arabidopsis. However, little is known about Fes1s in hexaploid wheat. In this study, we identified nine Fes1s in hexaploid wheat (TaFes1s) and found that they present as three triads. A phylogenetic relationship analysis revealed that these Fes1s grouped into Fes1A, Fes1B and Fes1C subclades, and Fes1As and Fes1Bs were divergent in monocots, but possibly not in dicots. The sequences, gene structures and protein motifs of TaFes1s homoeologues within a triad were highly conserved. Through cis-elements analysis including heat shock elements, and miRNA targets prediction, we found that regulation of three TaFes1s homoeologues may be different, while the expression patterns of three homoeologues were similar. The expression levels of TaFes1As were higher than those of TaFes1Bs and TaFes1Cs, and based on these expressions, TaFes1As were chosen for functional characterization. Intriguingly, neither TaFes1A-5A nor TaFes1A-5D could not rescue the thermotolerance defect of Arabidopsis fes1a mutants at seedling stage, but in the transgenic plants seed germination was accelerated under normal and heat stress condition. The functional characterization indicated that roles of Fes1As would be different in Arabidopsis and hexaploid wheat, and function retention of TaFes1As may occur during wheat evolution. In conclusion, our study comprehensively characterized the distribution and expression of Fes1s in hexaploid wheat and found that two TaFes1As could accelerate seed germination under normal and heat stress condition.

Single-molecule Real-time (SMRT) Sequencing Facilitates Transcriptome Research and Genome Annotation of the Fish Sillago sinica

As a newly described Sillaginidae species, Chinese sillago (Sillago sinica) needs a better understanding of gene annotation information. In this study, we reported the first full-length transcriptome data of S. sinica using the PacBio isoform sequencing Iso-seq and a description of transcriptome structure analysis. A total of 454,979 high-quality full-length transcripts were obtained by single-molecule real-time (SMRT) sequencing, which was corrected by Illumina sequencing data. After that, 66,948 non-redundant full-length transcripts were generated after mapping to the reference genome of S. sinica, including 49 fusion isoforms and 9,250 novel isoforms. 63,459 isoforms were successfully annotated by one of the Nr, Nt, SwissProt, Pfam, KOG, GO, and KEGG databases. Additionally, 30,987 alternative polyadenylation (APA) sites, 451,867 alternative splicing (AS) events, 21,928 long non-coding RNAs (lncRNAs) and 12,911 transcription factors (TFs) were identified. The full-length transcripts of S. sinica would provide a precious resource for characterizing the transcriptome of S. sinica and for the further study of gene function and regulatory mechanism of this species.

Wheat genomic study for genetic improvement of traits in China

Bread wheat (Triticum aestivum L.) is a major crop that feeds 40% of the world's population. Over the past several decades, advances in genomics have led to tremendous achievements in understanding the origin and domestication of wheat, and the genetic basis of agronomically important traits, which promote the breeding of elite varieties. In this review, we focus on progress that has been made in genomic research and genetic improvement of traits such as grain yield, end-use traits, flowering regulation, nutrient use efficiency, and biotic and abiotic stress responses, and various breeding strategies that contributed mainly by Chinese scientists. Functional genomic research in wheat is entering a new era with the availability of multiple reference wheat genome assemblies and the development of cutting-edge technologies such as precise genome editing tools, high-throughput phenotyping platforms, sequencing-based cloning strategies, high-efficiency genetic transformation systems, and speed-breeding facilities. These insights will further extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture in China and throughout the world.

Small HSPs play an important role in crosstalk between HSF-HSP and ROS pathways in heat stress response through transcriptomic analysis in lilies (Lilium longiflorum)

BACKGROUND

High temperature seriously limits the annual production of fresh cut lilies, which is one of the four major cut flowers in the global cut flower market. There were few transcriptomes focused on the gene expression of lilies under heat stress. In order to reveal the potential heat response patterns in bulbous plants and provide important genes for further genetic engineering techniques to improve thermotolerance of lily, RNA sequencing of lilies under heat treatments were conducted.

RESULTS

In this study, seedlings of Lilium longiflorum 'White Heaven' were heat-treated at 37 °C for different lengths of time (0 h, 0.5 h, 1 h, 3 h, 6 h, and 12 h with a 12 h-light/12 h-dark cycle). The leaves of these lily seedlings were immediately collected after heat treatments and quickly put into liquid nitrogen for RNA sequencing. 109,364,486-171,487,430 clean reads and 55,044 unigenes including 21,608 differentially expressed genes (DEGs) (fold change ≥2) were obtained after heat treatment. The number of DEGs increased sharply during the heat treatments of 0.5 h-1 h and 1 h-3 h compared to that of other periods. Genes of the heat stress transcription factor (HSF) family and the small heat shock proteins (small HSPs, also known as HSP20) family responded to heat stress early and quickly. Compared to that of the calcium signal and hormone pathways, DEGs of the HSF-HSP pathway and reactive oxygen species (ROS) pathway were significantly and highly induced. Moreover, they had the similar expression pattern in response to heat stress. Small HSPs family genes were the major components in the 50 most highly induced genes at each heat stress treatment and involved in ROS pathway in the rapid response to heat stress. Furthermore, the barley stripe mosaic virus induced gene silencing (BSMV-VIGS) of LlHsfA2 caused a significantly reduced thermotolerance phenotype in Lilium longiflorum 'White Heaven', meanwhile decreasing the expression of small HSPs family genes and increasing the ROS scavenging enzyme ascorbate peroxidase (APX) genes, indicating the potential interplay between these two pathways.

CONCLUSIONS

Based on our transcriptomic analysis, we provide a new finding that small HSPs play important roles in crosstalk between HSF-HSP and ROS pathways in heat stress response of lily, which also supply the groundwork for understanding the mechanism of heat stress in bulbous plants.

Global survey of alternative splicing and gene modules associated with fertility regulation in a thermosensitive genic male sterile wheat

Thermosensitive genic male sterile (TGMS) wheat lines are the core of two-line hybrid systems. Understanding the mechanism that regulates male sterility in TGMS wheat lines is helpful for promoting wheat breeding. Several studies have obtained information regarding the mechanisms associated with male sterility at the transcriptional level, but it is not clear how the post-transcriptional process of alternative splicing might contribute to controlling male sterility. In this study, we performed genome-wide analyses of alternative splicing during the meiosis stage in TGMS line BS366 using PacBio and RNA-Seq hybrid sequencing. Cytological observations indicated that cytoskeleton assembly in pollen cells, calcium deposition in pollen and tapetal cells, and vesicle transport in tapetal cells were deficient in BS366. According to our cytological findings, 49 differentially spliced genes were isolated. Moreover, 25 long non-coding RNA targets and three bHLH transcription factors were identified. Weighted gene co-expression network analysis detected four candidate differentially spliced genes that had strong co-relation with the seed setting percentage, which is the direct representation of male sterility in BS366. In this study, we obtained comprehensive data regarding the alternative splicing-mediated regulation of male sterility in TGMS wheat. The candidates identified may provide the molecular basis for an improved understanding of male sterility.

Combined Analysis of the Metabolome and Transcriptome to Explore Heat Stress Responses and Adaptation Mechanisms in Celery (Apium graveolens L.)

Celery is an important leafy vegetable that can grow during the cool season and does not tolerate high temperatures. Heat stress is widely acknowledged as one of the main abiotic stresses affecting the growth and yield of celery. The morphological and physiological indices of celery were investigated in the present study to explore the physiological mechanisms in response to high temperatures. Results showed that the antioxidant enzyme activity, proline, relative conductivity, and malondialdehyde were increased, while chlorophyll and the water content of leaves decreased under high-temperature conditions. Short-term heat treatment increased the stomatal conductance to cool off the leaves by transpiration; however, long-term heat treatment led to stomatal closure to prevent leaf dehydration. In addition, high temperature caused a disordered arrangement of palisade tissue and a loose arrangement of spongy tissue in celery leaves. Combined metabolomic and transcriptomic analyses were further used to reveal the regulatory mechanisms in response to heat stress at the molecular level in celery. A total of 1003 differential metabolites were identified and significantly enriched in amino acid metabolism and the tricarboxilic acid (TCA) cycle. Transcriptome sequencing detected 24,264 different genes, including multiple transcription factor families such as HSF, WRKY, MYB, AP2, bZIP, and bHLH family members that were significantly upregulated in response to heat stress, suggesting that these genes were involved in the response to heat stress. In addition, transcriptional and metabolic pathway analyses showed that heat stress inhibited the glycolysis pathway and delayed the TCA cycle but increased the expression of most amino acid synthesis pathways such as proline, arginine, and serine, consistent with the results of physiological indicators. qRT-PCR further showed that the expression pattern was similar to the expression abundance in the transcriptome. The important metabolites and genes in celery that significantly contributed to the response to high temperatures were identified in the present study, which provided the theoretical basis for breeding heat-resistant celery.

Hybrid RNA Sequencing Strategy for the Dynamic Transcriptomes of Winter Dormancy in an Evergreen Herbaceous Perennial, Iris japonica

Japanese iris (Iris japonica) is a popular perennial ornamental that originated in China; it has a long display period and remains green outdoors throughout the year. winter dormancy characteristics contribute greatly to the evergreenness of herbaceous perennials. Thus, it is crucial to explore the mechanism of winter dormancy in this evergreen herbaceous perennial. Here, we used the hybrid RNA-seq strategy including single-molecule real-time (SMRT) and next-generation sequencing (NGS) technologies to generate large-scale Full-length transcripts to examine the shoot apical meristems of Japanese iris. A total of 10.57 Gb clean data for SMRT and over 142 Gb clean data for NGS were generated. Using hybrid error correction, 58,654 full-length transcripts were acquired and comprehensively analysed, and their expression levels were validated by real-time qPCR. This is the first full-length RNA-seq study in the Iris genus; our results provide a valuable resource and improve understanding of RNA processing in this genus, for which little genomic information is available as yet. In addition, our data will facilitate in-depth analyses of winter dormancy mechanisms in herbaceous perennials, especially evergreen monocotyledons.

Transcriptomic and Metabolomic Analyses Reveal Key Metabolites, Pathways and Candidate Genes in Sophora davidii (Franch.) Skeels Seedlings Under Drought Stress

Soil aridification and desertification are particularly prominent in China's karst areas, severely limiting crop yields and vegetation restoration. Therefore, it is very important to identify naturally drought-tolerant plant species. Sophora davidii (Franch.) Skeels is resistant to drought and soil infertility, is deeply rooted and is an excellent plant material for soil and water conservation. We studied the transcriptomic and metabolomic changes in S. davidii in response to drought stress (CK, control; LD, mild drought stress; MD, moderate drought stress; and SD, severe drought stress). Sophora davidii grew normally under LD and MD stress but was inhibited under SD stress; the malondialdehyde (MDA), hydrogen peroxide (H₂O₂), soluble sugar, proline, chlorophyll a, chlorophyll b and carotenoid contents and ascorbate peroxidase (APX) activity significantly increased, while the superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT) activities and soluble protein content significantly decreased. In the LD/CK, MD/CK and SD/CK comparison groups, there were 318, 734 and 1779 DEGs, respectively, and 100, 168 and 281 differentially accumulated metabolites, respectively. Combined analysis of the transcriptomic and metabolomic data revealed the metabolic regulation of S. davidii in response to drought stress. First, key candidate genes such as PRR7, PRR5, GI, ELF3, PsbQ, PsaK, INV, AMY, E2.4.1.13, E3.2.1.2, NCED, PP2C, PYL, ABF, WRKY33, P5CS, PRODH, AOC3, HPD, GPX, GST, CAT and SOD1 may govern the drought resistance of S. davidii. Second, three metabolites (oxidised glutathione, abscisic acid and phenylalanine) were found to be related to drought tolerance. Third, several key candidate genes and metabolites involved in 10 metabolic pathways were identified, indicating that these metabolic pathways play an important role in the response to drought in S. davidii and possibly other plant species.

Genome-Wide Analysis of HSP70s in Hexaploid Wheat: Tandem Duplication, Heat Response, and Regulation

HSP70s play crucial roles in plant growth and development, as well as in stress response. Knowledge of the distribution and heat response of HSP70s is important to understand heat adaptation and facilitate thermotolerance improvement in wheat. In this study, we comprehensively analyzed the distribution of HSP70s in hexaploid wheat (TaHSP70s) and its relatives, and we found an obvious expansion of TaHSP70s in the D genome of hexaploid wheat. Meanwhile, a large portion of tandem duplication events occurred in hexaploid wheat. Among the 84 identified TaHSP70s, more than 64% were present as homeologs. The expression profiles of TaHSP70s in triads tended to be expressed more in non-stressful and heat stress conditions. Intriguingly, many TaHSP70s were especially heat responsive. Tandem duplicated TaHSP70s also participated in heat response and growth development. Further HSE analysis revealed divergent distribution of HSEs in the promoter regions of TaHSP70 homeologs, which suggested a distinct heat regulatory mechanism. Our results indicated that the heat response of TaHSP70s may experience a different regulation, and this regulation, together with the expression of tandem duplicated TaHSP70s, may help hexaploid wheat to adapt to heat conditions.

The activation of gene expression and alternative splicing in the formation and evolution of allopolyploid Brassica napus

Allopolyploids contain two or more sets of subgenomes. To establish a compatible relationship between subgenomes, a series of gene expression changes occurred in allopolyploids. What evolutionary changes of transcripts have taken place in Brassica napus during the early establishment and subsequent evolution was a fascinating scientific question. Here, we study this issue using a set of materials (natural, resynthesized B. napus and their progenitors/parents) by long-read RNA sequencing technology. The results showed that more genes were up-regulated in resynthesized B. napus compared with its two parents, and more up-regulated expressed genes were observed in natural B. napus compared with resynthesized B. napus. The presence of up-regulation genes in organism may help it adapt to the influence of "genomic shock" and cope with natural environment. Isoforms are produced from precursor mRNAs by alternative splicing (AS) events, and more than 60% of novel isoforms were identified in all materials, which could improve the reference genome information of B. napus. We found that the isoform numbers, the number of genes potentially involved in AS and alternative polyadenylation increased in B. napus after evolution, which may involve in the adaptation of plants to natural environment. In addition, all identified isoforms were functional annotated by searching 7 databases. In general, this study could improve our overall understanding of the full-length transcriptome of B. napus, and help us recognize the significant gene expression changes and isoform abundance changes occurred in allopolyploid B. napus during evolution.

Transcriptomic and Physiological Response of Durum Wheat Grain to Short-Term Heat Stress during Early Grain Filling

In a changing climate, extreme weather events such as heatwaves will be more frequent and could affect grain weight and the quality of crops such as wheat, one of the most significant crops in terms of global food security. In this work, we characterized the response of Triticum turgidum L. spp. durum wheat to short-term heat stress (HS) treatment at transcriptomic and physiological levels during early grain filling in glasshouse experiments. We found a significant reduction in grain weight (23.9%) and grain dimensions from HS treatment. Grain quality was also affected, showing a decrease in starch content (20.8%), in addition to increments in grain protein levels (14.6%), with respect to the control condition. Moreover, RNA-seq analysis of durum wheat grains allowed us to identify 1590 differentially expressed genes related to photosynthesis, response to heat, and carbohydrate metabolic process. A gene regulatory network analysis of HS-responsive genes uncovered novel transcription factors (TFs) controlling the expression of genes involved in abiotic stress response and grain quality, such as a member of the DOF family predicted to regulate glycogen and starch biosynthetic processes in response to HS in grains. In summary, our results provide new insights into the extensive transcriptome reprogramming that occurs during short-term HS in durum wheat grains.

Molecular Analysis Uncovers the Mechanism of Fertility Restoration in Temperature-Sensitive Polima Cytoplasmic Male-Sterile Brassica napus

Temperature-sensitive male sterility is a heritable agronomic trait affected by genotype-environment interactions. In rapeseed (Brassica napus), Polima (pol) temperature-sensitive cytoplasmic male sterility (TCMS) is commonly used for two-line breeding, as the fertility of pol TCMS lines can be partially restored at certain temperatures. However, little is known about the underlying molecular mechanism that controls fertility restoration. Therefore, we aimed to investigate the fertility conversion mechanism of the pol TCMS line at two different ambient temperatures (16 °C and 25 °C). Our results showed that the anthers developed and produced vigorous pollen at 16 °C but not at 25 °C. In addition, we identified a novel co-transcript of orf224-atp6 in the mitochondria that might lead to fertility conversion of the pol TCMS line. RNA-seq analysis showed that 1637 genes were significantly differentially expressed in the fertile flowers of 596-L when compared to the sterile flower of 1318 and 596-H. Detailed analysis revealed that differentially expressed genes were involved in temperature response, ROS accumulation, anther development, and mitochondrial function. Single-molecule long-read isoform sequencing combined with RNA sequencing revealed numerous genes produce alternative splicing transcripts at high temperatures. Here, we also found that alternative oxidase, type II NAD(P)H dehydrogenases, and transcription factor Hsfs might play a crucial role in male fertility under the low-temperature condition. RNA sequencing and bulked segregant analysis coupled with whole-genome sequencing identified the candidate genes involved in the post-transcriptional modification of orf224. Overall, our study described a putative mechanism of fertility restoration in a pol TCMS line controlled by ambient temperature that might help utilise TCMS in the two-line breeding of Brassica crops.

Single-molecule real-time transcript sequencing of developing cotton anthers facilitates genome annotation and fertility restoration candidate gene discovery

Heterosis refers to the superior phenotypes observed in hybrids. Cytoplasmic male sterility (CMS) system plays an important role in cotton heterosis utilization. However, the global gene expression patterns of CMS-D2 and its interaction with the restorer gene Rf1 remain unclear. Here, the full-length transcript sequencing was performed in anthers of the CMS-D2 restorer line using PacBio single-molecule real-time sequencing technology. Combining PacBio SMRT long-read isoforms and Illumina RNA-seq data, 107,066 isoforms from 44,338 loci were obtained, including 10,086 novel isoforms of novel genes and 66,419 new isoforms of known genes. Totally 56,572 alternative splicing (AS) events, 1146 lncRNAs, 61 fusion transcripts and 10,466 genes exhibited alternative polyadenylation (APA), and 60,995 novel isoforms with predicted open reading frames (ORFs) were further identified. Furthermore, the specifically expressed genes in restorer line were selected and confirmed by qRT-PCR. These findings provide a basis for upland cotton genome annotation and transcriptome research, and will help to reveal the molecular mechanism of interaction between Rf1 and CMS-D2 cytoplasm.

PacBio Single-Molecule Long-Read Sequencing Provides New Light on the Complexity of Full-Length Transcripts in Cattle

Cattle (Bos taurus) is one of the most widely distributed livestock species in the world, and provides us with high-quality milk and meat which have a huge impact on the quality of human life. Therefore, accurate and complete transcriptome and genome annotation are of great value to the research of cattle breeding. In this study, we used error-corrected PacBio single-molecule real-time (SMRT) data to perform whole-transcriptome profiling in cattle. Then, 22.5 Gb of subreads was generated, including 381,423 circular consensus sequences (CCSs), among which 276,295 full-length non-chimeric (FLNC) sequences were identified. After correction by Illumina short reads, we obtained 22,353 error-corrected isoforms. A total of 305 alternative splicing (AS) events and 3,795 alternative polyadenylation (APA) sites were detected by transcriptome structural analysis. Furthermore, we identified 457 novel genes, 120 putative transcription factors (TFs), and 569 novel long non-coding RNAs (lncRNAs). Taken together, this research improves our understanding and provides new insights into the complexity of full-length transcripts in cattle.

Comparative Analysis of Genomic and Transcriptome Sequences Reveals Divergent Patterns of Codon Bias in Wheat and Its Ancestor Species

The synonymous codons usage shows a characteristic pattern of preference in each organism. This codon usage bias is thought to have evolved for efficient protein synthesis. Synonymous codon usage was studied in genes of the hexaploid wheat Triticum aestivum (AABBDD) and its progenitor species, Triticum urartu (AA), Aegilops tauschii (DD), and Triticum turgidum (AABB). Triticum aestivum exhibited stronger usage bias for G/C-ending codons than did the three progenitor species, and this bias was especially higher compared to T. turgidum and Ae. tauschii. High GC content is a primary factor influencing codon usage in T. aestivum. Neutrality analysis showed a significant positive correlation (p<0.001) between GC12 and GC3 in the four species with regression line slopes near zero (0.16–0.20), suggesting that the effect of mutation on codon usage was only 16–20%. The GC3s values of genes were associated with gene length and distribution density within chromosomes. tRNA abundance data indicated that codon preference corresponded to the relative abundance of isoaccepting tRNAs in the four species. Both mutation and selection have affected synonymous codon usage in hexaploid wheat and its progenitor species. GO enrichment showed that GC biased genes were commonly enriched in physiological processes such as photosynthesis and response to acid chemical. In some certain gene families with important functions, the codon usage of small parts of genes has changed during the evolution process of T. aestivum.

Reading between the Lines: RNA-seq Data Mining Reveals the Alternative Message of the Rice Leaf Transcriptome in Response to Heat Stress

Rice (Oryza sativa L.) is a major food crop but heat stress affects its yield and grain quality. To identify mechanistic solutions to improve rice yield under rising temperatures, molecular responses of thermotolerance must be understood. Transcriptional and post-transcriptional controls are involved in a wide range of plant environmental responses. Alternative splicing (AS), in particular, is a widespread mechanism impacting the stress defence in plants but it has been completely overlooked in rice genome-wide heat stress studies. In this context, we carried out a robust data mining of publicly available RNA-seq datasets to investigate the extension of heat-induced AS in rice leaves. For this, datasets of interest were subjected to filtering and quality control, followed by accurate transcript-specific quantifications. Powerful differential gene expression (DE) and differential AS (DAS) identified 17,143 and 2162 heat response genes, respectively, many of which are novel. Detailed analysis of DAS genes coding for key regulators of gene expression suggests that AS helps shape transcriptome and proteome diversity in response to heat. The knowledge resulting from this study confirmed a widespread transcriptional and post-transcriptional response to heat stress in plants, and it provided novel candidates for rapidly advancing rice breeding in response to climate change.

Full-length transcriptome profiling reveals insight into the cold response of two kiwifruit genotypes (A. arguta) with contrasting freezing tolerances

BACKGROUND

Kiwifruit (Actinidia Lindl.) is considered an important fruit species worldwide. Due to its temperate origin, this species is highly vulnerable to freezing injury while under low-temperature stress. To obtain further knowledge of the mechanism underlying freezing tolerance, we carried out a hybrid transcriptome analysis of two A. arguta (Actinidi arguta) genotypes, KL and RB, whose freezing tolerance is high and low, respectively. Both genotypes were subjected to - 25 °C for 0 h, 1 h, and 4 h.

RESULTS

SMRT (single-molecule real-time) RNA-seq data were assembled using the de novo method, producing 24,306 unigenes with an N50 value of 1834 bp. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of DEGs showed that they were involved in the 'starch and sucrose metabolism', the 'mitogen-activated protein kinase (MAPK) signaling pathway', the 'phosphatidylinositol signaling system', the 'inositol phosphate metabolism', and the 'plant hormone signal transduction'. In particular, for 'starch and sucrose metabolism', we identified 3 key genes involved in cellulose degradation, trehalose synthesis, and starch degradation processes. Moreover, the activities of beta-GC (beta-glucosidase), TPS (trehalose-6-phosphate synthase), and BAM (beta-amylase), encoded by the abovementioned 3 key genes, were enhanced by cold stress. Three transcription factors (TFs) belonging to the AP2/ERF, bHLH (basic helix-loop-helix), and MYB families were involved in the low-temperature response. Furthermore, weighted gene coexpression network analysis (WGCNA) indicated that beta-GC, TPS5, and BAM3.1 were the key genes involved in the cold response and were highly coexpressed together with the CBF3, MYC2, and MYB44 genes.

CONCLUSIONS

Cold stress led various changes in kiwifruit, the 'phosphatidylinositol signaling system', 'inositol phosphate metabolism', 'MAPK signaling pathway', 'plant hormone signal transduction', and 'starch and sucrose metabolism' processes were significantly affected by low temperature. Moreover, starch and sucrose metabolism may be the key pathway for tolerant kiwifruit to resist low temperature damages. These results increase our understanding of the complex mechanisms involved in the freezing tolerance of kiwifruit under cold stress and reveal a series of candidate genes for use in breeding new cultivars with enhanced freezing tolerance.

Wheat heat tolerance is impaired by heightened deletions in the distal end of 4AL chromosomal arm

Heat stress (HS) causes substantial damages to worldwide crop production. As a cool season crop, wheat (Triticum aestivum) is sensitive to HS-induced damages. To support the genetic improvement of wheat HS tolerance (HST), we conducted fine mapping of TaHST1, a locus required for maintaining wheat vegetative and reproductive growth under elevated temperatures. TaHST1 was mapped to the distal terminus of 4AL chromosome arm using genetic populations derived from two BC₆ F₆ breeding lines showing tolerance (E6015-4T) or sensitivity (E6015-3S) to HS. The 4AL region carrying TaHST1 locus was approximately 0.949 Mbp and contained the last 19 high confidence genes of 4AL according to wheat reference genome sequence. Resequencing of E6015-3S and E6015-4T and haplotype analysis of 3087 worldwide wheat accessions revealed heightened deletion polymorphisms in the distal 0.949 Mbp region of 4AL, which was confirmed by the finding of frequent gene losses in this region in eight genome-sequenced hexaploid wheat cultivars. The great majority (86.36%) of the 3087 lines displayed different degrees of nucleotide sequence deletions, with only 13.64% of them resembling E6015-4T in this region. These deletions can impair the presence and/or function of TaHST1 and surrounding genes, thus rendering global wheat germplasm vulnerable to HS or other environmental adversities. Therefore, conscientious and urgent efforts are needed in global wheat breeding programmes to optimize the structure and function of 4AL distal terminus by ensuring the presence of TaHST1 and surrounding genes. The new information reported here will help to accelerate the ongoing global efforts in improving wheat HST.

Metabolite Analysis of Jerusalem Artichoke (Helianthus tuberosus L.) Seedlings in Response to Polyethylene Glycol-Simulated Drought Stress

Jerusalem artichokes are a perennial crop with high drought tolerance and high value as a raw material to produce biofuels, functional feed, and food. However, there are few comprehensive metabolomic studies on Jerusalem artichokes under drought conditions. Methods: Ultra-performance liquid chromatography and tandem mass spectrometry were used to identify differential metabolites in Jerusalem artichoke seedling leaves under polyethylene glycol (PEG) 6000-simulated drought stress at 0, 18, 24, and 36 h. Results: A total of 661 metabolites and 236 differential metabolites were identified at 0 vs. 18, 18 vs. 24, and 24 vs. 36 h. 146 differential metabolites and 56 common were identified and at 0 vs. 18, 24, and 36 h. Kyoto Encyclopedia of Genes and Genomes enrichment identified 236 differential metabolites involved in the biosynthesis of secondary metabolites and amino acids. Metabolites involved in glycolysis, phenolic metabolism, tricarboxylic cycle, glutamate-mediated proline biosynthesis, urea cycle, amino acid metabolism, unsaturated fatty acid biosynthesis, and the met salvage pathway responded to drought stress. Conclusion: A metabolic network in the leaves of Jerusalem artichokes under drought stress is proposed. These results will improve understanding of the metabolite response to drought stress in Jerusalem artichokes and develop a foundation for breeding drought-resistant varieties.

Phenotypic Characteristics and Transcriptome of Cucumber Male Flower Development Under Heat Stress

Cucumber (Cucumis sativus L.) is an important vegetable crop, which is thermophilic not heat resistant. High-temperature stress always results in sterility at reproductive stage. In the present study, we evaluate the male flower developmental changes under normal (CK) and heat stress (HS) condition. After HS, the activities of peroxidase (POD) and superoxide dismutase (SOD) and the contents of malondialdehyde (MDA) were increased. In addition, the pollen fertility was significantly decreased; and abnormal tapetum and microspore were observed by paraffin section. Transcriptome analysis results presented that total of 5828 differentially expressed genes (DEGs) were identified after HS. Among these DEGs, 20 DEGs were found at four stages, including DNA binding transcription factor, glycosyltransferase, and wound-responsive family protein. The gene ontology term of carbohydrate metabolic process was significantly enriched in all anther stages, and many saccharides and starch synthase-related genes, such as invertase, sucrose synthase, and starch branching enzyme, were significantly different expressed in HS compared with CK. Furthermore, co-expression network analysis showed a module (midnightblue) strongly consistent with HS, and two hub genes (CsaV3_6G004180 and CsaV3_5G034860) were found with a high degree of connectivity to other genes. Our results provide comprehensive understandings on male flower development in cucumber under HS.

Alternative Splicing Diversified the Heat Response and Evolutionary Strategy of Conserved Heat Shock Protein 90s in Hexaploid Wheat (Triticum aestivum L.)

Crops are challenged by the increasing high temperature. Heat shock protein 90 (HSP90), a molecular chaperone, plays a critical role in the heat response in plants. However, the evolutionary conservation and divergence of HSP90s homeologs in polyploidy crops are largely unknown. Using the newly released hexaploid wheat reference sequence, we identified 18 TaHSP90s that are evenly distributed as homeologous genes among three wheat subgenomes, and were highly conserved in terms of sequence identity and gene structure among homeologs. Intensive time-course transcriptomes showed uniform expression and transcriptional response profiles among the three TaHSP90 homeologs. Based on the comprehensive isoforms generated by combining full-length single-molecule sequencing and Illumina short read sequencing, 126 isoforms, including 90 newly identified isoforms of TaHSP90s, were identified, and each TaHSP90 generated one to three major isoforms. Intriguingly, the numbers and the splicing modes of the major isoforms generated by three TaHSP90 homeologs were obviously different. Furthermore, the quantified expression profiles of the major isoforms generated by three TaHSP90 homeologs are also distinctly varied, exhibiting differential alternative splicing (AS) responses of homeologs. Our results showed that the AS diversified the heat response of the conserved TaHSP90s and provided a new perspective for understanding about functional conservation and divergence of homologous genes.

Transcriptomic analysis of short-term heat stress response in Pinellia ternata provided novel insights into the improved thermotolerance by spermidine and melatonin

Heat stress has been a major environmental factor limiting the growth and development of Pinellia ternata which is an important Chinese traditional medicine. It has been reported that spermidine (SPD) and melatonin (MLT) play pivotal roles in modulating heat stress response (HSR). However, the roles of SPD and MLT in HSR of P. ternata, and the potential mechanism is still unknown. Here, exogenous SPD and MLT treatments alleviated heat-induced damages in P. ternata, which was supported by the increased chlorophyll content, OJIP curve, and relative water content, and the decreased malondialdehyde and electrolyte leakage. Then, RNA sequencing between CK (control) and Heat (1 h of heat treatment) was conducted to analyze how genes were in response to short-term heat stress in P. ternata. A total of 14,243 (7870 up- and 6373 down-regulated) unigenes were differentially expressed after 1 h of heat treatment. Bioinformatics analysis revealed heat-responsive genes mainly included heat shock proteins (HSPs), ribosomal proteins, ROS-scavenging enzymes, genes involved in calcium signaling, hormone signaling transduction, photosynthesis, pathogen resistance, and transcription factors such as heat stress transcription factors (HSFs), NACs, WRKYs, and bZIPs. Among them, PtABI5, PtNAC042, PtZIP17, PtSOD1, PtHSF30, PtHSFB2b, PtERF095, PtWRKY75, PtGST1, PtHSP23.2, PtHSP70, and PtLHC1 were significantly regulated by SPD or MLT treatment with same or different trends under heat stress condition, indicating that exogenous application of MLT and SPD might enhance heat tolerance in P. ternata through regulating these genes but may with different regulatory patterns. These findings contributed to the identification of potential genes involved in short-term HSR and the improved thermotolerance by MLT and SPD in P. ternata, which provided important clues for improving thermotolerance of P. ternata.

SMRT and Illumina RNA sequencing reveal novel insights into the heat stress response and crosstalk with leaf senescence in tall fescue

BACKGROUND

As a cool-season grass species, tall fescue (Festuca arundinacea) is challenged by increasing temperatures. Heat acclimation or activation of leaf senescence, are two main strategies when tall fescue is exposed to heat stress (HS). However, lacking a genome sequence, the complexity of hexaploidy nature, and the short read of second-generation sequencing hinder a comprehensive understanding of the mechanism. This study aims to characterize the molecular mechanism of heat adaptation and heat-induced senescence at transcriptional and post-transcriptional levels.

RESULTS

Transcriptome of heat-treated (1 h and 72 h) and senescent leaves of tall fescue were generated by combining single-molecular real-time and Illumina sequencing. In total, 4076; 6917, and 11,918 differentially expressed genes (DEGs) were induced by short- and long-term heat stress (HS), and senescence, respectively. Venn and bioinformatics analyses of DEGs showed that short-term HS strongly activated heat shock proteins (Hsps) and heat shock factors (Hsfs), as well as specifically activated FK506-binding proteins (FKBPs), calcium signaling genes, glutathione S-transferase genes, photosynthesis-related genes, and phytohormone signaling genes. By contrast, long-term HS shared most of DEGs with senescence, including the up-regulated chlorophyll catabolic genes, phytohormone synthesis/degradation genes, stress-related genes, and NACs, and the down-regulated photosynthesis-related genes, FKBPs, and catalases. Subsequently, transient overexpression in tobacco showed that FaHsfA2a (up-regulated specifically by short-term HS) reduced cell membrane damages caused by HS, but FaNAC029 and FaNAM-B1 (up-regulated by long-term HS and senescence) increased the damages. Besides, alternative splicing was widely observed in HS and senescence responsive genes, including Hsps, Hsfs, and phytohormone signaling/synthesis genes.

CONCLUSIONS

The short-term HS can stimulate gene responses and improve thermotolerance, but long-term HS is a damage and may accelerate leaf senescence. These results contribute to our understanding of the molecular mechanism underlying heat adaptation and heat-induced senescence.

Memory of 5-min heat stress in Arabidopsis thaliana

An ability of plants memorizing past heat exposure to modulate the expression of stress response transcripts during recovery is essential for efficient acquired thermotolerance. In this study, we demonstrated that expression of heat response transcripts spiked at 30 min or 1 h, but dramatically declined at 3 h during recoveries following exposure to 5-min heat stress in Arabidopsis. In contrast, expression of transcripts up-regulated by 45-min heat stress was sustained for 30 min or 1 h then declined during recovery. These results suggest that heat memory can be differently modulated depending on the duration of heat exposure, and indicate that plants can memorize even 5-min heat stress to regulate acclimatory responses during recovery. Later hypothesis can be supported by the finding that accumulation of heat response proteins was also modulated during recovery following 5-min heat stress. In addition, 5-min heat stress followed by 3 h recovery was efficient to activate acquired thermotolerance of plants, although spike of transcript expression was observed at 1 h during recovery. These results suggest that plants possess the ability to quickly memorize heat stress and reset cellular states during recovery to adapt to subsequent severe heat stress.

PacBio single molecule long-read sequencing provides insight into the complexity and diversity of the Pinctada fucata martensii transcriptome

BACKGROUND

The pearl oyster Pinctada fucata martensii is an economically valuable shellfish for seawater pearl production, and production of pearls depends on its growth. To date, the molecular mechanisms of the growth of this species remain poorly understood. The transcriptome sequencing has been considered to understanding of the complexity of mechanisms of the growth of P. f. martensii. The recently released genome sequences of P. f. martensii, as well as emerging Pacific Bioscience (PacBio) single-molecular sequencing technologies, provide an opportunity to thoroughly investigate these molecular mechanisms.

RESULTS

Herein, the full-length transcriptome was analysed by combining PacBio single-molecule long-read sequencing (PacBio sequencing) and Illumina sequencing. A total of 20.65 Gb of clean data were generated, including 574,561 circular consensus reads, among which 443,944 full-length non-chimeric (FLNC) sequences were identified. Through transcript clustering analysis of FLNC reads, 32,755 consensus isoforms were identified, including 32,095 high-quality consensus sequences. After removing redundant reads, 16,388 transcripts were obtained, and 641 fusion transcripts were derived by performing fusion transcript prediction of consensus sequences. Alternative splicing analysis of the 16,388 transcripts was performed after accounting for redundancy, and 9097 gene loci were detected, including 1607 new gene loci and 14,946 newly discovered transcripts. The original boundary of 11,235 genes on the chromosomes was corrected, 12,025 complete open reading frame sequences and 635 long non-coding RNAs (LncRNAs) were predicted, and functional annotation of 13,482 new transcripts was achieved. Two thousand three hundred eighteen alternative splicing events were detected. A total of 228 differentially expressed transcripts (DETs) were identified between the largest (L) and smallest (S) pearl oysters. Compared with the S, the L showed 99 and 129 significantly up-and down-regulated DETs, respectively. Six of these DETs were further confirmed by quantitative real-time RT-PCR (RT-qPCR) in independent experiment.

CONCLUSIONS

Our results significantly improve existing gene models and genome annotations, optimise the genome structure, and in-depth understanding of the complexity and diversity of the differential growth patterns of P. f. martensii.

Transcriptome analysis of heat stress and drought stress in pearl millet based on Pacbio full-length transcriptome sequencing

BACKGROUND

Heat and drought are serious threats for crop growth and development. As the sixth largest cereal crop in the world, pearl millet can not only be used for food and forage but also as a source of bioenergy. Pearl millet is highly tolerant to heat and drought. Given this, it is considered an ideal crop to study plant stress tolerance and can be used to identify heat-resistant genes.

RESULTS

In this study, we used Pacbio sequencing data as a reference sequence to analyze the Illumina data of pearl millet that had been subjected to heat and drought stress for 48 h. By summarizing previous studies, we found 26,299 new genes and 63,090 new transcripts, and the number of gene annotations increased by 20.18%. We identified 2792 transcription factors and 1223 transcriptional regulators. There were 318 TFs and 149 TRs differentially expressed under heat stress, and 315 TFs and 128 TRs were differentially expressed under drought stress. We used RNA sequencing to identify 6920 genes and 6484 genes differentially expressed under heat stress and drought stress, respectively.

CONCLUSIONS

Through Pacbio sequencing, we have identified more new genes and new transcripts. On the other hand, comparing the differentially expressed genes under heat tolerance with the DEGs under drought stress, we found that even in the same pathway, pearl millet responds with a different protein.

Isoform Sequencing Provides Insight Into Freezing Response of Common Wheat (Triticum aestivum L.)

The objective of the study is to reveal the freezing tolerance mechanisms of wheat by combining the emerging single-molecule real-time (SMRT) sequencing technology PacBio Sequel and Illumina sequencing. Commercial semiwinter wheat Zhoumai 18 was exposed to -6°C for 4 h at the four-leave stage. Leaves of the control group and freezing-treated group were used to perform cDNA library construction. PacBio SMRT sequencing yielded 51,570 high-quality isoforms from leaves of control sample of Zhoumai 18, encoded by 20,366 gene loci. In total, 73,695 transcript isoforms, corresponding to 23,039 genes, were identified from the freezing-treated leaves. Compared with transcripts from the International Wheat Genome Sequencing Consortium RefSeq v1.1, 57,667 novel isoforms were discovered, which were annotated 21,672 known gene loci, as well as 3,399 novel gene loci. Transcriptome characterization including alterative spliced events, alternative polydenylation sites, transcription factors, and fusion transcripts were also analyzed. Freezing-responsive genes and signals were uncovered and proved that the ICE-ERF-COR pathway and ABA signal transduction play a vital role in the freezing response of wheat. In this study, PacBio sequencing and Illumina sequencing were applied to investigate the freezing tolerance in common wheat, and the transcriptome results provide insights into the molecular regulation mechanisms under freezing treatment.

Varying Architecture of Heat Shock Elements Contributes to Distinct Magnitudes of Target Gene Expression and Diverged Biological Pathways in Heat Stress Response of Bread Wheat

The heat shock transcription factor (HSF) binds to cis-regulatory motifs known as heat shock elements (HSEs) to mediate the transcriptional response of HSF target genes. However, the HSF-HSEs interaction is not clearly understood. Using the newly released genome reference sequence of bread wheat, we identified 39,478 HSEs (95.6% of which were non-canonical HSEs) and collapsed them into 30,604 wheat genes, accounting for 27.6% wheat genes. Using the intensively heat-responsive transcriptomes of wheat, we demonstrated that canonical HSEs have a higher propensity to induce a response in the closest downstream genes than non-canonical HSEs. However, the response magnitude induced by non-canonical HSEs was comparable to that induced by canonical HSEs. Significantly, some non-canonical HSEs that contain mismatched nucleotides at specific positions within HSEs had a larger response magnitude than that of canonical HSEs. Consistently, most of the HSEs identified in the promoter regions of heat shock proteins were non-canonical HSEs, suggesting an important role for these non-canonical HSEs. Lastly, distinct diverged biological processes were observed between genes containing different HSE types, suggesting that sequence variation in HSEs plays a key role in the evolution of heat responses and adaptation. Our results provide a new perspective to understand the regulatory network underlying heat responses.

Full-Length Transcriptome Sequencing and Different Chemotype Expression Profile Analysis of Genes Related to Monoterpenoid Biosynthesis in Cinnamomum porrectum

Leaves of C. porrectum are rich in essential oils containing monoterpenes, sesquiterpenes and aromatic compounds, but the molecular mechanism of terpenoid biosynthesis in C. porrectum is still unclear. In this paper, the differences in the contents and compositions of terpenoids among three chemotypes were analyzed using gas chromatography mass spectrometry (GC/MS). Furthermore, the differential expression of gene transcripts in the leaf tissues of the three C. porrectum chemotypes were analyzed through a comparison of full-length transcriptomes and expression profiles. The essential oil of the three C. porrectum chemotypes leaves was mainly composed of monoterpenes. In the full-length transcriptome of C. porrectum, 104,062 transcripts with 306,337,921 total bp, an average length of 2944 bp, and an N50 length of 5449 bp, were obtained and 94025 transcripts were annotated. In the eucalyptol and linalool chemotype, the camphor and eucalyptol chemotype, and the camphor and linalool chemotype comparison groups, 21, 22 and 18 terpene synthase (TPS) unigenes were identified respectively. Three monoterpene synthase genes, CpTPS3, CpTPS5 and CpTPS9, were upregulated in the eucalyptol chemotype compared to the linalool chemotype and camphor chemotype. CpTPS1 was upregulated in the camphor chemotype compared to the linalool chemotype and the eucalyptol chemotype. CpTPS4 was upregulated in the linalool chemotype compared to the camphor chemotype and the eucalyptol chemotype. Different unigenes had different expression levels among the three chemotypes, but the unigene expression levels of the 2-C-methyl-D-erythritol 4phosphate (MEP) pathway were generally higher than those of the mevalonate acid (MVA) pathway. Quantitative reverse transcription PCR(qRT-PCR) further validated these expression levels. The present study provides new clues for the functional exploration of the terpenoid synthesis mechanism and key genes in different chemotypes of C. porrectum.

Integration of Transcriptomics and Metabolomics for Pepper (Capsicum annuum L.) in Response to Heat Stress

Heat stress (HS), caused by extremely high temperatures, is one of the most severe forms of abiotic stress in pepper. In the present study, we studied the transcriptome and metabolome of a heat-tolerant cultivar (17CL30) and a heat-sensitive cultivar (05S180) under HS. Briefly, we identified 5754 and 5756 differentially expressed genes (DEGs) in 17CL30 and 05S180, respectively. Moreover, we also identified 94 and 108 differentially accumulated metabolites (DAMs) in 17CL30 and 05S180, respectively. Interestingly, there were many common HS-responsive genes (approximately 30%) in both pepper cultivars, despite the expression patterns of these HS-responsive genes being different in both cultivars. Notably, the expression changes of the most common HS-responsive genes were typically much more significant in 17CL30, which might explain why 17CL30 was more heat tolerant. Similar results were also obtained from metabolome data, especially amino acids, organic acids, flavonoids, and sugars. The changes in numerous genes and metabolites emphasized the complex response mechanisms involved in HS in pepper. Collectively, our study suggested that the glutathione metabolic pathway played a critical role in pepper response to HS and the higher accumulation ability of related genes and metabolites might be one of the primary reasons contributing to the heat resistance.