Soybean transcription factor ORFeome associated with drought resistance: a valuable resource to accelerate research on abiotic stress resistance

Overexpression of PtoMYB99 diminishes poplar tolerance to osmotic stress by suppressing ABA and JA biosynthesis

Drought poses a serious challenge to sustained plant growth and crop yields in the context of global climate change. Drought tolerance in poplars and their underlying mechanisms still remain largely unknown. In this article, we investigated the overexpression of PtoMYB99 - both a drought and abscisic acid (ABA) induced gene constraining drought tolerance in poplars (as compared with wild type poplars). First, we found that PtoMYB99-OE lines exhibited increased stomatal opening and conductance, higher transpiration and photosynthetic rates, as well as reduced levels of ABA and jasmonic acid (JA). Second, PtoMYB99-OE lines accumulated more reactive oxygen species (ROS), including H2O2 and O2-, as well as malonaldehyde (MDA), proline, and soluble sugar under osmotic stress; conversely, the activity of antioxidant enzymes (SOD, POD, and CAT), was weakened in the PtoMYB99-OE lines. Third, the expression of ABA biosynthetic genes, PtoNCED3.1 and PtoNCED3.2, as well as JA biosynthetic genes, PtoOPR3.1 and PtoOPR3.2, was significantly reduced in the PtoMYB99-OE lines under both normal conditions and osmotic stress. Based on our results, we conclude that the overexpression of PtoMYB99 compromises tolerance to osmotic stress in poplar. These findings contribute to the understanding of the role of the MYB genes in drought stress and the biosynthesis of ABA and JA.

Overexpression of soybean GmNAC19 and GmGRAB1 enhances root growth and water-deficit stress tolerance in soybean

Soybean (Glycine max) is an important crop in agricultural production where water shortage limits yields in soybean. Root system plays important roles in water-limited environments, but the underlying mechanisms are largely unknown. In our previous study, we produced a RNA-seq dataset generated from roots of soybean at three different growth stages (20-, 30-, and 44-day-old plants). In the present study, we performed a transcriptome analysis of the RNA-seq data to select candidate genes with probable association with root growth and development. Candidate genes were functionally examined in soybean by overexpression of individual genes using intact soybean composite plants with transgenic hairy roots. Root growth and biomass in the transgenic composite plants were significantly increased by overexpression of the GmNAC19 and GmGRAB1 transcriptional factors, showing up to 1.8-fold increase in root length and/or 1.7-fold increase in root fresh/dry weight. Furthermore, greenhouse-grown transgenic composite plants had significantly higher seed yield by about 2-fold than control plants. Expression profiling in different developmental stages and tissues showed that GmNAC19 and GmGRAB1 were most highly expressed in roots, displaying a distinct root-preferential expression. Moreover, we found that under water-deficit conditions, overexpression of GmNAC19 enhanced water stress tolerance in transgenic composite plants. Taken together, these results provide further insights into the agricultural potential of these genes for development of soybean cultivars with improved root growth and enhanced tolerance to water-deficit conditions.

CRISPR/Cas9 mediated gene-editing of GmHdz4 transcription factor enhances drought tolerance in soybean (Glycine max [L.] Merr.)

The HD-Zip transcription factors play a crucial role in plant development, secondary metabolism, and abiotic stress responses, but little is known about HD-Zip I genes in soybean. Here, a homeodomain-leucine zipper gene designated GmHdz4 was isolated. Chimeric soybean plants, GmHdz4 overexpressing (GmHdz4-oe), and gene-editing via CRISPR/Cas9 (gmhdz4) in hairy roots, were generated to examine the GmHdz4 gene response to polyethylene glycol (PEG)-simulated drought stress. Bioinformatic analysis showed GmHdz4 belonged to clade δ, and was closely related to other drought tolerance-related HD-Zip I family genes such as AtHB12, Oshox12, and Gshdz4. The GmHdz4 was located in the plant nucleus and showed transcriptional activation activity by yeast hybrid assay. Quantitative real-time PCR analysis revealed that GmHdz4 expression varied in tissues and was induced by PEG-simulated drought stress. The gmhdz4 showed promoted growth of aboveground parts, and its root system architecture, including the total root length, the root superficial area, and the number of root tips were significantly higher than those of GmHdz4-oe even the non-transgenic line (NT) on root tips number. The better maintenance of turgor pressure by osmolyte accumulation, and the higher activity of antioxidant enzymes to scavenge reactive oxygen species, ultimately suppressed the accumulation of hydrogen peroxide (H₂O₂), superoxide anion (O^2-), and malondialdehyde (MDA), conferring higher drought tolerance in gmhdz4 compared with both GmHdz4-oe and NT. Together, our results provide new insights for future research on the mechanisms by which GmHdz4 gene-editing via CRISPR/Cas9 system could promote drought stress and provide a potential target for molecular breeding in soybean.

The PIN gene family in relic plant L. chinense: Genome-wide identification and gene expression profiling in different organizations and abiotic stress responses

The auxin efflux carrier PIN-FORMED (PIN) proteins are required for the polar transport of auxin between cells through their asymmetric distribution on the plasma membrane, thus mediating the differential distribution of auxin in plants, finally, affecting plant growth and developmental processes. In this study, 11 LcPIN genes were identified. The structural characteristics and evolutionary status of LcPIN genes were thoroughly investigated and interpreted combining physicochemical property analysis, evolutionary analysis, gene structure analysis, chromosomal localization, etc. Multi-species protein sequence analysis showed that angiosperm PIN genes have strong purification options and some functional sites were predicted about PIN protein polarity, trafficking and activity in L. chinense. Further qRT-PCR and transcriptome data analysis indicated that the long LcPINs have highly expressed from globular embryo to plantlet, and the LcPIN6a started upregulated in cotyledon embryo. The LcPIN3 and LcPIN6a are both highly expressed during the development of stamens and petals and the expression of LcPIN2 is related to root elongation, suggesting that they may play an important role in these processes. Experiment data indicates that LcPIN5 and LcPIN8 might play a key role in auxin transport in Liriodendron stems and leaves under abiotic stress. Analyzed the response of LcPIN genes to abiotic stress and as a basis for uncovering the biological role of LcPIN genes in development and adaption to adverse environments. This study provides a foundation for further genetic and functional analyses.

Auxin Profiling and GmPIN Expression in Phytophthora sojae-Soybean Root Interactions

Auxin (indole-3-acetic acid, IAA) has been implicated as a susceptibility factor in both beneficial and pathogenic molecular plant-microbe interactions. Previous studies have identified a large number of auxin-related genes underlying quantitative disease resistance loci (QDRLs) for Phytophthora sojae. Thus, we hypothesized that auxin may be involved the P. sojae-soybean interaction. The levels of IAA and related metabolites were measured in mycelia and media supernatant as well as in mock and inoculated soybean roots in a time course assay. The expression of 11 soybean Pin-formed (GmPIN) auxin efflux transporter genes was also examined. Tryptophan, an auxin precursor, was detected in the P. sojae mycelia and media supernatant. During colonization of roots, levels of IAA and related metabolites were significantly higher in both moderately resistant Conrad and moderately susceptible Sloan inoculated roots compared with mock controls at 48 h postinoculation (hpi) in one experiment and at 72 hpi in a second, with Sloan accumulating higher levels of the auxin catabolite IAA-Ala than Conrad. Additionally, one GmPIN at 24 hpi, one at 48 hpi, and three at 72 hpi had higher expression in inoculated compared with the mock control roots in Conrad. The ability of resistant cultivars to cope with auxin accumulation may play an important role in quantitative disease resistance. Levels of jasmonic acid (JA), another plant hormone associated with defense responses, were also higher in inoculated roots at these same time points, suggesting that JA also plays a role during the later stages of infection.

Salt-responsive transcriptome analysis of triticale reveals candidate genes involved in the key metabolic pathway in response to salt stress

Triticale is tolerant of many environmental stresses, especially highly resistant to salt stress. However, the molecular regulatory mechanism of triticale seedlings under salt stress conditions is still unclear so far. In this study, a salt-responsive transcriptome analysis was conducted to identify candidate genes or transcription factors related to salt tolerance in triticale. The root of salt-tolerant triticale cultivars TW004 with salt-treated and non-salt stress at different time points were sampled and subjected to de novo transcriptome sequencing. Total 877,858 uniquely assembled transcripts were identified and most contigs were annotated in public databases including nr, GO, KEGG, eggNOG, Swiss-Prot and Pfam. 59,280, 49,345, and 85,922 differentially expressed uniquely assembled transcripts between salt treated and control triticale root samples at three different time points (C12_vs_T12, C24_vs_T24, and C48_vs_T48) were identified, respectively. Expression profile and functional enrichment analysis of DEGs found that some DEGs were significantly enriched in metabolic pathways related to salt tolerance, such as reduction–oxidation pathways, starch and sucrose metabolism. In addition, several transcription factor families that may be associated with salt tolerance were also identified, including AP2/ERF, NAC, bHLH, WRKY and MYB. Furthermore, 14 DEGs were selected to validate the transcriptome profiles via quantitative RT-PCR. In conclusion, these results provide a foundation for further researches on the regulatory mechanism of triticale seedlings adaptation to salt stress in the future.

Transcriptomic profile analysis of the halophyte Suaeda rigida response and tolerance under NaCl stress

Suaeda rigida is a lignified, true haplotype that predominantly grows in the Tarim basin, China. It has significant economic and ecological value. Herein, with aim to determine the genes associated with salt tolerance, transcriptome sequencing was performed on its stem, leaves and root over three set NaCl gradients regimens at treatment intervals of 3 h and 5 days. From our findings, we identified 829,095 unigenes, with 331,394 being successfully matched to at least one annotation database. In roots, under 3 h treatment, no up-regulated DEGs were identified in 100 and 500 mM NaCl treated samples. Under 5 days treatment, 97, 60 and 242 up-regulated DEGs were identified in 100, 300, 500 mM NaCl treated samples, respectively. We identified 50, 22 and 255 down-regulated DEGs in 100, 300, 500 mM NaCl treated samples, respectively. GO biological process enrichment analysis established that down-regulated DEGs were associated with nitrogen compound transport, organic substance transport and intracellular protein transport while the up-regulated genes were enriched in cell wall biogenesis, such as plant-type cell wall biogenesis, cell wall assembly, extracellular matrix organization and plant-type cell wall organization. These findings provide valuable knowledge on genes associated with salt tolerance of Suaeda rigida, and can be applied in other downstream haplotype studies.

ABAS1 from soybean is a 1R-subtype MYB transcriptional repressor that enhances ABA sensitivity

Transcription factors (TFs) help plants respond to environmental stresses by regulating gene expression. Up till now, studies on the MYB family of TFs have mainly focused on the highly abundant R2R3-subtype. While the less well-known 1R-subtype has been generally shown to enhance abscisic acid (ABA) sensitivity by acting as transcriptional activators, the mechanisms of their functions are unclear. Here we identified an ABA sensitivity-associated gene from soybean, ABA-Sensitive 1 (GmABAS1), of the 1R-subtype of MYB. Using the GFP-GmABAS1 fusion protein, we demonstrated that GmABAS1 is localized in the nucleus, and with yeast reporter systems, we showed that it is a transcriptional repressor. We then identified the target gene of GmABAS1 to be Glyma.01G060300, an annotated ABI five-binding protein 3 and showed that GmABAS1 binds to the promoter of Glyma.01G060300 both in vitro and in vivo. Furthermore, Glyma.01G060300 and GmABAS1 exhibited reciprocal expression patterns under osmotic stress, inferring that GmABAS1 is a transcriptional repressor of Glyma.01G060300. As a further confirmation, AtAFP2, an orthologue of Glyma.01G060300, was down-regulated in GmABAS1-transgenic Arabidopsis thaliana, enhancing the plant's sensitivity to ABA. This is the first time a 1R-subtype of MYB from soybean has been reported to enhance ABA sensitivity by acting as a transcriptional repressor.

Diel pattern of circadian clock and storage protein gene expression in leaves and during seed filling in cowpea (Vigna unguiculata)

BACKGROUND

Cowpea (Vigna unguiculata) is an important source of protein supply for animal and human nutrition. The major storage globulins VICILIN and LEGUMIN (LEG) are synthesized from several genes including LEGA, LEGB, LEGJ and CVC (CONVICILIN). The current hypothesis is that the plant circadian core clock genes are conserved in a wide array of species and that primary metabolism is to a large extent controlled by the plant circadian clock. Our aim was to investigate a possible link between gene expression of storage proteins and the circadian clock.

RESULTS

We identified cowpea orthologues of the core clock genes VunLHY, VunTOC1, VunGI and VunELF3, the protein storage genes VunLEG, VunLEGJ, and VunCVC as well as nine candidate reference genes used in RT-PCR. ELONGATION FACTOR 1-A (ELF1A) resulted the most suitable reference gene. The clock genes VunELF3, VunGI, VunTOC1 and VunLHY showed a rhythmic expression profile in leaves with a typical evening/night and morning/midday phased expression. The diel patterns were not completely robust and only VungGI and VungELF3 retained a rhythmic pattern under free running conditions of darkness. Under field conditions, rhythmicity and phasing apparently faded during early pod and seed development and was regained in ripening pods for VunTOC1 and VunLHY. Mature seeds showed a rhythmic expression of VunGI resembling leaf tissue under controlled growth chamber conditions. Comparing time windows during developmental stages we found that VunCVC and VunLEG were significantly down regulated during the night in mature pods as compared to intermediate ripe pods, while changes in seeds were non-significant due to high variance. The rhythmic expression under field conditions was lost under growth chamber conditions.

CONCLUSIONS

The core clock gene network is conserved in cowpea leaves showing a robust diel expression pattern except VunELF3 under growth chamber conditions. There appears to be a clock transcriptional reprogramming in pods and seeds compared to leaves. Storage protein deposition may be circadian regulated under field conditions but the strong environmental signals are not met under artificial growth conditions. Diel expression pattern in field conditions may result in better usage of energy for protein storage.

The PIN gene family in cotton (Gossypium hirsutum): genome-wide identification and gene expression analyses during root development and abiotic stress responses

Background

Cell elongation and expansion are significant contributors to plant growth and morphogenesis, and are often regulated by environmental cues and endogenous hormones. Auxin is one of the most important phytohormones involved in the regulation of plant growth and development and plays key roles in plant cell expansion and elongation. Cotton fiber cells are a model system for studying cell elongation due to their large size. Cotton is also the world’s most utilized crop for the production of natural fibers for textile and garment industries, and targeted expression of the IAA biosynthetic gene iaaM increased cotton fiber initiation. Polar auxin transport, mediated by PIN and AUX/LAX proteins, plays a central role in the control of auxin distribution. However, very limited information about PIN-FORMED (PIN) efflux carriers in cotton is known.

Results

In this study, 17 PIN-FORMED (PIN) efflux carrier family members were identified in the Gossypium hirsutum (G. hirsutum) genome. We found that PIN1–3 and PIN2 genes originated from the At subgenome were highly expressed in roots. Additionally, evaluation of gene expression patterns indicated that PIN genes are differentially induced by various abiotic stresses. Furthermore, we found that the majority of cotton PIN genes contained auxin (AuxREs) and salicylic acid (SA) responsive elements in their promoter regions were significantly up-regulated by exogenous hormone treatment.

Conclusions

Our results provide a comprehensive analysis of the PIN gene family in G. hirsutum, including phylogenetic relationships, chromosomal locations, and gene expression and gene duplication analyses. This study sheds light on the precise roles of PIN genes in cotton root development and in adaption to stress responses.

Electronic supplementary material

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

Comprehensive Analysis of the Soybean (Glycine max) GmLAX Auxin Transporter Gene Family

The phytohormone auxin plays a critical role in regulation of plant growth and development as well as plant responses to abiotic stresses. This is mainly achieved through its uneven distribution in plant via a polar auxin transport process. Auxin transporters are major players in polar auxin transport. The AUXIN RESISTENT 1/LIKE AUX1 (AUX/LAX) auxin influx carriers belong to the amino acid permease family of proton-driven transporters and function in the uptake of indole-3-acetic acid (IAA). In this study, genome-wide comprehensive analysis of the soybean AUX/LAX (GmLAX) gene family, including phylogenic relationships, chromosome localization, and gene structure, was carried out. A total of 15 GmLAX genes, including seven duplicated gene pairs, were identified in the soybean genome. They were distributed on 10 chromosomes. Despite their higher percentage identities at the protein level, GmLAXs exhibited versatile tissue-specific expression patterns, indicating coordinated functioning during plant growth and development. Most GmLAXs were responsive to drought and dehydration stresses and auxin and abscisic acid (ABA) stimuli, in a tissue- and/or time point- sensitive mode. Several GmLAX members were involved in responding to salt stress. Sequence analysis revealed that promoters of GmLAXs contained different combinations of stress-related cis-regulatory elements. These studies suggest that the soybean GmLAXs were under control of a very complex regulatory network, responding to various internal and external signals. This study helps to identity candidate GmLAXs for further analysis of their roles in soybean development and adaption to adverse environments.

Salt-Responsive Transcriptome Profiling of Suaeda glauca via RNA Sequencing

Background

Suaeda glauca, a succulent halophyte of the Chenopodiaceae family, is widely distributed in coastal areas of China. Suaeda glauca is highly resistant to salt and alkali stresses. In the present study, the salt-responsive transcriptome of Suaeda glauca was analyzed to identify genes involved in salt tolerance and study halophilic mechanisms in this halophyte.

Results

Illumina HiSeq 2500 was used to sequence cDNA libraries from salt-treated and control samples with three replicates each treatment. De novo assembly of the six transcriptomes identified 75,445 unigenes. A total of 23,901 (31.68%) unigenes were annotated. Compared with transcriptomes from the three salt-treated and three salt-free samples, 231 differentially expressed genes (DEGs) were detected (including 130 up-regulated genes and 101 down-regulated genes), and 195 unigenes were functionally annotated. Based on the Gene Ontology (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) classifications of the DEGs, more attention should be paid to transcripts associated with signal transduction, transporters, the cell wall and growth, defense metabolism and transcription factors involved in salt tolerance.

Conclusions

This report provides a genome-wide transcriptional analysis of a halophyte, Suaeda glauca, under salt stress. Further studies of the genetic basis of salt tolerance in halophytes are warranted.

Comprehensive Analysis and Expression Profiling of the OsLAX and OsABCB Auxin Transporter Gene Families in Rice (Oryza sativa) under Phytohormone Stimuli and Abiotic Stresses

The plant hormone auxin regulates many aspects of plant growth and developmental processes. Auxin gradient is formed in plant as a result of polar auxin transportation by three types of auxin transporters such as OsLAX, OsPIN, and OsABCB. We report here the analysis of two rice auxin transporter gene families, OsLAX and OsABCB, using bioinformatics tools, publicly accessible microarray data, and quantitative RT-PCR. There are 5 putative OsLAXs and 22 putative OsABCBs in rice genome, which were mapped on 8 chromosomes. The exon-intron structure of OsLAX genes and properties of deduced proteins were relatively conserved within grass family, while that of OsABCB genes varied greatly. Both constitutive and organ/tissue specific expression patterns were observed in OsLAXs and OsABCBs. Analysis of evolutionarily closely related "gene pairs" together with organ/tissue specific expression revealed possible "function gaining" and "function losing" events during rice evolution. Most OsLAX and OsABCB genes were regulated by drought and salt stress, as well as hormonal stimuli [auxin and Abscisic Acid (ABA)], which suggests extensive crosstalk between abiotic stresses and hormone signaling pathways. The existence of large number of auxin and stress related cis-regulatory elements in promoter regions might account for their massive responsiveness of these genes to these environmental stimuli, indicating complexity of regulatory networks involved in various developmental and physiological processes. The comprehensive analysis of OsLAX and OsABCB auxin transporter genes in this study would be helpful for understanding the biological significance of these gene families in hormone signaling and adaptation of rice plants to unfavorable environments.

Genome-wide analysis and expression profiling of the PIN auxin transporter gene family in soybean (Glycine max)

BACKGROUND

The plant phytohormone auxin controls many aspects of plant growth and development, which largely depends on its uneven distribution in plant tissues. Transmembrane proteins of the PIN family are auxin efflux facilitators. They play a key role in polar auxin transport and are associated with auxin asymmetrical distribution in plants. PIN genes have been characterized in several plant species, while comprehensive analysis of this gene family in soybean has not been reported yet.

RESULTS

In this study, twenty-three members of the PIN gene family were identified in the soybean genome through homology searches. Analysis of chromosome distribution and phylogenetic relationships of the soybean PIN genes indicated nine pairs of duplicated genes and a legume specific subfamily. Organ/tissue expression patterns and promoter activity assays of the soybean PINs suggested redundant functions for most duplicated genes and complementary and tissue-specific functions during development for non-duplicated genes. The soybean PIN genes were differentially regulated by various abiotic stresses and phytohormone stimuli, implying crosstalk between auxin and abiotic stress signaling pathways. This was further supported by the altered auxin distribution under these conditions as revealed by DR5::GUS transgenic soybean hairy root. Our data indicates that GmPIN9, a legume-specific PIN gene, which was responsive to several abiotic stresses, might play a role in auxin re-distribution in soybean root under abiotic stress conditions.

CONCLUSIONS

This study provided the first comprehensive analysis of the soybean PIN gene family. Information on phylogenetic relationships, gene structure, protein profiles and expression profiles of the soybean PIN genes in different tissues and under various abiotic stress treatments helps to identity candidates with potential roles in specific developmental processes and/or environmental stress conditions. Our study advances our understanding of plant responses to abiotic stresses and serves as a basis for uncovering the biological role of PIN genes in soybean development and adaption to adverse environments.

Soybean transcription factor ORFeome associated with drought resistance: a valuable resource to accelerate research on abiotic stress resistance

BACKGROUND

Whole genome sequencing provides the most comprehensive collection of an organism's genetic information. The availability of complete genome sequences is expected to dramatically deliver a high impact on biology. However, to achieve this impact in the area of crop improvement, significant efforts are still required on functional genomics, including the areas of gene annotation, cloning, expression profiling, and functional validation.

RESULTS

Here we report our efforts in generating the first transcription factor (TF) open reading frame (ORF)eome resource associated with drought resistance in soybean (Glycine max), a major oil/protein crop grown worldwide. This study provides a highly annotated soybean TF-ORFeome associated with drought resistance. It contains information from experimentally verified protein-coding sequences (CDS), expression profiling under several abiotic stresses (drought, salinity, dehydration and ABA), and computationally predicted protein subcellular localization and cis-regulatory elements (CREs) analysis. All the information is available to plant researchers through a freely accessible and user-friendly database, Soybean Knowledge Base (SoyKB).

CONCLUSIONS

The soybean TF-ORFeome provides a valuable public resource for functional genomics studies, especially in the area of plant abiotic stresses. It will accelerate findings in the areas of abiotic stresses and lead to the generation of crops with enhanced resistance to multiple stresses.