De novo transcriptome sequencing and comprehensive analysis of the drought-responsive genes in the desert plant Cynanchum komarovii

Physiology and transcriptomics highlight the underlying mechanism of sunflower responses to drought stress and rehydration

Drought can adversely influence the crop growth and production. Accordingly, sunflowers have strong adaptability to drought; hence, we conducted analyses for sunflower seedlings with drought stress and rehydration drought acclimation through physiological measurements and transcriptomics. It showed that drought can cause the accumulation of ROS and enhance the activity of antioxidant enzymes and the content of osmolytes. After rehydration, the contents of ROS and MDA were significantly reduced concomitant with increased antioxidant activity and osmotic adjustment. Totally, 2,589 DEGs were identified among treatments. Functional enrichment analysis showed that DEGs were mainly involved in plant hormone signal transduction, MAPK signaling, and biosynthesis of secondary metabolites. Comparison between differentially spliced genes and DEGs indicated that bHLH025, NAC53, and SINAT3 may be pivotal genes involved in sunflower drought resistance. Our results not only highlight the underlying mechanism of drought stress and rehydration in sunflower but also provide a theoretical basis for crop genetic breeding.

Functional characterization of alkaline phosphatases involved alarm pheromone in the vetch aphid Megoura viciae

The alkaline phosphatases (ALPs) are highly promiscuous enzymes and have been extensively investigated in mammals for their medical significance, but their functional promiscuity is relatively poorly understood in insects. Here, we first identified four ALP genes (designated as MvALP1-4) in the vetch aphid Megoura viciae that contained one alkaline phosphatase site, three metal-binding sites, and varied other functional sites. Phylogenetic analysis, molecular docking and the spatiotemporal expression profiling of MvALP1-4 were very different, indicating a promiscuous functionality. We also found that MvALP4 involved the biosynthesis of aphid alarm pheromones (EβF) in vitro and in vivo. Finally, transcriptome analysis in the stimulated and unstimulated aphids supported the involvement of MvALPs in the biosynthesis of aphid alarm pheromones. Our study identified a multifunctional ALP involved terpene synthase enzyme activity in the aphid, which contributes to the understanding of the functional plasticity of ALPs in insects.

Study of the Metabolite Changes in Ganoderma lucidum under Pineapple Leaf Residue Stress via LC-MS/MS Coupled with a Non-Targeted Metabolomics Approach

The effects of fermentation metabolites of G. lucidum under different pineapple leaf residue additions were separated and identified using liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS). The mass spectra showed that the metabolites had good response values only in the positive ion mode, and 3019 metabolites with significant differences, mainly distributed in 95 metabolic pathways, were identified. The multivariate analyses, including the principal component analysis (PCA), orthogonal least squares discriminant analysis (OPLS-DA), and volcano plots (VP), revealed that the G. lucidum metabolites exhibited significant differences (p < 0.05) and were well clustered under various pineapple leaf residue additions, featuring 494–545 upregulated and 998–1043 downregulated metabolites. The differential metabolic pathway analysis proved that two metabolic pathways related to the biosynthesis of amino acids and ABC transporters were particularly significant under the addition of pineapple leaf residue, where amino acids such as histidine and lysine were upregulated in contrast to downregulated tyrosine, valine, L-alanine, and L-asparagine. These study results are considered instrumental in substantiating the application of pineapple leaf residue in the cultivation of G. lucidum and improving its utilization rate and added value.

Transcriptomic and proteomic analyses uncover the drought adaption landscape of Phoebe zhennan

BACKGROUND

Phoebe zhennan S.Lee (nanmu) is listed as a threatened tree species in China, whose growth and development, especially during the seedling stage, can be severely limited by drought. Previous studies on nanmu responses to drought stress involved physiological and biochemical analyses, while the molecular mechanisms remained unclear. Therefore, it is of great significance to carry out molecular biology research on the drought resistance of nanmu and reveal the genetic background and molecular regulation mechanism of nanmu drought resistance.

RESULTS

Drought stress enhanced the soluble sugar (SS), free proline(PRO), superoxide anion (O2·-), and hydrogen peroxide (H₂O₂) contents as well as the peroxidase (POD) and monodehydroascorbate reductase (MDHAR) activities of nanmu. However, glutathione S-transferase (GST) activity was sensitive to drought stress. Further transcriptomic and proteomic analyses revealed the abundant members of the differentially expressed genes(DEGs) and differentially expressed proteins(DEPs) that were related to phenylpropanoid and flavonoid biosynthesis, hormone biosynthesis and signal transduction, chlorophyll metabolism, photosynthesis, and oxidation-reduction reaction, which suggested their involvement in the drought response of nanmu. These enhanced the osmotic regulation, detoxification, and enzyme-induced and non-enzyme-induced antioxidant ability of nanmu. Moreover, 52% (447/867) of proteins that were up-regulated and 34% (307/892) down-regulated ones were attributed to the increase and decrease of transcription abundance. Transcript up (T^U) and protein up (P^U) groups had 447 overlaps, while transcript down (T^D) and protein down (P^D) groups had 307 overlaps, accounting for 54% of up and 35% of down-regulated proteins. The lack of overlap between DEGs and DEPs also suggested that post-transcriptional regulation has a critical role in nanmu response to drought.

CONCLUSIONS

Our research results provide significant insights into the regulatory mechanisms of drought stress in nanmu.

Transcriptional analysis reveals sodium nitroprusside affects alfalfa in response to PEG-induced osmotic stress at germination stage

Drought is one of the most common environmental factors that affect alfalfa germination and development. Nitric oxide (NO) could mediate stress tolerance in plants. The goal of this study was to determine exogenous NO donor-mediated drought adaption molecular mechanisms during the alfalfa germination stage. In this study, physiological and transcriptome analyses were performed on 7 days of the growth period seedlings by sodium nitroprusside (SNP) and polyethylene glycol (PEG) treatment. The results showed that SNP supplementation alleviated malondialdehyde accumulation, increased levels of proline and soluble sugars, and enhanced antioxidant enzyme activity under osmotic stress conditions. RNA-Seq experiments identified 5828 genes exhibiting differential expression in seedlings treated with PEG, SNP, or SNP+PEG relative to seedlings treated with distilled water. Of these DEGs, 3235 were upregulated, and 2593 were downregulated relative to the controls. Fifteen DEGs were amplified by qRT-PCR to verify the changes in expression determined by RNA-Seq, revealing that PIF3, glnA, PLCG1, and RP-S11e exhibited enhanced expression under the SNP+PEG treatment. SNP was found to modulate redox homeostasis-related genes such as GSTs, SOD2, GPX, and RBOH, and triggered calcium signaling transduction. It also induced some key genes relating to the abscisic acid, ethylene, and auxin signaling transduction in response to PEG stress. Conversely, genes associated with secondary metabolite biosynthesis and the metabolism of starch and sucrose during osmotic stress were downregulated by SNP. These results provide new insights into SNP-mediated drought adaption mechanisms at transcriptome-wide in alfalfa and reveal key drought tolerance pathways in this species.

Kurup S, Cheruth AJ, Subramaniam S, Hassan Tawfik N, M.A. Amiri K. Transcriptome Profiling of Haloxylon persicum (Bunge ex Boiss and Buhse) an Endangered Plant Species under PEG-Induced Drought Stress

Haloxylon persicum is an endangered western Asiatic desert plant species, which survives under extreme environmental conditions. In this study, we focused on transcriptome analysis of H. persicum to understand the molecular mechanisms associated with drought tolerance. Two different periods of polyethylene glycol (PEG)-induced drought stress (48 h and 72 h) were imposed on H. persicum under in vitro conditions, which resulted in 18 million reads, subsequently assembled by de novo method with more than 8000 transcripts in each treatment. The N50 values were 1437, 1467, and 1524 for the control sample, 48 h samples, and 72 h samples, respectively. The gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis resulted in enrichment of mitogen-activated protein kinase (MAPK) and plant hormone signal transduction pathways under PEG-induced drought conditions. The differential gene expression analysis (DGEs) revealed significant changes in the expression pattern between the control and the treated samples. The KEGG analysis resulted in mapping transcripts with 138 different pathways reported in plants. The differential expression of drought-responsive transcription factors depicts the possible signaling cascades involved in drought tolerance. The present study provides greater insight into the fundamental transcriptome reprogramming of desert plants under drought.

Transcriptomic response to cold of thermophilous medicinal plant Marsdenia tenacissima

Extracts from Marsdenia tenacissima, involving tenacissoside H, I and G, have been used as remedies of cancer, inflammation and asthma. Low temperature serves as one of the main factors constrain the planting expansion and quality of M. tenacissima, but its functional mechanism has been known scarcely for the lack of genomic information and transcriptional profile. Here we investigated the transcriptomic responses of M. tenacissima under cold stress to gain insight into the molecular mechanism of low temperature sensitivity. Total RNAs were collected from samples obtained at 4-time points (after 0, 3, 6 and 48 h cold treatments with 4 °C, respectively), then used for library construction and sequenced on the Illumina Hiseq™ 4000 platform. Passing quality assessments, 500794 transcripts, and 206137 unigenes were de novo assembly out in Trinity v2.4.0, holding contig N50 of 2566 bp and unigene mean length of 754 bp. 44.20% of assembled unigenes were annotated to the well-known public protein database on a basis of sequence similarity. Using statistical comparison of the fragments per kilo base of transcript per million reads mapped (FPKM) values between conditions, 6082 group-specific differentially expressed genes (DEGs) were identified and considered as cold-responsive genes, which contained copious transcription factors and active secondary metabolism. Among them, 43 unigenes were constantly up-regulated expression along with cold time, which mainly implicated in the biosynthesis of secondary metabolites, carbon metabolism, RNA and DNA metabolism. Conversely, 21 unigenes involved in photosynthesis, cell wall, protein degradation, and transporters were downregulated continually with cold timescale. Experimentally, MtEF1α was chosen as the best housekeeping gene. Functional enrichments found that damaging of cold stress on M. tenacissima may be ascribed to inability of photosynthesis, ribsome processing, flavonoid biosynthesis and terpenoids degradation. Correlation analysis between cold induced transcription factors and tenacissoside biosynthesis-related genes indicated that 3β-HSD significant positively correlated with bHLH51, and 4-MSO with NF-YB, GRAS3, Trihelix, FAR1, MYB60, MYBS1, bZIP43. Further promoter clone found MYB-binding site in the promoter of 4-MSO. In view of the reported cold tolerance of MYB60, it is recommended as a potential candidate suitable for future molecular design of exaptation cultivation with high bioactive constituents.

Identification of candidate genes involved in steroidal alkaloids biosynthesis in organ-specific transcriptomes of Veratrum nigrum L

Veratrum nigrum is protected plant of Melanthiaceae family, able to synthetize unique steroidal alkaloids important for pharmacy. Transcriptomes from leaves, stems and rhizomes of in vitro maintained V. nigrum plants were sequenced and annotated for genes and markers discovery. Sequencing of samples derived from the different organs resulted in a total of 108,511 contigs with a mean length of 596 bp. Transcripts derived from leaf and stalk were annotated at 28%, and 38% in Nr nucleotide database, respectively. The sequencing revealed 949 unigenes related with lipid metabolism, including 73 transcripts involved in steroids and genus-specific steroid alkaloids biosynthesis. Additionally, 3203 candidate SSRs markers we identified in unigenes with average density of one SSR locus every 6.2 kb sequence. Unraveling of biochemical machinery of the pathway responsible for steroidal alkaloids will open possibility to design and optimize biotechnological process. The transcriptomic data provide valuable resources for biochemical, molecular genetics, comparative transcriptomics, functional genomics, ecological and evolutionary studies of V. nigrum.

De novo transcriptome assembly and co-expression network analysis of Cynanchum thesioides: Identification of genes involved in resistance to drought stress

Cynanchum thesioides are upright, xerophytic shrubs that are widely distributed in arid and semi-arid areas of China, North Korea, Mongolia and Siberia. To date, little is known about the molecular mechanisms of drought resistance in C. thesioides. To better understand drought resistance, we used transcriptome analysis and Illumina sequencing technology on C. thesioides, to identify drought-responsive genes. Using de novo assembly 55,268 unigenes were identified from 207.58 Gb of clean data. Amongst these, 36,265 were annotated with gene descriptions, conserved domains, gene ontology terms and metabolic pathways. The sequencing results showed that genes that were differentially expressed (DEGs) under drought stress were enriched in pathways such as carbon metabolism, starch and sucrose metabolism, amino acid biosynthesis, phenylpropanoid biosynthesis and plant hormone signal transduction. Moreover, many functional genes were up-regulated under severe drought stress to enhance tolerance. Weighted gene co-expression network analysis showed that there were key hub genes related to drought stress. Hundreds of candidate genes were identified under severe drought stress, including transcriptional factors such as MYB, G2-like, ERF, C2H2, NAC, NF-X1, GRF, HD-ZIP, HB-other, HSF, C3H, GRAS, WRKY, bHLH and Trihelix. These data are a valuable resource for further investigation into the molecular mechanism for drought stress in C. thesioides and will facilitate exploration of drought resistance genes.

De novo assembly of Agave sisalana transcriptome in response to drought stress provides insight into the tolerance mechanisms

Agave, monocotyledonous succulent plants, is endemic to arid regions of North America, exhibiting exceptional tolerance to their xeric environments. They employ various strategies to overcome environmental constraints, such as crassulacean acid metabolism, wax depositions, and protective leaf morphology. Genomic resources of Agave species have received little attention irrespective of their cultural, economic and ecological importance, which so far prevented the understanding of the molecular bases underlying their adaptations to the arid environment. In this study, we aimed to elucidate molecular mechanism(s) using transcriptome sequencing of A. sisalana. A de novo approach was applied to assemble paired-end reads. The expression study unveiled 3,095 differentially expressed unigenes between well-irrigated and drought-stressed leaf samples. Gene ontology and KEGG analysis specified a significant number of abiotic stress responsive genes and pathways involved in processes like hormonal responses, antioxidant activity, response to stress stimuli, wax biosynthesis, and ROS metabolism. We also identified transcripts belonging to several families harboring important drought-responsive genes. Our study provides the first insight into the genomic structure of A. sisalana underlying adaptations to drought stress, thus providing diverse genetic resources for drought tolerance breeding research.

Plant Life in Extreme Environments: How Do You Improve Drought Tolerance?

Systems studies of drought stress in resurrection plants and other xerophytes are rapidly identifying a large number of genes, proteins and metabolites that respond to severe drought stress or desiccation. This has provided insight into drought resistance mechanisms, which allow xerophytes to persist under such extreme environmental conditions. Some of the mechanisms that ensure cellular protection during severe dehydration appear to be unique to desert species, while many other stress signaling pathways are in common with well-studied model and crop species. However, despite the identification of many desiccation inducible genes, there are few "gene-to-field" examples that have led to improved drought tolerance and yield stability derived from resurrection plants, and only few examples have emerged from model species. This has led to many critical reviews on the merit of the experimental approaches and the type of plants used to study drought resistance mechanisms. This article discusses the long-standing arguments between the ecophysiology and molecular biology communities, on how to "drought-proof" future crop varieties. It concludes that a more positive and inclusive dialogue between the different disciplines is needed, to allow us to move forward in a much more constructive way.

Transcriptomic studies reveal a key metabolic pathway contributing to a well-maintained photosynthetic system under drought stress in foxtail millet (Setaria italica L.)

Drought stress is one of the most important abiotic factors limiting crop productivity. A better understanding of the effects of drought on millet (Setaria italica L.) production, a model crop for studying drought tolerance, and the underlying molecular mechanisms responsible for drought stress responses is vital to improvement of agricultural production. In this study, we exposed the drought resistant F1 hybrid, M79, and its parental lines E1 and H1 to drought stress. Subsequent physiological analysis demonstrated that M79 showed higher photosynthetic energy conversion efficiency and drought tolerance than its parents. A transcriptomic study using leaves collected six days after drought treatment, when the soil water content was about ∼20%, identified 3066, 1895, and 2148 differentially expressed genes (DEGs) in M79, E1 and H1 compared to the respective untreated controls, respectively. Further analysis revealed 17 Gene Ontology (GO) enrichments and 14 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in M79, including photosystem II (PSII) oxygen-evolving complex, peroxidase (POD) activity, plant hormone signal transduction, and chlorophyll biosynthesis. Co-regulation analysis suggested that these DEGs in M79 contributed to the formation of a regulatory network involving multiple biological processes and pathways including photosynthesis, signal transduction, transcriptional regulation, redox regulation, hormonal signaling, and osmotic regulation. RNA-seq analysis also showed that some photosynthesis-related DEGs were highly expressed in M79 compared to its parental lines under drought stress. These results indicate that various molecular pathways, including photosynthesis, respond to drought stress in M79, and provide abundant molecular information for further analysis of the underlying mechanism responding to this stress.

De novo transcriptome sequencing and comprehensive analysis of the heat stress response genes in the basidiomycetes fungus Ganoderma lucidum

Ganoderma lucidum is a valuable basidiomycete with numerous pharmacological compounds, which is widely consumed throughout China. We previously found that the polysaccharide content of Ganoderma lucidum fruiting bodies could be significantly improved by 45.63% with treatment of 42 °C heat stress (HS) for 2 h. To further investigate genes involved in HS response and explore the mechanisms of HS regulating the carbohydrate metabolism in Ganoderma lucidum, high-throughput RNA-Seq was conducted to analyse the difference between control and heat-treated mycelia at transcriptome level. We sequenced six cDNA libraries with three from control group (mycelia cultivated at 28 °C) and three from heat-treated group (mycelia subjected to 42 °C for 2 h). A total of 99,899 transcripts were generated using Trinity method and 59,136 unigenes were annotated by seven public databases. Among them, 2790 genes were identified to be differential expressed genes (DEGs) under HS condition, which included 1991 up-regulated and 799 down-regulated. 176 DEGs were then manually classified into five main responsive-related categories according to their putative functions and possible metabolic pathways. These groups include stress resistance-related factors; protein assembly, transportation and degradation; signal transduction; carbohydrate metabolism and energy provision-related process; other related functions, suggesting that a series of metabolic pathways in Ganoderma lucidum are activated by HS and the response mechanism involves a complex molecular network which needs further study. Remarkably, 48 DEGs were found to regulate carbohydrate metabolism, both in carbohydrate hydrolysis for energy provision and polysaccharide synthesis. In summary, this comprehensive transcriptome analysis will provide enlarged resource for further investigation into the molecular mechanisms of basidiomycete under HS condition.

Transcriptome sequence analysis and mining of SSRs in Jhar Ber (Ziziphus nummularia (Burm.f.) Wight & Arn) under drought stress

Ziziphus nummularia (Burm.f.) Wight & Arn., a perennial shrub that thrives in the arid regions, is naturally tolerant to drought. However, there are limited studies on the genomics of drought tolerance in Ziziphus sp. In this study, RNA-sequencing of one month old seedlings treated with PEG 6000 was performed using Roche GS-FLX454 Titanium pyrosequencing. A total of 367,176 raw sequence reads were generated, and upon adapter trimming and quality filtration 351,872 reads were assembled de novo into 32,739 unigenes. Further characterization of the unigenes indicated that 73.25% had significant hits in the protein database. Kyoto encyclopedia of genes and genomes database (KEGG) identified 113 metabolic pathways from the obtained unigenes. A large number of drought-responsive genes were obtained and among them differential gene expression of 16 highly induced genes was validated by qRT-PCR analysis. To develop genic-markers, 3,425 simple sequence repeats (SSRs) were identified in 2,813 unigene sequences. The data generated shall serve as an important reservoir for the identification and characterization of drought stress responsive genes for development of drought tolerant crops.

Transcriptome sequencing of Paeonia suffruticosa 'Shima Nishiki' to identify differentially expressed genes mediating double-color formation

Paeonia suffruticosa 'Shima Nishiki' is one of extremely rare double-color cultivars in the world. It usually shows the two beautiful colors of red and white in the same flower, and this trait undoubtedly makes the flowers more charming for the ornamental market. However, few studies have been done to unravel the molecular mechanisms of double-color formation in P. suffruticosa 'Shima Nishiki'. In this study, we measured the anthocyanin composition and concentration, and sequenced the transcriptomes of the red and white petals. We found that the total content of Pg-based glycosides was at a significantly higher level in the red petals. Furthermore, we assembled and annotated 92,671 unigenes. Comparative analyses of the two transcriptomes showed 227 differentially expressed genes (DEGs), among which 57 were up-regulated, and 170 were down-regulated in the red petals. Subsequently, we identified 3 DEGs and the other 6 structural genes in the anthocyanin biosynthetic pathway including PsCHS, PsCHI, PsF3H, PsF3'H, PsDFR, PsANS, PsAOMT, PsMYB, and PsWD40. Among them, PsDFR and PsMYB expressed at a significantly higher level and showed positive correlations between their expression and anthocyanin concentration in the red petals. However, PsWD40 expressed at a significantly lower level and exhibited an inverse relationship in the red petals. Furthermore, we further confirmed the relative expression of the 9 candidate genes using quantitative real-time PCR. Based on the above results, we concluded that the significant differential expression of PsDFR, PsMYB and PsWD40 may play a key role in anthocyanin concentration in the red and white petals, thereby mediating double-color formation. These data will provide a valuable resource to better understand the molecular mechanisms of double-color formation of P. suffruticosa 'Shima Nishiki'.

Transcriptomic studies reveal a key metabolic pathway contributing to a well-maintained photosynthetic system under drought stress in foxtail millet (Setaria italica L.)

Drought stress is one of the most important abiotic factors limiting crop productivity. A better understanding of the effects of drought on millet (Setaria italica L.) production, a model crop for studying drought tolerance, and the underlying molecular mechanisms responsible for drought stress responses is vital to improvement of agricultural production. In this study, we exposed the drought resistant F₁ hybrid, M79, and its parental lines E1 and H1 to drought stress. Subsequent physiological analysis demonstrated that M79 showed higher photosynthetic energy conversion efficiency and drought tolerance than its parents. A transcriptomic study using leaves collected six days after drought treatment, when the soil water content was about ∼20%, identified 3066, 1895, and 2148 differentially expressed genes (DEGs) in M79, E1 and H1 compared to the respective untreated controls, respectively. Further analysis revealed 17 Gene Ontology (GO) enrichments and 14 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in M79, including photosystem II (PSII) oxygen-evolving complex, peroxidase (POD) activity, plant hormone signal transduction, and chlorophyll biosynthesis. Co-regulation analysis suggested that these DEGs in M79 contributed to the formation of a regulatory network involving multiple biological processes and pathways including photosynthesis, signal transduction, transcriptional regulation, redox regulation, hormonal signaling, and osmotic regulation. RNA-seq analysis also showed that some photosynthesis-related DEGs were highly expressed in M79 compared to its parental lines under drought stress. These results indicate that various molecular pathways, including photosynthesis, respond to drought stress in M79, and provide abundant molecular information for further analysis of the underlying mechanism responding to this stress.

A genome-wide transcriptome map of pistachio (Pistacia vera L.) provides novel insights into salinity-related genes and marker discovery

Background

Pistachio (Pistacia vera L.) is one of the most important commercial nut crops worldwide. It is a salt-tolerant and long-lived tree, with the largest cultivation area in Iran. Climate change and subsequent increased soil salt content have adversely affected the pistachio yield in recent years. However, the lack of genomic/global transcriptomic sequences on P. vera impedes comprehensive researches at the molecular level. Hence, whole transcriptome sequencing is required to gain insight into functional genes and pathways in response to salt stress.

Results

RNA sequencing of a pooled sample representing 24 different tissues of two pistachio cultivars with contrasting salinity tolerance under control and salt treatment by Illumina Hiseq 2000 platform resulted in 368,953,262 clean 100 bp paired-ends reads (90 Gb). Following creating several assemblies and assessing their quality from multiple perspectives, we found that using the annotation-based metrics together with the length-based parameters allows an improved assessment of the transcriptome assembly quality, compared to the solely use of the length-based parameters. The generated assembly by Trinity was adopted for functional annotation and subsequent analyses. In total, 29,119 contigs annotated against all of five public databases, including NR, UniProt, TAIR10, KOG and InterProScan. Among 279 KEGG pathways supported by our assembly, we further examined the pathways involved in the plant hormone biosynthesis and signaling as well as those to be contributed to secondary metabolite biosynthesis due to their importance under salinity stress. In total, 11,337 SSRs were also identified, which the most abundant being dinucleotide repeats. Besides, 13,097 transcripts as candidate stress-responsive genes were identified. Expression of some of these genes experimentally validated through quantitative real-time PCR (qRT-PCR) that further confirmed the accuracy of the assembly. From this analysis, the contrasting expression pattern of NCED3 and SOS1 genes were observed between salt-sensitive and salt-tolerant cultivars.

Conclusion

This study, as the first report on the whole transcriptome survey of P. vera, provides important resources and paves the way for functional and comparative genomic studies on this major tree to discover the salinity tolerance-related markers and stress response mechanisms for breeding of new pistachio cultivars with more salinity tolerance.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3989-7) contains supplementary material, which is available to authorized users.

De novo assembly and analysis of the Pugionium cornutum (L.) Gaertn. transcriptome and identification of genes involved in the drought response

Pugionium cornutum (L.) Gaertn. is a xerophytic plant species widely distributed in sandy and desert habitats in northwest China. However, the molecular mechanism of drought tolerance in P. cornutum has received little attention. At present, there is limited available transcriptome information for P. cormutum in public databases. Illumina sequencing was used to identify drought-responsive genes and to further characterize the molecular basis of drought tolerance in P. cornutum. In total, 51,385 unigenes with an average length of 825.32bp were obtained by de novo transcriptome assembly. Among these unigenes, 35,276 were annotated with gene descriptions, conserved domains, gene ontology terms, and metabolic pathways. In addition, the results showed that differentially expressed genes (DEGs) were mainly involved in photosynthesis, nitrogen metabolism, and plant hormone signal transduction pathways, notably ascorbate and aldarate metabolism, which could be an alternative pathway to enhance antioxidant capacity in P. cornutum in response to drought stress. These results provide an important clue about the effects of accumulation of ROS on ascorbic acid biosynthesis in P. cornutum. In addition, we found that transcription of most genes involved in ascorbic acid metabolism was altered under drought stress. Additionally, 93 drought-inducible transcription factor genes were identified in the DEGs under drought conditions; these included DREB, AP2/EREBP, B-2a, ERF2, MYB and Zinc finger family. The results of this study provide further insight into the molecular mechanisms of stress tolerance in P. cornutum, and also identify some attractive candidate genes and valuable information for improving drought stress tolerance in other species through genetic engineering.

ABNORMAL POLLEN VACUOLATION1 (APV1) is required for male fertility by contributing to anther cuticle and pollen exine formation in maize

Anther cuticle and pollen exine are the major protective barriers against various stresses. The proper functioning of genes expressed in the tapetum is vital for the development of pollen exine and anther cuticle. In this study, we report a tapetum-specific gene, Abnormal Pollen Vacuolation1 (APV1), in maize that affects anther cuticle and pollen exine formation. The apv1 mutant was completely male sterile. Its microspores were swollen, less vacuolated, with a flat and empty anther locule. In the mutant, the anther epidermal surface was smooth, shiny, and plate-shaped compared with the three-dimensional crowded ridges and randomly formed wax crystals on the epidermal surface of the wild-type. The wild-type mature pollen had elaborate exine patterning, whereas the apv1 pollen surface was smooth. Only a few unevenly distributed Ubisch bodies were formed on the apv1 mutant, leading to a more apparent inner surface. A significant reduction in the cutin monomers was observed in the mutant. APV1 encodes a member of the P450 subfamily, CYP703A2-Zm, which contains 530 amino acids. APV1 appeared to be widely expressed in the tapetum at the vacuolation stage, and its protein signal co-localized with the endoplasmic reticulum (ER) signal. RNA-Seq data revealed that most of the genes in the fatty acid metabolism pathway were differentially expressed in the apv1 mutant. Altogether, we suggest that APV1 functions in the fatty acid hydroxylation pathway which is involved in forming sporopollenin precursors and cutin monomers that are essential for the development of pollen exine and anther cuticle in maize.

De Novo Assembly, Annotation, and Characterization of Root Transcriptomes of Three Caladium Cultivars with a Focus on Necrotrophic Pathogen Resistance/Defense-Related Genes

Roots are vital to plant survival and crop yield, yet few efforts have been made to characterize the expressed genes in the roots of non-model plants (root transcriptomes). This study was conducted to sequence, assemble, annotate, and characterize the root transcriptomes of three caladium cultivars (Caladium × hortulanum) using RNA-Seq. The caladium cultivars used in this study have different levels of resistance to Pythiummyriotylum, the most damaging necrotrophic pathogen to caladium roots. Forty-six to 61 million clean reads were obtained for each caladium root transcriptome. De novo assembly of the reads resulted in approximately 130,000 unigenes. Based on bioinformatic analysis, 71,825 (52.3%) caladium unigenes were annotated for putative functions, 48,417 (67.4%) and 31,417 (72.7%) were assigned to Gene Ontology (GO) and Clusters of Orthologous Groups (COG), respectively, and 46,406 (64.6%) unigenes were assigned to 128 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. A total of 4518 distinct unigenes were observed only in Pythium-resistant "Candidum" roots, of which 98 seemed to be involved in disease resistance and defense responses. In addition, 28,837 simple sequence repeat sites and 44,628 single nucleotide polymorphism sites were identified among the three caladium cultivars. These root transcriptome data will be valuable for further genetic improvement of caladium and related aroids.

Integrated mRNA and microRNA analysis identifies genes and small miRNA molecules associated with transcriptional and post-transcriptional-level responses to both drought stress and re-watering treatment in tobacco

BACKGROUND

Drought stress is one of the most severe problem limited agricultural productivity worldwide. It has been reported that plants response to drought-stress by sophisticated mechanisms at both transcriptional and post-transcriptional levels. However, the precise molecular mechanisms governing the responses of tobacco leaves to drought stress and water status are not well understood. To identify genes and miRNAs involved in drought-stress responses in tobacco, we performed both mRNA and small RNA sequencing on tobacco leaf samples from the following three treatments: untreated-control (CL), drought stress (DL), and re-watering (WL).

RESULTS

In total, we identified 798 differentially expressed genes (DEGs) between the DL and CL (DL vs. CL) treatments and identified 571 DEGs between the WL and DL (WL vs. DL) treatments. Further analysis revealed 443 overlapping DEGs between the DL vs. CL and WL vs. DL comparisons, and, strikingly, all of these genes exhibited opposing expression trends between these two comparisons, strongly suggesting that these overlapping DEGs are somehow involved in the responses of tobacco leaves to drought stress. Functional annotation analysis showed significant up-regulation of genes annotated to be involved in responses to stimulus and stress, (e.g., late embryogenesis abundant proteins and heat-shock proteins) antioxidant defense (e.g., peroxidases and glutathione S-transferases), down regulation of genes related to the cell cycle pathway, and photosynthesis processes. We also found 69 and 56 transcription factors (TFs) among the DEGs in, respectively, the DL vs. CL and the WL vs. DL comparisons. In addition, small RNA sequencing revealed 63 known microRNAs (miRNA) from 32 families and 368 novel miRNA candidates in tobacco. We also found that five known miRNA families (miR398, miR390, miR162, miR166, and miR168) showed differential regulation under drought conditions. Analysis to identify negative correlations between the differentially expressed miRNAs (DEMs) and DEGs revealed 92 mRNA-miRNA interactions between CL and DL plants, and 32 mRNA-miRNA interactions between DL and WL plants.

CONCLUSIONS

This study provides a global view of the transcriptional and the post-transcriptional responses of tobacco under drought stress and re-watering conditions. Our results establish an empirical foundation that should prove valuable for further investigations into the molecular mechanisms through which tobacco, and plants more generally, respond to drought stress at multiple molecular genetic levels.

De novo assembly and transcriptome characterization of spruce dwarf mistletoe Arceuthobium sichuanense uncovers gene expression profiling associated with plant development

Background

The parasitic flowering plant dwarf mistletoe (Arceuthobium spp., Viscaceae) is one of the most destructive forest pests, posing a major threat to numerous conifer species worldwide. Arceuthobium sichuanense (spruce dwarf mistletoe, SDM) infects Qinghai spruce (Picea crassifolia) and causes severe damage to spruce forests in Northwest China. SDM is a Chinese native parasitic plant and acquires carbohydrates and mineral nutrition from its hosts. However, underlying molecular basis of the physiological development is largely unknown. Investigations of these physiological traits have been hampered by the lack of genomic resources for this species.

Results

In this study, to investigate the transcriptomic processes underlying physiological traits and development in SDM, we used RNA from four major tissues (i.e., shoots, flowers, fruits, and seeds) for de novo assembly and to annotate the transcriptome of this species. We uncovered the annotated transcriptome and performed whole genome expression profiling to uncover transcriptional dynamics during physiological development, and we identified key gene categories involved in the process of sexual development. The assembled SDM transcriptome reported in this work contains 331,347 assembled transcripts; 226,687 unigenes were functionally annotated by Gene Ontology analysis. RNA-Seq analysis using this reference transcriptome identified 22,641 differentially expressed genes from shoots, flowers, fruits, and seeds. These genes are enriched in processes including organic substance metabolism, cellular metabolism, biosynthesis, and cellular component. In addition, genes related to transport, transcription, hormone biosynthesis and signaling, carbohydrate metabolism, and photosynthesis were differentially expressed between tissues.

Conclusion

This work reveals tissue-specific gene expression patterns and pathways of SDM and implied to a difference between photosynthetic and non-photosynthetic tissues in plants. The data can potentially be used for future investigations on endophytic parasitism and SDM-spruce interaction, and it dramatically increases the available genomic resources for Arceuthobium and dwarf mistletoe communities. This preliminary study of the Arceuthobium transcriptome provides excellent opportunities for characterizing plant parasitic genes with unknown functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3127-y) contains supplementary material, which is available to authorized users.