Monitoring of gene expression profiles and identification of candidate genes involved in drought responses in Festuca mairei

Comparative transcriptome analysis of differentially expressed genes of Medicago falcata L. breeding lines response to saline-alkaline stress

BACKGROUND

Salt-alkali stress is an abiotic stress that inhibits crop growth and reduces yield. It significantly affects various physiological processes in plants, including photosynthesis, osmotic regulation, and antioxidant defense. However, studies on the transcriptional response mechanisms of Medicago falcata L. under salt-alkali stress are limited. In this study, RNA-seq technology was used to analyze differentially expressed genes (DEGs) in salt-alkali tolerant M.falcata breeding lines (LM18) and the salt-alkali sensitive Hulunbeier (HL) under salt-alkali stress. Furthermore, physiological indicators such as chlorophyll content, proline accumulation, and superoxide dismutase (SOD) activity were assessed to compare the responses of LM18 and HL to salt-alkali stress. By integrating transcriptomic and physiological analyses, this study provides new insights into the physiological and molecular regulatory mechanisms of M. falcata in response to salt-alkali stress.

RESULTS

The results showed that compared to the untreated controls, 10,289 and 2,478 DEGs were detected in LM18 and HL M.falcata seedlings, with 788 shared DEGs detected in both. GO functional analysis classified these DEGs into three categories: Biological Process, Cellular Components, and Molecular Functions, with significant enrichment in GO terms such as "response to osmotic stress", "intramolecular oxidoreductase activity" and "antioxidant activity". Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed the involvement of these DEGs in key metabolic pathways, including "Phenylpropanoid biosynthesis", "Plant hormone signal transduction", "Plant-pathogen interaction", "Isoflavonoid biosynthesis", "Circadian rhythm-plant" and "Photosynthesis-antenna proteins". Physiological indicators and membership function analysis confirmed that LM18 has greater salt-alkali tolerance than HL. Transcription factor analysis identified 42 transcription factor families, with the ERF family being the most abundant, followed by MYB-related, WRKY, bHLH, and MYB families. Weighted Gene Co-expression Network Analysis (WGCNA) showed that the MEturquoise module exhibited a significant positive correlation with salt-alkali stress and several physiological indicators. Module gene network analysis and GO enrichment revealed that MS.gene64536(MYBP), MS.gene76249(SRM1) and MS.gene049843 (MPK3) have functions related to "response to salt stress" and "positive regulation of response to salt stress", suggesting their key roles in salt-alkali tolerance in M.falcata. All three genes were upregulated in the salt-alkali tolerant LM18.

CONCLUSIONS

The GO terms and KEGG pathways significantly enriched in LM18 involved a significantly higher number of DEGs compared to HL, suggesting a more robust and effective mechanism in LM18. These findings highlight the robust molecular and physiological adaptations of LM18 in response to salt-alkali stress.

A comparison of shared patterns of differential gene expression and gene ontologies in response to water-stress in roots and leaves of four diverse genotypes of Lolium and Festuca spp. temperate pasture grasses

Ryegrasses (Lolium spp.) and fescues (Festuca spp.) are closely related and widely cultivated perennial forage grasses. As such, resilience in the face of abiotic stresses is an important component of their traits. We have compared patterns of differentially expressed genes (DEGs) in roots and leaves of two perennial ryegrass genotypes and a single genotype of each of a festulolium (predominantly Italian ryegrass) and meadow fescue with the onset of water stress, focussing on overall patterns of DEGs and gene ontology terms (GOs) shared by all four genotypes. Plants were established in a growing medium of vermiculite watered with nutrient solution. Leaf and root material were sampled at 35% (saturation) and, as the medium dried, at 15%, 5% and 1% estimated water contents (EWCs) and RNA extracted. Differential gene expression was evaluated comparing the EWC sampling points from RNAseq data using a combination of analysis methods. For all genotypes, the greatest numbers of DEGs were identified in the 35/1 and 5/1 comparisons in both leaves and roots. In total, 566 leaf and 643 root DEGs were common to all 4 genotypes, though a third of these leaf DEGs were not regulated in the same up/down direction in all 4 genotypes. For roots, the equivalent figure was 1% of the DEGs. GO terms shared by all four genotypes were often enriched by both up- and down-regulated DEGs in the leaf, whereas generally, only by either up- or down-regulated DEGs in the root. Overall, up-regulated leaf DEGs tended to be more genotype-specific than down-regulated leaf DEGs or root DEGs and were also associated with fewer GOs. On average, only 5–15% of the DEGs enriching common GO terms were shared by all 4 genotypes, suggesting considerable variation in DEGs between related genotypes in enacting similar biological processes.

Differential gene expression and gene ontologies associated with increasing water-stress in leaf and root transcriptomes of perennial ryegrass (Lolium perenne)

Perennial ryegrass (Lolium perenne) is a forage and amenity grass species widely cultivated in temperate regions worldwide. As such, perennial ryegrass populations are exposed to a range of environmental conditions and stresses on a seasonal basis and from year to year. One source of potential stress is limitation on water availability. The ability of these perennial grasses to be able to withstand and recover after periods of water limitation or drought can be a key component of grassland performance. Thus, we were interested in looking at changes in patterns of gene expression associated with increasing water stress. Clones of a single genotype of perennial ryegrass were grown under non-flowering growth room conditions in vermiculite supplemented with nutrient solution. Leaf and root tissue was sampled at 4 times in quadruplicate relating to estimated water contents of 35%, 15%, 5% and 1%. RNA was extracted and RNAseq used to generate transcriptome profiles at each sampling point. Transcriptomes were assembled using the published reference genome sequence and differential gene expression analysed using 3 different programmes, DESeq2, edgeR and limma (with the voom transformation), individually and in combination, deriving Early, Middle and Late stage comparisons. Identified differentially expressed genes were then associated with enriched GO terms using BLAST2GO. For the leaf, up-regulated differentially expressed genes were strongly associated with GO terms only during the Early stage and the majority of GO terms were associated with only down-regulated genes at the Middle or Late stages. For the roots, few differentially expressed genes were identified at either Early or Middle stages. Only one replicate at 1% estimated water content produced high quality data for the root, however, this indicated a high level of differential expression. Again the majority of enriched GO terms were associated with down-regulated genes. The performance of the different analysis programmes and the annotations associated with identified differentially expressed genes is discussed.

Modulations in primary and secondary metabolic pathways and adjustment in physiological behaviour of Withania somnifera under drought stress

In general medicinal plants grown under water limiting conditions show much higher concentrations of secondary metabolites in comparison to control plants. In the present study, Withania somnifera plants were subjected to water stress and data related to drought tolerance phenomenon was collected and a putative mechanistic concept considering growth responses, physiological behaviour, and metabolite content and gene expression aspects is presented. Drought induced metabolic and physiological responses as well as drastic decrease in CO₂ uptake due to stomatal limitations. As a result, the consumption of reduction equivalents (NADPH²⁺) for CO₂ assimilation via the calvin cycle declines significantly resulting in the generation of a large oxidative stress and an oversupply of antioxidant enzymes. Drought also results in the shifting of metabolic processes towards biosynthetic activities that consume reduction equivalents. Thus, biosynthesis of reduced compounds (isoprenoids, phenols and alkaloids) is enhanced. The dynamics of various metabolites have been discussed in the light of gene expression analysis of control and drought treated leaves. Gene encoding enzymes of pathways leading to glucose, fructose and fructan production, conversion of triose phosphates to hexoses and hexose phosphorylation were up-regulated in the drought stressed leaves. The down-regulated Calvin cycle genes were co-ordinately regulated with the down-regulation of chloroplast triosephosphate/phosphate translocator, cytoplasmic fructose-1,6-bisphosphate aldolase and fructose bisphosphatase. Expression of gene encoding Squalene Synthase (SQS) was highly upregulated under drought stress which is responsible for the diversion of carbon flux towards withanolides biosynthesis from isoprenoid pathway.

Identification of Genes Involved in the Responses of Tangor (C. reticulata × C. sinensis) to Drought Stress

Drought is the major abiotic stress with adverse effects on citrus, decreasing the agronomical yield and influencing the fruit quality. In this study, cDNA-amplified fragment length polymorphism (cDNA-AFLP) technique was used to investigate the transcriptional profile changes and identify drought-responsive genes in “Amakusa” tangor (C. reticulata × C. sinensis), a hybrid citrus sensitive to water stress. The 255 out of 6,245 transcript-derived fragments (TDFs) displayed altered expression patterns including (A) induction, (B) repression, (C) upregulation, and (D) downregulation. With BLAST search, the gene products of differentially expressed fragments (DEFs) could be classified into several categories: cellular processes, transcription, transport, metabolism, stress/stimuli response, and developmental processes. Downregulated genes were highly represented by photosynthesis and basic metabolism, while upregulated ones were enriched in genes that were involved in transcription regulation, defense, energy, and transport. Present result also revealed some transient and up- and then downregulated genes such as aquaporin protein and photosystem enzyme. Expression patterns of 17 TDFs among 18 homologous to function-known genes were confirmed by qRT-PCR analysis. The present results revealed potential mechanism of drought tolerance in fruit crop and also provided candidate genes for future experiments in citrus.

Green leaf volatiles, fire and nonanoic acid activate MAPkinases in the model grass species Lolium temulentum

Background

Previously it has been shown that mechanical wounding, salinity and heat activated a 46 kDa and 44 kDa mitogen-activated protein kinases (MAPKs) in forage related grasses. Forage and turf related grasses are utilized in diverse environments where they are routinely subjected to herbicides and exposed to fire and volatiles after cutting, however very little is known concerning the perception or molecular responses to these different stresses or compounds.

Results

In the model grass species Lolium temulentum (Lt), a 46 kDa mitogen-activated protein kinase (MAPK) was activated in the leaves within 5 min and a 44 kDa MAPK 15 min after exposure to green leaf volatiles released from grass clippings. When the tips of leaves of Lt plants were scorched by fire, the 46 kDa MAPK and 44 kDa MAPK were rapidly activated within 5 min and 20 min respectively in the treated leaf, and 15 min systemically in an adjacent untreated tiller after exposure to fire. Nonanoic acid (pelargonic acid), a component in herbicides used on grasses, activated a 46 kDa MAPK in the treated leaves within 5 min of exposure and 15 min in systemic tissues. At concentrations normally used in the herbicides, nonanoic acid was found to only weakly activate the 44 kDa MAPK after an hour in treated leaves, but strongly activated it in the systemic tillers 30 min after treatment. Acetic acid, HCl and NaOH also were found to activate these MAPKs in treated tillers.

Conclusion

The rapid activation of these MAPKs to a wide range of stress stimuli, suggest that these MAPKs play a role in the perception and response to these stresses and compounds. The activation of the MAPK by green leaf volatiles indicates a role for these compounds in wound signaling in grasses.

Electronic supplementary material

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

Comparative analysis of latex transcriptome reveals putative molecular mechanisms underlying super productivity of Hevea brasiliensis

Increasing demand for natural rubber prompts studies into the mechanisms governing the productivity of rubber tree (Heveabrasiliensis). It is very interesting to notice that a rubber tree of clone PR107 in Yunnan, China is reported to yield more than 20 times higher than the average rubber tree. This super-high-yielding (SHY) rubber tree (designated as SY107), produced 4.12 kg of latex (cytoplasm of rubber producing laticifers, containing about 30% of rubber) per tapping, more than 7-fold higher than that of the control. This rubber tree is therefore a good material to study how the rubber production is regulated at a molecular aspect. A comprehensive cDNA-AFLP transcript profiling was performed on the latex of SY107 and its average counterparts by using the 384 selective primer pairs for two restriction enzyme combinations (ApoI/MseI and TaqI/MseI). A total of 746 differentially expressed (DE) transcript-derived fragments (TDFs) were identified, of which the expression patterns of 453 TDFs were further confirmed by RT-PCR. These RT-PCR confirmed TDFs represented 352 non-redundant genes, of which 215 had known or partially known functions and were grouped into 10 functional categories. The top three largest categories were transcription and protein synthesis (representing 24.7% of the total genes), defense and stress (15.3%), and primary and secondary metabolism (14.0%). Detailed analysis of the DE-genes suggests notable characteristics of SHY phenotype in improved sucrose loading capability, rubber biosynthesis-preferred sugar utilization, enhanced general metabolism and timely stress alleviation. However, the SHY phenotype has little correlation with rubber-biosynthesis pathway genes.

Evidence of current impact of climate change on life: a walk from genes to the biosphere

We review the evidence of how organisms and populations are currently responding to climate change through phenotypic plasticity, genotypic evolution, changes in distribution and, in some cases, local extinction. Organisms alter their gene expression and metabolism to increase the concentrations of several antistress compounds and to change their physiology, phenology, growth and reproduction in response to climate change. Rapid adaptation and microevolution occur at the population level. Together with these phenotypic and genotypic adaptations, the movement of organisms and the turnover of populations can lead to migration toward habitats with better conditions unless hindered by barriers. Both migration and local extinction of populations have occurred. However, many unknowns for all these processes remain. The roles of phenotypic plasticity and genotypic evolution and their possible trade-offs and links with population structure warrant further research. The application of omic techniques to ecological studies will greatly favor this research. It remains poorly understood how climate change will result in asymmetrical responses of species and how it will interact with other increasing global impacts, such as N eutrophication, changes in environmental N : P ratios and species invasion, among many others. The biogeochemical and biophysical feedbacks on climate of all these changes in vegetation are also poorly understood. We here review the evidence of responses to climate change and discuss the perspectives for increasing our knowledge of the interactions between climate change and life.

A silver-staining cDNA-AFLP protocol suitable for transcript profiling in the latex of Hevea brasiliensis (para rubber tree)

cDNA amplified fragment length polymorphism (cDNA-AFLP) is a powerful transcript-profiling tool widely used in diverse plant species. When applied to a new biological system, however, existing protocols usually require substantial modifications. Furthermore, the usage of radioactive isotope in typical protocols excludes their application in many labs. Latex, as the cytoplasm of rubber-producing cells sees a critical role in elucidating rubber biosynthesis and its regulation in rubber tree (Hevea brasiliensis). This paper describes a detailed step-by-step silver-staining cDNA-AFLP procedure, which is suitable to latex transcript profiling analysis. Theoretical analysis revealed that with the combination of two restriction enzyme pairs (ApoI/MseI and TaqI/MseI), approximately 94% of latex whole transcriptome could be visualized. After varying multiple parameters, including the amounts of primary and secondary template usage, pre-amplification cycle number and gel development, we obtained a high-quality silver-staining fingerprint. In the ApoI/MseI system, an average of 88.6 discernable bands (100-1,000 bp) was produced for each selective primer pair, and 97.2 bands for another system (TaqI/MseI). TaqI/MseI was the first pair of 4-bp cutters used in cDNA-AFLP analysis and proved to be efficient and reliable. The sensitivity and reliability of our method were further verified by an application example in detecting differential gene expression in the latex of Hevea tree.