Identification and Comparative Analysis of Differential Gene Expression in Soybean Leaf Tissue under Drought and Flooding Stress Revealed by RNA-Seq

Differential Regulation of Drought Responses in Two Phaseolus vulgaris Genotypes

Drought is probably the most harmful stress affecting common bean crops. Domestication, worldwide spread and local farming practices has entailed the development of a wide variety of common bean genotypes with different degrees of resistance to water stress. In this work, physiological and molecular responses to water stress have been compared in two common bean accessions, PHA-0683 and PMB-0220, previously identified as highly and moderately resistant to water stress, respectively. Our hypothesis was that only quantitative differences in the expression patterns of key genes should be found if molecular mechanisms regulating drought resistance are similar in the two accessions. However, results presented here indicate that the resistance to drought in PMB-0220 and PHA-0683 common bean accessions is regulated by different molecular mechanisms. Differential regulation of ABA synthesis and ABA signaling related genes among the two genotypes, and the control of the drought-induced senescence have a relevant contribution to the higher resistance level of PHA-0683 accession. Our results also suggest that expression patterns of key senescence-related transcription factors could be considered in the screening for drought resistance in common bean germplasm collections.

Review: Proteomic Techniques for the Development of Flood-Tolerant Soybean

Soybean, which is rich in protein and oil as well as phytochemicals, is cultivated in several climatic zones. However, its growth is markedly decreased by flooding stress, which is caused by climate change. Proteomic techniques were used for understanding the flood-response and -tolerant mechanisms in soybean. Subcellular proteomics has potential to elucidate localized cellular responses and investigate communications among subcellular components during plant growth and under stress stimuli. Furthermore, post-translational modifications play important roles in stress response and tolerance to flooding stress. Although many flood-response mechanisms have been reported, flood-tolerant mechanisms have not been fully clarified for soybean because of limitations in germplasm with flooding tolerance. This review provides an update on current biochemical and molecular networks involved in soybean tolerance against flooding stress, as well as recent developments in the area of functional genomics in terms of developing flood-tolerant soybeans. This work will expedite marker-assisted genetic enhancement studies in crops for developing high-yielding stress-tolerant lines or varieties under abiotic stress.

Integration of sRNA, degradome, transcriptome analysis and functional investigation reveals gma-miR398c negatively regulates drought tolerance via GmCSDs and GmCCS in transgenic Arabidopsis and soybean

BACKGROUND

Drought conditions adversely affect soybean growth, resulting in severe yield losses worldwide. Increasing experimental evidence indicates miRNAs are important post-transcriptional regulators of gene expression. However, the drought-responsive molecular mechanism underlying miRNA-mRNA interactions remains largely uncharacterized in soybean. Meanwhile, the miRNA-regulated drought response pathways based on multi-omics approaches remain elusive.

RESULTS

We combined sRNA, transcriptome and degradome sequencing to elucidate the complex regulatory mechanism mediating soybean drought resistance. One-thousand transcripts from 384 target genes of 365 miRNAs, which were enriched in the peroxisome, were validated by degradome-seq. An integrated analysis showed 42 miRNA-target pairs exhibited inversely related expression profiles. Among these pairs, a strong induction of gma-miR398c as a major gene negatively regulates multiple peroxisome-related genes (GmCSD1a/b, GmCSD2a/b/c and GmCCS). Meanwhile, we detected that alternative splicing of GmCSD1a/b might affect soybean drought tolerance by bypassing gma-miR398c regulation. Overexpressing gma-miR398c in Arabidopsis thaliana L. resulted in decreased percentage germination, increased leaf water loss, and reduced survival under water deficiency, which displayed sensitivity to drought during seed germination and seedling growth. Furthermore, overexpressing gma-miR398c in soybean decreased GmCSD1a/b, GmCSD2a/b/c and GmCCS expression, which weakened the ability to scavenge O2.-, resulting in increased relative electrolyte leakage and stomatal opening compared with knockout miR398c and wild-type soybean under drought conditions.

CONCLUSION

The study indicates that gma-miR398c negatively regulates soybean drought tolerance, and provides novel insights useful for breeding programs to improve drought resistance by CRISPR technology.

Effects of drought stress on physiological responses and gene expression changes in herbaceous peony (Paeonia lactiflora Pall.)

Herbaceous peony (Paeonia lactiflora Pall.) is known as the flower phase. This phase is somewhat resistant to drought, but long-term drought and severe water shortage will affect its normal growth and development. In this study, physiological indices and the transcriptome of P. lactiflora were determined to clarify its physiological responses and gene expression changes under drought stress. The results showed that under drought stress, soluble sugar content, peroxidase (POD), catalase (CAT) and ascorbate peroxidase (APX) activities, and chlorophyll, carotenoid and flavonoid contents were significantly increased, and soluble protein content, superoxide dismutase (SOD), glutathione reductase (GR), dehydroascorbate reductase (DHAR), monodehydroascorbate reductase (MDHAR), ascorbic acid (AsA) and glutathione (GSH) activity first increased and then decreased after day 14. Moreover, drought stress also significantly reduced chlorophyll content, photosynthesis and chlorophyll fluorescence parameters. Transcriptomic analysis revealed that compared with the Control, 10,747 differentially expressed genes (DEGs) were upregulated and 11,835 downregulated under drought stress. These DEGs were classified into three categories and 46 functional groups by GO function classification. The 3,179 DEGs were enriched into 128 pathways by KEGG pathway enrichment. The ROS system, chlorophyll degradation and photosynthetic capacity, as well as secondary pathways of biosynthesis and sugar metabolism were included. Additionally, relevant genes expressed in some metabolic pathways were discovered. These results provide a theoretical basis for understanding the responses of P. lactiflora to drought stress.

Genome-wide characterization and expression analysis of soybean trihelix gene family

Trihelix transcription factors play multiple roles in plant growth, development and various stress responses. In this study, we identified 71 trihelix family genes in the soybean genome. These trihelix genes were located at 19 out of 20 soybean chromosomes unevenly and were classified into six distinct subfamilies: GT-1, GT-2, GTγ, SIP1, SH4 and GTδ. The gene structure and conserved functional domain of these trihelix genes were similar in the same subfamily but diverged between different subfamilies. Thirteen segmental duplicated gene pairs were identified and all of them experienced a strong purifying selective pressure during evolution. Various stress-responsive cis-elements presented in the promoters of soybean trihelix genes, suggesting that the trihelix genes might respond to the environmental stresses in soybean. The expression analysis suggests that trihelix genes are involved in diverse functions during soybean development, flood or salinity tolerance, and plant immunity. Our results provide genomic information of the soybean trihelix genes and a basis for further characterizing their roles in response to environmental stresses.

High-Resolution Mapping in Two RIL Populations Refines Major "QTL Hotspot" Regions for Seed Size and Shape in Soybean (Glycine max L.)

Seed size and shape are important traits determining yield and quality in soybean. However, the genetic mechanism and genes underlying these traits remain largely unexplored. In this regard, this study used two related recombinant inbred line (RIL) populations (ZY and K3N) evaluated in multiple environments to identify main and epistatic-effect quantitative trait loci (QTLs) for six seed size and shape traits in soybean. A total of 88 and 48 QTLs were detected through composite interval mapping (CIM) and mixed-model-based composite interval mapping (MCIM), respectively, and 15 QTLs were common among both methods; two of them were major (R2 > 10%) and novel QTLs (viz., qSW-1-1ZN and qSLT-20-1K3N). Additionally, 51 and 27 QTLs were identified for the first time through CIM and MCIM methods, respectively. Colocalization of QTLs occurred in four major QTL hotspots/clusters, viz., "QTL Hotspot A", "QTL Hotspot B", "QTL Hotspot C", and "QTL Hotspot D" located on Chr06, Chr10, Chr13, and Chr20, respectively. Based on gene annotation, gene ontology (GO) enrichment, and RNA-Seq analysis, 23 genes within four "QTL Hotspots" were predicted as possible candidates, regulating soybean seed size and shape. Network analyses demonstrated that 15 QTLs showed significant additive x environment (AE) effects, and 16 pairs of QTLs showing epistatic effects were also detected. However, except three epistatic QTLs, viz., qSL-13-3ZY, qSL-13-4ZY, and qSW-13-4ZY, all the remaining QTLs depicted no main effects. Hence, the present study is a detailed and comprehensive investigation uncovering the genetic basis of seed size and shape in soybeans. The use of a high-density map identified new genomic regions providing valuable information and could be the primary target for further fine mapping, candidate gene identification, and marker-assisted breeding (MAB).

Comparative transcriptome, physiological and biochemical analyses reveal response mechanism mediated by CBF4 and ICE2 in enhancing cold stress tolerance in Gossypium thurberi

Low temperature is one of the key environmental stresses that impair plant growth and significantly restricts the productivity and spatial distribution of crop plants. Gossypium thurberi, a wild diploid cotton species, has adapted to a wide range of temperatures and exhibits a better tolerance to chilling stress. Here, we compared phenotypes and physiochemical changes in G. thurberi under cold stress and found this species indeed showed better cold tolerance. Therefore, to understand the molecular mechanisms of the cold tolerance in G. thurberi, we compared transcription changes in leaves of G. thurberi under cold stress by high-throughput transcriptome sequencing. In total, 35 617 unigenes were identified in the whole-genome transcription profile, and 4226 differentially expressed genes (DEGs) were discovered in the leaves upon cold treatment. Gene Ontology (GO) classification analyses showed that the majority of DEGs belonged to categories of signal transduction, transcription factors (TFs) and carbohydrate transport and metabolism. The expression of several cold-responsive genes such as ICE1, CBF4, RAP2-7 and abscisic acid (ABA) biosynthesis genes involved in different signalling pathways were induced after G. thurberi seedlings were exposed to cold stress. Furthermore, cold sensitivity was increased in CBF4 and ICE2 virus-induced gene silencing (VIGS) plants, and high level of malondialdehyde (MDA) showed that the CBF4 and ICE2 silenced plants were under oxidative stress compared to their wild types, which relatively had higher levels of antioxidant enzyme activity, as evident by high levels of proline and superoxide dismutase (SOD) content. In conclusion, our findings reveal a new regulatory network of cold stress response in G. thurberi and broaden our understanding of the cold tolerance mechanism in cotton, which might accelerate functional genomics studies and genetic improvement for cold stress tolerance in cultivated cotton.

Advances in Omics Approaches for Abiotic Stress Tolerance in Tomato

Tomato, one of the most important crops worldwide, has a high demand in the fresh fruit market and processed food industries. Despite having considerably high productivity, continuous supply as per the market demand is hard to achieve, mostly because of periodic losses occurring due to biotic as well as abiotic stresses. Although tomato is a temperate crop, it is grown in almost all the climatic zones because of widespread demand, which makes it challenge to adapt in diverse conditions. Development of tomato cultivars with enhanced abiotic stress tolerance is one of the most sustainable approaches for its successful production. In this regard, efforts are being made to understand the stress tolerance mechanism, gene discovery, and interaction of genetic and environmental factors. Several omics approaches, tools, and resources have already been developed for tomato growing. Modern sequencing technologies have greatly accelerated genomics and transcriptomics studies in tomato. These advancements facilitate Quantitative trait loci (QTL) mapping, genome-wide association studies (GWAS), and genomic selection (GS). However, limited efforts have been made in other omics branches like proteomics, metabolomics, and ionomics. Extensive cataloging of omics resources made here has highlighted the need for integration of omics approaches for efficient utilization of resources and a better understanding of the molecular mechanism. The information provided here will be helpful to understand the plant responses and the genetic regulatory networks involved in abiotic stress tolerance and efficient utilization of omics resources for tomato crop improvement.

Genome-wide transcriptional profiling for elucidating the effects of brassinosteroids on Glycine max during early vegetative development

Soybean is a widely grown grain legume and one of the most important economic crop species. Brassinosteroids play a crucial role in plant vegetative growth and reproductive development. However, it remains unclear how BRs regulate the developmental processes in soybean, and the molecular mechanism underlying soybean early development is largely unexplored. In this study, we first characterized how soybean early vegetative growth was specifically regulated by the BR biosynthesis inhibitor propiconazole; this characterization included shortened root and shoot lengths, reduced leaf area, and decreased chlorophyll content. In addition, the growth inhibition induced by Pcz could be rescued by exogenous brassinolide application. The RNA-seq technique was employed to investigate the BR regulatory networks during soybean early vegetative development. Identification and analysis of differentially expressed genes indicated that BRs orchestrate a wide range of cellular activities and biological processes in soybean under various BR concentrations. The regulatory networks between BRs and multiple hormones or stress-related pathways were investigated. The results provide a comprehensive view of the physiological functions of BRs and new insights into the molecular mechanisms at the transcriptional level of BR regulation of soybean early development.

Identification of Genes/Proteins Related to Submergence Tolerance by Transcriptome and Proteome Analyses in Soybean

Flooding can lead to yield reduction of soybean. Therefore, identification of flooding tolerance genes has great significance in production practice. In this study, Qihuang 34, a highly-resistant variety to flooding stress, was selected for submergence treatments. Transcriptome and proteome analyses were conducted, by which twenty-two up-regulated differentially expressed genes (DEGs)/differentially expressed proteins (DEPs) associated with five KEGG pathways were isolated. The number of the DEGs/DEPs enriched in glycolysis/gluconeogenesis was the highest. Four of these genes were confirmed by RT-qPCR, suggesting that glycolysis/gluconeogenesis may be activated to generate energy for plant survival under anaerobic conditions. Thirty-eight down-regulated DEGs/DEPs associated with six KEGG pathways were identified under submergence stress. Eight DEGs/DEPs enriched in phenylpropanoid biosynthesis were assigned to peroxidase, which catalyzes the conversion of coumaryl alcohol to hydroxy-phenyl lignin in the final step of lignin biosynthesis. Three of these genes were confirmed by RT-qPCR. The decreased expression of these genes led to the inhibition of lignin biosynthesis, which may be the cause of plant softening under submergence stress for a long period of time. This study revealed a number of up-/down-regulated pathways and the corresponding DEGs/DEPs, by which, a better understanding of the mechanisms of submergence tolerance in soybean may be achieved.

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.

Transcriptome profiling illustrates expression signatures of dehydration tolerance in developing grasspea seedlings

This study highlights dehydration-mediated temporal changes in physicochemical, transcriptome and metabolome profiles indicating altered gene expression and metabolic shifts, underlying endurance and adaptation to stress tolerance in the marginalized crop, grasspea. Grasspea, often regarded as an orphan legume, is recognized to be fairly tolerant to water-deficit stress. In the present study, 3-week-old grasspea seedlings were subjected to dehydration by withholding water over a period of 144 h. While there were no detectable phenotypic changes in the seedlings till 48 h, the symptoms appeared during 72 h and aggravated upon prolonged dehydration. The physiological responses to water-deficit stress during 72-96 h displayed a decrease in pigments, disruption in membrane integrity and osmotic imbalance. We evaluated the temporal effects of dehydration at the transcriptome and metabolome levels. In total, 5201 genes of various functional classes including transcription factors, cytoplasmic enzymes and structural cell wall proteins, among others, were found to be dehydration-responsive. Further, metabolome profiling revealed 59 dehydration-responsive metabolites including sugar alcohols and amino acids. Despite the lack of genome information of grasspea, the time course of physicochemical and molecular responses suggest a synchronized dehydration response. The cross-species comparison of the transcriptomes and metabolomes with other legumes provides evidence for marked molecular diversity. We propose a hypothetical model that highlights novel biomarkers and explain their relevance in dehydration-response, which would facilitate targeted breeding and aid in commencing crop improvement efforts.

SlNAC2 overexpression in Arabidopsis results in enhanced abiotic stress tolerance with alteration in glutathione metabolism

Plant NAC (NAM, ATAF, and CUC) transcription factors (TF) have important roles to play in abiotic stress responses through activation of a battery of functional genes/transcriptional regulators responsible for stress tolerance. Here we report the cloning of a novel Solanum lycopersicum L., NAC2 TF having 960 nucleotides long CDS (GenBank: KT740994.1). Phylogenetic analysis depicted the similarity of SlNAC2 to other orthologs. SlNAC2 was overexpressed in Arabidopsis thaliana to assess and characterize its role in plant abiotic stress responses. The transgenic events were first confirmed by genomic DNA PCR and qRT PCR; then the T3 generation plants were used for stress assays. Soil stress assay depicted better survivability of the transgenic plants under both salt (NaCl) and drought (PEG) stress. The transgenic plants showed enhanced endurance; with better antioxidative response, reduced accumulation of reactive oxygen species (ROS) molecules and better retention of water in tissue. This study for the very first time analyzed the different stakeholders of the glutathione metabolism in SlNAC2 overexpressing transgenic lines on exposure to both salinity and PEG stress. The expression of the two genes (ɤ-ECS, GS) responsible for glutathione biosynthesis increased with SlNAC2 overexpression. Further glutathione reductase responsible for reduction of glutathione disulfide (GSSG) to glutathione (GSH) also increased significantly which suggested the regulation of glutathione metabolism as a mechanism for the osmotic stress tolerance conferred to plants upon NAC overexpression.

Effect of wet storage conditions on potato tuber transcriptome, phytohormones and growth

BACKGROUND

Stored potato (Solanum tuberosum L.) tubers are sensitive to wet conditions that can cause rotting in long-term storage. To study the effect of water on the tuber surface during storage, microarray analysis, RNA-Seq profiling, qRT-PCR and phytohormone measurements were performed to study gene expression and hormone content in wet tubers incubated at two temperatures: 4 °C and 15 °C. The growth of the plants was also observed in a greenhouse after the incubation of tubers in wet conditions.

RESULTS

Wet conditions induced a low-oxygen response, suggesting reduced oxygen availability in wet tubers at both temperatures when compared to that in the corresponding dry samples. Wet conditions induced genes coding for heat shock proteins, as well as proteins involved in fermentative energy production and defense against reactive oxygen species (ROS), which are transcripts that have been previously associated with low-oxygen stress in hypoxic or anoxic conditions. Wet treatment also induced senescence-related gene expression and genes involved in cell wall loosening, but downregulated genes encoding protease inhibitors and proteins involved in chloroplast functions and in the biosynthesis of secondary metabolites. Many genes involved in the production of phytohormones and signaling were also affected by wet conditions, suggesting altered regulation of growth by wet conditions. Hormone measurements after incubation showed increased salicylic acid (SA), abscisic acid (ABA) and auxin (IAA) concentrations as well as reduced production of jasmonate 12-oxo-phytodienoic acid (OPDA) in wet tubers. After incubation in wet conditions, the tubers produced fewer stems and more roots compared to controls incubated in dry conditions.

CONCLUSIONS

In wet conditions, tubers invest in ROS protection and defense against the abiotic stress caused by reduced oxygen due to excessive water. Changes in ABA, SA and IAA that are antagonistic to jasmonates affect growth and defenses, causing induction of root growth and rendering tubers susceptible to necrotrophic pathogens. Water on the tuber surface may function as a signal for growth, similar to germination of seeds.

Comparative Transcriptome Combined with Proteome Analyses Revealed Key Factors Involved in Alfalfa (Medicago sativa) Response to Waterlogging Stress

Alfalfa (Medicago sativa) is the most widely grown and most important forage crop in the world. However, alfalfa is susceptible to waterlogging stress, which is the major constraint for its cultivation area and crop production. So far, the molecular mechanism of alfalfa response to the waterlogging is largely unknown. Here, comparative transcriptome combined with proteomic analyses of two cultivars (M12, tolerant; M25, sensitive) of alfalfa showing contrasting tolerance to waterlogging were performed to understand the mechanism of alfalfa in response to waterlogging stress. Totally, 748 (581 up- and 167 down-regulated) genes were differentially expressed in leaves of waterlogging-stressed alfalfa compared with the control (M12_W vs. M12_CK), whereas 1193 (740 up- and 453 down-regulated) differentially abundant transcripts (DATs) were detected in the leaves of waterlogging-stressed plants in comparison with the control plants (M25_W vs. M25_CK). Furthermore, a total of 187 (122 up- and 65 down-regulated) and 190 (105 up- and 85 down-regulated) differentially abundant proteins (DAPs) were identified via isobaric tags for relative and absolute quantification (iTRAQ) method in M12_W vs. M12_CK and M25_W vs. M25_CK comparison, respectively. Compared dataset analysis of proteomics and transcriptomics revealed that 27 and eight genes displayed jointly up-regulated or down-regulated expression profiles at both mRNA and protein levels in M12_W vs. M12_CK comparison, whereas 30 and 27 genes were found to be co-up-regulated or co-down-regulated in M25_W vs. M25_CK comparison, respectively. The strongly enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways for co-up-regulated genes at mRNA and protein levels in M12_W vs. M12_CK comparison were 'Amino sugar and nucleotide sugar metabolism', 'Arginine and proline metabolism' and 'Starch and sucrose metabolism', whereas co-up-regulated protein-related pathways including 'Arginine and proline metabolism' and 'Valine, leucine and isoleucine degradation' were largely enriched in M25_W vs. M25_CK comparison. Importantly, the identified genes related to beta-amylase, Ethylene response Factor (ERF), Calcineurin B-like (CBL) interacting protein kinases (CIPKs), Glutathione peroxidase (GPX), and Glutathione-S-transferase (GST) may play key roles in conferring alfalfa tolerance to waterlogging stress. The present study may contribute to our understanding the molecular mechanism underlying the responses of alfalfa to waterlogging stress, and also provide important clues for further study and in-depth characterization of waterlogging-resistance breeding candidate genes in alfalfa.

Physiological and Transcriptomic Analysis of Tree Peony (Paeonia section Moutan DC.) in Response to Drought Stress

Tree peony (Paeonia section Moutan DC.) is a famous ornamental plant, and P. ostii has been used for seed oil production in China because it is rich in α-linolenic acid. P. ostii has some resistance to drought, but lack of water can severely hinder its growth and development in arid areas. In order to clarify drought stress induced physiological and molecular changes of P. ostia, its physiological and transcriptomic analyses were performed under drought stress, and we found that P. ostii leaves drooped significantly 12 days after treatment and observed a significant increase in all detected physiological indices in response to drought treatment except leaf water content, chlorophyll, and carotenoid content. Meanwhile, the activity of three antioxidant enzymes basically increased under drought treatment. Moreover, drought treatment significantly reduced photosynthetic and chlorophyll fluorescence parameters except non-photochemical quenching (qN), and maintained more intact mesophyll cell structures. Additionally, many differentially expressed genes (DEGs) were found by transcriptome sequencing, which play an important role in P. ostia drought tolerance by controlling a variety of biological processes, including the reactive oxygen species (ROS) system, chlorophyll degradation and photosynthetic competency, fatty acid metabolism, proline metabolism, biosynthesis of secondary metabolism, and plant hormone metabolism. These results provide a better understanding of P. ostii responses to drought stress.

Physiological and transcriptomic responses of reproductive stage soybean to drought stress

KEY MESSAGE

The dynamic alterations of the physiological and molecular processes in reproductive stage soybean indicated the dramatic impact caused by drought. Drought is a major abiotic stress that limits soybean (Glycine max) production. Most prior studies were focused on either model species or crops that are at their vegetative stages. It is known that the reproductive stage of soybean is more susceptible to drought. Therefore, an understanding on the responsive mechanisms during this stage will not only be important for basic plant physiology, but the knowledge can also be used for crop improvement via either genetic engineering or molecular breeding. In this study, physiological measurements and RNA-Seq analysis were used to dissect the metabolic alterations and molecular responses in the leaves of soybean grown at drought condition. Photosynthesis rate, stomata conductance, transpiration, and water potential were reduced. The activities of SOD and CAT were increased, while the activity of POD stayed unchanged. A total of 2771 annotated genes with at least twofold changes were found to be differentially expressed in the drought-stressed plants in which 1798 genes were upregulated and 973 were downregulated. Via KEGG analysis, these genes were assigned to multiple molecular pathways, including ABA biogenesis, compatible compound accumulation, secondary metabolite synthesis, fatty acid desaturation, plant transcription factors, etc. The large number of differentially expressed genes and the diverse pathways indicated that soybean employs complicated mechanisms to cope with drought. Some of the identified genes and pathways can be used as targets for genetic engineering or molecular breeding to improve drought resistance in soybean.

Transcriptome profilling analysis characterized the gene expression patterns responded to combined drought and heat stresses in soybean

Heat and drought are the two major abiotic stress limiting soybean growth and output worldwide. Knowledge of the molecular mechanisms underlying the responses to heat, drought, and combined stress is essential for soybean molecular breeding. In this study, RNA-sequencing was used to determine the transcriptional responses of soybean to heat, drought and combined stress. RNA-sequencing analysis demonstrated that many genes involved in the defense response, photosynthesis, metabolic process, etc. are differentially expressed in response to drought and heat. However, 1468 and 1220 up-regulated and 1146 and 686 down-regulated genes were confirmed as overlapping differentially expressed genes at 8 h and 24 h after treatment, and these genes are mainly involved in transport, binding and defense response. Furthermore, we compared the heat, drought and the combined stress-responsive genes and identified potential new targets for enhancing stress tolerance of soybean. Comparison of single and combined stress suggests the combined stress did not result in a simple additive response, and that there may be a synergistic response to the combination of drought and heat in soybean.

Comparative transcriptomic analysis reveals common molecular factors responsive to heat and drought stress in Agrostis stolonifera

Heat and drought stress are primary abiotic stresses confining growth of cool-season grass species during summer. The objective of this study was to identify common molecular factors and metabolic pathways associated with heat and drought responses in creeping bentgrass (Agrostis stolonifera) by comparative analysis of transcriptomic profiles between plants exposed to heat and drought stress. Plants were exposed to heat stress (35/30 °C day/night temperature) or drought stress by withholding irrigation for 21 d in growth chambers. Transcriptomic profiling by RNA-seq in A. stolonifera (cv. 'Penncross') found 670 commonly up-regulated and 812 commonly down-regulated genes by heat and drought stress. Transcriptional up-regulations of differentially expressed genes (DEGs) due to heat and drought stress include genes that were highly enriched in oxylipin biosynthetic process and proline biosynthetic process. Transcriptional down-regulations of genes under heat and drought stress were highly enriched and involved in thiamine metabolic process and calcium sensing receptor. These commonly-regulated genes by heat and drought stress identified in A. stolonifera suggested that drought and heat responses shared such common molecular factors and pathways, which could be potential candidate genes for genetic modification of improving plant tolerance to the combined heat and drought stress.

Characterization of the XTH Gene Family: New Insight to the Roles in Soybean Flooding Tolerance

Xyloglucan endotransglycosylases/hydrolases (XTHs) are a class of enzymes involved in the construction and remodeling of cellulose/xyloglucan crosslinks and play an important role in regulating cell wall extensibility. However, little is known about this class of enzymes in soybean. Here, 61 soybean XTH genes (GmXTHs) were identified and classified into three subgroups through comparative phylogenetic analysis. Genome duplication greatly contributed to the expansion of GmXTH genes in soybean. A conserved amino acid motif responsible for the catalytic activity was identified in all GmXTHs. Further expression analysis revealed that most GmXTHs exhibited a distinct organ-specific expression pattern, and the expression level of many GmXTH genes was significantly associated with ethylene and flooding stress. To illustrate a possible role of XTH genes in regulating stress responses, the ArabidopsisAtXTH31 gene was overexpressed in soybean. The generated transgenic plants exhibited improved tolerance to flooding stress, with a higher germination rate and longer roots/hypocotyls during the seedling stage and vegetative growth stages. In summary, our combined bioinformatics and gene expression pattern analyses suggest that GmXTH genes play a role in regulating soybean stress responses. The enhanced soybean flooding tolerance resulting from the expression of an Arabidopsis XTH also supports the role of XTH genes in regulating plant flooding stress responses.

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.

Molecular and physiological responses in roots of two full-sib poplars uncover mechanisms that contribute to differences in partial submergence tolerance

Poplar is a major afforestation tree species in flood-prone areas. Here, we compared molecular and physiological responses in the roots of two full-sib poplar clones, LS1 (flood-tolerant) and LS2 (flood-susceptive), subjected to stagnant flooding using transcript and metabolite profiling. LS1 displayed less phenotypic damage and superior leaf gas exchange and plant growth compared with those of LS2. We concluded that three characteristics might contribute to the differences in flood tolerance between LS1 and LS2. First, fermentation was initiated through lactic dehydrogenation in LS1 roots under flooding and subsequently dominated by alcohol fermentation. However, lactic dehydrogenase was persistently active in flooded LS2. Second, 13 differentially expressed genes associated with energy and O2 consumption processes under soil flooding had lower transcript levels in LS1 than those in LS2, which might contribute to better energy-/O2-saving abilities and behaviours in flood-tolerant LS1 than those in flood-susceptible LS2 under hypoxic stress. Third, LS1 possessed increased reactive oxygen species scavenging abilities compared with those of LS2 under edaphic flooding. Our data are a valuable contribution to understanding the mechanisms involved in the flood tolerance of poplar.

Improving Flooding Tolerance of Crop Plants

A major problem of climate change is the increasing duration and frequency of heavy rainfall events. This leads to soil flooding that negatively affects plant growth, eventually leading to death of plants if the flooding persists for several days. Most crop plants are very sensitive to flooding, and dramatic yield losses occur due to flooding each year. This review summarizes recent progress and approaches to enhance crop resistance to flooding. Most experiments have been done on maize, barley, and soybean. Work on other crops such as wheat and rape has only started. The most promising traits that might enhance crop flooding tolerance are anatomical adaptations such as aerenchyma formation, the formation of a barrier against radial oxygen loss, and the growth of adventitious roots. Metabolic adaptations might be able to improve waterlogging tolerance as well, but more studies are needed in this direction. Reasonable approaches for future studies are quantitative trait locus (QTL) analyses or genome-wide association (GWA) studies in combination with specific tolerance traits that can be easily assessed. The usage of flooding-tolerant relatives or ancestral cultivars of the crop of interest in these experiments might enhance the chances of finding useful tolerance traits to be used in breeding.

Transcriptomic analysis of drought stress responses of sea buckthorn (Hippophae rhamnoidessubsp. sinensis) by RNA-Seq

Sea buckthorn is one of the most important eco-economic tree species in China due to its ability to grow and produce acceptable yields under limited water and fertilizer availability. In this study, the differentially expressed genes under drought stress (DS) of sea buckthorn were identified and compared with control (CK) by RNA-Seq. A total of 122,803 unigenes were identified in sea buckthorn, and 70,025 unigenes significantly matched a sequence in at least one of the seven databases. A total of 24,060 (19.59%) unigenes can be assigned to 19 KEGG pathways, and 1,644 unigenes were differentially expressed between DS and CK, of which 519 unigenes were up-regulated and 1,125 unigenes down-regulated. Of the 47 significantly enriched GO terms, 14, 7 and 26 items were related to BP, CC and MF, respectively. KEGG enrichment analysis showed 398 DEGs involved in 97 different pathways, of which 119 DEGs were up-regulated and 279 DEGs were down-regulated under drought stress. In addition, we found 4438 transcriptor factors (TFs) in sea buckthorn, of which 100 were differentially expressed between DS and CK. These results lay a first foundation for further investigations of the very specific functions of these unigenes in sea buckthorn in response to drought stress.

Plant-associated bacteria mitigate drought stress in soybean

Agriculture accounts for ~ 70% of all water use and the world population is increasing annually; soon more people will need to be fed, while also using less water. The use of plant-associated bacteria (PAB) is an eco-friendly alternative that can increase crop water use efficiency. This work aimed to study the effect of some PAB on increasing soybean tolerance to drought stress, the mechanisms of the drought tolerance process, and the effect of the PAB on promoting plant growth and on the biocontrol of Sclerotinia sclerotiorum. PAB were isolated from soybean rhizosphere and S. sclerotiorum sclerotia. The strains identified as UFGS1 (Bacillus subtilis), UFGS2 (Bacillus thuringiensis), UFGRB2 and UFGRB3 (Bacillus cereus) were selected on their ability to grow in media with reduced water activity. Soybean plants were inoculated with the PAB and evaluated for growth promotion, physiological and molecular parameters, after drought stress. Under drought stress, UFGS2 and UFGRB2 sustained potential quantum efficiency of PSII (Fv/Fm), while a decrease was found in the control plants. Moreover, UFGS2 and UFGRB3 maintained the photosynthetic rates in non-stressed conditions compared to the control. UFGS2-treated plants showed a higher stomatal conductance and higher transpiration than the control, after drought stress. Some PAB-treated plants also had other beneficial phenotypes, such as increases in fresh and dried biomass relative to the control. Differential gene expression analysis of genes involved in plant stress pathways shows changes in expression in PAB-treated plants. Results from this study suggest that PAB can mitigate drought stress in soybean and may improve water efficiency under certain conditions.

Subcellular Proteomics: Application to Elucidation of Flooding-Response Mechanisms in Soybean

Soybean, which is rich in protein and oil, is cultivated in several climatic zones; however, its growth is markedly decreased by flooding. Proteomics is a useful tool for understanding the flooding-response mechanism in soybean. Subcellular proteomics has the potential to elucidate localized cellular responses and investigate communications among subcellular components during plant growth and during stress. Under flooding, proteins related to signaling, stress and the antioxidative system are increased in the plasma membrane; scavenging enzymes for reactive-oxygen species are suppressed in the cell wall; protein translation is suppressed through inhibition of proteins related to preribosome biogenesis and mRNA processing in the nucleus; levels of proteins involved in the electron transport chain are reduced in the mitochondrion; and levels of proteins related to protein folding are decreased in the endoplasmic reticulum. This review discusses the advantages of a gel-free/label-free proteomic technique and methods of plant subcellular purification. It also summarizes cellular events in soybean under flooding and discusses future prospects for generation of flooding-tolerant soybean.

Evolutionary and expression analyses of soybean basic Leucine zipper transcription factor family

BACKGROUND

Soybean, a major legume crop native to East Asia, presents a wealth of resources for utilization. The basic leucine zipper (bZIP) transcription factors play important roles in various biological processes including developmental regulation and responses to environmental stress stimuli. Currently, little information is available regarding the bZIP family in the legume crop soybean.

RESULTS

Using a genome-wide domain analysis, we identified 160 GmbZIP genes in soybean genome, named from GmbZIP1 to GmbZIP160. These 160GmbZIP genes, distributed unevenly across 20 chromosomes, were grouped into 12 subfamilies based on phylogenetic analysis. Gene structure and conserved motif analyses showed that GmbZIP within the same subfamily shared similar intron-exon organizations and motif composition. Syntenic and phylogenetic analyses identified 40 Arabidopsis bZIP genes and 83 soybean bZIP genes as orthologs. By investigating the expression profiling of GmbZIP in different tissues and under drought and flooding stresses, we showed that a majority of GmbZIP (83.44%) exhibited transcript abundance in all examined tissues and 75.6% displayed transcript changes after drought and flooding treatment, suggesting that GmbZIP may play a broad role in soybean development and response to water stress.

CONCLUSIONS

One hundred sixty GmbZIP genes were identified in soybean genome. Our results provide insights for the evolutionary history of bZIP family in soybean and shed light on future studies on the function of bZIP genes in response to water stress in soybean.

NAC transcription factor JUNGBRUNNEN1 enhances drought tolerance in tomato

Water deficit (drought stress) massively restricts plant growth and the yield of crops; reducing the deleterious effects of drought is therefore of high agricultural relevance. Drought triggers diverse cellular processes including the inhibition of photosynthesis, the accumulation of cell-damaging reactive oxygen species and gene expression reprogramming, besides others. Transcription factors (TF) are central regulators of transcriptional reprogramming and expression of many TF genes is affected by drought, including members of the NAC family. Here, we identify the NAC factor JUNGBRUNNEN1 (JUB1) as a regulator of drought tolerance in tomato (Solanum lycopersicum). Expression of tomato JUB1 (SlJUB1) is enhanced by various abiotic stresses, including drought. Inhibiting SlJUB1 by virus-induced gene silencing drastically lowers drought tolerance concomitant with an increase in ion leakage, an elevation of hydrogen peroxide (H₂ O₂ ) levels and a decrease in the expression of various drought-responsive genes. In contrast, overexpression of AtJUB1 from Arabidopsis thaliana increases drought tolerance in tomato, alongside with a higher relative leaf water content during drought and reduced H₂ O₂ levels. AtJUB1 was previously shown to stimulate expression of DREB2A, a TF involved in drought responses, and of the DELLA genes GAI and RGL1. We show here that SlJUB1 similarly controls the expression of the tomato orthologs SlDREB1, SlDREB2 and SlDELLA. Furthermore, AtJUB1 directly binds to the promoters of SlDREB1, SlDREB2 and SlDELLA in tomato. Our study highlights JUB1 as a transcriptional regulator of drought tolerance and suggests considerable conservation of the abiotic stress-related gene regulatory networks controlled by this NAC factor between Arabidopsis and tomato.

Insights into soybean transcriptome reconfiguration under hypoxic stress: Functional, regulatory, structural, and compositional characterization

Soybean (Glycine max) is one of the major crops worldwide and flooding stress affects the production and expansion of cultivated areas. Oxygen is essential for mitochondrial aerobic respiration to supply the energy demand of plant cells. Because oxygen diffusion in water is 10,000 times lower than in air, partial (hypoxic) or total (anoxic) oxygen deficiency is important component of flooding. Even when oxygen is externally available, oxygen deficiency frequently occurs in bulky, dense or metabolically active tissues such as phloem, meristems, seeds, and fruits. In this study, we analyzed conserved and divergent root transcriptional responses between flood-tolerant Embrapa 45 and flood-sensitive BR 4 soybean cultivars under hypoxic stress conditions with RNA-seq. To understand how soybean genes evolve and respond to hypoxia, stable and differentially expressed genes were characterized structurally and compositionally comparing its mechanistic relationship. Between cultivars, Embrapa 45 showed less up- and more down-regulated genes, and stronger induction of phosphoglucomutase (Glyma05g34790), unknown protein related to N-terminal protein myristoylation (Glyma06g03430), protein suppressor of phyA-105 (Glyma06g37080), and fibrillin (Glyma10g32620). RNA-seq and qRT-PCR analysis of non-symbiotic hemoglobin (Glyma11g12980) indicated divergence in gene structure between cultivars. Transcriptional changes for genes in amino acids and derivative metabolic process suggest involvement of amino acids metabolism in tRNA modifications, translation accuracy/efficiency, and endoplasmic reticulum stress in both cultivars under hypoxia. Gene groups differed in promoter TATA box, ABREs (ABA-responsive elements), and CRT/DREs (C-repeat/dehydration-responsive elements) frequency. Gene groups also differed in structure, composition, and codon usage, indicating biological significances. Additional data suggests that cis-acting ABRE elements can mediate gene expression independent of ABA in soybean roots under hypoxia.

Molecular mapping and genomics of soybean seed protein: a review and perspective for the future

KEY MESSAGE

Genetic improvement of soybean protein meal is a complex process because of negative correlation with oil, yield, and temperature. This review describes the progress in mapping and genomics, identifies knowledge gaps, and highlights the need of integrated approaches. Meal protein derived from soybean [Glycine max (L) Merr.] seed is the primary source of protein in poultry and livestock feed. Protein is a key factor that determines the nutritional and economical value of soybean. Genetic improvement of soybean seed protein content is highly desirable, and major quantitative trait loci (QTL) for soybean protein have been detected and repeatedly mapped on chromosomes (Chr.) 20 (LG-I), and 15 (LG-E). However, practical breeding progress is challenging because of seed protein content's negative genetic correlation with seed yield, other seed components such as oil and sucrose, and interaction with environmental effects such as temperature during seed development. In this review, we discuss rate-limiting factors related to soybean protein content and nutritional quality, and potential control factors regulating seed storage protein. In addition, we describe advances in next-generation sequencing technologies for precise detection of natural variants and their integration with conventional and high-throughput genotyping technologies. A syntenic analysis of QTL on Chr. 15 and 20 was performed. Finally, we discuss comprehensive approaches for integrating protein and amino acid QTL, genome-wide association studies, whole-genome resequencing, and transcriptome data to accelerate identification of genomic hot spots for allele introgression and soybean meal protein improvement.

Global gene expression in cotton (Gossypium hirsutum L.) leaves to waterlogging stress

Cotton is sensitive to waterlogging stress, which usually results in stunted growth and yield loss. To date, the molecular mechanisms underlying the responses to waterlogging in cotton remain elusive. Cotton was grown in a rain-shelter and subjected to 0 (control)-, 10-, 15- and 20-d waterlogging at flowering stage. The fourth-leaves on the main-stem from the top were sampled and immediately frozen in liquid nitrogen for physiological measurement. Global gene transcription in the leaves of 15-d waterlogged plants was analyzed by RNA-Seq. Seven hundred and ninety four genes were up-regulated and 1018 genes were down-regulated in waterlogged cotton leaves compared with non-waterlogged control. The differentially expressed genes were mainly related to photosynthesis, nitrogen metabolism, starch and sucrose metabolism, glycolysis and plant hormone signal transduction. KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis indicated that most genes related to flavonoid biosynthesis, oxidative phosphorylation, amino acid metabolism and biosynthesis as well as circadian rhythm pathways were differently expressed. Waterlogging increased the expression of anaerobic fermentation related genes, such as alcohol dehydrogenase (ADH), but decreased the leaf chlorophyll concentration and photosynthesis by down-regulating the expression of photosynthesis related genes. Many genes related to plant hormones and transcription factors were differently expressed under waterlogging stress. Most of the ethylene related genes and ethylene-responsive factor-type transcription factors were up-regulated under water-logging stress, suggesting that ethylene may play key roles in the survival of cotton under waterlogging stress.

Elucidation of the molecular responses to waterlogging in Sesbania cannabina roots by transcriptome profiling

Sesbania cannabina, a multipurpose leguminous crop, is highly resistant to waterlogging stress. However, the scant genomic resources in the genus Sesbania have greatly hindered further exploration of the mechanisms underlying its waterlogging tolerance. Here, the genetic basis of flooding tolerance in S. cannabina was examined by transcriptome-wide gene expression changes using RNA-Seq in seedlings exposed to short-term (3 h) and long-term (27 h) waterlogging. After de- novo assembly, 213990 unigenes were identified, of which 145162 (79.6%) were annotated. Gene Ontology and pathway enrichment analyses revealed that the glycolysis and fermentation pathways were stimulated to produce ATP under hypoxic stress conditions. Energy-consuming biosynthetic processes were dramatically repressed by short and long term waterlogging, while amino acid metabolism was greatly induced to maintain ATP levels. The expression pattern of 10 unigenes involved in phenylpropanoid biosynthesis, glycolysis, and amino acid metabolism revealed by qRT-PCR confirmed the RNA-Seq data. The present study is a large-scale assessment of genomic resources of Sesbania and provides guidelines for probing the molecular mechanisms underlying S. cannabina waterlogging tolerance.

Metabolic profiles of flooding-tolerant mechanism in early-stage soybean responding to initial stress

KEY MESSAGE

Metabolomic analysis of flooding-tolerant mutant and abscisic acid-treated soybeans suggests that accumulated fructose might play a role in initial flooding tolerance through regulation of hexokinase and phosphofructokinase. Soybean is sensitive to flooding stress, which markedly reduces plant growth. To explore the mechanism underlying initial-flooding tolerance in soybean, mass spectrometry-based metabolomic analysis was performed using flooding-tolerant mutant and abscisic-acid treated soybeans. Among the commonly-identified metabolites in both flooding-tolerant materials, metabolites involved in carbohydrate and organic acid displayed same profile at initial-flooding stress. Sugar metabolism was highlighted in both flooding-tolerant materials with the decreased and increased accumulation of sucrose and fructose, respectively, compared to flooded soybeans. Gene expression of hexokinase 1 was upregulated in flooded soybean; however, it was downregulated in both flooding-tolerant materials. Metabolites involved in carbohydrate/organic acid and proteins related to glycolysis/tricarboxylic acid cycle were integrated. Increased protein abundance of phosphofructokinase was identified in both flooding-tolerant materials, which was in agreement with its enzyme activity. Furthermore, sugar metabolism was pointed out as the tolerant-responsive process at initial-flooding stress with the integration of metabolomics, proteomics, and transcriptomics. Moreover, application of fructose declined the increased fresh weight of plant induced by flooding stress. These results suggest that fructose might be the critical metabolite through regulation of hexokinase and phosphofructokinase to confer initial-flooding stress in soybean.

Genome-wide identification of gene expression in contrasting maize inbred lines under field drought conditions reveals the significance of transcription factors in drought tolerance

Drought is a major threat to maize growth and production. Understanding the molecular regulation network of drought tolerance in maize is of great importance. In this study, two maize inbred lines with contrasting drought tolerance were tested in the field under natural soil drought and well-watered conditions. In addition, the transcriptomes of their leaves was analyzed by RNA-Seq. In total, 555 and 2,558 genes were detected to specifically respond to drought in the tolerant and the sensitive line, respectively, with a more positive regulation tendency in the tolerant genotype. Furthermore, 4,700, 4,748, 4,403 and 4,288 genes showed differential expression between the two lines under moderate drought, severe drought and their well-watered controls, respectively. Transcription factors were enriched in both genotypic differentially expressed genes and specifically responsive genes of the tolerant line. It was speculated that the genotype-specific response of 20 transcription factors in the tolerance line and the sustained genotypically differential expression of 22 transcription factors might enhance tolerance to drought in maize. Our results provide new insight into maize drought tolerance-related regulation systems and provide gene resources for subsequent studies and drought tolerance improvement.

Metabolomic Profiling of Soybeans (Glycine max L.) Reveals the Importance of Sugar and Nitrogen Metabolism under Drought and Heat Stress

Soybean is an important crop that is continually threatened by abiotic stresses, especially drought and heat stress. At molecular levels, reduced yields due to drought and heat stress can be seen as a result of alterations in metabolic homeostasis of vegetative tissues. At present an incomplete understanding of abiotic stress-associated metabolism and identification of associated metabolites remains a major gap in soybean stress research. A study with a goal to profile leaf metabolites under control conditions (28/24 °C), drought [28/24 °C, 10% volumetric water content (VWC)], and heat stress (43/35 °C) was conducted in a controlled environment. Analyses of non-targeted metabolomic data showed that in response to drought and heat stress, key metabolites (carbohydrates, amino acids, lipids, cofactors, nucleotides, peptides and secondary metabolites) were differentially accumulated in soybean leaves. The metabolites for various cellular processes, such as glycolysis, the tricarboxylic acid (TCA) cycle, the pentose phosphate pathway, and starch biosynthesis, that regulate carbohydrate metabolism, amino acid metabolism, peptide metabolism, and purine and pyrimidine biosynthesis, were found to be affected by drought as well as heat stress. Computationally based regulatory networks predicted additional compounds that address the possibility of other metabolites and metabolic pathways that could also be important for soybean under drought and heat stress conditions. Metabolomic profiling demonstrated that in soybeans, keeping up with sugar and nitrogen metabolism is of prime significance, along with phytochemical metabolism under drought and heat stress conditions.

Analysis of EF-Hand Proteins in Soybean Genome Suggests Their Potential Roles in Environmental and Nutritional Stress Signaling

Calcium ion (Ca2+) is a universal second messenger that plays a critical role in plant responses to diverse physiological and environmental stimuli. The stimulus-specific signals are perceived and decoded by a series of Ca2+ binding proteins serving as Ca2+ sensors. The majority of Ca2+ sensors possess the EF-hand motif, a helix-loop-helix structure which forms a turn-loop structure. Although EF-hand proteins in model plant such as Arabidopsis have been well described, the identification, classification, and the physiological functions of EF-hand-containing proteins from soybean are not systemically reported. In this study, a total of at least 262 genes possibly encoding proteins containing one to six EF-hand motifs were identified in soybean genome. These genes include 6 calmodulins (CaMs), 144 calmodulin-like proteins (CMLs), 15 calcineurin B-like proteins, 50 calcium-dependent protein kinases (CDPKs), 13 CDPK-related protein kinases, 2 Ca2+- and CaM-dependent protein kinases, 17 respiratory burst oxidase homologs, and 15 unclassified EF-hand proteins. Most of these genes (87.8%) contain at least one kind of hormonal signaling- and/or stress response-related cis-elements in their -1500 bp promoter regions. Expression analyses by exploring the published microarray and Illumina transcriptome sequencing data revealed that the expression of these EF-hand genes were widely detected in different organs of soybean, and nearly half of the total EF-hand genes were responsive to various environmental or nutritional stresses. Quantitative RT-PCR was used to confirm their responsiveness to several stress treatments. To confirm the Ca2+-binding ability of these EF-hand proteins, four CMLs (CML1, CML13, CML39, and CML95) were randomly selected for SDS-PAGE mobility-shift assay in the presence and absence of Ca2+. Results showed that all of them have the ability to bind Ca2+. This study provided the first comprehensive analyses of genes encoding for EF-hand proteins in soybean. Information on the classification, phylogenetic relationships and expression profiles of soybean EF-hand genes in different tissues and under various environmental and nutritional stresses will be helpful for identifying candidates with potential roles in Ca2+ signal-mediated physiological processes including growth and development, plant-microbe interactions and responses to biotic and abiotic stresses.

Association of Candidate Genes With Submergence Response in Perennial Ryegrass

Perennial ryegrass is a popular cool-season grass species due to its high quality for forage and turf. The objective of this study was to identify associations of candidate genes with growth and physiological traits to submergence stress and recovery after de-submergence in a global collection of 94 perennial ryegrass accessions. Accessions varied largely in leaf color, plant height (HT), leaf fresh weight (LFW), leaf dry weight (LDW), and chlorophyll fluorescence (Fv/Fm) at 7 days of submergence and in HT, LFW and LDW at 7 days of recovery in two experiments. Among 26 candidate genes tested by various models, single nucleotide polymorphisms (SNPs) in 10 genes showed significant associations with traits including 16 associations for control, 10 for submergence, and 8 for recovery. Under submergence, Lp1-SST encoding sucrose:sucrose 1-fructosyltransferase and LpGA20ox encoding gibberellin 20-oxidase were associated with LFW and LDW, and LpACO1 encoding 1-aminocyclopropane-1-carboxylic acid oxidase was associated with LFW. Associations between Lp1-SST and HT, Lp6G-FFT encoding fructan:fructan 6G-fructosyltransferase and Fv/Fm, LpCAT encoding catalase and HT were also detected under submergence stress. Upon de-submergence, Lp1-SST, Lp6G-FFT, and LpPIP1 encoding plasma membrane intrinsic protein type 1 were associated with LFW or LDW, while LpCBF1b encoding C-repeat binding factor were associated with HT. Nine significant SNPs in Lp1-SST, Lp6G-FFT, LpCAT, and LpACO1 resulted in amino acid changes with five substitutions found in Lp1-SST under submergence or recovery. The results indicated that allelic diversity in genes involved in carbohydrate and antioxidant metabolism, ethylene and gibberellin biosynthesis, and transcript factor could contribute to growth variations in perennial ryegrass under submergence stress and recovery after de-submergence.

Metagenome-Wide Association Study and Machine Learning Prediction of Bulk Soil Microbiome and Crop Productivity

Areas within an agricultural field in the same season often differ in crop productivity despite having the same cropping history, crop genotype, and management practices. One hypothesis is that abiotic or biotic factors in the soils differ between areas resulting in these productivity differences. In this study, bulk soil samples collected from a high and a low productivity area from within six agronomic fields in Illinois were quantified for abiotic and biotic characteristics. Extracted DNA from these bulk soil samples were shotgun sequenced. While logistic regression analyses resulted in no significant association between crop productivity and the 26 soil characteristics, principal coordinate analysis and constrained correspondence analysis showed crop productivity explained a major proportion of the taxa variance in the bulk soil microbiome. Metagenome-wide association studies (MWAS) identified more Bradyrhizodium and Gammaproteobacteria in higher productivity areas and more Actinobacteria, Ascomycota, Planctomycetales, and Streptophyta in lower productivity areas. Machine learning using a random forest method successfully predicted productivity based on the microbiome composition with the best accuracy of 0.79 at the order level. Our study showed that crop productivity differences were associated with bulk soil microbiome composition and highlighted several nitrogen utility-related taxa. We demonstrated the merit of MWAS and machine learning for the first time in a plant-microbiome study.

Overexpression of the ascorbate peroxidase gene from eggplant and sponge gourd enhances flood tolerance in transgenic Arabidopsis

Previously, we found that the flood resistance of eggplant (Solanum melongena) and sponge gourd (Luffa cylindrica) enhanced ascorbate peroxidase (APX) activity under flooding, and consequently, both the SmAPX and LcAPX genes were cloned. In this study, the SmAPX and LcAPX genes were transferred under a ubiquitin promoter to Arabidopsis (At) via Agrobacterium tumefaciens. The expression and amount of APX and APX activities of the SmAPX and LcAPX transgenic lines were significantly higher than those of non-transgenic (NT) plants under a waterlogged condition. Furthermore, the SmAPX, LcAPX, At-sucrose synthases (SUS)-1, phosphoenolpyruvate carboxylase (PEPC), and lactate dehydrogenase (LDH) genes were overexpressed in all transgenic Arabidopsis lines after flooding treatment. Compared to NT plants, the malondialdehyde (MDA) contents and H₂O₂ accumulation were significantly lower, but germination rates were significantly higher in all transgenic lines with higher APX activity, indicating that the overexpression of SmAPX and LcAPX in Arabidopsis could enhance flood tolerance by eliminating H₂O₂. Moreover, Arabidopsis seedlings overexpressing SmAPX and LcAPX also displayed greater resistance to flooding and less oxidative injury than NT plants subjected to flooding condition.

Identification and Analysis of NaHCO3 Stress Responsive Genes in Wild Soybean (Glycine soja) Roots by RNA-seq

Soil alkalinity is a major abiotic constraint to crop productivity and quality. Wild soybean (Glycine soja) is considered to be more stress-tolerant than cultivated soybean (G. max), and has considerable genetic variation for increasing alkalinity tolerance of soybean. In this study, we analyzed the transcriptome profile in the roots of an alkalinity tolerant wild soybean variety N24852 at 12 and 24 h after 90 mM NaHCO₃ stress by RNA-sequencing. Compared with the controls, a total of 449 differentially expressed genes (DEGs) were identified, including 95 and 140 up-regulated genes, and 108 and 135 down-regulated genes at 12 and 24 h after NaHCO₃ treatment, respectively. Quantitative RT-PCR analysis of 14 DEGs showed a high consistency with their expression profiles by RNA-sequencing. Gene Ontology (GO) terms related to transcription factors and transporters were significantly enriched in the up-regulated genes at 12 and 24 h after NaHCO₃ stress, respectively. Nuclear factor Y subunit A transcription factors were enriched at 12 h after NaHCO₃ stress, and high percentages of basic helix-loop-helix, ethylene-responsive factor, Trihelix, and zinc finger (C2H2, C3H) transcription factors were found at both 12 and 24 h after NaHCO₃ stress. Genes related to ion transporters such as ABC transporter, aluminum activated malate transporter, glutamate receptor, nitrate transporter/proton dependent oligopeptide family, and S-type anion channel were enriched in up-regulated DEGs at 24 h after NaHCO₃ treatment, implying their roles in maintaining ion homeostasis in soybean roots under alkalinity. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis showed "phenylpropanoid biosynthesis" and "phenylalanine metabolism" pathways might participate in soybean response to alkalinity. This study provides a foundation to further investigate the functions of NaHCO₃ stress-responsive genes and the molecular basis of soybean tolerance to alkalinity.

Plant Aquaporins: Genome-Wide Identification, Transcriptomics, Proteomics, and Advanced Analytical Tools

Aquaporins (AQPs) are channel-forming integral membrane proteins that facilitate the movement of water and many other small molecules. Compared to animals, plants contain a much higher number of AQPs in their genome. Homology-based identification of AQPs in sequenced species is feasible because of the high level of conservation of protein sequences across plant species. Genome-wide characterization of AQPs has highlighted several important aspects such as distribution, genetic organization, evolution and conserved features governing solute specificity. From a functional point of view, the understanding of AQP transport system has expanded rapidly with the help of transcriptomics and proteomics data. The efficient analysis of enormous amounts of data generated through omic scale studies has been facilitated through computational advancements. Prediction of protein tertiary structures, pore architecture, cavities, phosphorylation sites, heterodimerization, and co-expression networks has become more sophisticated and accurate with increasing computational tools and pipelines. However, the effectiveness of computational approaches is based on the understanding of physiological and biochemical properties, transport kinetics, solute specificity, molecular interactions, sequence variations, phylogeny and evolution of aquaporins. For this purpose, tools like Xenopus oocyte assays, yeast expression systems, artificial proteoliposomes, and lipid membranes have been efficiently exploited to study the many facets that influence solute transport by AQPs. In the present review, we discuss genome-wide identification of AQPs in plants in relation with recent advancements in analytical tools, and their availability and technological challenges as they apply to AQPs. An exhaustive review of omics resources available for AQP research is also provided in order to optimize their efficient utilization. Finally, a detailed catalog of computational tools and analytical pipelines is offered as a resource for AQP research.