Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes

Co-orthologues of ribosome biogenesis factors in A. thaliana are differentially regulated by transcription factors

Different genes coding for one ribosome biogenesis factor are differentially expressed and are likely under the control of distinct transcription factors, which contributes to the regulatory space for ribosome maturation. Maturation of ribosomes including rRNA processing and modification, rRNA folding and ribosome protein association requires the function of many ribosome biogenesis factors (RBFs). Recent studies document plant-specific variations of the generally conserved process of ribosome biogenesis. For instance, distinct rRNA maturation pathways and intermediates have been identified, the existence of plant specific RBFs has been proposed and several RBFs are encoded by multiple genes. The latter in combination with the discussed ribosome heterogeneity points to a possible function of the different proteins representing one RBF in diversification of ribosomal compositions. Such factor-based regulation would require a differential regulation of their expression, may be even controlled by different transcription factors. We analyzed the expression profiles of genes coding for putative RBFs and transcription factors. Most of the genes coding for RBFs are expressed in a comparable manner, while different genes coding for a single RBF are often differentially expressed. Based on a selected set of genes we document a function of the transcription factors AtMYC1, AtMYC2, AtbHLH105 and AtMYB26 on the regulation of different RBFs. Moreover, on the example of the RBFs LSG1 and BRX1, both encoded by two genes, we give a first hint on a differential transcription factor dependence of expression. Consistent with this observation, the phenotypic analysis of RBF mutants suggests a relation between LSG1-1 and BRX1-1 expression and the transcription factor MYC1. In summary, we propose that the multiple genes coding for one RBF are required to enlarge the regulatory space for ribosome biogenesis.

Identification of a Novel Melon Transcription Factor CmNAC60 as a Potential Regulator of Leaf Senescence

NAC transcription factors (TFs) play important roles in plants’ responses to abiotic stresses and developmental processes, including leaf senescence. Oriental melon (Cucumis melo var. makuwa Makino) is an important vegetable crop in China and eastern Asia countries. However, little is known about the functions of the melon NAC family members. In this study, a phylogenetic tree was constructed to show that CmNAC60 and the senescence regulator AtNAP were in the same cluster, which implied that CmNAC60 might be a NAC related to leaf senescence. The expression analysis of CmNAC60 in different melon organs showed that the expression of CmNAC60 was highest in the male flowers and lowest in the hypocotyl. In addition, the expression level of CmNAC60 in the senescing leaves was significantly higher than in the non-senescing leaves. Similarly, the expression level of CmNAC60 in the dark-treated leaves was significantly higher than in the untreated leaves. Furthermore, the subcellular localization and transcriptional activation assays indicated that CmNAC60 was a nucleus localized NAC transcription factor with a C-terminal transactivation domain. An analysis of the tissue specific expression showed that the promoter of CmNAC60 may contain cis-acting regulatory elements responsive to leaf senescence. CmNAC60 overexpressing lines of Arabidopsis showed a precocious senescence compared with the wild type (WT). Collectively, our results showed that CmNAC60 was associated with leaf senescence, and could be potentially utilized in molecular breeding to improve melon yield or to extend the postharvest shelf life by delaying leaf senescence.

Comprehensive analyses of ZFP gene family and characterization of expression profiles during plant hormone response in cotton

BACKGROUND

Zinc finger proteins (ZFPs) containing only a single zinc finger domain play important roles in the regulation of plant growth and development, as well as in biotic and abiotic stress responses. To date, the evolutionary history and functions of the ZFP gene family have not been identified in cotton.

RESULTS

In this paper, we identified 29 ZFP genes in Gossypium hirsutum. This gene family was divided into seven subfamilies, 22 of which were distributed over 17 chromosomes. Bioinformatic analysis revealed that 20 GhZFP genes originated from whole genome duplications and two originated from dispersed duplication events, indicating that whole genome duplication is the main force in the expansion of the GhZFP gene family. Most GhZFP8 subfamily genes, except for GhZFP8-3, were highly expressed during fiber cell growth, and were induced by brassinosteroids in vitro. Furthermore, we found that a large number of GhZFP genes contained gibberellic acid responsive elements, auxin responsive elements, and E-box elements in their promoter regions. Exogenous application of these hormones significantly stimulated the expression of these genes.

CONCLUSIONS

Our findings reveal that GhZFP8 genes are involved in cotton fiber development and widely induced by auxin, gibberellin and BR, which provides a foundation for the identification of more downstream genes with potential roles in phytohormone stimuli, and a basis for breeding better cotton varieties in the future.

The NAC Protein from Tamarix hispida, ThNAC7, Confers Salt and Osmotic Stress Tolerance by Increasing Reactive Oxygen Species Scavenging Capability

Plant specific NAC (NAM, ATAF1/2 and CUC2) transcription factors (TFs) play important roles in response to abiotic stress. In this study, we identified and characterized a NAC protein, ThNAC7, from Tamarix hispida. ThNAC7 is a nuclear localized protein and has transcriptional activation activity. ThNAC7 expression was markedly induced by salt and osmotic stresses. Transiently transformed T. hispida seedlings overexpressing ThNAC7 (OE) or with RNA interference (RNAi) silenced ThNAC7 were generated to investigate abiotic stress tolerance via the gain- and loss- of function. Overexpressing ThNAC7 showed an increased reactive oxygen species (ROS) scavenging capabilities and proline content, which was accomplished by enhancing the activities of superoxide dismutase (SOD) and peroxidase (POD) in transiently transformed T. hispida and stably transformed Arabidopsis plants. Additionally, ThNAC7 activated these physiological changes by regulating the transcription level of P5CS, SOD and PODgenes. RNA-sequencing (RNA-seq) comparison between wild-type and ThNAC7-transformed Arabidopsis showed that more than 40 known salt tolerance genes might regulated by ThNAC7, including stress tolerance-related genes and TF genes. The results indicated that ThNAC7 induces the transcription level of genes associated with stress tolerance to enhance salt and osmotic stress tolerance via an increase in osmotic potential and enhanced ROS scavenging.

Large-Scale Genome Comparison Based on Cumulative Fourier Power and Phase Spectra: Central Moment and Covariance Vector

Genome comparison is a vital research area of bioinformatics. For large-scale genome comparisons, the Multiple Sequence Alignment (MSA) methods have been impractical to use due to its algorithmic complexity. In this study, we propose a novel alignment-free method based on the one-to-one correspondence between a DNA sequence and its complete central moment vector of the cumulative Fourier power and phase spectra. In addition, the covariance between the four nucleotides in the power and phase spectra is included. We use the cumulative Fourier power and phase spectra to define a 28-dimensional vector for each DNA sequence. Euclidean distances between the vectors can measure the dissimilarity between DNA sequences. We perform testing with datasets of different sizes and types including simulated DNA sequences, exon-intron and complete genomes. The results show that our method is more accurate and efficient for performing hierarchical clustering than other alignment-free methods and MSA methods.

Understanding the Role of the WRKY Gene Family under Stress Conditions in Pigeonpea (Cajanus Cajan L.)

Pigeonpea (Cajanus cajan L.), a protein-rich legume, is a major food component of the daily diet for residents in semi-arid tropical regions of the word. Pigeonpea is also known for its high level of tolerance against biotic and abiotic stresses. In this regard, understanding the genes involved in stress tolerance has great importance. In the present study, identification, and characterization of WRKY, a large transcription factor gene family involved in numerous biological processes like seed germination, metabolism, plant growth, biotic and abiotic stress responses was performed in pigeonpea. A total of 94 WRKY genes identified in the pigeonpea genome were extensively characterized for gene-structures, localizations, phylogenetic distribution, conserved motif organizations, and functional annotation. Phylogenetic analysis revealed three major groups (I, II, and III) of pigeonpea WRKY genes. Subsequently, expression profiling of 94 CcWRKY genes across different tissues like root, nodule, stem, petiole, petal, sepal, shoot apical meristem (SAM), mature pod, and mature seed retrieved from the available RNAseq data identified tissue-specific WRKY genes with preferential expression in the vegetative and reproductive stages. Gene co-expression networks identified four WRKY genes at the center of maximum interaction which may play a key role in the entire WRKY regulations. Furthermore, quantitative real-time polymerase chain reaction (qRT-PCR) expression analysis of WRKY genes in root and leaf tissue samples from plants under drought and salinity stress identified differentially expressed WRKY genes. The study will be helpful to understand the evolution, regulation, and distribution of the WRKY gene family, and additional exploration for the development of stress tolerance cultivars in pigeonpea and other legumes crops.

Exploring regulatory networks in plants: transcription factors of starch metabolism

Biological networks are complex (non-linear), redundant (cyclic) and compartmentalized at the subcellular level. Rational manipulation of plant metabolism may have failed due to inherent difficulties of a comprehensive understanding of regulatory loops. We first need to identify key factors controlling the regulatory loops of primary metabolism. The paradigms of plant networks are revised in order to highlight the differences between metabolic and transcriptional networks. Comparison between animal and plant transcription factors (TFs) reveal some important differences. Plant transcriptional networks function at a lower hierarchy compared to animal regulatory networks. Plant genomes contain more TFs than animal genomes, but plant proteins are smaller and have less domains as animal proteins which are often multifunctional. We briefly summarize mutant analysis and co-expression results pinpointing some TFs regulating starch enzymes in plants. Detailed information is provided about biochemical reactions, TFs and cis regulatory motifs involved in sucrose-starch metabolism, in both source and sink tissues. Examples about coordinated responses to hormones and environmental cues in different tissues and species are listed. Further advancements require combined data from single-cell transcriptomic and metabolomic approaches. Cell fractionation and subcellular inspection may provide valuable insights. We propose that shuffling of promoter elements might be a promising strategy to improve in the near future starch content, crop yield or food quality.

Transcription factors involved in abiotic stress responses in Maize (Zea mays L.) and their roles in enhanced productivity in the post genomics era

BACKGROUND

Maize (Zea mays L.) is a principal cereal crop cultivated worldwide for human food, animal feed, and more recently as a source of biofuel. However, as a direct consequence of water insufficiency and climate change, frequent occurrences of both biotic and abiotic stresses have been reported in various regions around the world, and recently, this has become a constant threat in increasing global maize yields. Plants respond to abiotic stresses by utilizing the activities of transcription factors (TFs), which are families of genes coding for specific TF proteins. TF target genes form a regulon that is involved in the repression/activation of genes associated with abiotic stress responses. Therefore, it is of utmost importance to have a systematic study on each TF family, the downstream target genes they regulate, and the specific TF genes involved in multiple abiotic stress responses in maize and other staple crops.

METHOD

In this review, the main TF families, the specific TF genes and their regulons that are involved in abiotic stress regulation will be briefly discussed. Great emphasis will be given on maize abiotic stress improvement throughout this review, although other examples from different plants like rice, Arabidopsis, wheat, and barley will be used.

RESULTS

We have described in detail the main TF families in maize that take part in abiotic stress responses together with their regulons. Furthermore, we have also briefly described the utilization of high-efficiency technologies in the study and characterization of TFs involved in the abiotic stress regulatory networks in plants with an emphasis on increasing maize production. Examples of these technologies include next-generation sequencing, microarray analysis, machine learning, and RNA-Seq.

CONCLUSION

In conclusion, it is expected that all the information provided in this review will in time contribute to the use of TF genes in the research, breeding, and development of new abiotic stress tolerant maize cultivars.

Functional characterization of LkERF-B2 for improved salt tolerance ability in Arabidopsis thaliana

The ethylene response factors have been reported to play critical roles in developmental and environmental responses in plants. In the present study, an ERF transcription factor gene was aimed to be identified from Larix kaempferi. Molecular characteristics and function of this gene were further explored. The result showed that a 1344 bp ERF transcription factor gene containing initiation and termination codon was obtained by RT-PCR and named LkERF-B2. LkERF-B2 gene encoded 447 amino acids containing a typical AP2/ERF domain. Alignment of predicted amino acid sequence of LkERF-B2 in various plant species showed that this ERF transcription factor was highly homologous (79.0%) with that of Picea sitchensi. To elucidate the function of LkERF-B2, LkERF-B2 overexpression vector was successfully constructed and transformed to Arabidopsis thaliana via dip flower. Compared with control plant, LkERF-B2 overexpressed transgenic A. thaliana showed a significantly higher survival rate under cold, heat, NaCl and drought stresses. NaCl stress analysis revealed that control and transgenic Arabidopsis were both flowering earlier under 100 and 150 mM/L NaCl treatment. While under 200-300 mM/L NaCl treatment, the growth of control plant was significantly inhibited compared with transgenic A. thaliana. Salt injury rate and salt injury index of transgenic Arabidopsis were lower than those of the control. Further investigation showed that transgenic Arabidopsis exhibited much higher content of chloroplast pigments under different NaCl concentration. Meanwhile, the activity of SOD and POD was also enhanced in transgenic A. thaliana. These results suggested that LkERF-B2 was a key transcription factor and could lead to enhanced salt stress tolerance.

Comparative genome analysis of the SPL gene family reveals novel evolutionary features in maize

SPLs are plant-specific transcription factors that play important regulatory roles in plant growth and development. Systematic analysis of the SPL family has been performed in numerous plants, such as Arabidopsis, rice, and Populus. However, no comparative analysis has been performed across different species to examine evolutionary features. In this study, we present a comparative analysis of SPLs in different species. The results showed that 84 SPLs of the four species can be divided into six groups according to phylogeny. We found that most of the SPL-containing regions in maize showed extensive conservation with duplicated regions of rice and sorghum. A gene duplication analysis in maize indicated that ZmSPLs showed a significant excess of segmental duplication. The Ka/Ks analysis indicated that 9 out of 18 duplicated pairs in maize experienced positive selection, while SPL gene pairs of rice and sorghum mainly evolved under purifying selection, suggesting novel evolutionary features for ZmSPLs. The 31 ZmSPLs were further analyzed by describing their gene structure, phylogenetic relationships, chromosomal location, and expression, Among the ZmSPLs, 13 were predicated to be targeted by miR156s and involved in drought stress response. These results provide the foundation for future functional analyses of ZmSPLs.

Genome-wide characterization of aspartic protease (AP) gene family in Populus trichocarpa and identification of the potential PtAPs involved in wood formation

BACKGROUND

Aspartic protease (AP) is one of four large proteolytic enzyme families that are involved in plant growth and development. Little is known about the AP gene family in tree species, although it has been characterized in Arabidopsis, rice and grape. The AP genes that are involved in tree wood formation remain to be determined.

RESULTS

A total of 67 AP genes were identified in Populus trichocarpa (PtAP) and classified into three categories (A, B and C). Chromosome mapping analysis revealed that two-thirds of the PtAP genes were located in genome duplication blocks, indicating the expansion of the AP family by segmental duplications in Populus. The microarray data from the Populus eFP browser demonstrated that PtAP genes had diversified tissue expression patterns. Semi-qRT-PCR analysis further determined that more than 10 PtAPs were highly or preferentially expressed in the developing xylem. When the involvement of the PtAPs in wood formation became the focus, many SCW-related cis-elements were found in the promoters of these PtAPs. Based on PtAP_promoter::GUS techniques, the activities of PtAP66 promoters were observed only in fiber cells, not in the vessels of stems as the xylem and leaf veins developed in the transgenic Populus tree, and strong GUS signals were detected in interfascicular fiber cells, roots, anthers and sepals of PtAP17_promoter::GUS transgenic plants. Intensive GUS activities in various secondary tissues implied that PtAP66 and PtAP17 could function in wood formation. In addition, most of the PtAP proteins were predicted to contain N- and (or) O-glycosylation sites, and the integration of PNGase F digestion and western blotting revealed that the PtAP17 and PtAP66 proteins were N-glycosylated in Populus.

CONCLUSIONS

Comprehensive characterization of the PtAP genes suggests their functional diversity during Populus growth and development. Our findings provide an overall understanding of the AP gene family in trees and establish a better foundation to further describe the roles of PtAPs in wood formation.

Genome-wide identification, phylogeny, evolutionary expansion and expression analyses of bZIP transcription factor family in tartaty buckwheat

Background

In reported plants, the bZIP family is one of the largest transcription factor families. bZIP genes play roles in the light signal, seed maturation, flower development, cell elongation, seed accumulation protein, abiotic and biological stress and other biological processes. While, no detailed identification and genome-wide analysis of bZIP family genes in Fagopyum talaricum (tartary buckwheat) has previously been published. The recently reported genome sequence of tartary buckwheat provides theoretical basis for us to study and discuss the characteristics and expression of bZIP genes in tartary buckwheat based on the whole genome.

Results

In this study, 96 FtbZIP genes named from FtbZIP1 to FtbZIP96 were identified and divided into 11 subfamilies according to their genetic relationship with 70 bZIPs of A. thaliana. FtbZIP genes are not evenly distributed on the chromosomes, and we found tandem and segmental duplication events of FtbZIP genes on 8 tartary buckwheat chromosomes. According to the results of gene and motif composition, FtbZIP located in the same group contained analogous intron/exon organizations and motif composition. By qRT-PCR, we quantified the expression of FtbZIP members in stem, root, leaf, fruit, and flower and during fruit development. Exogenous ABA treatment increased the weight of tartary buckwheat fruit and changed the expressions of FtbZIP genes in group A.

Conclusions

Through our study, we identified 96 FtbZIP genes in tartary buckwheat and synthetically further analyzed the structure composition, evolution analysis and expression pattern of FtbZIP proteins. The expression pattern indicates that FtbZIP is important in the course of plant growth and development of tartary buckwheat. Through comprehensively analyzing fruit weight and FtbZIP genes expression after ABA treatment and endogenous ABA content of tartary buckwheat fruit, ABA may regulate downstream gene expression by regulating the expression of FtPinG0003523300.01 and FtPinG0003196200.01, thus indirectly affecting the fruit development of tartary buckwheat. This will help us to further study the function of FtbZIP genes in the tartary buckwheat growth and improve the fruit of tartary buckwheat.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5882-z) contains supplementary material, which is available to authorized users.

A survey of transcriptome complexity using PacBio single-molecule real-time analysis combined with Illumina RNA sequencing for a better understanding of ricinoleic acid biosynthesis in Ricinus communis

Background

Ricinus communis is a highly economically valuable oil crop plant from the spurge family, Euphorbiaceae. However, the available reference genomes are incomplete and to date studies on ricinoleic acid biosynthesis at the transcriptional level are limited.

Results

In this study, we combined PacBio single-molecule long read isoform and Illumina RNA sequencing to identify the alternative splicing (AS) events, novel isoforms, fusion genes, long non-coding RNAs (lncRNAs) and alternative polyadenylation (APA) sites to unveil the transcriptomic complexity of castor beans and identify critical genes related to ricinoleic acid biosynthesis. Here, we identified 11,285 AS-variants distributed in 21,448 novel genes and detected 520 fusion genes, 320 lncRNAs and 9511 (APA-sites). Furthermore, a total of 6067, 5983 and 4058 differentially expressed genes between developing beans of the R. communis lines 349 and 1115 with extremely different oil content were identified at 7, 14 and 21 days after flowering, respectively. Specifically, 14, 18 and 11 DEGs were annotated encoding key enzymes related to ricinoleic acid biosynthesis reflecting the higher castor oil content of 1115 compared than 349. Quantitative real-time RT-PCR further validated fifteen of these DEGs at three-time points.

Conclusion

Our results significantly improved the existed gene models of R. communis, and a putative model of key genes was built to show the differences between strains 349 and 1115, illustrating the molecular mechanism of castor oil biosynthesis. A multi-transcriptome database and candidate genes were provided to further improve the level of ricinoleic acid in transgenic crops.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5832-9) contains supplementary material, which is available to authorized users.

Transcriptional Regulation of Abscission Zones

Precise and timely regulation of organ separation from the parent plant (abscission) is consequential to improvement of crop productivity as it influences both the timing of harvest and fruit quality. Abscission is tightly associated with plant fitness as unwanted organs (petals, sepals, filaments) are shed after fertilization while seeds, fruits, and leaves are cast off as means of reproductive success or in response to abiotic/biotic stresses. Floral organ abscission in Arabidopsis has been a useful model to elucidate the molecular mechanisms that underlie the separation processes, and multiple abscission signals associated with the activation and downstream pathways have been uncovered. Concomitantly, large-scale analyses of omics studies in diverse abscission systems of various plants have added valuable insights into the abscission process. The results suggest that there are common molecular events linked to the biosynthesis of a new extracellular matrix as well as cell wall disassembly. Comparative analysis between Arabidopsis and soybean abscission systems has revealed shared and yet disparate regulatory modules that affect the separation processes. In this review, we discuss our current understanding of the transcriptional regulation of abscission in several different plants that has improved on the previously proposed four-phased model of organ separation.

Genome-Wide Analysis of TCP Family Genes in Zea mays L. Identified a Role for ZmTCP42 in Drought Tolerance

The Teosinte-branched 1/Cycloidea/Proliferating (TCP) plant-specific transcription factors (TFs) have been demonstrated to play a fundamental role in plant development and organ patterning. However, it remains unknown whether or not the TCP gene family plays a role in conferring a tolerance to drought stress in maize, which is a major constraint to maize production. In this study, we identified 46 ZmTCP genes in the maize genome and systematically analyzed their phylogenetic relationships and synteny with rice, sorghum, and ArabidopsisTCP genes. Expression analysis of the 46 ZmTCP genes in different tissues and under drought conditions, suggests their involvement in maize response to drought stress. Importantly, genetic variations in ZmTCP32 and ZmTCP42 are significantly associated with drought tolerance at the seedling stage. RT-qPCR results suggest that ZmTCP32 and ZmTCP42 RNA levels are both induced by ABA, drought, and polyethylene glycol treatments. Based on the significant association between the genetic variation of ZmTCP42 and drought tolerance, and the inducible expression of ZmTCP42 by drought stress, we selected ZmTCP42, to investigate its function in drought response. We found that overexpression of ZmTCP42 in Arabidopsis led to a hypersensitivity to ABA in seed germination and enhanced drought tolerance, validating its function in drought tolerance. These results suggested that ZmTCP42 functions as an important TCP TF in maize, which plays a positive role in drought tolerance.

Molecular cloning, expression profiling, and functional analysis of a broad-complex isoform 2/3 (Br-Z2/Z3) transcription factor in the diamondback moth, Plutella xylostella (L.)

The diamondback moth, Plutella xylostella (L.) (Lepidoptera: Plutellidae), is a widespread and destructive pest of cruciferous crops. New strategies for controlling it are needed because it is rapidly developing resistance to conventional pesticides. In insects, transcription factors (TFs) including broad-complex (Br-C) are thought to be useful for insecticide development because they are able to regulate the transcription of functional genes involved in responses to external stimuli including insecticides. In the present study, we cloned and sequenced the open reading frames (ORFs) of three BTB-ZF encoding genes from the diamondback moth deposited in the National Center for Biotechnology Information (NCBI) database under accessions MG753773, MG288674, and MG753772. The lengths of these ORFs were 1,680, 1,428, and 1,647 bp, respectively. The phylogenetic analysis based on the predicted amino acid sequences of ZF domains showed that MG753773 and MG288674 belonged to Z2/Z3 and Z7 of Br-C while MG753772 belonged to Ttk types. In the agreement, the highest expression level of MG753773 occurred during the prepupal stage, MG288674 and MG753772 were expressed during all stages and peaked in the adult and egg stages, respectively. RNA interference silencing of MG753773 in the late third instar larvae significantly decreased survival and pupation of the insects. With precocene II, transcription of MG753773 increased (4×) in the fourth instar larva 24 hr later; 48 hr later the rate of prepupation and pupation was significantly higher. These findings will contribute to the development of new regulators of the growth and development for diamondback moth control.

Transcriptome analysis of harvested bell peppers (Capsicum annuum L.) in response to cold stress

Bell peppers are valued for their plentiful vitamin C and nutritional content. Pepper fruits are susceptible to cold storage, which leads to chilling injury (CI); however, the crucial metabolic product and molecular basis response to cold stress have not been elucidated definitely yet. To comprehensively understand the gene regulation network and CI mechanisms in response to cold stress on a molecular level, we performed high-throughput RNA-Seq analysis to investigate genome-wide expression profiles in bell peppers at different storage temperatures (4 °C and 10 °C). A total of 61.55 Gb of clean data were produced; 3863 differentially expressed genes (DEGs) including 1669 up-regulated and 2194 down-regulated were annotated and classified between the CI group and control. Together, a total of 41 cold-induced transcription factor families comprising 250 transcription factors (TFs) were identified. Notably, numerous DEGs involved in biomembrane stability, dehydration and osmoregulation, and plant hormone signal transduction processes were discovered. The transcriptional level of 20 DEGs was verified by reverse transcription quantitative polymerase chain reaction (RT-qPCR). Our results present transcriptome profiles of bell peppers in response to cold stress; the data obtained may be useful for the identification of key candidate genes and elucidation of the mechanisms underlying membrane damage during chilling injury.

DRL1, Encoding A NAC Transcription Factor, Is Involved in Leaf Senescence in Grapevine

The NAC (for NAM, ATAF1,2, and CUC2) proteins family are plant-specific transcription factors, which play important roles in leaf development and response to environmental stresses. In this study, an NAC gene, DRL1, isolated from grapevine Vitis vinifera L. "Yatomi Rose", was shown to be involved in leaf senescence. The quantity of DRL1 transcripts decreased with advancing leaf senescence in grapevine. Overexpressing the DRL1 gene in tobacco plants significantly delayed leaf senescence with respect to chlorophyll concentration, potential quantum efficiency of photosystem II (Fv/Fm), and ion leakage. Moreover, exogenous abscisic acid (ABA) markedly reduced the expression of DRL1, and the ABA and salicylic acid (SA) concentration was lower in the DRL1-overexpressing transgenic plants than in the wild-type plants. The DRL1 transgenic plants exhibited reduced sensitivity to ABA-induced senescence but no significant change in the sensitivity to jasmonic acid-, SA- or ethylene-induced senescence. Transcriptomic analysis and RNA expression studies also indicated that the transcript abundance of genes associated with ABA biosynthesis and regulation, including 9-cis-epoxycarotenoid dioxygenase (NCED1), NCED5, zeaxanthin epoxidase1 (ZEP1), ABA DEFICIENT2 (ABA2), ABA4, and ABA INSENSITIVE 2 (ABI2), was markedly reduced in the DRL1-overexpressing plants. These results suggested that DRL1 plays a role as a negative regulator of leaf senescence by regulating ABA synthesis.

Characterization of Dof family in Pyrus bretschneideri and role of PbDof9.2 in flowering time regulation

DNA binding with One Finger (Dof) proteins are plant-specific transcription factors with highly conserved Dof domain, including C2-C2 type zinc finger motifs. In this study, we identified 45 PbDofs in pear (Pyrusbretschneideri). PbDofs were classified into eight subfamilies by phylogenetic analysis. Conserved motifs of PbDof proteins were analyzed by MEME. PbDofs in subfamily D₁ werehomologous to CDFs in Arabidopsis. In this study, we showed that PbDof9.2 was regulated by both the circadian clock and photoperiod. PbDof9.2-GFP proteinwas localized in the nucleus. Overexpression of PbDof9.2 in Arabidopsis caused delayed flowering time. PbDof9.2 suppressed the flowering time regulator FT and could repress flowering time by promoting activity of PbTFL1a and PbTFL1b promoter. These results suggest that Doftranscription factors have conserved functions in plant development.

Revisiting the Role of Transcription Factors in Coordinating the Defense Response Against Citrus Bark Cracking Viroid Infection in Commercial Hop (Humulus Lupulus L.)

Transcription factors (TFs) play a major role in controlling gene expression by intricately regulating diverse biological processes such as growth and development, the response to external stimuli and the activation of defense responses. The systematic identification and classification of TF genes are essential to gain insight into their evolutionary history, biological roles, and regulatory networks. In this study, we performed a global mining and characterization of hop TFs and their involvement in Citrus bark cracking viroid CBCVd infection by employing a digital gene expression analysis. Our systematic analysis resulted in the identification of a total of 3,818 putative hop TFs that were classified into 99 families based on their conserved domains. A phylogenetic analysis classified the hop TFs into several subgroups based on a phylogenetic comparison with reference TF proteins from Arabidopsis thaliana providing glimpses of their evolutionary history. Members of the same subfamily and subgroup shared conserved motif compositions. The putative functions of the CBCVd-responsive hop TFs were predicted using their orthologous counterparts in A. thaliana. The analysis of the expression profiling of the CBCVd-responsive hop TFs revealed a massive differential modulation, and the expression of the selected TFs was validated using qRT-PCR. Together, the comprehensive integrated analysis in this study provides better insights into the TF regulatory networks associated with CBCVd infections in the hop, and also offers candidate TF genes for improving the resistance in hop against viroids.

Dynamics of Gene Expression in Single Root Cells of Arabidopsis thaliana

Single cell RNA sequencing can yield high-resolution cell-type-specific expression signatures that reveal new cell types and the developmental trajectories of cell lineages. Here, we apply this approach to Arabidopsis (Arabidopsis thaliana) root cells to capture gene expression in 3,121 root cells. We analyze these data with Monocle 3, which orders single cell transcriptomes in an unsupervised manner and uses machine learning to reconstruct single cell developmental trajectories along pseudotime. We identify hundreds of genes with cell-type-specific expression, with pseudotime analysis of several cell lineages revealing both known and novel genes that are expressed along a developmental trajectory. We identify transcription factor motifs that are enriched in early and late cells, together with the corresponding candidate transcription factors that likely drive the observed expression patterns. We assess and interpret changes in total RNA expression along developmental trajectories and show that trajectory branch points mark developmental decisions. Finally, by applying heat stress to whole seedlings, we address the longstanding question of possible heterogeneity among cell types in the response to an abiotic stress. Although the response of canonical heat-shock genes dominates expression across cell types, subtle but significant differences in other genes can be detected among cell types. Taken together, our results demonstrate that single cell transcriptomics holds promise for studying plant development and plant physiology with unprecedented resolution.

Genome-wide analyses and expression patterns under abiotic stress of NAC transcription factors in white pear (Pyrus bretschneideri)

BACKGROUND

Although the genome of Chinese white pear ('Dangshansuli') has been released, little is known about the functions, evolutionary history and expression patterns of NAC families in this species to date.

RESULTS

In this study, we identified a total of 183 NAC transcription factors (TFs) in the pear genome, among which 146 pear NAC (PbNAC) members were mapped onto 16 chromosomes, and 37 PbNAC genes were located on scaffold contigs. No PbNAC genes were mapped to chromosome 2. Based on gene structure, protein motif analysis, and topology of the phylogenetic tree, the pear PbNAC family was classified into 33 groups. By comparing and analyzing the unique NAC subgroups in Rosaceae, we identified 19 NAC subgroups specific to pear. We also found that whole-genome duplication (WGD)/segmental duplication played critical roles in the expansion of the NAC family in pear, such as the 83 PbNAC duplicated gene pairs dated back to the two WGD events. Further, we found that purifying selection was the primary force driving the evolution of PbNAC family genes. Next, we used transcriptomic data to study responses to drought and cold stresses in pear, and we found that genes in groups C2f, C72b, and C100a were related to drought and cold stress response.

CONCLUSIONS

Through the phylogenetic, evolutionary, and expression analyses of the NAC gene family in Chinese white pear, we indentified 11 PbNAC TFs associated with abiotic stress in pear.

MYB-CC transcription factor, TaMYBsm3, cloned from wheat is involved in drought tolerance

BACKGROUND

MYB-CC transcription factors (TFs) genes have been demonstrated to be involved in the response to inorganic phosphate (Pi) starvation and regulate some Pi-starvation-inducible genes. However, their role in drought stress has not been investigated in bread wheat. In this study, the TaMYBsm3 genes, including TaMYBsm3-A, TaMYBsm3-B, and TaMYBsm3-D, encoding MYB-CC TF proteins in bread wheat, were isolated to investigate the possible molecular mechanisms related to drought-tolerance in plants.

RESULTS

TaMYBsm3-A, TaMYBsm3-B, and TaMYBsm3-D were mapped on chromosomes 6A, 6B, and 6D in wheat, respectively. TaMYBsm3 genes belonged to MYB-CC TFs, containing a conserved MYB DNA-binding domain and a conserved coiled-coil domain. TaMYBsm3-D was localized in the nucleus, and the N-terminal region was a transcriptional activation domain. TaMYBsm3 genes were ubiquitously expressed in different tissues of wheat, and especially highly expressed in the stamen and pistil. Under drought stress, transgenic plants exhibited milder wilting symptoms, higher germination rates, higher proline content, and lower MDA content comparing with the wild type plants. P5CS1, DREB2A, and RD29A had significantly higher expression in transgenic plants than in wild type plants.

CONCLUSION

TaMYBsm3-A, TaMYBsm3-B, and TaMYBsm3-D were associated with enhanced drought tolerance in bread wheat. Overexpression of TaMYBsm3-D increases the drought tolerance of transgenic Arabidopsis through up-regulating P5CS1, DREB2A, and RD29A.

Genome-wide identification and expression profiling of trihelix gene family under abiotic stresses in wheat

BACKGROUND

The trihelix gene family is a plant-specific transcription factor family that plays important roles in plant growth, development, and responses to abiotic stresses. However, to date, no systemic characterization of the trihelix genes has yet been conducted in wheat and its close relatives.

RESULTS

We identified a total of 94 trihelix genes in wheat, as well as 22 trihelix genes in Triticum urartu, 29 in Aegilops tauschii, and 31 in Brachypodium distachyon. We analyzed the chromosomal locations and orthology relations of the identified trihelix genes, and no trihelix gene was found to be located on chromosome 7A, 7B, or 7D of wheat, thereby reflecting the uneven distributions of wheat trihelix genes. Phylogenetic analysis indicated that the 186 identified trihelix proteins in wheat, rice, B. distachyon, and Arabidopsis were clustered into five major clades. The trihelix genes belonging to the same clades usually shared similar motif compositions and exon/intron structural patterns. Five pairs of tandem duplication genes and three pairs of segmental duplication genes were identified in the wheat trihelix gene family, thereby validating the supposition that more intrachromosomal gene duplication events occur in the genome of wheat than in that of other grass species. The tissue-specific expression and differential expression profiling of the identified genes under cold and drought stresses were analyzed by using RNA-seq data. qRT-PCR was also used to confirm the expression profiles of ten selected wheat trihelix genes under multiple abiotic stresses, and we found that these genes mainly responded to salt and cold stresses.

CONCLUSIONS

In this study, we identified trihelix genes in wheat and its close relatives and found that gene duplication events are the main driving force for trihelix gene evolution in wheat. Our expression profiling analysis demonstrated that wheat trihelix genes responded to multiple abiotic stresses, especially salt and cold stresses. The results of our study built a basis for further investigation of the functions of wheat trihelix genes and provided candidate genes for stress-resistant wheat breeding programs.

Cell-Free Protein Synthesis of Plant Transcription Factors

Biotic and abiotic stimuli induce profound transcriptional reprograming in plants through sophisticated regulation of transcription factors (TFs). Recombinant proteins of TFs play an important role in unveiling their molecular functions. Cell-free protein synthesis (CFPS) system from wheat germ has been developed as one of the most efficient protein synthesis platforms. However, preparation of linear DNA templates for in vitro transcription is time-consuming and laborious. Here, we describe a versatile method for in vitro transcription and translation of the wheat germ CFPS system. Our two-step PCR method enables researchers to generate a variety of transcription templates from a single plasmid including fusion proteins of an N- or C-terminal tag and truncated proteins. Thus, this method supports a rapid and high-throughput expression of proteins for a large-scale proteomics analysis.

The Arabidopsis thaliana transcription factor MYB59 regulates calcium signalling during plant growth and stress response

Transcription factor MYB59 is involved in plant growth and stress responses by acting as negative regulator of Ca signalling and homeostasis. The Arabidopsis thaliana transcription factor MYB59 is induced by cadmium (Cd) and plays a key role in the regulation of cell cycle progression and root elongation, but its mechanism of action is poorly understood. We investigated the expression of MYB59 and differences between wild-type plants, the myb59 mutant and MYB59-overexpressing lines (obtained by transformation in the mutant genotype) during plant growth and in response to various forms of stress. We also compared the transcriptomes of wild-type and myb59 mutant plants to determine putative MYB59 targets. The myb59 mutant has longer roots, smaller leaves and smaller cells than wild-type plants and responds differently to stress in germination assay. Transcriptomic analysis revealed the upregulation in the myb59 mutant of multiple genes involved in calcium (Ca) homeostasis and signalling, including those encoding calmodulin-like proteins and Ca transporters. Notably, MYB59 was strongly induced by Ca deficiency, and the myb59 mutant was characterized by higher levels of cytosolic Ca in root cells and showed a modest alteration of Ca transient frequency in guard cells, associated with the absence of Ca-induced stomatal closure. These results indicate that MYB59 negatively regulates Ca homeostasis and signalling during Ca deficiency, thus controlling plant growth and stress responses.

Genome-wide analysis reveals the evolution and structural features of WRINKLED1 in plants

WRINKLED1 (WRI1), an AP2/ERE transcription factor, is one of the most important regulators of oil accumulation. It has been extensively studied in angiosperms, but its evolution and overview features in plants remain unknown. In this study, WRI1s, as well as WRI1-likes in non-WRI1 species, were investigated in 64 genome-sequenced plants. Their origin, distribution, duplication, evolution, functional domains, motifs, properties, and cis-elements were analyzed. Results suggest that WRI1 and WRI1-like may originate from Chlorophyta, and WRI1-likes in angiosperms resemble phylogenetically and structurally WRI1s from Chlorophyta and non-vascular plants. WRI1 or WRI1-like may be essential to vascular plants but not to non-vascular plants. Two YRG elements and two RAYD elements, as well as their phosphorylation sites and the 14-3-3 binding motif, are relatively conserved from Chlorophyta to angiosperm. The predicted DNA-binding domains are slightly shorter than the combination of one YRG element and one RAYD element. WRI1 gradually evolves from alkalinity to acidity. More motifs were developed in N-terminuses and C-terminuses in vascular plants. A short acidic amino-acid-enriched domain in the C-terminal region is predicted to be the putative transactivation domain. The VYL exon appears randomly in different WRI1 transcripts and it is not important for the function of WRI1. In addition, more cis-elements developed during WRI1 evolution may suggest its more complicated regulation and physiological functions. These results will assist future function studies of WRI1 and evolution studies of fatty acid biosynthesis regulation in plants.

Analysis of Rice Proteins with DLN Repressor Motif/S

Transcriptional regulation includes both activation and repression of downstream genes. In plants, a well-established class of repressors are proteins with an ERF-associated amphiphilic repression/EAR domain. They contain either DLNxxP or LxLxL as the identifying hexapeptide motif. In rice (Oryza sativa), we have identified a total of 266 DLN repressor proteins, with the former motif and its modifications thereof comprising 227 transcription factors and 39 transcriptional regulators. Apart from DLNxxP motif conservation, DLNxP and DLNxxxP motifs with variable numbers/positions of proline and those without any proline conservation have been identified. Most of the DLN repressome proteins have a single DLN motif, with higher relative percentage in the C-terminal region. We have designed a simple yeast-based experiment wherein a DLN motif can successfully cause strong repression of downstream reporter genes, when fused to a transcriptional activator of rice or yeast. The DLN hexapeptide motif is essential for repression, and at least two “DLN” residues cause maximal repression. Comparatively, rice has more DLN repressor encoding genes than Arabidopsis, and DLNSPP motif from rice is 40% stronger than the known Arabidopsis SRDX motif. The study reports a straightforward assay to analyze repressor activity, along with the identification of a strong DLN repressor from rice.

Genome-wide Analysis of bZIP Transcription Factors in wheat and Functional Characterization of a TabZIP under Abiotic Stress

The basic leucine zipper (bZIP) represents one of the largest as well as most diverse transcription factor (TFs) families. They are known to play role in both stress as well as in various plant developmental processes. In the present study, a total of 191 bZIP transcription factors have been identified from Triticum aestivum. Expression analysis during various stress conditions, developmental stages, different varieties and gene ontology enrichment analysis suggest their possible roles in abiotic stress as well as in developmental responses. In the current analysis, one of the members named as TabZIP (Traes_7AL_25850F96F.1) was selected for detailed analysis to understand its role under different abiotic stress conditions. Gene expression studies revealed differential expression of TabZIP in various abiotic stress conditions like heat, salinity and dehydration suggesting the possible role of bZIP in various stress mitigation mechanism. Arabidopsis transgenics overexpressing TabZIP showed enhanced tolerance to salinity, drought, heat and oxidative stress. Thus TabZIP (Traes_7AL_25850F96F.1) can serve as a candidate gene for improving heat as well as other abiotic stress tolerance and can be helpful in enhancing the crop productivity under stress conditions.

Field transcriptome revealed a novel relationship between nitrate transport and flowering in Japanese beech

Recent advances in molecular and genetic studies about flowering time control have been increasingly available to elucidate the physiological mechanism underlying masting, the intermittent and synchronized production of a large amount of flowers and seeds in plant populations. To identify unexplored developmental and physiological processes associated with masting, genome-wide transcriptome analysis is a promising tool, but such analyses have yet to be performed. We established a field transcriptome using a typical masting species, Japanese beech (Fagus crenata Blume), over two years, and analyzed the data using a nonlinear time-series analysis called convergent cross mapping. Our field transcriptome was found to undergo numerous changes depending on the status of floral induction and season. An integrated approach of high-throughput transcriptomics and causal inference was successful at detecting novel causal regulatory relationships between nitrate transport and florigen synthesis/transport in a forest tree species. The synergistic activation of nitrate transport and floral transition could be adaptive to simultaneously satisfy floral transition at the appropriate timing and the nitrogen demand needed for flower formation.

Characterization of the Rosa roxburghii Tratt transcriptome and analysis of MYB genes

Rosa roxburghii Tratt (Rosaceae) has a fruit that is flavorful, economically valuable, and highly nutritious, providing numerous health benefits. Myeloblastosis (MYB) proteins play key roles in the development and fruit quality of R. roxburghii. However, there is little available genomic and transcriptomic information for R. roxburghii. Here, a normalized cDNA library was constructed from five tissues, including the stem, leaf, flower, young fruit, and mature fruit, using the Illumina HiSeq 3000 platform. De novo assembly was performed, and 470.66 million clean reads were obtained. In total, 63,727 unigenes, with an average GC content of 42.08%, were discovered, 60,406 of which were annotated. In addition, 9,354 unigenes were assigned to Gene Ontology categories, and 20,202 unigenes were assigned to 25 Eukaryotic Ortholog Groups. Additionally, 19,508 unigenes were classified into 140 pathways of the Kyoto Encyclopedia of Genes and Genomes database. Based on the transcriptome, 163 unigenes associated with MYBs were detected. Among these genes, 75 genes were significantly expressed in the various tissues, including 10 R1 MYB, 42 R2R3 MYB, one R1R2R3 MYB, three R4 MYB and 19 atypical MYB-like proteins. The expression levels of the 12 MYB genes randomly selected for quantitative real-time PCR analysis corroborated the RNA sequencing results. A total of 37,545 microsatellites were detected, with an average expressed sequence tag–simple sequence repeat frequency of 0.59 (37,545/63,727). This transcriptome data improves our understanding of the role of MYB in R. roxburghii and will be valuable for identifying genes of interest.

Overexpression of RcLEC1-B, a HAP3 transcription factor of PLB from Rosa canina, increases the level of endogenous gibberellin and alters the development of cuticle and floral organs in Arabidopsis

The HAP3 subfamily gene RcLEC1-B, was isolated from protocorm-like body (PLB) of Rosa canina, encodes 213 amino acid residues. It was shown that RcLEC1-B was specifically expressed in PLB of R. canina and its subcellular localization is in the nucleus. Overexpression of RcLEC1-B in Arabidopsis resulted in a decrease in endogenous ABA level, an increase in GA, IAA and CTK contents, and an increased number of branches. RcLEC1-B promotes the formation of spontaneous embryoids, suggesting that it may be a homolog of the Arabidopsis LEC1 gene. RcLEC1-B-OE changed the number and morphology of flower organs and resulted in open carpels and exposed ovules, along with a reduced percentage of fertile fruit. This is the first observation that overexpression of a homolog of LEC1 in Arabidopsis can lead to morphological changes in floral organs, cuticle defects, and adhesions between organs; this may result from the increased level of gibberellin in the transgenic plants.

Genome-wide analysis, expansion and expression of the NAC family under drought and heat stresses in bread wheat (T. aestivum L.)

The NAC family is one of the largest plant-specific transcription factor families, and some of its members are known to play major roles in plant development and response to biotic and abiotic stresses. Here, we inventoried 488 NAC members in bread wheat (Triticum aestivum). Using the recent release of the wheat genome (IWGS RefSeq v1.0), we studied duplication events focusing on genomic regions from 4B-4D-5A chromosomes as an example of the family expansion and neofunctionalization of TaNAC members. Differentially expressed TaNAC genes in organs and in response to abiotic stresses were identified using publicly available RNAseq data. Expression profiling of 23 selected candidate TaNAC genes was studied in leaf and grain from two bread wheat genotypes at two developmental stages in field drought conditions and revealed insights into their specific and/or overlapping expression patterns. This study showed that, of the 23 TaNAC genes, seven have a leaf-specific expression and five have a grain-specific expression. In addition, the grain-specific genes profiles in response to drought depend on the genotype. These genes may be considered as potential candidates for further functional validation and could present an interest for crop improvement programs in response to climate change. Globally, the present study provides new insights into evolution, divergence and functional analysis of NAC gene family in bread wheat.

Insights into the complex role of GRAS transcription factors in the arbuscular mycorrhiza symbiosis

To improve access to limiting nutrients, the vast majority of land plants forms arbuscular mycorrhizal (AM) symbioses with Glomeromycota fungi. We show here that AM-related GRAS transcription factors from different subgroups are upregulated during a time course of mycorrhization. Based on expression studies in mutants defective in arbuscule branching (ram1-1, with a deleted MtRam1 GRAS transcription factor gene) or in the formation of functional arbuscules (pt4-2, mutated in the phosphate transporter gene MtPt4), we demonstrate that the five AM-related GRAS transcription factor genes MtGras1, MtGras4, MtGras6, MtGras7, and MtRad1 can be differentiated by their dependency on MtRAM1 and MtPT4, indicating that the network of AM-related GRAS transcription factors consists of at least two regulatory modules. One module involves the MtRAM1- and MtPT4-independent transcription factor MtGRAS4 that activates MtGras7. Another module is controlled by the MtRAM1- and MtPT4-dependent transcription factor MtGRAS1. Genome-wide expression profiles of mycorrhized MtGras1 knockdown and ram1-1 roots differ substantially, indicating different targets. Although an MtGras1 knockdown reduces transcription of AM-related GRAS transcription factor genes including MtRam1 and MtGras7, MtGras1 overexpression alone is not sufficient to activate MtGras genes. MtGras1 knockdown roots display normal fungal colonization, with a trend towards the formation of smaller arbuscules.

Systematic analysis of NAC transcription factors' gene family and identification of post-flowering drought stress responsive members in sorghum

SbNAC genes (131) encoding 183 proteins were identified from the sorghum genome and characterized. The expression patterns of SbSNACs were evaluated at three sampling time points under post-flowering drought stress. NAC proteins are specific transcription factors in plants, playing vital roles in development and response to various environmental stresses. Despite the fact that Sorghum bicolor is well-known for its drought-tolerance, it suffers from grain yield loss due to pre and post-flowering drought stress. In the present study, 131 SbNAC genes encoding 183 proteins were identified from the sorghum genome. The phylogenetic trees were constructed based on the NAC domains of sorghum, and also based on sorghum with Arabidopsis and 8 known NAC domains of other plants, which classified the family into 15 and 19 subfamilies, respectively. Based on the obtained results, 13 SbNAC proteins joined the SNAC subfamily, and these proteins are expected to be involved in response to abiotic stresses. Promoter analysis revealed that all SbNAC genes comprise different stress-associated cis-elements in their promoters. UTRs analysis indicated that 101 SbNAC transcripts had upstream open reading frames, while 39 of the transcripts had internal ribosome entry sites in their 5'UTR. Moreover, 298 miRNA target sites were predicted to exist in the UTRs of SbNAC transcripts. The expression patterns of SbSNACs were evaluated in three genotypes at three sampling time points under post-flowering drought stress. Based on the results, it could be suggested that some gene members are involved in response to drought stress at the post-flowering stage since they act as positive or negative transcriptional regulators. Following further functional analyses, some of these genes might be perceived to be promising candidates for breeding programs to enhance drought tolerance in crops.

Molecular cloning and expression analysis of GT-2-like genes in strawberry

GT-2 factors are the members of trihelix transcription factors (TFs) which can function in regulating plant development and responding to different abiotic stress. These proteins contain two structural domains composed by three tandem repeats helix-loop-helix-loop-helix. The strawberry (Fragaria × ananassa Duch.) is one of the most prevalent fruit crops due to its high economic and nutritional value. Nevertheless, strawberry production is limited by a range of biotic and abiotic stresses (such as drought, extreme temperature) that cause significant losses every year. Despite the potential roles of GT-2 transcription factor in plants, the functional and systematic analysis of the strawberry GT-2 subfamily has not been reported yet. In this research, we identified six GT-2 factors in 'Benihoppe' strawberry (Fragaria × ananassa) and all the FaGT-2-like proteins contain two trihelix domains. In addition, bioinformatics analysis showed that FaGT-2-like proteins might participate in transcription or transcription regulation. Compared with other reported GT-2 proteins, the similarity between FaGT-2-like and FvGT-2-like amino acid sequences was the highest, which can reach to 100%. Expression of these TFs indicated all of the FaGT-2-like genes could express in different tissues: root, stem, and leaf within distinct expression patterns. Furthermore, quantitative real-time PCR (qRT-PCR) analysis provided us with cues that all the FaGT-2-like genes were downregulated in response to various abiotic stress and hormone treatment. All the gene expressions can be inhibited by salt, drought, cold and ABA treatments, indicating that all the FaGT-2-like genes in 'Benihoppe' strawberry might act as the negative regulatory factors to respond to the abiotic stress. In summary, these results would lay a useful foundation for FaGT-2-like genes on functional study.

Insight into the B3Transcription Factor Superfamily and Expression Profiling of B3 Genes in Axillary Buds after Topping in Tobacco(Nicotiana tabacum L.)

Members of the plant-specific B3 transcription factor superfamily play important roles in various growth and developmental processes in plants. Even though there are many valuable studies on B3 genes in other species, little is known about the B3 superfamily in tobacco. We identified 114 B3 proteins from tobacco using comparative genome analysis. These proteins were classified into four subfamilies based on their phylogenetic relationships, and include the ARF, RAV, LAV, and REM subfamilies. The chromosomal locations, gene structures, conserved protein motifs, and sub-cellular localizations of the tobacco B3 proteins were analyzed. The patterns of exon-intron numbers and arrangement and the protein structures of the tobacco B3 proteins were in general agreement with their phylogenetic relationships. The expression patterns of 114 B3 genes revealed that many B3 genes show tissue-specific expression. The expression levels of B3 genes in axillary buds after topping showed that the REM genes are mainly up-regulated in response to topping, while the ARF genes are down-regulated after topping.

Transcriptome Profiling Reveals Effects of Drought Stress on Gene Expression in Diploid Potato Genotype P3-198

Potato (Solanum tuberosum L.) is one of the three most important food crops worldwide; however, it is strongly affected by drought stress. The precise molecular mechanisms of drought stress response in potato are not very well understood. The diploid potato genotype P3-198 has been verified to be highly resistant to drought stress. Here, a time-course experiment was performed to identify drought resistance response genes in P3-198 under polyethylene glycol (PEG)-induced stress using RNA-sequencing. A total of 1665 differentially expressed genes (DEGs) were specifically identified, and based on gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the transcription factor activity, protein kinase activity, and the plant hormone signal transduction process were significantly enriched. Annotation revealed that these DEGs mainly encode transcription factors, protein kinases, and proteins related to redox regulation, carbohydrate metabolism, and osmotic adjustment. In particular, genes encoding abscisic acid (ABA)-dependent signaling molecules were significantly differentially expressed, which revealed the important roles of the ABA-dependent signaling pathway in the early response of P3-198 to drought stress. Quantitative real-time PCR experimental verification confirmed the differential expression of genes in the drought resistance signaling pathway. Our results provide valuable information for understanding potato drought-resistance mechanisms, and also enrich the gene resources available for drought-resistant potato breeding.

Transcriptional regulators involved in responses to volatile organic compounds in plants

Field studies have shown that plants growing next to herbivore-infested plants acquire higher resistance to herbivore damage. This increased resistance is partly due to regulation of plant gene expression by volatile organic compounds (VOCs) released by plants that sense environmental challenges such as herbivores. The molecular basis for VOC sensing in plants, however, is poorly understood. Here, we report the identification of TOPLESS-like proteins (TPLs) that have VOC-binding activity and are involved in VOC sensing in tobacco. While screening for volatiles that induce stress-responsive gene expression in tobacco BY-2 cells and tobacco plants, we found that some sesquiterpenes induce the expression of stress-responsive genes. These results provided evidence that plants sense these VOCs and motivated us to analyze the mechanisms underlying volatile sensing using tobacco as a model system. Using a pulldown assay with caryophyllene derivative-linked beads, we identified TPLs as transcriptional co-repressors that bind volatile caryophyllene analogs. Overexpression of TPLs in cultured BY-2 cells or tobacco leaves reduced caryophyllene-induced gene expression, indicating that TPLs are involved in the responses to caryophyllene analogs in tobacco. We propose that unlike animals, which use membrane receptors for sensing odorants, a transcriptional co-repressor plays a role in sensing and mediating VOC signals in plant cells.

Signaling Overview of Plant Somatic Embryogenesis

Somatic embryogenesis (SE) is a means by which plants can regenerate bipolar structures from a somatic cell. During the process of cell differentiation, the explant responds to endogenous stimuli, which trigger the induction of a signaling response and, consequently, modify the gene program of the cell. SE is probably the most studied plant regeneration model, but to date it is the least understood due to the unclear mechanisms that occur at a cellular level. In this review, the authors seek to emphasize the importance of signaling on plant SE, highlighting the interactions between the different plant growth regulators (PGR), mainly auxins, cytokinins (CKs), ethylene and abscisic acid (ABA), during the induction of SE. The role of signaling is examined from the start of cell differentiation through the early steps on the embryogenic pathway, as well as its relation to a plant's tolerance of different types of stress. Furthermore, the role of genes encoded to transcription factors (TFs) during the embryogenic process such as the LEAFY COTYLEDON (LEC), WUSCHEL (WUS), BABY BOOM (BBM) and CLAVATA (CLV) genes, Arabinogalactan-proteins (AGPs), APETALA 2 (AP2) and epigenetic factors is discussed.

Genome-wide analysis of the NAC transcription factor family in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

The NAC (NAM, ATAF1/2, and CUC2) transcription factor family represents a group of large plant-specific transcriptional regulators, participating in plant development and response to external stress. However, there is no comprehensive study on the NAC genes of Tartary buckwheat (Fagopyrum tataricum), a large group of extensively cultivated medicinal and edible plants. The recently published Tartary buckwheat genome permits us to explore all the FtNAC genes on a genome-wide basis.

RESULTS

In the present study, 80 NAC (FtNAC) genes of Tartary buckwheat were obtained and named uniformly according to their distribution on chromosomes. Phylogenetic analysis of NAC proteins in both Tartary buckwheat and Arabidopsis showed that the FtNAC proteins are widely distributed in 15 subgroups with one subgroup unclassified. Gene structure analysis found that multitudinous FtNAC genes contained three exons, indicating that the structural diversity in Tartary buckwheat NAC genes is relatively low. Some duplication genes of FtNAC have a conserved structure that was different from others, indicating that these genes may have a variety of functions. By observing gene expression, we found that FtNAC genes showed abundant differences in expression levels in various tissues and at different stages of fruit development.

CONCLUSIONS

In this research, 80 NAC genes were identified in Tartary buckwheat, and their phylogenetic relationships, gene structures, duplication, global expression and potential roles in Tartary buckwheat development were studied. Comprehensive analysis will be useful for a follow-up study of functional characteristics of FtNAC genes and for the development of high-quality Tartary buckwheat varieties.

The rice ethylene response factor OsERF83 positively regulates disease resistance to Magnaporthe oryzae

Rice blast caused by Magnaporthe oryzae is one of the most destructive diseases of rice (Oryza sativa) worldwide. Here, we report the identification and functional characterization of a novel ethylene response factor (ERF) gene, OsERF83, which was expressed in rice leaves in response to rice blast fungus infection. OsERF83 expression was also induced by treatments with methyl jasmonate, ethephon, and salicylic acid, indicating that multiple phytohormones could be involved in the regulation of OsERF83 expression under biotic stress. Subcellular localization and transactivation analyses demonstrated that OsERF83 is a nucleus-localized transcriptional activator. A gel-shift assay using recombinant OsERF83 protein indicated that, like other ERFs, it binds to the GCC box. Transgenic rice plants overexpressing OsERF83 exhibited significantly suppressed lesion formation after rice blast infection, indicating that OsERF83 positively regulates disease resistance in rice. Genes encoding several classes of pathogenesis-related (PR) proteins, including PR1, PR2, PR3, PR5, and PR10, were upregulated in the OsERF83ox plants. Taken together, our findings show that OsERF83 is a novel ERF transcription factor that confers blast resistance by regulating the expression of defense-related genes in rice.

Impact of Alternatively Polyadenylated Isoforms of ETHYLENE RESPONSE FACTOR4 with Activator and Repressor Function on Senescence in Arabidopsis thaliana L

Leaf senescence is highly regulated by transcriptional reprogramming, implying an important role for transcriptional regulators. ETHYLENE RESPONSE FACTOR4 (ERF4) was shown to be involved in senescence regulation and to exist in two different isoforms due to alternative polyadenylation of its pre-mRNA. One of these isoforms, ERF4-R, contains an ERF-associated amphiphilic repression (EAR) motif and acts as repressor, whereas the other form, ERF4-A, is lacking this motif and acts as activator. Here, we analyzed the impact of these isoforms on senescence. Both isoforms were able to complement the delayed senescence phenotype of the erf4 mutant with a tendency of ERF4-A for a slightly better complementation. However, overexpression led to accelerated senescence of 35S:ERF4-R plants but not of 35S:ERF4-A plants. We identified CATALASE3 (CAT3) as direct target gene of ERF4 in a yeast-one-hybrid screen. Both isoforms directly bind to the CAT3 promoter but have antagonistic effects on gene expression. The ratio of ERF4-A to ERF4-R mRNA changed during development, leading to a complex age-dependent regulation of CAT3 activity. The RNA-binding protein FPA shifted the R/A-ratio and fpa mutants are pointing towards a role of alternative polyadenylation regulators in senescence.

Genome-wide characterization and expression analyses of the MYB superfamily genes during developmental stages in Chinese jujube

The MYB transcription factor (TF) superfamily, one of the largest gene superfamilies, regulates a variety of physiological processes in plants. Although many MYB superfamily genes have been identified in plants, the MYB TFs in Chinese jujube (Ziziphus jujuba Mill.) have not been fully identified and characterized. Additionally, the functions of these genes remain unclear. In total, we identified 171 MYB superfamily genes in jujube and divided them into five subfamilies containing 99 genes of the R2R3-MYB subfamily, 58 genes of the MYB-related subfamily, four genes of the R1R2R3-MYB subfamily, one gene of the 4R-MYB subfamily, and nine genes of the atypical MYB subfamily. The 99 R2R3-MYB genes of jujube were divided into 35 groups, C1–C35, and the 58 MYB-related genes were divided into the following groups: the R-R-type, CCA1-like, I-box-binding-like, TBP-like, CPC-like, and Chinese jujube-specific groups. ZjMYB genes in jujube were well supported by additional highly conserved motifs and exon/intron structures. Most R1 repeats of MYB-related proteins comprised the R2 repeat and had highly conserved EED and EEE residue groups in jujube. Three tandem duplicated gene pairs were found on 12 chromosomes in jujube. According to an expression analysis of 126 ZjMYB genes, MYB-related genes played important roles in jujube development and fruit-related biological processes. The total flavonoid content of jujube fruit decreased as ripening progressed. A total of 93 expressed genes were identified in the RNA-sequencing data from jujube fruit, and 56 ZjMYB members presented significant correlations with total flavonoid contents by correlation analysis. Five pairs of paralogous MYB genes within jujube were composed of nine jujube MYB genes. A total of 14 ZjMYB genes had the same homology to the MYB genes of Arabidopsis and peach, indicating that these 14 MYB genes and their orthologs probably existed before the ancestral divergence of the MYB superfamily. We used a synteny analysis of MYB genes between jujube and Arabidopsis to predict that the functions of the ZjMYBs involve flavonoid/phenylpropanoid metabolism, the light signaling pathway, auxin signal transduction, and responses to various abiotic stresses (cold, drought, and salt stresses). Additionally, we speculate that ZjMYB108 is an important TF involved in the flavonoid metabolic pathway. This comprehensive analysis of MYB superfamily genes in jujube lay a solid foundation for future comprehensive analyses of ZjMYB gene functions.

Genome-Wide Characterization and Identification of Trihelix Transcription Factor and Expression Profiling in Response to Abiotic Stresses in Rice (Oryza sativa L.)

Trihelix transcription factors play a role in plant growth, development and various stress responses. Here, we identified 41 trihelix family genes in the rice genome. These OsMSLs (Myb/SANT-LIKE) were located on twelve chromosomes. Synteny analysis indicated only six duplicated gene pairs in the rice trihelix family. Phylogenetic analysis of these OsMSLs and the trihelix genes from other species divided them into five clusters. OsMSLs from different groups significantly diverged in terms of gene structure and conserved functional domains. However, all OsMSLs contained the same five cis-elements. Some of these were responsive to light and dehydration stress. All OsMSLs expressed in four tissues and six developmental stages of rice but with different expression patterns. Quantitative real-time PCR analysis revealed that the OsMSLs responded to abiotic stresses including drought and high salt stress and stress signal molecule including ABA (abscisic acid), hydrogen peroxide. OsMSL39 were simultaneously expressed under all treatments, while OsMSL28 showed high expression under hydrogen peroxide, drought, and high salt treatments. Moreover, OsMSL16/27/33 displayed significant expression under ABA and drought treatments. Nevertheless, their responses were regulated by light. The expression levels of the 12 chosen OsMSLs differed between light and dark conditions. In conclusion, our results helped elucidate the biological functions of rice trihelix genes and provided a theoretical basis for further characterizing their biological roles in responding to abiotic stresses.

Genotype-specific suppression of multiple defense pathways in apple root during infection by Pythium ultimum

The genotype-specific defense activation in the roots of perennial tree crops to soilborne necrotrophic pathogens remains largely unknown. A recent phenotyping study indicated that the apple rootstock genotypes B.9 and G.935 have contrasting resistance responses to infection by Pythium ultimum. In the current study, a comparative transcriptome analysis by Illumina Solexa HiSeq 3000 platform was carried out to identify the global transcriptional regulation networks between the susceptible B.9 and the resistant G.935 to P. ultimum infection. Thirty-six libraries were sequenced to cover three timepoints after pathogen inoculation, with three biological replicates for each sample. The transcriptomes in the roots of the susceptible genotype B.9 were reflected by overrepresented differentially expressed genes (DEGs) with downregulated patterns and systematic suppression of cellular processes at 48 h post inoculation (hpi). In contrast, DEGs with annotated functions, such as kinase receptors, MAPK signaling, JA biosynthesis enzymes, transcription factors, and transporters, were readily induced at 24 hpi and continued up-regulation at 48 hpi in G.935 roots. The earlier and stronger defense activation is likely associated with an effective inhibition of necrosis progression in G.935 roots. Lack of effector-triggered immunity or existence of a susceptibility gene could contribute to the severely disturbed transcriptome and susceptibility in B.9 roots. The identified DEGs constitute a valuable resource for hypothesis-driven studies to elucidate the resistance/tolerance mechanisms in apple roots and validating their potential association with resistance traits.

A Novel Sweetpotato Transcription Factor Gene IbMYB116 Enhances Drought Tolerance in Transgenic Arabidopsis

Several members of the MYB transcription factor family have been found to regulate growth, developmental processes, metabolism, and biotic and abiotic stress responses in plants. However, the role of MYB116 in plants is still unclear. In this study, a MYB transcription factor gene IbMYB116 was cloned and characterized from the sweetpotato [Ipomoea batatas (L.) Lam.] line Xushu55-2, a line that is considered to be drought resistant. We show here that IbMYB116 is a nuclear protein and that it possesses a transactivation domain at the C terminus. This gene exhibited a high expression level in the leaf tissues of Xushu55-2 and was strongly induced by PEG6000 and methyl-jasmonate (MeJA). The IbMYB116-overexpressing Arabidopsis plants showed significantly enhanced drought tolerance, increased MeJA content, and a decreased H₂O₂ level under drought stress. The overexpression of IbMYB116 in Arabidopsis systematically upregulated jasmonic acid (JA) biosynthesis genes and activated the JA signaling pathway as well as reactive oxygen species (ROS)-scavenging system genes under drought stress conditions. The overall results suggest that the IbMYB116 gene might enhance drought tolerance by activating a ROS-scavenging system through the JA signaling pathway in transgenic Arabidopsis. These findings reveal, for the first time, the crucial role of IbMYB116 in the drought tolerance of plants.