De novo assembly of red clover transcriptome based on RNA-Seq data provides insight into drought response, gene discovery and marker identification

Multiomics Analyses Reveal MsC3H29 Positively Regulates Flavonoid Biosynthesis to Improve Drought Resistance of Autotetraploid Cultivated Alfalfa (Medicago sativa L.)

Alfalfa (Medicago sativa subsp. sativa), the "queen of forage," is the most important perennial legume, with high productivity and an excellent nutritional profile. Medicago sativa subsp. falcata is a subspecies of the alfalfa complex and exhibits better drought tolerance. However, drought stress significantly hampers their development and yield. The molecular mechanisms underlying the aboveground and underground tissues of sativa and falcata responding to drought stress remain obscure. Here, we performed a comprehensive comparative analysis of the physiological and transcriptomic responses of sativa and falcata under drought stress. The results showed that photosynthesis was inhibited, and antioxidant enzymes were activated under drought stress. MsC3H29, a CCCH-type zinc finger protein, was identified as a hub gene through weighted gene coexpression network analysis (WGCNA) and was significantly induced by drought in underground tissue. The MsC3H29 protein was localized in the nucleus. Overexpression (OE) of MsC3H29 can increase the primary root length and fresh weight of transgenic alfalfa hairy roots, while RNA interference (RNAi) decreases them under drought stress. The 2',7'-dichlorodihydrofluorescein diacetate (H2DCFDA) staining revealed that MsC3H29 promoted drought tolerance of alfalfa hairy roots through decreasing ROS accumulation. The targeted metabolome analysis showed that the overexpression of MsC3H29 resulted in higher levels of accumulation for flavonoid monomers, including vicenin, daidzein, apigenin, isorhamnetin, quercetin, and tricin, in transgenic alfalfa hairy roots before and after drought stress, while RNAi led to a reduction. Our study provided a key candidate gene for molecular breeding to improve drought resistance in alfalfa.

The MYB Transcription Factor GmMYB78 Negatively Regulates Phytophthora sojae Resistance in Soybean

Phytophthora root rot is a devastating disease of soybean caused by Phytophthora sojae. However, the resistance mechanism is not yet clear. Our previous studies have shown that GmAP2 enhances sensitivity to P. sojae in soybean, and GmMYB78 is downregulated in the transcriptome analysis of GmAP2-overexpressing transgenic hairy roots. Here, GmMYB78 was significantly induced by P. sojae in susceptible soybean, and the overexpressing of GmMYB78 enhanced sensitivity to the pathogen, while silencing GmMYB78 enhances resistance to P. sojae, indicating that GmMYB78 is a negative regulator of P. sojae. Moreover, the jasmonic acid (JA) content and JA synthesis gene GmAOS1 was highly upregulated in GmMYB78-silencing roots and highly downregulated in overexpressing ones, suggesting that GmMYB78 could respond to P. sojae through the JA signaling pathway. Furthermore, the expression of several pathogenesis-related genes was significantly lower in GmMYB78-overexpressing roots and higher in GmMYB78-silencing ones. Additionally, we screened and identified the upstream regulator GmbHLH122 and downstream target gene GmbZIP25 of GmMYB78. GmbHLH122 was highly induced by P. sojae and could inhibit GmMYB78 expression in resistant soybean, and GmMYB78 was highly expressed to activate downstream target gene GmbZIP25 transcription in susceptible soybean. In conclusion, our data reveal that GmMYB78 triggers soybean sensitivity to P. sojae by inhibiting the JA signaling pathway and the expression of pathogenesis-related genes or through the effects of the GmbHLH122-GmMYB78-GmbZIP25 cascade pathway.

Phenotypic characterization of drought responses in red clover (Trifolium pratense L.)

Introduction

Red clover (Trifolium pratense) is a protein-rich, short-lived perennial forage crop that can achieve high yields, but suffers increasingly from drought in different cultivation areas. Breeding for increased adaptation to drought is becoming essential, but at this stage it is unclear which traits breeders should target to phenotype responses to drought that allow them to identify the most promising red clover genotypes. In this study, we assessed how prolonged periods of drought affected plant growth in field conditions, and which traits could be used to distinguish better adapted plant material.

Methods

A diverse panel of 395 red clover accessions was evaluated during two growing seasons. We simulated 6-to-8-week drought periods during two consecutive summers, using mobile rain-out shelters, while an irrigated control field was established in an adjacent parcel. Plant growth was monitored throughout both growing seasons using multiple flights with a drone equipped with RGB and thermal sensors. At various observation moments throughout both growing seasons, we measured canopy cover (CC) and canopy height (CH). The crop water stress index (CWSI) was determined at two moments, during or shortly after the drought event.

Results

Manual and UAV-derived measurements for CH were well correlated, indicating that UAV-derived measurements can be reliably used in red clover. In both years, CC, CH and CWSI were affected by drought, with measurable growth reductions by the end of the drought periods, and during the recovery phase. We found that the end of the drought treatment and the recovery phase of approximately 20 days after drought were suitable periods to phenotype drought responses and to distinguish among genotypes.

Discussion

Multifactorial analysis of accession responses revealed interactions of the maturity type with drought responses, which suggests the presence of two independent strategies in red clover: 'drought tolerance' and 'drought recovery'. We further found that a large proportion of the accessions able to perform well under well-watered conditions were also the ones that were less affected by drought. The results of this investigation are interpreted in view of the development of breeding for adaptation to drought in red clover.

Two common, often coexisting grassland plant species differ in their evolutionary potential in response to experimental drought

For terrestrial plant communities, the increase in frequency and intensity of drought events is considered as one of the most severe consequences of climate change. While single-species studies demonstrate that drought can lead to relatively rapid adaptive genetic changes, the evolutionary potential and constraints to selection need to be assessed in comparative approaches to draw more general conclusions. In a greenhouse experiment, we compare the phenotypic response and evolutionary potential of two co-occurring grassland plant species, Bromus erectus and Trifolium pratense, in two environments differing in water availability. We quantified variation in functional traits and reproductive fitness in response to drought and compared multivariate genetic variance-covariance matrices and predicted evolutionary responses between species. Species showed different drought adaptation strategies, reflected in both their species-specific phenotypic plasticity and predicted responses to selection indicating contrasting evolutionary potential under drought. In T. pratense we found evidence for stronger genetic constraints under drought compared to more favourable conditions, and for some traits plastic and predicted evolutionary responses to drought had opposing directions, likely limiting the potential for adaptive change. Our study contributes to a more detailed understanding of the evolutionary potential of species with different adaptive strategies in response to climate change and may help to inform future scenarios for semi-natural grassland ecosystems.

Genome-wide discovery of di-nucleotide SSR markers based on whole genome re-sequencing data of Cicer arietinum L. and Cicer reticulatum Ladiz

Simple sequence repeats (SSRs) are valuable genetic markers due to their co-dominant inheritance, multi-allelic and reproducible nature. They have been largely used for exploiting genetic architecture of plant germplasms, phylogenetic analysis, and mapping studies. Among the SSRs, di-nucleotide repeats are the most frequent of the simple repeats distributed throughout the plant genomes. In present study, we aimed to discover and develop di-nucleotide SSR markers by using the whole genome re-sequencing (WGRS) data from Cicer arietinum L. and C. reticulatum Ladiz. A total of 35,329 InDels were obtained in C. arietinum, whereas 44,331 InDels in C. reticulatum. 3387 InDels with 2 bp length were detected in C. arietinum, there were 4704 in C. reticulatum. Among 8091 InDels, 58 di-nucleotide regions that were polymorphic between two species were selected and used for validation. We tested primers for evaluation of genetic diversity in 30 chickpea genotypes including C. arietinum, C. reticulatum, C. echinospermum P.H. Davis, C. anatolicum Alef., C. canariense A. Santos & G.P. Lewis, C. microphyllum Benth., C. multijugum Maesen, C. oxyodon Boiss. & Hohen. and C. songaricum Steph ex DC. A total of 244 alleles were obtained for 58 SSR markers giving an average of 2.36 alleles per locus. The observed heterozygosity was 0.08 while the expected heterozygosity was 0.345. Polymorphism information content was found to be 0.73 across all loci. Phylogenetic tree and principal coordinate analysis clearly divided the accessions into four groups. The SSR markers were also evaluated in 30 genotypes of a RIL population obtained from an interspecific cross between C. arietinum and C. reticulatum. Chi-square (χ²) test revealed an expected 1:1 segregation ratio in the population. These results demonstrated the success of SSR identification and marker development for chickpea with the use of WGRS data. The newly developed 58 SSR markers are expected to be useful for chickpea breeders.

A genome-wide association study of freezing tolerance in red clover (Trifolium pratense L.) germplasm of European origin

Improvement of persistency is an important breeding goal in red clover (Trifolium pratense L.). In areas with cold winters, lack of persistency is often due to poor winter survival, of which low freezing tolerance (FT) is an important component. We conducted a genome wide association study (GWAS) to identify loci associated with freezing tolerance in a collection of 393 red clover accessions, mostly of European origin, and performed analyses of linkage disequilibrium and inbreeding. Accessions were genotyped as pools of individuals using genotyping-by-sequencing (pool-GBS), generating both single nucleotide polymorphism (SNP) and haplotype allele frequency data at accession level. Linkage disequilibrium was determined as a squared partial correlation between the allele frequencies of pairs of SNPs and found to decay at extremely short distances (< 1 kb). The level of inbreeding, inferred from the diagonal elements of a genomic relationship matrix, varied considerably between different groups of accessions, with the strongest inbreeding found among ecotypes from Iberia and Great Britain, and the least found among landraces. Considerable variation in FT was found, with LT50-values (temperature at which 50% of the plants are killed) ranging from -6.0°C to -11.5°C. SNP and haplotype-based GWAS identified eight and six loci significantly associated with FT (of which only one was shared), explaining 30% and 26% of the phenotypic variation, respectively. Ten of the loci were found within or at a short distance (<0.5 kb) from genes possibly involved in mechanisms affecting FT. These include a caffeoyl shikimate esterase, an inositol transporter, and other genes involved in signaling, transport, lignin synthesis and amino acid or carbohydrate metabolism. This study paves the way for a better understanding of the genetic control of FT and for the development of molecular tools for the improvement of this trait in red clover through genomics assisted breeding.

Genotyping by sequencing and a newly developed mRNA-GBS approach to link population genetic and transcriptome analyses reveal pattern differences between sites and treatments in red clover (Trifolium pratense L.)

The important worldwide forage crop red clover (Trifolium pratense L.) is widely cultivated as cattle feed and for soil improvement. Wild populations and landraces have great natural diversity that could be used to improve cultivated red clover. However, to date, there is still insufficient knowledge about the natural genetic and phenotypic diversity of the species. Here, we developed a low-cost complexity reduced mRNA analysis (mRNA-GBS) and compared the results with population genetic (GBS) and previously published mRNA-Seq data, to assess whether analysis of intraspecific variation within and between populations and transcriptome responses is possible simultaneously. The mRNA-GBS approach was successful. SNP analyses from the mRNA-GBS approach revealed comparable patterns to the GBS results, but due to site-specific multifactorial influences of environmental responses as well as conceptual and methodological limitations of mRNA-GBS, it was not possible to link transcriptome analyses with reduced complexity and sequencing depth to previously published greenhouse and field expression studies. Nevertheless, the use of short sequences upstream of the poly(A) tail of mRNA to reduce complexity are promising approaches that combine population genetics and expression profiling to analyze many individuals with trait differences simultaneously and cost-effectively, even in non-model species. Nevertheless, our study design across different regions in Germany was also challenging. The use of reduced complexity differential expression analyses most likely overlays site-specific patterns due to highly complex plant responses under natural conditions.

Genes Associated with Biological Nitrogen Fixation Efficiency Identified Using RNA Sequencing in Red Clover (Trifolium pratense L.)

Commonly studied in the context of legume-rhizobia symbiosis, biological nitrogen fixation (BNF) is a key component of the nitrogen cycle in nature. Despite its potential in plant breeding and many years of research, information is still lacking as to the regulation of hundreds of genes connected with plant-bacteria interaction, nodulation, and nitrogen fixation. Here, we compared root nodule transcriptomes of red clover (Trifolium pratense L.) genotypes with contrasting nitrogen fixation efficiency, and we found 491 differentially expressed genes (DEGs) between plants with high and low BNF efficiency. The annotation of genes expressed in nodules revealed more than 800 genes not yet experimentally confirmed. Among genes mediating nodule development, four nod-ule-specific cysteine-rich (NCR) peptides were confirmed in the nodule transcriptome. Gene duplication analyses revealed that genes originating from tandem and dispersed duplication are significantly over-represented among DEGs. Weighted correlation network analysis (WGCNA) organized expression profiles of the transcripts into 16 modules linked to the analyzed traits, such as nitrogen fixation efficiency or sample-specific modules. Overall, the results obtained broaden our knowledge about transcriptomic landscapes of red clover's root nodules and shift the phenotypic description of BNF efficiency on the level of gene expression in situ.

Combined analysis of transcriptome and metabolome reveals the molecular mechanism and candidate genes of Haloxylon drought tolerance

Haloxylon ammodendron and Haloxylon persicum, as typical desert plants, show strong drought tolerance and environmental adaptability. They are ideal model plants for studying the molecular mechanisms of drought tolerance. Transcriptomic and metabolomic analyses were performed to reveal the response mechanisms of H. ammodendron and H. persicum to a drought environment at the levels of transcription and physiological metabolism. The results showed that the morphological structures of H. ammodendron and H. persicum showed adaptability to drought stress. Under drought conditions, the peroxidase activity, abscisic acid content, auxin content, and gibberellin content of H. ammodendron increased, while the contents of proline and malondialdehyde decreased. The amino acid content of H. persicum was increased, while the contents of proline, malondialdehyde, auxin, and gibberellin were decreased. Under drought conditions, 12,233 and 17,953 differentially expressed genes (DEGs) were identified in H. ammodendron and H. persicum , respectively, including members of multiple transcription factor families such as FAR1, AP2/ERF, C2H2, bHLH, MYB, C2C2, and WRKY that were significantly up-regulated under drought stress. In the positive ion mode, 296 and 452 differential metabolites (DEMs) were identified in H. ammodendron and H. persicum, respectively; in the negative ion mode, 252 and 354 DEMs were identified, primarily in carbohydrate and lipid metabolism. A combined transcriptome and metabolome analysis showed that drought stress promoted the glycolysis/gluconeogenesis pathways of H. ammodendron and H. persicum and increased the expression of amino acid synthesis pathways, consistent with the physiological results. In addition, transcriptome and metabolome were jointly used to analyze the expression changes of the genes/metabolites of H. ammodendron and H. persicum that were associated with drought tolerance but were regulated differently in the two plants. This study identified drought-tolerance genes and metabolites in H. ammodendron and H. persicum and has provided new ideas for studying the drought stress response of Haloxylon.

Transcriptional Memory in Taraxacum mongolicum in Response to Long-Term Different Grazing Intensities

Grazing, as an important land use method in grassland, has a significant impact on the morphological and physiological traits of plants. However, little is known about how the molecular mechanism of plant responds to different grazing intensities. Here, we investigated the response of Taraxacum mongolicum to light grazing and heavy grazing intensities in comparison with a non-grazing control. Using de novo transcriptome assembly, T. mongolicum leaves were compared for the expression of the different genes under different grazing intensities in natural grassland. In total, 194,253 transcripts were de novo assembled and comprised in nine leaf tissues. Among them, 11,134 and 9058 genes were differentially expressed in light grazing and heavy grazing grassland separately, with 5867 genes that were identified as co-expression genes in two grazing treatments. The Nr, SwissProt, String, GO, KEGG, and COG analyses by BLASTx searches were performed to determine and further understand the biological functions of those differentially expressed genes (DEGs). Analysis of the expression patterns of 10 DEGs by quantitative real-time RT-PCR (qRT-PCR) confirmed the accuracy of the RNA-Seq results. Based on a comparative transcriptome analysis, the most significant transcriptomic changes that were observed under grazing intensity were related to plant hormone and signal transduction pathways, carbohydrate and secondary metabolism, and photosynthesis. In addition, heavy grazing resulted in a stronger transcriptomic response compared with light grazing through increasing the of the secondary metabolism- and photosynthesis-related genes. These changes in key pathways and related genes suggest that they may synergistically respond to grazing to increase the resilience and stress tolerance of T. mongolicum. Our findings provide important clues for improving grassland use and protection and understanding the molecular mechanisms of plant response to grazing.

Physiological and transcriptome analyses highlight multiple pathways involved in drought stress in Medicago falcata

Medicago falcata is one of the leguminous forage crops, which grows well in arid and semiarid region. To fully investigate the mechanism of drought resistance response in M. falcata, we challenged the M. falcata plants with 30% PEG-6000, and performed physiological and transcriptome analyses. It was found that, the activities of antioxidant enzymes (eg. SOD, POD, and CAT) and soluble sugar content were all increased in the PEG-treated group, as compared to the control group. Transcriptome results showed that a total of 706 genes were differentially expressed in the PEG-treated plants in comparison with the control. Gene enrichment analyses on differentially expressed genes revealed that a number of genes in various pathway were significantly enriched, including the phenylpropanoid biosynthesis (ko00940) and glycolysis/gluconeogenesis (ko00010), indicating the involvement of these key pathways in drought response. Furthermore, the expression levels of seven differentially expressed genes were verified to be involved in drought response in M. falcata by qPCR. Taken together, these results will provide valuable information related to drought response in M. falcata and lay a foundation for molecular studies and genetic breeding of legume crops in future research.

Progress of Genomics-Driven Approaches for Sustaining Underutilized Legume Crops in the Post-Genomic Era

Legume crops, belonging to the Fabaceae family, are of immense importance for sustaining global food security. Many legumes are profitable crops for smallholder farmers due to their unique ability to fix atmospheric nitrogen and their intrinsic ability to thrive on marginal land with minimum inputs and low cultivation costs. Recent progress in genomics shows promise for future genetic gains in major grain legumes. Still it remains limited in minor legumes/underutilized legumes, including adzuki bean, cluster bean, horse gram, lathyrus, red clover, urd bean, and winged bean. In the last decade, unprecedented progress in completing genome assemblies of various legume crops and resequencing efforts of large germplasm collections has helped to identify the underlying gene(s) for various traits of breeding importance for enhancing genetic gain and contributing to developing climate-resilient cultivars. This review discusses the progress of genomic resource development, including genome-wide molecular markers, key breakthroughs in genome sequencing, genetic linkage maps, and trait mapping for facilitating yield improvement in underutilized legumes. We focus on 1) the progress in genomic-assisted breeding, 2) the role of whole-genome resequencing, pangenomes for underpinning the novel genomic variants underlying trait gene(s), 3) how adaptive traits of wild underutilized legumes could be harnessed to develop climate-resilient cultivars, 4) the progress and status of functional genomics resources, deciphering the underlying trait candidate genes with putative function in underutilized legumes 5) and prospects of novel breeding technologies, such as speed breeding, genomic selection, and genome editing. We conclude the review by discussing the scope for genomic resources developed in underutilized legumes to enhance their production and play a critical role in achieving the "zero hunger" sustainable development goal by 2030 set by the United Nations.

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.

Non-Coding RNAs in Response to Drought Stress

Drought stress causes changes in the morphological, physiological, biochemical and molecular characteristics of plants. The response to drought in different plants may vary from avoidance, tolerance and escape to recovery from stress. This response is genetically programmed and regulated in a very complex yet synchronized manner. The crucial genetic regulations mediated by non-coding RNAs (ncRNAs) have emerged as game-changers in modulating the plant responses to drought and other abiotic stresses. The ncRNAs interact with their targets to form potentially subtle regulatory networks that control multiple genes to determine the overall response of plants. Many long and small drought-responsive ncRNAs have been identified and characterized in different plant varieties. The miRNA-based research is better documented, while lncRNA and transposon-derived RNAs are relatively new, and their cellular role is beginning to be understood. In this review, we have compiled the information on the categorization of non-coding RNAs based on their biogenesis and function. We also discuss the available literature on the role of long and small non-coding RNAs in mitigating drought stress in plants.

Time-course transcriptomic analysis of Petunia ×hybrida leaves under water deficit stress using RNA sequencing

Water deficit limits plant growth and development, resulting in quality loss of horticultural crops. However, there is limited information on gene regulation and signaling pathways related to water deficit stress response at multiple time points. The objective of this research was to investigate global gene expression patterns under water deficit stress to provide an insight into how petunia (Petunia ×hybrida 'Mitchell Diploid') responded in the process of stress. Nine-week-old petunias were irrigated daily or placed under water stress by withholding water. Stressed plants reduced stomatal conductance after five days of water deficit, indicating they perceived stress and initiated stress response mechanisms. To analyze transcriptomic changes at the early stage of water deficit, leaf tissue samples were collected 1, 3, and 5 days after water was withheld for RNA sequencing. Under water deficit stress, 154, 3611, and 980 genes were upregulated and 41, 2806, and 253 genes were downregulated on day 1, 3, and 5, respectively. Gene Ontology analysis revealed that redox homeostasis processes through sulfur and glutathione metabolism pathways, and hormone signal transduction, especially abscisic acid and ethylene, were enriched under water deficit stress. Thirty-four transcription factor families were identified, including members of AP2/ERF, NAC, MYB-related, C2H2, and bZIP families, and TFs in AP2/ERF family was the most abundant in petunia. Interestingly, only one member of GRFs was upregulated on day 1, while most of TFs were differentially expressed on day 3 and/or 5. The transcriptome data from this research will provide valuable molecular resources for understanding the early stages of water stress-responsive networks as well as engineering petunia with enhanced water stress tolerance.

Isoflavone levels, nodulation and gene expression profiles of a CRISPR/Cas9 deletion mutant in the isoflavone synthase gene of red clover

Isoflavones are not involved in rhizobial signaling in red clover, but likely play a role in defense in the rhizosphere. Red clover (Trifolium pratense) is a high-quality forage legume, well suited for grazing and hay production in the temperate regions of the world. Like many legumes, red clover produces a number of phenylpropanoid compounds including anthocyanidins, flavan-3-ols, flavanols, flavanones, flavones, and isoflavones. The study of isoflavone biosynthesis and accumulation in legumes has come into the forefront of biomedical and agricultural research due to potential for medicinal, antimicrobial, and environmental implications. CRISPR/Cas9 was used to knock out the function of a key enzyme in the biosynthesis of isoflavones, isoflavone synthase (IFS1). A hemizygous plant carrying a 9-bp deletion in the IFS1 gene was recovered and was intercrossed to obtain homozygous mutant plants. Levels of the isoflavones formononetin, biochanin A and genistein were significantly reduced in the mutant plants. Wild-type and mutant plants were inoculated with rhizobia to test the effect of the mutation on nodulation, but no significant differences were observed, suggesting that these isoflavones do not play important roles in nodulation. Gene expression profiling revealed an increase in expression of the upstream genes producing the precursors for IFS1, namely, phenylalanine ammonium lyase and chalcone synthase, but there were no significant differences in IFS1 gene expression or in the downstream genes in the production of specific isoflavones. Higher expression in genes involved in ethylene response was observed in the mutant plants. This response is normally associated with biotic stress, suggesting that the plants may have been responding to cues in the surrounding rhizosphere due to lower levels of isoflavones.

Remodeling of the cell wall as a drought-tolerance mechanism of a soybean genotype revealed by global gene expression analysis

Drought stress is major abiotic stress that affects soybean production. Therefore, it is widely desirable that soybean becomes more tolerant to stress. To provide insights into regulatory mechanisms of the stress response, we compared the global gene expression profiles from leaves of two soybean genotypes that display different responses to water-deficit (BR 16 and Embrapa 48, drought-sensitive and drought-tolerant, respectively). After the RNA-seq analysis, a total of 5335 down-regulated and 3170 up-regulated genes were identified in the BR16. On the other hand, the number of genes differentially expressed was markedly lower in the Embrapa 48, 355 up-regulated and 471 down-regulated genes. However, induction and expression of protein kinases and transcription factors indicated signaling cascades involved in the drought tolerance. Overall, the results suggest that the metabolism of pectin is differently modulated in response to drought stress and may play a role in the soybean defense mechanism against drought. This occurs via an increase of the cell wall plasticity and crosslink, which contributed to a higher hydraulic conductance (K _f) and relative water content (RWC%). The drought-tolerance mechanism of the Embrapa 48 genotype involves remodeling of the cell wall and increase of the hydraulic conductance to the maintenance of cell turgor and metabolic processes, resulting in the highest leaf RWC, photosynthetic rate (A), transpiration (E) and carboxylation (A/C _i). Thus, we concluded that the cell wall adjustment under drought is important for a more efficient water use which promoted a more active photosynthetic metabolism, maintaining higher plant growth under drought stress.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-021-00043-4.

Transcriptome analysis reveals major transcriptional changes during regrowth after mowing of red clover (Trifolium pratense)

BACKGROUND

Red clover (Trifolium pratense) is globally used as a fodder plant due its high nutritional value and soil improving qualities. In response to mowing, red clover exhibits specific morphological traits to compensate the loss of biomass. The morphological reaction is well described, but the underlying molecular mechanisms and its role for plants grown in the field are unclear.

RESULTS

Here, we characterize the global transcriptional response to mowing of red clover by comparing plants grown under greenhouse conditions with plants growing on agriculturally used fields. Unexpectedly, we found that biotic and abiotic stress related changes of plants grown in the field overlay their regrowth related transcriptional changes and characterized transcription related protein families involved in these processes. Further, we can show that gibberellins, among other phytohormones, also contribute to the developmental processes related to regrowth after biomass-loss.

CONCLUSIONS

Our findings show that massive biomass loss triggers less transcriptional changes in field grown plants than their struggle with biotic and abiotic stresses and that gibberellins also play a role in the developmental program related to regrowth after mowing in red clover. Our results provide first insights into the physiological and developmental processes of mowing on red clover and may serve as a base for red clover yield improvement.

Identification of Differentially Expressed Drought-Responsive Genes in Guar [Cyamopsis tetragonoloba (L.) Taub]

Drought remains one of the most serious environmental stresses because of the continuous reduction in soil moisture, which requires the improvement of crops with features such as drought tolerance. Guar [Cyamopsis tetragonoloba (L.) Taub], a forage and industrial crop, is a nonthirsty plant. However, the information on the transcriptome changes that occur under drought stress in guar is very limited; therefore, a gene expression analysis is necessary in this context. Here, we studied the differentially expressed genes (DEGs) in response to drought stress and their metabolic pathways. RNA-Seq via an expectation-maximization algorithm was used to estimate gene abundance. Subsequently, an Empirical Analysis of Digital Gene Expression Data in the R Bioconductor package was used to identify DEGs. Blast2GO, InterProScan, and the Kyoto Encyclopedia of Genes and Genomes were used to explore functional annotation, protein analysis, enzymes, and metabolic pathways. Transcription factors were identified using the PlantTFDB database. Our study identified 499 upregulated and 191 downregulated genes in response to drought stress. Of those, 32 upregulated and six downregulated genes were deemed as novel genes exclusive to guar. An aggregate of 137 protein families, 306 domains, 12 repeats, and two sites were upregulated. The proton-dependent oligopeptide transporter family and transferase, aquaporin transporter, calcium/calmodulin-dependent/calcium-dependent protein kinase, aspartic peptidase A1 family, UDP-glucuronosyl/UDP-glucosyltransferase, and major intrinsic protein were the most upregulated protein families. The upregulated unigenes were associated with 88 enzymes and 77 KEGG pathways. Finally, the MYB-related, MYB, and ERF transcription factor families were upregulated. These data may be useful for understanding the plant molecular response to drought stress.

Quick and efficient approach to develop genomic resources in orphan species: Application in Lavandula angustifolia

Next-Generation Sequencing (NGS) technologies, by reducing the cost and increasing the throughput of sequencing, have opened doors to generate genomic data in a range of previously poorly studied species. In this study, we propose a method for the rapid development of a large-scale molecular resources for orphan species. We studied as an example the true lavender (Lavandula angustifolia Mill.), a perennial sub-shrub plant native from the Mediterranean region and whose essential oil have numerous applications in cosmetics, pharmaceuticals, and alternative medicines. The heterozygous clone “Maillette” was used as a reference for DNA and RNA sequencing. We first built a reference Unigene, compound of coding sequences, thanks to de novo RNA-seq assembly. Then, we reconstructed the complete genes sequences (with introns and exons) using an Unigene-guided DNA-seq assembly approach. This aimed to maximize the possibilities of finding polymorphism between genetically close individuals despite the lack of a reference genome. Finally, we used these resources for SNP mining within a collection of 16 commercial lavender clones and tested the SNP within the scope of a genetic distance analysis. We obtained a cleaned reference of 8, 030 functionally in silico annotated genes. We found 359K polymorphic sites and observed a high SNP frequency (mean of 1 SNP per 90 bp) and a high level of heterozygosity (more than 60% of heterozygous SNP per genotype). On overall, we found similar genetic distances between pairs of clones, which is probably related to the out-crossing nature of the species and the restricted area of cultivation. The proposed method is transferable to other orphan species, requires little bioinformatics resources and can be realized within a year. This is also the first reported large-scale SNP development on Lavandula angustifolia. All the genomics resources developed herein are publicly available and provide a rich pool of molecular resources to explore and exploit lavender genetic diversity in breeding programs.

Molecular study of drought response in the Mediterranean conifer Pinus pinaster Ait.: Differential transcriptomic profiling reveals constitutive water deficit-independent drought tolerance mechanisms

Adaptation of long-living forest trees to respond to environmental changes is essential to secure their performance under adverse conditions. Water deficit is one of the most significant stress factors determining tree growth and survival. Maritime pine (Pinus pinaster Ait.), the main source of softwood in southwestern Europe, is subjected to recurrent drought periods which, according to climate change predictions for the years to come, will progressively increase in the Mediterranean region. The mechanisms regulating pine adaptive responses to environment are still largely unknown. The aim of this work was to go a step further in understanding the molecular mechanisms underlying maritime pine response to water stress and drought tolerance at the whole plant level. A global transcriptomic profiling of roots, stems, and needles was conducted to analyze the performance of siblings showing contrasted responses to water deficit from an ad hoc designed full-sib family. Although P. pinaster is considered a recalcitrant species for vegetative propagation in adult phase, the analysis was conducted using vegetatively propagated trees exposed to two treatments: well-watered and moderate water stress. The comparative analyses led us to identify organ-specific genes, constitutively expressed as well as differentially expressed when comparing control versus water stress conditions, in drought-sensitive and drought-tolerant genotypes. Different response strategies can point out, with tolerant individuals being pre-adapted for coping with drought by constitutively expressing stress-related genes that are detected only in latter stages on sensitive individuals subjected to drought.

Sex-Specific Transcriptome Sequencing of Zoea I Larvae and Identification of Sex-Linked Genes Using Bulked Segregant Analysis in Pacific White Shrimp Litopenaeus vannamei

In order to learn the sex determination and differentiation process in shrimp, we performed sex-specific transcriptome sequencing of the Zoea I larvae in Litopenaeus vannamei (L. vannamei) using a DNA/RNA co-extraction method. The sex-specific genes and the sex-linked SNPs were identified. De novo assembly of all the clean reads generated 41,270 unigenes with an average length of 1026 bp and an N50 of 1922 bp. A total of 60 differentially expressed genes (DEGs) between females and males at Zoea I stage were obtained, in which 41 unigenes had annotations. DEGs that might be related to sex development in L. vannamei were analyzed. Besides, paralogous genes of some known sex-related genes in model organisms such as sex-lethal (Sxl) showed no sex-biased difference, which indicated they may not play roles in sex development at Zoea I stage. Candidate sex-linked SNPs were validated in six populations of L. vannamei, and SNPs in Unigene0020898 and Unigene0020336 were fully linked to genders, suggesting the two genes located in the sex-determining region of L. vannamei. Our report provided the sex-specific gene expression in early developmental stage of L. vannamei, and found some sex-biased genes that may participate in the sex development process. We also detected two sex-linked genes which may be located in the sex determination region of L. vannamei through bulked segregant analysis, and the SNPs on these genes also provided a method for genetic sex identification at transcriptional level. Our findings will facilitate further researches on molecular mechanism of sex determination and differentiation in shrimps.

Insights into phylogeny, age and evolution of Allium (Amaryllidaceae) based on the whole plastome sequences

BACKGROUND AND AIMS

The genus Allium L., one of the largest monocotyledonous genera and one that includes many economically important crops with nutritional and medicinal value, has been the focus of classification or phylogeny studies for centuries. Recent studies suggested that the genus can be divided into 15 subgenera and 72 sections, which were further classified into three evolutionary lineages. However, the phylogenetic relationships reconstructed by one or two loci showed weaker support, especially for the third evolutionary lineage, which might not show the species relationships very clearly and could hinder further adaptive and evolutionary study.

METHODS

In this study, a total of 39 complete chloroplast genomes of Allium (covering 12 Allium subgenera) were collected, and combining these with 125 species of plastomes from 19 other families of monocots, we reconstructed the phylogeny of the genus Allium, estimated the origin and divergence time of the three evolutionary lineages and investigated the adaptive evolution in this genus and related families.

RESULTS

Our phylogenetic analysis confirmed the monophyly and three evolutionary lineages of Allium, while new species relationships were detected within the third evolutionary lineage. The divergence time of the three evolutionary lineages was estimated to be in the early Eocene to the middle Miocene, and numerous positive selected genes (PSGs) and PSGs with high average Ka/Ks values were found in Allium species.

CONCLUSIONS

Our results detected a well-supported phylogenetic relationship of Allium. The PSGs and PSGs with high Ka/Ks values, as well as diversified morphologies, complicated chromosome characteristics and unique reproductive modes may play important roles in the adaptation and evolution of Allium species. This is the first study that conducted phylogenetic and evolutionary analyses on the genus Allium combined with the plastome and morphological and cytological data. We hope that this study can contribute to further analysis of Allium for other researchers.

Phylogeny and highland adaptation of Chinese species in Allium section Daghestanica (Amaryllidaceae) revealed by transcriptome sequencing

Allium L. is one of the largest monocotyledonous genera with extensive distribution in the Northern Hemisphere. The fundamental phylogenies of Allium have been investigated using many morphological and molecular characters. However, the morphological characters may not agree with the molecular results in some Allium groups or sections (such as the Chinese Allium section Daghestanica), which may result in ambiguous species relationships and hinder further evolutionary and adaptive researches. Here, transcriptome sequences of the six Chinese endemics from Allium section Daghestanica were collected, with their single-copy genes (SCGs) were extracted. The interspecies relationships were analyzed using concatenation and coalescent methods. The branch-site model (BSM) was conducted to detect the positively selected genes (PSGs) in five highland species of this section. Based on 1644, 1281 and 1580 SCGs in flowers, leaves, and flowers-leaves combination respectively, a robust consistent and well-resolved phylogeny was generated from the concatenation method. Strong conflicts among individual gene trees were detected in the coalescent method, and morphological characters were incongruent with molecular relationships to some degree. Many PSGs were involved in responses of various stresses and stimuli (e.g. hypoxia, low temperature, aridity), DNA repair, metabolism, nutrient or energy intake, photosynthesis, and signal transduction. Our study revealed a clear interspecies relationship of Chinese endemics in Allium section Daghestanica and suggested that the discordance between morphological characters and molecular relationships might result from that the former are more susceptible to convergence compared with the latter. PSGs detected in our study may provide some insights into highland adaptation in Allium species.

De novo transcriptome characterisation of two auxin-related genes associated with plant growth habit in Astragalus adsurgens Pall

Astragalus adsurgens Pall., a perennial legume native to China, is commonly used as a forage crop. And it has great value for sustainable development of grasslands in arid and semi arid regions. However, to date, little is known regarding the A. adsurgens genome, and no studies have determined whether it would be possible to improve the germplasm of A. adsurgens through genetic modification. In this study, we used an RNA-seq protocol to generate a de novo transcriptome including 151,516 unigenes of A. adsurgens. We compared the transcriptomes of A. adsurgens having different growth habits (prostrate/erect) and identified 14,133 single nucleotide polymorphism sites (SNP) in 8,139 unigenes. Differential expression gene (DEG) analysis suggested that 10,982 unigenes were up-regulated in the prostrate plant relative to the erect plant, while 10,607 unigenes were down-regulated. Of the 21,589 DEG, Unigene72782_All (LAX4) and CL12494.Contig3_All (TIR1), an auxin transporter gene and an auxin transport inhibitor gene, respectively, were predicted to influence the growth habit of A. adsurgens, which were verified by qRT-PCR in these phenotypes. These results suggest that auxin transport was more active in the prostrate plant than in the erect plant, resulting in asymmetric distribution of auxin that affects the growth habit of A. adsurgens. Overall, this study may provide a basis for future research on key genes in A. adsurgens and may deepen our understanding of the molecular mechanisms regulating plant growth habit.

Dynamic Transcriptional Profiles of Arabidopsis thaliana Infected by Tomato spotted wilt virus

Tomato spotted wilt virus (TSWV) is a negative-stranded RNA virus that infects hundreds of plant species, causing great economic loss. Infected Arabidopsis thaliana plants develop symptoms including chlorosis and wilt, which can lead to cell death. From 9 to 15 days after TSWV infection, symptoms progress through a three-stage process of appearance, severity, and death. In this study, deep sequencing technology was first used to explore gene expression in response to TSWV infection in model plant A. thaliana at different symptom development stages. We found that plant immune defense and protein degradation are induced by TSWV infection and that both inductions became stronger over time. The photosynthesis pathway was attenuated with TSWV infection. Cell wall metabolism had a large extent of downregulation while some genes were upregulated. These results illustrate the dynamic nature of TSWV infection in A. thaliana at the whole-transcriptome level. The link between biological processes and subpathway metabolism was further analyzed. Our study provides new insight into host regulatory networks and dynamic processes in response to TSWV infection.

Seedling response to water stress in valley oak (Quercus lobata) is shaped by different gene networks across populations

Drought is a major stress for plants, creating a strong selection pressure for traits that enable plant growth and survival in dry environments. Many drought responses are conserved species-wide responses, while others vary among populations distributed across heterogeneous environments. We tested how six populations of the widely distributed California valley oak (Quercus lobata) sampled from contrasting climates would differ in their response to soil drying relative to well-watered controls in a common environment by measuring ecophysiological traits in 93 individuals and gene expression (RNA-seq) in 42 individuals. Populations did not differ in their adjustment of turgor loss point during soil drying, suggesting a generalized species-wide response. Differential expression analysis identified 689 genes with a common response to treatment across populations and 470 genes with population-specific responses. Weighted gene co-expression network analysis (WGCNA) identified groups of genes with similar expression patterns that may be regulated together (gene modules). Several gene modules responded differently to water stress among populations, suggesting regional differences in gene network regulation. Populations from sites with a high mean annual temperature responded to the imposed water stress with significantly greater changes in gene module expression, indicating that these populations may be locally adapted to respond to drought. We propose that this variation among valley oak populations provides a mechanism for differential tolerance to the increasingly frequent and severe droughts in California.

Identification of Single Nucleotide Polymorphism in Red Clover (Trifolium pratense L.) Using Targeted Genomic Amplicon Sequencing and RNA-seq

Red clover (Trifolium pratense L.) is a diploid, naturally cross-pollinated, cool-season species. As a perennial forage legume, red clover is mostly cultivated in temperate regions worldwide. Being a non-model crop species, genomic resources for red clover have been underdeveloped. Thus far, genomic analysis used in red clover has mainly relied on simple sequence repeat (SSR) markers. However, SSR markers are sparse in the genome and it is often difficult to unambiguously map them using short reads generated by next generation sequencing technology. Single nucleotide polymorphisms (SNPs) have been successfully applied in genomics assisted breeding in several agriculturally important species. Due to increasing importance of legumes in forage production, there is a clear need to develop SNP based markers for red clover that can be applied in breeding applications. In this study, we first developed an analytical pipeline that can confidently identify SNPs in a set of 72 different red clover genotypes using sequences generated by targeted amplicon sequencing. Then, with the same filtering stringency used in this pipeline, we used sequences from publicly available RNA-seq data to identify confident SNPs in different red clover varieties. Using this strategy, we have identified a total of 69,975 SNPs across red clover varieties. Among these, 28% (19,116) of them are missense mutations. Using Medicago truncatula as the reference, we annotated the regions affected by these missense mutations. We identified 2,909 protein coding regions with missense mutations. Pathway analysis of these coding regions indicated several biological processes impacted by these mutations. Specifically, three domains (homeobox domain, pentatricopeptide repeat containing plant-like, and regulator of Vps4 activity) were identified with five or more missense SNPs. These domain might also be a functional contributor in the molecular mechanisms of self-incompatibility in red clover. Future in-depth sequence diversity analysis of these three genes may yield valuable insights into the molecular mechanism involved in self-incompatibility in red clover.

RNAseq analysis reveals drought-responsive molecular pathways with candidate genes and putative molecular markers in root tissue of wheat

Drought is one of the major impediments in wheat productivity. Traditional breeding and marker assisted QTL introgression had limited success. Available wheat genomic and RNA-seq data can decipher novel drought tolerance mechanisms with putative candidate gene and marker discovery. Drought is first sensed by root tissue but limited information is available about how roots respond to drought stress. In this view, two contrasting genotypes, namely, NI5439 41 (drought tolerant) and WL711 (drought susceptible) were used to generate ~78.2 GB data for the responses of wheat roots to drought. A total of 45139 DEGs, 13820 TF, 288 miRNAs, 640 pathways and 435829 putative markers were obtained. Study reveals use of such data in QTL to QTN refinement by analysis on two model drought-responsive QTLs on chromosome 3B in wheat roots possessing 18 differentially regulated genes with 190 sequence variants (173 SNPs and 17 InDels). Gene regulatory networks showed 69 hub-genes integrating ABA dependent and independent pathways controlling sensing of drought, root growth, uptake regulation, purine metabolism, thiamine metabolism and antibiotics pathways, stomatal closure and senescence. Eleven SSR markers were validated in a panel of 18 diverse wheat varieties. For effective future use of findings, web genomic resources were developed. We report RNA-Seq approach on wheat roots describing the drought response mechanisms under field drought conditions along with genomic resources, warranted in endeavour of wheat productivity.

Comparative analysis of proso millet (Panicum miliaceum L.) leaf transcriptomes for insight into drought tolerance mechanisms

BACKGROUND

Drought stress is a major abiotic stress that causes huge losses in agricultural production. Proso millet (Panicum miliaceum L.) can efficiently adapt to drought stress and provides important information and gene resources to improve drought tolerance. However, its complex drought-responsive mechanisms remain unclear.

RESULTS

Among 37 core Chinese proso millet cultivars, Jinshu 6 (JS6) was selected as the drought-sensitive test material, whereas Neimi 5 (NM5) was selected as the drought-tolerant test material under PEG-induced water stress. After sequencing, 1695 differentially expressed genes (DEGs) were observed in JS6 and NM5 without PEG-induced water stress (JS6CK and NM5CK). A total of 833 and 2166 DEGs were found in the two cultivars under simulated drought by using 20% PEG-6000 for 6 (JS6T6 and NM5T6) and 24 h (JS6T24 and NM5T24), respectively. The DEGs in JS6T6 and JS6T24 treatments were approximately 0.298- and 0.754-fold higher than those in NM5T6 and NM5T24, respectively. Compared with the respective controls, more DEGs were found in T6 treatments than in T24 treatments. A delay in the transcriptional responses of the ROS scavenging system to simulated drought treatment and relatively easy recovery of the expression of photosynthesis-associated genes were observed in NM5. Compared with JS6, different regulation strategies were observed in the jasmonic acid (JA) signal transduction pathway of NM5.

CONCLUSION

Under PEG-induced water stress, NM5 maintained highly stable gene expression levels. Compared with drought-sensitive cultivars, the different regulation strategies in the JA signal transduction pathway in drought-tolerant cultivars may be one of the driving forces underlying drought stress tolerance.

Transcriptional profiling of contrasting genotypes revealed key candidates and nucleotide variations for drought dissection in Camellia sinensis (L.) O. Kuntze

Tea is popular health beverage consumed by millions of people worldwide. Drought is among the acute abiotic stress severely affecting tea cultivation, globally. In current study, transcriptome sequencing of four diverse tea genotypes with inherent contrasting genetic response to drought (tolerant & sensitive) generated more than 140 million reads. De novo and reference-based assembly and functional annotation of 67,093 transcripts with multifarious public protein databases yielded 54,484 (78.2%) transcripts with significant enrichment of GO and KEGG drought responsive pathways in tolerant genotypes. Comparative DGE and qRT analysis revealed key role of ABA dependent & independent pathways, potassium & ABC membrane transporters (AtABCG22, AtABCG11, AtABCC5 & AtABCC4) and antioxidant defence system against oxidative stress in tolerant genotypes, while seems to be failed in sensitive genotypes. Additionally, highly expressed UPL3HECT E3 ligases and RING E3 ligases possibly enhance drought tolerance by actively regulating functional modification of stress related genes. Further, ascertainment of, 80803 high quality putative SNPs with functional validation of key non-synonymous SNPs suggested their implications for developing high-throughput genotyping platform in tea. Futuristically, functionally relevant genomic resources can be potentially utilized for gene discovery, genetic engineering and marker-assisted genetic improvement for better yield and quality in tea under drought conditions.

Comparative analysis of two sister Erythrophleum species (Leguminosae) reveal contrasting transcriptome-wide responses to early drought stress

With the ongoing climate change, African rainforests are expected to experience severe drought events in the future. In Africa, the tropical genus Erythrophleum (Fabaceae) includes two forest sister timber tree species displaying contrasting geographical distributions. Erythrophleum ivorense is adapted to wet evergreen Guineo-Congolian forests, whereas E. suaveolens occurs in a wider range of climates, being found in moist dense forests but also in gallery forests under a relatively drier climate. This geographical distribution pattern suggests that the two species might cope differently to drought at the genomic level. Yet, the genetic basis of tolerance response to drought stress in both species is still uncharacterized. To bridge this gap, we performed an RNA-seq approach on seedlings from each species to monitor their transcriptional responses at different levels of drought stress (0, 2 and 6 weeks after stopping watering seedlings). Monitoring of wilting symptoms revealed that E. suaveolens displayed an earlier phenotypic response to drought stress than E. ivorense. At the transcriptomic level, results revealed 2020 (1204 down-regulated/816 up-regulated) and 1495 differentially expressed genes (DEGs) in response to drought stress from a total of 67,432 and 66,605 contigs assembled in E. ivorense and E. suaveolens, respectively. After identifying 30,374 orthologs between species, we found that only 7 of them were DEGs shared between species, while 587 and 458 were differentially expressed only in E. ivorense or E. suaveolens, respectively. GO and KEGG enrichment analysis revealed that the two species differ in terms of significantly regulated pathways as well as the number and expression profile of DEGs (Up/Down) associated with each pathway in the two stress stages. Our results suggested that the two studied species react differently to drought. E. suaveolens seems displaying a prompt response to drought at its early stage strengthened by the down-regulation of many DEGs encoding for signaling and metabolism-related pathways. A considerable up-regulation of these pathways was also found in E. ivorense at the late stage of drought, suggesting this species may be a late responder. Overall, our data may serve as basis for further understanding the genetic control of drought tolerance in tropical trees and favor the selection of crucial genes for genetically enhancing drought resistance.

Development of polymorphic EST-SSR markers and characterization of the autotetraploid genome of sainfoin (Onobrychis viciifolia)

Background

Sainfoin (Onobrychis viciifolia) is a highly nutritious, tannin-containing, and tetraploid forage legume. Due to the lack of detailed transcriptomic and genomic information on this species, genetic and breeding projects for sainfoin improvement have been significantly hindered.

Methods

In this study, a total of 24,630,711 clean reads were generated from 14 different sainfoin tissues using Illumina paired-end sequencing technology and deposited in the NCBI SRA database (SRX3763386). From these clean reads, 77,764 unigene sequences were obtained and 6,752 EST-SSRs were identified using de novo assembly. A total of 2,469 primer pairs were designed, and 200 primer pairs were randomly selected to analyze the polymorphism in five sainfoin wild accessions.

Results

Further analysis of 40 sainfoin individuals from the five wild populations using 61 EST-SSR loci showed that the number of alleles per locus ranged from 4 to 15, and the expected heterozygosity varied from 0.55 to 0.91. Additionally, by counting the EST-SSR band number and sequencing the three or four bands in one sainfoin individual, sainfoin was confirmed to be autotetraploid. This finding provides a high level of information about this plant.

Discussion

Through this study, 61 EST-SSR markers were successfully developed and shown to be useful for genetic studies and investigations of population genetic structures and variabilities among different sainfoin accessions.

Dicer functions transcriptionally and posttranscriptionally in a multilayer antiviral defense

In antiviral RNA interference (RNAi), Dicer plays a primary role in processing double-stranded RNA (dsRNA) molecules into small-interfering RNAs (siRNAs) that guide Argonaute effectors to posttranscriptional suppression of target viral genes. Here, we show a distinct role for Dicer in the siRNA-independent transcriptional induction of certain host genes upon viral infection in a filamentous fungus. Previous studies have shown that the two key players, dicer-like 2 (dcl2) and argonaute-like 2 (agl2), of antiviral RNAi in a phytopathogenic ascomycete, Cryphonectria parasitica, are highly transcriptionally induced upon infection with certain RNA mycoviruses, including the positive-stranded RNA hypovirus mutant lacking the RNAi suppressor (Cryphonectria hypovirus 1-Δp69, CHV1-Δp69). This induction is regulated by the Spt-Ada-Gcn5 acetyltransferase (SAGA) complex, a well-known transcriptional coactivator. The present study shows that diverse host genes, in addition to dcl2 and agl2, were up-regulated more than 10-fold by SAGA upon infection with CHV1-Δp69. Interestingly, DCL2, but not AGL2, was essential for SAGA-mediated global gene up-regulation. Moreover, deletion of certain virus-induced genes enhanced a CHV1-Δp69 symptom (growth rate) but not its accumulation. Constitutive, modest levels of dcl2 expression drastically reduced viral siRNA accumulation but were sufficient for full-scale up-regulation of host genes, suggesting that high induction of dcl2 and siRNA production are not essential for the transcriptional up-regulation function of DCL2. These data clearly demonstrate the dual functionality of DCL2: as a dsRNA-specific nuclease in posttranscriptional antiviral RNA silencing and as a key player in SAGA-mediated host gene induction, which independently represses viral replication and alleviates virus-induced symptom expression.

De Novo Assembly and Discovery of Genes That Involved in Drought Tolerance in the Common Vetch

The common vetch (Vicia sativa) is often used as feed for livestock because of its high nutritional value. However, drought stress reduces forage production through plant damage. Here, we studied the transcriptional profiles of common vetch exposed to drought in order to understand the molecular mechanisms of drought tolerance in this species. The genome of the common vetch has not been sequenced, therefore we used Illumina sequencing to generate de novo transcriptomes. Nearly 500 million clean reads were used to generate 174,636 transcripts, including 122,299 unigenes. In addition, 5313 transcription factors were identified and these transcription factors were classified into 79 different gene families. We also identified 11,181 SSR loci from di- to hexa-nucleotides whose repeat number was greater than five. On the basis of differentially expressed genes, Gene Ontology analysis identified many drought-relevant categories, including “oxidation-reduction process”, “lipid metabolic process” and “oxidoreductase activity”. In addition to these, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis identified pathways, such as “Plant hormone signal transduction”, “Glycolysis/Gluconeogenesis” and “Phenylpropanoid biosynthesis”, as differentially expressed in the plants exposed to drought. The expression results in this study will be useful for further extending our knowledge on the drought tolerance of common vetch.

Genome-wide identification of drought-responsive microRNAs in two sets of Malus from interspecific hybrid progenies

Drought stress can negatively impact apple fruit quality and yield. Apple microRNAs (miRNAs) participate in apple tree and fruit development, as well as in biotic stress tolerance; however, it is largely unknown whether these molecules are involved in the drought response. To identify drought-responsive miRNAs in Malus, we first examined the drought stress tolerance of ten F₁ progenies of R3 (M. × domestica) × M. sieversii. We performed Illumina sequencing on pooled total RNA from both drought-tolerant and drought-sensitive plants. The sequencing results identified a total of 206 known miRNAs and 253 candidate novel miRNAs from drought-tolerant plants and drought-sensitive plants under control or drought conditions. We identified 67 miRNAs that were differentially expressed in drought-tolerant plants compared with drought-sensitive plants under drought conditions. Under drought stress, 61 and 35 miRNAs were differentially expressed in drought-tolerant and drought-sensitive plants, respectively. We determined the expression levels of seven out of eight miRNAs by stem-loop qPCR analysis. We also predicted the target genes of all differentially expressed miRNAs and identified the expression of some genes. Gene Ontology analyses indicated that the target genes were mainly involved in stimulus response and cellular and metabolic processes. Finally, we confirmed roles of two miRNAs in apple response to mannitol. Our results reveal candidate miRNAs and their associated mRNAs that could be targeted for improving drought tolerance in Malus species, thus providing a foundation for understanding the molecular networks involved in the response of apple trees to drought stress.

Transcriptome Profiling and Characterization of Drought-Tolerant Potato Plant (Solanum tuberosum L.)

Potato (Solanum tuberosum L.) is the third most important food crop, and breeding drought-tolerant varieties is vital research goal. However, detailed molecular mechanisms in response to drought stress in potatoes are not well known. In this study, we developed EMS-mutagenized potatoes that showed significant tolerance to drought stress compared to the wild-type (WT) 'Desiree' cultivar. In addition, changes to transcripts as a result of drought stress in WT and drought-tolerant (DR) plants were investigated by de novo assembly using the Illumina platform. One-week-old WT and DR plants were treated with -1.8 Mpa polyethylene glycol-8000, and total RNA was prepared from plants harvested at 0, 6, 12, 24, and 48 h for subsequent RNA sequencing. In total, 61,100 transcripts and 5,118 differentially expressed genes (DEGs) displaying up- or down-regulation were identified in pairwise comparisons of WT and DR plants following drought conditions. Transcriptome profiling showed the number of DEGs with up-regulation and down-regulation at 909, 977, 1181, 1225 and 826 between WT and DR plants at 0, 6, 12, 24, and 48 h, respectively. Results of KEGG enrichment showed that the drought tolerance mechanism of the DR plant can mainly be explained by two aspects, the 'photosynthetic-antenna protein' and 'protein processing of the endoplasmic reticulum'. We also divided eight expression patterns in four pairwise comparisons of DR plants (DR0 vs DR6, DR12, DR24, DR48) under PEG treatment. Our comprehensive transcriptome data will further enhance our understanding of the mechanisms regulating drought tolerance in tetraploid potato cultivars.

Unigene-based RNA-seq provides insights on drought stress responses in Marsdenia tenacissima

Marsdenia tenacissima is a well-known anti-cancer medicinal plant used in traditional Chinese medicine, which often grows on the karst landform and the water conservation capacity of land is very poorly and drought occurrences frequently. We found M. tenacissima has strong drought resistance because of continuousdrought16 d, the leaves of M. tenacissima were fully curly and dying. But the leaves were fully almost recovering after re-watering 24h. The activity of SOD and POD were almost doubled under drought stress. The content of osmotic regulating substance proline and soluble sugar were three times than control group. But after re-watering, these indexes were declined rapidly. Three cDNA libraries of control, drought stress, and re-watering treatments were constructed. There were 43,129,228, 47,116,844, and 42,815,454 clean reads with Q20 values of 98.06, 98.04, and 97.88respectively.SRA accession number of raw data was PRJNA498187 on NCBI. A total of 8672, 6043, and 6537 differentially expressed genes (DEGs) were identified in control vs drought stress, control vs re-watering, and drought stress vs re-watering, respectively. In addition, 1039, 1016, and 980 transcription factors (TFs) were identified, respectively. Among them, 363, 267, and 299 TFs were identified as DEGs in drought stress, re-watering, and drought stress and re-watering, respectively. These differentially expressed TFs mainly belonged to the bHLH, bZIP, C2H2, ERF, MYB, MYB-related, and NAC families. A comparative analysis found that 1174 genes were up-regulated and 2344 were down-regulated under drought stress and this pattern was the opposite to that found after re-watering. Among the up-regulated genes, 64 genes were homologous to known functional genes that directly protect plants against drought stress. Furthermore, 44 protein kinases and 38 TFs with opposite expression patterns under drought stress and re-watering were identified, which are possibly candidate regulators for drought stress resistance in M. tenacissima. Our study is the first to characterize the M. tenacissima transcriptome in response to drought stress, and will serve as a useful resource for future studies on the functions of candidate protein kinases and TFs involved in M. tenacissima drought stress resistance.

Analysis of transcripts and splice isoforms in red clover (Trifolium pratense L.) by single-molecule long-read sequencing

BACKGROUND

Red clover (Trifolium pratense L.) is an important cool-season legume plant, which is the most widely planted forage legume after alfalfa. Although a draft genome sequence was published already, the sequences and completed structure of mRNA transcripts remain unclear, which limit further explore on red clover.

RESULTS

In this study, the red clover transcriptome was sequenced using single-molecule long-read sequencing to identify full-length splice isoforms, and 29,730 novel isoforms from known genes and 2194 novel isoforms from novel genes were identified. A total of 5492 alternative splicing events was identified and the majority of alter spliced events in red clover was corrected as intron retention. In addition, of the 15,229 genes detected by SMRT, 8719 including 186,517 transcripts have at least one poly(A) site. Furthermore, we identified 4333 long non-coding RNAs and 3762 fusion transcripts.

CONCLUSIONS

We analyzed full-length transcriptome of red clover with PacBio SMRT. Those new findings provided important information for improving red clover draft genome annotation and fully characterization of red clover transcriptome.

Legume genetic resources and transcriptome dynamics under abiotic stress conditions

Grain legumes are an important source of nutrition and income for billions of consumers and farmers around the world. However, the low productivity of new legume varieties, due to the limited genetic diversity available for legume breeding programmes and poor policymaker support, combined with an increasingly unpredictable global climate is resulting in a large gap between current yields and the increasing demand for legumes as food. Hence, there is a need for novel approaches to develop new high-yielding legume cultivars that are able to cope with a range of environmental stressors. Next-generation technologies are providing the tools that could enable the more rapid and cost-effective genomic and transcriptomic studies for most major crops, allowing the identification of key functional and regulatory genes involved in abiotic stress resistance. In this review, we provide an overview of the recent achievements regarding abiotic stress resistance in a wide range of legume crops and highlight the transcriptomic and miRNA approaches that have been used. In addition, we critically evaluate the availability and importance of legume genetic resources with desirable abiotic stress resistance traits.

Transcriptome analysis of leaf senescence in red clover (Trifolium pratense L.)

Red clover (Trifolium pratense L.) is an important cool-season legume plant, which is used as forage. Leaf senescence is a critical developmental process that negatively affects plant quality and yield. The regulatory mechanism of leaf senescence has been studied, and genes involved in leaf senescence have been cloned and characterized in many plants. However, those works mainly focused on model plants. Information about regulatory pathways and the genes involved in leaf senescence in red clover is very sparse. In this study, to better understand leaf senescence in red clover, transcriptome analysis of mature and senescent leaves was investigated using RNA-Seq. A total of about 35,067 genes were identified, and 481 genes were differentially expressed in mature and senescent leaves. Some identified differentially expressed genes showed similar expression patterns as those involved in leaf senescence in other species, such as Arabidopsis, Medicago truncatula and rice. Differentially expressed genes were confirmed by quantitative real-time PCR (qRT-PCR). Genes involved in signal transduction, transportation and metabolism of plant hormones, transcription factors and plant senescence were upregulated, while the downregulated genes were primarily involved in nutrient cycling, lipid/carbohydrate metabolism, hormone response and other processes. There were 64 differentially expressed transcription factor genes identified by RNA-Seq, including ERF, WRKY, bHLH, MYB and NAC. A total of 90 genes involved in biosynthesis, metabolism and transduction of plant hormones, including abscisic acid, jasmonic acid, cyokinin, brassinosteroid, salicylic acid and ethylene, were identified. Furthermore, 207 genes with direct roles in leaf senescence were demonstrated, such as senescence-associated genes. These genes were associated with senescence in other plants. Transcriptome analysis of mature and senescent leaves in red clover provides a large number of differentially expressed genes. Further analysis and identification of senescence-associated genes can provide new insight into the regulatory mechanisms of leaf development and senescence in legume plant and red clover.

High-throughput SSR marker development and its application in a centipedegrass (Eremochloa ophiuroides (Munro) Hack.) genetic diversity analysis

Centipedegrass (Eremochloa ophiuroides (Munro) Hack.) is a perennial, warm-season C4 grass species that shows great potential for use as a low-maintenance turfgrass species in tropical and subtropical regions. However, limited genetic and genomic information is available for this species, which has impeded systematic studies on the enhancement of its turf quality and resistance against biotic and abiotic stress. In this study, Illumina HiSeq high-throughput sequencing technology was performed to generate centipedegrass transcriptome sequences. A total of 352,513 assembled sequences were used to search for simple sequence repeat (SSR) loci, and 64,470 SSR loci were detected in 47,638 SSR containing sequences. The tri-nucleotides were the most frequent repeat motif, followed by di-nucleotides, tetra-nucleotides hexnucleotides, and pentanucleotides. A total of 48,061 primer pairs were successfully designed in the flanking sequences of the SSRs, and 100 sets of primers were randomly selected for the initial validation in four centipedegrass accessions. In total, 56 (56.0%) of the 100 primer pairs tested successfully amplified alleles from all four centipedegrass accessions, while 50 were identified as polymorphic markers and were then used to assess the level of genetic diversity among 43 centipedegrass core collections. The genetic diversity analysis exhibited that the number of alleles (Na) per locus ranged from 3 to 13, and the observed heterozygosity (Ho) ranged from 0.17 to 0.83. The polymorphism information content (PIC) value of the markers ranged from 0.15 to 0.78, and the genetic distances (coefficient Nei72) between the accessions varied from 0.07 to 0.48. The UPGMA-based dendrogram clustered all 43 core collections into two main groups and six subgroups, which further validated the effectiveness of these newly developed SSR markers. Hence, these newly developed SSR markers will be valuable and potentially useful for future genetic and genomic studies of E. ophiuroides.

Genome-wide atlas of alternative polyadenylation in the forage legume red clover

Studies on prevalence and significance of alternative polyadenylation (APA) in plants have been so far limited mostly to the model plants. Here, a genome-wide analysis of APA was carried out in different tissue types in the non-model forage legume red clover (Trifolium pratense L). A profile of poly(A) sites in different tissue types was generated using so-called 'poly(A)-tag sequencing' (PATseq) approach. Our analysis revealed tissue-wise dynamics of usage of poly(A) sites located at different genomic locations. We also identified poly(A) sites and underlying genes displaying APA in different tissues. Functional categories enriched in groups of genes manifesting APA between tissue types were determined. Analysis of spatial expression of genes encoding different poly(A) factors showed significant differential expression of genes encoding orthologs of FIP1(V) and PCFS4, suggesting that these two factors may play a role in regulating spatial APA in red clover. Our analysis also revealed a high degree of conservation in diverse plant species of APA events in mRNAs encoding two key polyadenylation factors, CPSF30 and FIP1(V). Together with our previously reported study of spatial gene expression in red clover, this study will provide a comprehensive account of transcriptome dynamics in this non-model forage legume.

Improving the Yield and Nutritional Quality of Forage Crops

Despite being some of the most important crops globally, there has been limited research on forages when compared with cereals, fruits, and vegetables. This review summarizes the literature highlighting the significance of forage crops, the current improvements and some of future directions for improving yield and nutritional quality. We make the point that the knowledge obtained from model plant and grain crops can be applied to forage crops. The timely development of genomics and bioinformatics together with genome editing techniques offer great scope to improve forage crops. Given the social, environmental and economic importance of forage across the globe and especially in poorer countries, this opportunity has enormous potential to improve food security and political stability.

Transcriptomic and proteomic feature of salt stress-regulated network in Jerusalem artichoke (Helianthus tuberosus L.) root based on de novo assembly sequencing analysis

Ribosome activation and sugar metabolic process mainly act on the regulation of salt tolerance in the bioenergy crop Helianthus tuberosus L. as dissected by integrated transcriptomic and proteomic analyses. Helianthus tuberosus L. is an important halophyte plant that can survive in saline-alkali soil. It is vitally necessary to build an available genomic resource to investigate the molecular mechanisms underlying salt tolerance in H. tuberosus. De novo assembly and annotation of transcriptomes were built for H. tuberosus using a HiSeq 4000 platform. 293,823 transcripts were identified and annotated into 190,567 unigenes. In addition, iTRAQ-labeled quantitative proteomics was carried out to detect global protein profiling as a response to salt stress. Comparative omics analysis showed that 5432 genes and 43 proteins were differentially expressed in H. tuberosus under salt stress, which were enriched in the following processes: carbohydrate metabolism, ribosome activation and translation, oxidation-reduction and ion binding. The reprogramming of transcript and protein works suggested that the induced activity of ribosome and sugar signaling may endue H. tuberosus with salt tolerance. With high-quality sequencing and annotation, the obtained transcriptomics and proteomics provide a robust genomic resource for dissecting the regulatory molecular mechanism of H. tuberosus in response to salt stress.

Mining and Development of Novel SSR Markers Using Next Generation Sequencing (NGS) Data in Plants

Microsatellites, or simple sequence repeats (SSRs), are one of the most informative and multi-purpose genetic markers exploited in plant functional genomics. However, the discovery of SSRs and development using traditional methods are laborious, time-consuming, and costly. Recently, the availability of high-throughput sequencing technologies has enabled researchers to identify a substantial number of microsatellites at less cost and effort than traditional approaches. Illumina is a noteworthy transcriptome sequencing technology that is currently used in SSR marker development. Although 454 pyrosequencing datasets can be used for SSR development, this type of sequencing is no longer supported. This review aims to present an overview of the next generation sequencing, with a focus on the efficient use of de novo transcriptome sequencing (RNA-Seq) and related tools for mining and development of microsatellites in plants.

Enriching Genomic Resources and Transcriptional Profile Analysis of Miscanthus sinensis under Drought Stress Based on RNA Sequencing

Miscanthus × giganteus is wildly cultivated as a potential biofuel feedstock around the world; however, the narrow genetic basis and sterile characteristics have become a limitation for its utilization. As a progenitor of M. × giganteus, M. sinensis is widely distributed around East Asia providing well abiotic stress tolerance. To enrich the M. sinensis genomic databases and resources, we sequenced and annotated the transcriptome of M. sinensis by using an Illumina HiSeq 2000 platform. Approximately 316 million high-quality trimmed reads were generated from 349 million raw reads, and a total of 114,747 unigenes were obtained after de novo assembly. Furthermore, 95,897 (83.57%) unigenes were annotated to at least one database including NR, Swiss-Prot, KEGG, COG, GO, and NT, supporting that the sequences obtained were annotated properly. Differentially expressed gene analysis indicates that drought stress 15 days could be a critical period for M. sinensis response to drought stress. The high-throughput transcriptome sequencing of M. sinensis under drought stress has greatly enriched the current genomic available resources. The comparison of DEGs under different periods of drought stress identified a wealth of candidate genes involved in drought tolerance regulatory networks, which will facilitate further genetic improvement and molecular studies of the M. sinensis.

Transcriptome Analysis of Two Species of Jute in Response to Polyethylene Glycol (PEG)- induced Drought Stress

Drought stress results in significant crop yield losses. Comparative transcriptome analysis between tolerant and sensitive species can provide insights into drought tolerance mechanisms in jute. We present a comprehensive study on drought tolerance in two jute species-a drought tolerant species (Corchorus olitorius L., GF) and a drought sensitive species (Corchorus capsularis L., YY). In total, 45,831 non-redundant unigenes with average sequence length of 1421 bp were identified. Higher numbers of differentially expressed genes (DEGs) were discovered in YY (794) than in GF (39), implying that YY was relatively more vulnerable or hyper-responsive to drought stress at the molecular level; the two main pathways, phenylpropanoid biosynthesis and peroxisome pathway, significantly involved in scavenging of reactive oxygen species (ROS) and 14 unigenes in the two pathways presented a significant differential expression in response to increase of superoxide. Our classification analysis showed that 1769 transcription factors can be grouped into 81 families and 948 protein kinases (PKs) into 122 families. In YY, we identified 34 TF DEGs from and 23 PK DEGs, including 19 receptor-like kinases (RLKs). Most of these RLKs were downregulated during drought stress, implying their role as negative regulators of the drought tolerance mechanism in jute.

Dehydration induced transcriptomic responses in two Tibetan hulless barley (Hordeum vulgare var. nudum) accessions distinguished by drought tolerance

Background

The harsh environment on the Qinghai-Tibetan Plateau gives Tibetan hulless barley (Hordeum vulgare var. nudum) great ability to resist adversities such as drought, salinity, and low temperature, and makes it a good subject for the analysis of drought tolerance mechanism. To elucidate the specific gene networks and pathways that contribute to its drought tolerance, and for identifying new candidate genes for breeding purposes, we performed a transcriptomic analysis using two accessions of Tibetan hulless barley, namely Z772 (drought-tolerant) and Z013 (drought-sensitive).

Results

There were more up-regulated genes of Z772 than Z013 under both mild (5439-VS-2604) and severe (7203-VS-3359) dehydration treatments. Under mild dehydration stress, the pathways exclusively enriched in drought-tolerance genotype Z772 included Protein processing in endoplasmic reticulum, tricarboxylic acid (TCA) cycle, Wax biosynthesis, and Spliceosome. Under severe dehydration stress, the pathways that were mainly enriched in Z772 included Carbon fixation in photosynthetic organisms, Pyruvate metabolism, Porphyrin and chlorophyll metabolism. The main differentially expressed genes (DEGs) in response to dehydration stress and genes whose expression was different between tolerant and sensitive genotypes were presented in this study, respectively. The candidate genes for drought tolerance were selected based on their expression patterns.

Conclusions

The RNA-Seq data obtained in this study provided an initial overview on global gene expression patterns and networks that related to dehydration shock in Tibetan hulless barley. Furthermore, these data provided pathways and a targeted set of candidate genes that might be essential for deep analyzing the molecular mechanisms of plant tolerance to drought stress.

Electronic supplementary material

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