Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes

Transcriptomic analysis of grapevine Dof transcription factor gene family in response to cold stress and functional analyses of the VaDof17d gene

Dof genes enhance cold tolerance in grapevine and VaDof17d is tightly associated with the cold-responsive pathway and with the raffinose family oligosaccharides. DNA-binding with one finger (Dof) proteins comprise a large family that plays important roles in the regulation of abiotic stresses. No in-depth analysis of Dof genes has been performed in the grapevine. In this study, we analyzed a total of 25 putative Dof genes in grapevine at genomic and transcriptomic levels, compiled expression profiles of 11 selected VaDof genes under cold stress and studied the potential function of the VaDof17d gene in grapevine calli. The 25 Dof proteins can be classified into four phylogenetic groups. RNA-seq and qRT-PCR results demonstrated that a total of 11 VaDof genes responded to cold stress. Comparative mRNA sequencing of 35S::VaDof17d grape calli showed that VaDof17d was tightly associated with the cold-responsive pathway and with the raffinose family oligosaccharides (RFOs), as observed by the up-regulation of galactinol synthase (GolS) and raffinose synthase genes. We found that the Dof17d-ED (CRISPR/Cas9-mediated mutagenesis of Dof17d-ED) mutant had low cold tolerance with a decreased RFOs level during cold stress. These results formed the fundamental knowledge for further analysis of the biological roles of Dof genes in the grapevine's adaption to cold stresses.

The Arabidopsis sucrose non-fermenting-1-related protein kinase AtSnRK2.4 interacts with a transcription factor, AtMYB21, that is involved in salt tolerance

Sucrose non-fermenting-1-related protein kinase 2 s (SnRK2 s) are important stress-related plant protein kinases in plants. The interaction partners and phosphorylation substrates of group II and III SnRK2 s in Arabidopsis thaliana have been identified, but similar data for group I SnRK2 s are very limited. Here, we used a yeast two-hybrid (Y2H) screen to find proteins that interact with Arabidopsis AtSnRK2.4, a group I SnRK2. The transcription factor AtMYB21 was identified as an AtSnRK2.4 interaction partner, and its interaction with AtSnRK2.4 was confirmed by an in vitro pull-down assay and a bimolecular fluorescence complementation (BiFC) assay. A subcellular localization assay demonstrated that AtSnRK2.4 and AtMYB21 were located in the cytoplasm and nucleus of onion epidermal cells. AtSnRK2.4 and AtMYB21 were expressed in many tissues and upregulated in response to NaCl stress. Transgenic plants that overexpressed AtSnRK2.4 or AtMYB21 gene exhibited enhanced tolerance to salt stress at germination and post-germination stages. Moreover, the expression of downstream stress-responsive genes was upregulated in salt-stressed AtSnRK2.4 and AtMYB21 transgenic Arabidopsis. These results suggest that AtSnRK2.4 may act synergistically with AtMYB21 to mediate the response to salt stress through the upregulation of downstream stress-responsive genes.

Genome-wide analysis of AP2/ERF transcription factors in pineapple reveals functional divergence during flowering induction mediated by ethylene and floral organ development

The APETALA2/ethylene-responsive factor (AP2/ERF) has important roles in regulating developmental processes and hormone signaling transduction in plants. Pineapple demonstrates a special sensitivity to ethylene, and AP2/ERFs may contribute to this distinct sensitivity of pineapples to ethylene. However, little information is available on the AP2/ERF of pineapple. In this study, 97 AP2/ERF family members were identified from the pineapple genome. The AcAP2/ERF superfamily could be further divided into five subfamilies, and different subfamily existed functional divergence in multifarious biological processes. ERF and RAV subfamily genes might play important roles in the process of ethylene response of pineapple; ERF and DREB subfamily genes had particular functions in the floral organ development. This study is the first to provide detailed information on the features of AP2/ERFs in pineapple, provide new insights into the potential functional roles of the AP2/ERF superfamily members, and will facilitate a better understanding of the molecular mechanism of flower in pineapple.

The Barley HvWRKY6 Transcription Factor Is Required for Resistance Against Pyrenophora teres f. teres

Barley is an important cereal crop worldwide because of its use in the brewing and distilling industry. However, adequate supplies of quality malting barley are threatened by global climate change due to drought in some regions and excess precipitation in others, which facilitates epidemics caused by fungal pathogens. The disease net form net blotch caused by the necrotrophic fungal pathogen Pyrenophora teres f. teres (Ptt) has emerged as a global threat to barley production and diverse populations of Ptt have shown a capacity to overcome deployed genetic resistances. The barley line CI5791 exhibits remarkably effective resistance to diverse Ptt isolates from around the world that maps to two major QTL on chromosomes 3H and 6H. To identify genes involved in this effective resistance, CI5791 seed were γ-irradiated and two mutants, designated CI5791-γ3 and CI5791-γ8, with compromised Ptt resistance were identified from an M₂ population. Phenotyping of CI5791-γ3 and -γ8 × Heartland F₂ populations showed three resistant to one susceptible segregation ratios and CI5791-γ3 × -γ8 F₁ individuals were susceptible, thus these independent mutants are in a single allelic gene. Thirty-four homozygous mutant (susceptible) CI5791-γ3 × Heartland F₂ individuals, representing 68 recombinant gametes, were genotyped via PCR genotype by sequencing. The data were used for single marker regression mapping placing the mutation on chromosome 3H within an approximate 75 cM interval encompassing the 3H CI5791 resistance QTL. Sequencing of the mutants and wild-type (WT) CI5791 genomic DNA following exome capture identified independent mutations of the HvWRKY6 transcription factor located on chromosome 3H at ∼50.7 cM, within the genetically delimited region. Post transcriptional gene silencing of HvWRKY6 in barley line CI5791 resulted in Ptt susceptibility, confirming that it functions in NFNB resistance, validating it as the gene underlying the mutant phenotypes. Allele analysis and transcript regulation of HvWRKY6 from resistant and susceptible lines revealed sequence identity and upregulation upon pathogen challenge in all genotypes analyzed, suggesting a conserved transcription factor is involved in the defense against the necrotrophic pathogen. We hypothesize that HvWRKY6 functions as a conserved signaling component of defense mechanisms that restricts Ptt growth in barley.

Inositol pyrophosphates promote the interaction of SPX domains with the coiled-coil motif of PHR transcription factors to regulate plant phosphate homeostasis

Phosphorus is an essential nutrient taken up by organisms in the form of inorganic phosphate (Pi). Eukaryotes have evolved sophisticated Pi sensing and signaling cascades, enabling them to stably maintain cellular Pi concentrations. Pi homeostasis is regulated by inositol pyrophosphate signaling molecules (PP-InsPs), which are sensed by SPX domain-containing proteins. In plants, PP-InsP-bound SPX receptors inactivate Myb coiled-coil (MYB-CC) Pi starvation response transcription factors (PHRs) by an unknown mechanism. Here we report that a InsP₈–SPX complex targets the plant-unique CC domain of PHRs. Crystal structures of the CC domain reveal an unusual four-stranded anti-parallel arrangement. Interface mutations in the CC domain yield monomeric PHR1, which is no longer able to bind DNA with high affinity. Mutation of conserved basic residues located at the surface of the CC domain disrupt interaction with the SPX receptor in vitro and in planta, resulting in constitutive Pi starvation responses. Together, our findings suggest that InsP₈ regulates plant Pi homeostasis by controlling the oligomeric state and hence the promoter binding capability of PHRs via their SPX receptors.

Plants regulate phosphate homeostasis via the interaction of PHR transcription factors with SPX receptors bound to inositol pyrophosphate signaling molecules. Here the authors show that inositol pyrophosphate-bound SPX interacts with the coiled-coil domain of PHR, which regulates the oligomerization and activity of the transcription factor.

Orchid Bsister gene PeMADS28 displays conserved function in ovule integument development

The ovules and egg cells are well developed to be fertilized at anthesis in many flowering plants. However, ovule development is triggered by pollination in most orchids. In this study, we characterized the function of a B_sister gene, named PeMADS28, isolated from Phalaenopsis equestris, the genome-sequenced orchid. Spatial and temporal expression analysis showed PeMADS28 predominantly expressed in ovules between 32 and 48 days after pollination, which synchronizes with integument development. Subcellular localization and protein–protein interaction analyses revealed that PeMADS28 could form a homodimer as well as heterodimers with D-class and E-class MADS-box proteins. In addition, ectopic expression of PeMADS28 in Arabidopsis thaliana induced small curled rosette leaves, short silique length and few seeds, similar to that with overexpression of other species’ B_sister genes in Arabidopsis. Furthermore, complementation test revealed that PeMADS28 could rescue the phenotype of the ABS/TT16 mutant. Together, these results indicate the conserved function of B_sisterPeMADS28 associated with ovule integument development in orchid.

Systematic investigation and expression profiles of the GbR2R3-MYB transcription factor family in ginkgo (Ginkgo biloba L.)

As one of the largest families of transcription factors, the R2R3-MYB family plays a significant role in plant growth, development, and response to hormone and environmental stress. To explore its evolutionary mechanism and potential function in Ginkgo biloba, a gymnosperm of great economic and ecological value, we presented a comprehensive analysis of the R2R3-MYB genes in ginkgo. Sixty-nine GbR2R3-MYB genes were identified and these genes could be classified into 33 groups based on the characteristics of the amino acid sequence of the R2R3-MYB domain and gene structure. Syntenic analyses indicated that few tandem and segmental duplications possibly resulted in the contraction of the GbR2R3-MYB gene family. Based on the transcriptome data, expression profiles of eight different tissues and different developmental stages of leaf and kernel showed that GbR2R3-MYB genes had distinct temporal and spatial expression characteristics. Specific expression patterns of the sixteen GbR2R3-MYB genes were also identified in response to different abiotic stresses and hormonal exposures. Further investigation revealed that GbR2R3-MYB19 was located in the nucleus and possessed transcriptional activity, implying its potential roles in the regulation of multiple biological processes. Our findings provide a robust basis for future comprehensive evolutionary and functional analyses of GbR2R3-MYB genes in ginkgo.

Transcriptome analysis of salt stress responsiveness in the seedlings of wild and cultivated Ricinus communis L

Soil salinity is one of the major environmental factors, influencing agricultural productivity of crops. As a non-edible and ideal oilseed crop, castor (Ricinus communis L.) has great industrial value in biofuel, but molecular mechanisms of salt stress regulation are still unknown. In this study, the differentially expressed genes (DEGs) for differential salt tolerance in two castor cultivar (wild castor : Y, cultivated castor 'Tongbi 5': Z) were identified. 12 libraries were sampled for Illumina high-throughput sequencing to consider 132,426 nonredundant unigenes and 31,221 gene loci. Multiple phytohormones and transcription factors (TFs) were correlated with salt-tolerance and differently enriched in these two genotypes. The type 2C protein phosphatases (PP2C) homologs were all upregulated under salt stress. Importantly, IAA (1), DELLA (1) and Jasmonate zim domain (JAZ) (1) were also identified and found to be differentially expressed. Based on the co-expressed module by regulatory networks and heatmap analysis, ERF/AP2, WRKY and bHLH families were prominently participate in high salt stress response of wild and cultivated castor. Finally, these results highlight that the hub DEGs and families were more accumulated in cultivated castor than those in wild castor, providing novel insights into the salinity adaptive mechanisms and genetic improvement in castor.

Genome-Wide Analysis Reveals the Potential Role of MYB Transcription Factors in Floral Scent Formation in Hedychium coronarium

The MYB gene family is one of the largest groups of transcription factors (TFs) playing diverse roles in several biological processes. Hedychium coronarium (white ginger lily) is a renowned ornamental plant both in tropical and subtropical regions due to its flower shape and strong floral scent mainly composed of terpenes and benzenoids. However, there is no information available regarding the role of the MYB gene family in H. coronarium. In the current study, the MYB gene family was identified and extensively analyzed. The identified 253 HcMYB genes were unevenly mapped on 17 chromosomes at a different density. Promoter sequence analysis showed numerous phytohormones related to cis-regulatory elements. The majority of HcMYB genes contain two to three introns and motif composition analysis showed their functional conservation. Phylogenetic analysis revealed that HcMYBs could be classified into 15 distinct clades, and the segmental duplication events played an essential role in the expansion of the HcMYB gene family. Tissue-specific expression patterns of HcMYB genes displayed spatial and temporal expression. Furthermore, seven HcMYB (HcMYB7/8/75/79/145/238/248) were selected for further investigation. Through RT-qPCR, the response of candidates HcMYB genes toward jasmonic acid methyl ester (MeJA), abscisic acid (ABA), ethylene, and auxin was examined. Yeast one-hybrid (Y1H) assays revealed that candidate genes directly bind to the promoter of bottom structural volatile synthesis genes (HcTPS1, HcTPS3, HcTPS10, and HcBSMT2). Moreover, yeast two-hybrid (Y2H) assay showed that HcMYB7/8/75/145/248 interact with HcJAZ1 protein. In HcMYB7/8/79/145/248-silenced flowers, the floral volatile contents were decreased and downregulated the expression of key structural genes, suggesting that these genes might play crucial roles in floral scent formation in H. coronarium by regulating the expression of floral scent biosynthesis genes. Collectively, these findings indicate that HcMYB genes might be involved in the regulatory mechanism of terpenoids and benzenoid biosynthesis in H. coronarium.

IbERF71, with IbMYB340 and IbbHLH2, coregulates anthocyanin accumulation by binding to the IbANS1 promoter in purple-fleshed sweet potato (Ipomoea batatas L.)

KEY MESSAGE

The transcription factor (TF) IbERF71 forms a novel complex, IbERF71-IbMYB340-IbbHLH2, to coregulate anthocyanin biosynthesis by binding to the IbANS1 promoter in purple-fleshed sweet potatoes. Purple-fleshed sweet potato (Ipomoea batatas L.) is very popular because of its abundant anthocyanins, which are natural pigments with multiple physiological functions. TFs involved in regulating anthocyanin biosynthesis have been identified in many plants. However, the molecular mechanism of anthocyanin biosynthesis in purple-fleshed sweet potatoes has rarely been examined. In this study, TF IbERF71 and its partners were screened by bioinformatics and RT-qPCR analysis. The results showed that the expression levels of IbERF71 and partners IbMYB340 and IbbHLH2 were higher in purple-fleshed sweet potatoes than in other colors and that the expression levels positively correlated with anthocyanin contents. Moreover, transient expression assays showed that cotransformation of IbMYB340+IbbHLH2 resulted in anthocyanin accumulation in tobacco leaves and strawberry receptacles, and additional IbERF71 significantly increased visual aspects. Furthermore, the combination of the three TFs significantly increased the expression levels of FvANS and FvGST, which are involved in anthocyanin biosynthesis and transport of strawberry receptacles. The dual-luciferase reporter system verified that cotransformation of the three TFs enhanced the transcription activity of IbANS1. In addition, yeast two-hybrid and firefly luciferase complementation assays revealed that IbMYB340 interacted with IbbHLH2 and IbERF71 but IbERF71 could not interact with IbbHLH2 in vitro. In summary, our findings provide novel evidence that IbERF71 and IbMYB340-IbbHLH2 form the regulatory complex IbERF71-IbMYB340-IbbHLH2 that coregulates anthocyanin accumulation by binding to the IbANS1 promoter in purple-fleshed sweet potatoes. Thus, the present study provides a new regulatory network of anthocyanin biosynthesis and strong insight into the color development of purple-fleshed sweet potatoes.

Genome-Wide Analysis of the Catharanthus roseus RLK1-Like in Soybean and GmCrRLK1L20 Responds to Drought and Salt Stresses

Abiotic stresses, such as drought and salinity, severely affects the growth, development and productivity of the plants. The Catharanthus roseus RLK1-like (CrRLK1L) protein kinase family is involved in several processes in the plant life cycle. However, there have been few studies addressing the functions of CrRLK1L proteins in soybean. In this study, 38 CrRLK1L genes were identified in the soybean genome (Glycine max Wm82.a2.v1). Phylogenetic analysis demonstrated that soybean CrRLK1L genes were grouped into clusters, cluster I, II, III. The chromosomal mapping demonstrated that 38 CrRLK1L genes were located in 14 of 20 soybean chromosomes. None were discovered on chromosomes 1, 4, 6, 7, 11, and 14. Gene structure analysis indicated that 73.6% soybean CrRLK1L genes were characterized by a lack of introns.15.7% soybean CrRLK1L genes only had one intron and 10.5% soybean CrRLK1L genes had more than one intron. Five genes were obtained from soybean drought- and salt-induced transcriptome databases and were found to be highly up-regulated. GmCrRLK1L20 was notably up-regulated under drought and salinity stresses, and was therefore studied further. Subcellular localization analysis revealed that the GmCrRLK1L20 protein was located in the cell membrane. The overexpression of the GmCrRLK1L20 gene in soybean hairy roots improved both drought tolerance and salt stresses and enhanced the expression of the stress-responsive genes GmMYB84, GmWRKY40, GmDREB-like, GmGST15, GmNAC29, and GmbZIP78. These results indicated that GmCrRLK1L20 could play a vital role in defending against drought and salinity stresses in soybean.

AtMPK6-induced phosphorylation of AtERF72 enhances its DNA binding activity and interaction with TGA4/OBF4 in Arabidopsis

The ethylene-responsive element binding factor (ERF) family is a large family of transcription factors involved in plant development and environmental stress responses. We previously reported the identification of 29 putative substrates of Mitogen-activated Protein Kinase3 (AtMPK3), AtMPK4 and AtMPK6, based on a solid-phase phosphorylation screening using a lambda phage expression library in Arabidopsis thaliana. In this study, a putative MPK substrate, AtERF72 (At3g16770), was strongly phosphorylated by AtMPK6 on the serine residue at position 151 (Ser151). AtERF72 binds to the GCC box (AGCCGCC) in the promoters of several pathogenesis-related (PR) genes and activates their transcription. We also show that the DNA-binding activity of AtERF72 is enhanced upon phosphorylation by AtMPK6 in vitro. In addition, transient co-expression experiments in Arabidopsis protoplasts revealed that effector constructs expressing a mutant variant of AtERF72, AtERF72^S151D (carrying a Ser to aspartic acid [Asp] substitution at amino acid position 151) showed higher expression of the β-glucuronidase (GUS) reporter gene driven by the GCC box element than effector constructs expressing the wild-type AtERF72. Furthermore, yeast two-hybrid assays revealed that the interaction between AtERF72^S151D and TGA4/OBF4 was stronger than that between wild-type AtERF72 and TGA4/OBF4. Since AtERF72^S151D is equivalent to AtERF72 phosphorylated by AtMPK6 at Ser151, these results suggest that the phosphorylation of AtERF72 by AtMPK6 triggers an event of transcriptional regulation from defence signalling in Arabidopsis.

Identification and Expression Analysis of the CsMYB Gene Family in Root Knot Nematode-Resistant and Susceptible Cucumbers

MYB (myeloblastosis) transcription factors (TFs) play important roles in controlling various physiological processes in plants, such as responses to biotic and abiotic stress, metabolism, and defense. A previous study identified a gene, Csa6G410090, encoding a plant lipid transfer protein (LTP), as a possible regulator in cucumber (Cucumis sativus L.) of the resistance response to root-knot nematode (RKN) [Meloidogyne incognita Kofoid and White (Chitwood)]. Myb-type DNA-binding TFs were presumed to regulate downstream genes expression, including LTPs, however, the regulation mechanism remained unclear. To elucidate whether and which MYB TFs may be involved in regulation of the resistance response, this study identified 112 genes as candidate members of the CsMYB gene family by combining CDD and SMART databases, using the Hidden Markov Model (HMM) and manual calibration. Within this group, ten phylogenetic subgroups were resolved according to sequence-based classification, consistent with results from comprehensive investigation of gene structure, conserved motifs, chromosome locations, and cis-element analysis. Distribution and collinearity analysis indicated that amplification of the CsMYB gene family in cucumber has occurred mainly through tandem repeat events. Spatial gene expression analysis showed that 8 CsMYB genes were highly expressed at differing levels in ten different tissues or organs. The roots of RKN-resistant and susceptible cucumbers were inoculated with M. incognita, finding that CsMYB (Csa6G538700, Csa1G021940, and Csa5G641610) genes showed up-regulation coincident with upregulation of the "hub" gene LTP (Csa6G410090) previously implicated as a major gene in the resistance response to RKN in cucumber. Results of this study suggest hypotheses regarding the elements and regulation of the resistant response as well as possible RKN resistance-enhancing strategies in cucumber and perhaps more broadly in plants.

Intrinsic Disorder in Plant Transcription Factor Systems: Functional Implications

Eukaryotic cells are complex biological systems that depend on highly connected molecular interaction networks with intrinsically disordered proteins as essential components. Through specific examples, we relate the conformational ensemble nature of intrinsic disorder (ID) in transcription factors to functions in plants. Transcription factors contain large regulatory ID-regions with numerous orphan sequence motifs, representing potential important interaction sites. ID-regions may affect DNA-binding through electrostatic interactions or allosterically as for the bZIP transcription factors, in which the DNA-binding domains also populate ensembles of dynamic transient structures. The flexibility of ID is well-suited for interaction networks requiring efficient molecular adjustments. For example, Radical Induced Cell Death1 depends on ID in transcription factors for its numerous, structurally heterogeneous interactions, and the JAZ:MYC:MED15 regulatory unit depends on protein dynamics, including binding-associated unfolding, for regulation of jasmonate-signaling. Flexibility makes ID-regions excellent targets of posttranslational modifications. For example, the extent of phosphorylation of the NAC transcription factor SOG1 regulates target gene expression and the DNA-damage response, and phosphorylation of the AP2/ERF transcription factor DREB2A acts as a switch enabling heat-regulated degradation. ID-related phase separation is emerging as being important to transcriptional regulation with condensates functioning in storage and inactivation of transcription factors. The applicative potential of ID-regions is apparent, as removal of an ID-region of the AP2/ERF transcription factor WRI1 affects its stability and consequently oil biosynthesis. The highlighted examples show that ID plays essential functional roles in plant biology and has a promising potential in engineering.

Genome-wide identification of NF-YA gene family in cotton and the positive role of GhNF-YA10 and GhNF-YA23 in salt tolerance

Nuclear factor YA (NF-YA) genes play important roles in many biological processes, such as leaf growth, nitrogen nutrition, drought resistance, and salt stress. The functions of NF-YA genes in cotton have not been elucidated. The current study identified a total of 16, 16, 31, and 29 genes from Gossypium raimondii, G. arboretum, G. barbadense, and G. hirsutum, respectively. The NF-YA genes in cotton were phylogenetically classified into 4 groups. Analysis of gene structure, conserved motifs and multiple sequence alignments supported the evolutionary conservation of NF-YA family genes in cotton. Analysis of the expression patterns of GhNF-YAs in cotton suggested that GhNF-YAs play important roles in plant growth, development, and stress responses. The quantitative real-time PCR (qRT-PCR) validation of selected genes suggested that GhNF-YA genes are induced in response to salt, drought, ABA, and MeJA treatments. GhNF-YA genes may regulate salt and drought stress via the ABA or MeJA pathway. Silencing of GhNF-YA10 and GhNF-YA23 significantly reduced the salt tolerance of cotton seedlings, indicating that these genes participate in the regulation of the response of cotton to salt stress. These results establish a foundation for subsequent functional studies of the NF-YA gene family in cotton.

Co-expression network analysis of protein phosphatase 2A (PP2A) genes with stress-responsive genes in Arabidopsis thaliana reveals 13 key regulators

Abiotic and biotic stresses adversely affect plant growth and development and eventually result in less yield and threaten food security worldwide. In plants, several studies have been carried out to understand molecular responses to abiotic and biotic stresses. However, the complete circuitry of stress-responsive genes that plants utilise in response to those environmental stresses are still unknown. The protein phosphatase 2A (PP2A) gene has been known to have a crucial role in abiotic and biotic stresses; but how it regulates the stress response in plants is still not known completely. In this study, we constructed gene co-expression networks of PP2A genes with stress-responsive gene datasets from cold, drought, heat, osmotic, genotoxic, salt, and wounding stresses to unveil their relationships with the PP2A under different conditions of stress. The graph analysis identified 13 hub genes and several influential genes based on closeness centrality score (CCS). Our findings also revealed the count of unique genes present in different settings of stresses and subunits. We also formed clusters of influential genes based on the stress, CCS, and co-expression value. Analysis of cis-regulatory elements (CREs), recurring in promoters of these genes was also performed. Our study has led to the identification of 16 conserved CREs.

Gene co-expression analysis reveals transcriptome divergence between wild and cultivated chickpea under drought stress

Ancestral adaptations in crop wild relatives can provide a genetic reservoir for crop improvement. Here we document physiological changes to mild and severe drought stress, and the associated transcriptome dynamics in both wild and cultivated chickpea. Over 60% of transcriptional changes were related to metabolism, indicating that metabolic plasticity is a core and conserved drought response. In addition, changes in RNA processing and protein turnover were predominant in the data, suggestive of broad restructuring of the chickpea proteome in response to drought. While 12% of the drought-responsive transcripts have similar dynamics in cultivated and wild accessions, numerous transcripts had expression patterns unique to particular genotypes, or that distinguished wild from cultivated genotypes and whose divergence may be a consequence of domestication. These and other comparisons provide a transcriptional correlate of previously described species' genetic diversity, with wild accessions well differentiated from each other and from cultivars, and cultivars essentially indistinguishable at the broad transcriptome level. We identified metabolic pathways such as phenylpropanoid metabolism, and biological processes such as stomatal development, which are differentially regulated across genotypes with potential consequences on drought tolerance. These data indicate that wild Cicer reticulatum may provide both conserved and divergent mechanisms as a resource in breeding for drought tolerance in cultivated chickpea.

Overexpression of SlMBP22 in Tomato Affects Plant Growth and Enhances Tolerance to Drought Stress

MADS-box transcription factors play crucial and diverse roles in plant growth and development, and the responses to biotic and abiotic stresses. However, the implementation of MADS-box transcription factors in regulating plant architecture and stress responses has not been fully explored in tomato. Here, we found that a novel MADS-box transcription factor, SlMBP22, participated in the control of agronomical traits, tolerance to abiotic stress, and regulation of auxin and gibberellin signalling. Transgenic plants overexpressing SlMBP22 (SlMBP22-OE) displayed pleiotropic phenotypes, including reduced plant height and leaf size, by affecting auxin and/or gibberellin signalling. SlMBP22 was induced by dehydration treatment, and SlMBP22-OE plants were more tolerant to drought stress than wild-type (WT). Furthermore, SlMBP22 overexpression plants accumulated more chlorophyll, starch and soluble sugar than WT, indicating that the darker green leaves might be attributed to increased chlorophyll levels in the transgenic plants. RNA-Seq results showed that the transcript levels of a series of genes related to chloroplast development, chlorophyll metabolism, starch and sucrose metabolism, hormone signalling, and stress responses were altered. Collectively, our data demonstrate that SlMBP22 plays an important role in both regulating tomato growth and resisting drought stress.

AtHAP5A regulates iron translocation in iron-deficient Arabidopsis thaliana

Iron (Fe) deficient plants employ multiple strategies to increase root uptake and root-to-shoot translocation of Fe. The identification of genes that are responsible for these processes, and a comprehensive understanding of the regulatory effects of transcriptional networks on their expression, including transcription factors (TFs), is underway in Arabidopsis thaliana. Here, we show that a Histone- or heme-associated proteins (HAP) transcription factor (TF), HAP5A, is necessary for the response to Fe deficiency in Arabidopsis. Its expression was induced under Fe deficiency, and the lack of HAP5A significantly decreased Fe translocation from the root to the shoot, resulting in substantial chlorosis of the newly expanded leaves, compared with the wild-type (WT, Col-0). Further analysis found that the expression of a gene encoding nicotianamine (NA) synthase (NAS1) was dramatically decreased in the hap5a mutant, regardless of the Fe status. Yeast-one-hybrid and ChIP analyses suggested that HAP5A directly binds to the promoter region of NAS1. Moreover, overexpression of NAS1 could rescue the chlorosis phenotype of hap5a in Fe deficient conditions. In summary, a novel pathway was elucidated, showing that NAS1-dependent translocation of Fe from the root to the shoot is controlled by HAP5A in Fe-deficient Arabidopsis thaliana.

Identification and Characterization of PLATZ Transcription Factors in Wheat

The PLATZ (plant AT-rich protein and zinc-binding protein) transcription factor family is a class of plant-specific zinc-dependent DNA-binding proteins. PLATZ has essential roles in seed endosperm development, as well as promoting cell proliferation duration in the earlier stages of the crops. In the present study, 62 TaPLATZ genes were identified from the wheat genome, and they were unequally distributed on 15 chromosomes. According to the phylogenetic analysis, 62 TaPLATZ genes were classified into six groups, including two groups that were unique in wheat. Members in the same groups shared similar exon-intron structures. The polyploidization, together with genome duplication of wheat, plays a crucial role in the expansion of the TaPLATZs family. Transcriptome data indicated a distinct divergence expression pattern of TaPLATZ genes that could be clustered into four modules. The TaPLATZs in Module b possessed a seed-specific expression pattern and displayed obvious high expression in the earlier development stage of seeds. Subcellular localization data of TaPLATZs suggesting that they likely perform a function as a conventional transcription factor. This study provides insight into understanding the structure divergence, evolutionary features, expression profiles, and potential function of PLATZ in wheat.

Genome-Wide Identification and Expression Profiling Analysis of the Trihelix Gene Family Under Abiotic Stresses in Medicago truncatula

The trihelix transcription factor (GT) family is widely involved in regulating plant growth and development, and most importantly, responding to various abiotic stresses. Our study first reported the genome-wide identification and analysis of GT family genes in Medicago truncatula. Overall, 38 trihelix genes were identified in the M. truncatula genome and were classified into five subfamilies (GT-1, GT-2, SH4, GTγ and SIP1). We systematically analyzed the phylogenetic relationship, chromosomal distribution, tandem and segmental duplication events, gene structures and conserved motifs of MtGTs. Syntenic analysis revealed that trihelix family genes in M. truncatula had the most collinearity relationship with those in soybean followed by alfalfa, but very little collinearity with those in the maize and rice. Additionally, tissue-specific expression analysis of trihelix family genes suggested that they played various roles in the growth and development of specific tissues in M. truncatula. Moreover, the expression of some MtGT genes, such as MtGT19, MtGT20, MtGT22, and MtGT33, was dramatically induced by drought, salt, and ABA treatments, illustrating their vital roles in response to abiotic stresses. These findings are helpful for improving the comprehensive understanding of trihelix family; additionally, the study provides candidate genes for achieving the genetic improvement of stress resistance in legumes.

The Aureochrome Photoreceptor PtAUREO1a Is a Highly Effective Blue Light Switch in Diatoms

Aureochromes represent a unique type of blue light photoreceptors that possess a blue light sensing flavin-binding LOV-domain and a DNA-binding bZIP domain, thus being light-driven transcription factors. The diatom Phaeodactylum tricornutum, a member of the essential marine primary producers, possesses four aureochromes (PtAUREO1a, 1b, 1c, 2). Here we show a dramatic change in the global gene expression pattern of P. tricornutum wild-type cells after a shift from red to blue light. About 75% of the genes show significantly changed transcript levels already after 10 and 60 min of blue light exposure, which includes genes of major transcription factors as well as other photoreceptors. Very surprisingly, this light-induced regulation of gene expression is almost completely inhibited in independent PtAureo1a knockout lines. Such a massive and fast transcriptional change depending on one single photoreceptor is so far unprecedented. We conclude that PtAUREO1a plays a key role in diatoms upon blue light exposure.

•

Blue light induces a very fast transcriptional response in the diatom P. tricornutum

•

This strong response is almost completely inhibited when Aureochrome 1a is absent

•

The results imply a key role of PtAureo1a in blue light-induced responses in diatoms

In Silico Identification of MYB and bHLH Families Reveals Candidate Transcription Factors for Secondary Metabolic Pathways in Cannabis sativa L

Plant secondary metabolic pathways are finely regulated by the activity of transcription factors, among which members of the bHLH and MYB subfamilies play a main role. Cannabis sativa L. is a unique officinal plant species with over 600 synthesized phytochemicals having diverse scale-up industrial and pharmaceutical usage. Despite comprehensive knowledge of cannabinoids' metabolic pathways, very little is known about their regulation, while the literature on flavonoids' metabolic pathways is still scarce. In this study, we provide the first genome-wide analysis of bHLH and MYB families in C. sativa reference cultivar CBDRx and identification of candidate coding sequences for these transcription factors. Cannabis sativa bHLHs and MYBs were then classified into functional subfamilies through comparative phylogenetic analysis with A. thaliana transcription factors. Analyses of gene structure and motif distribution confirmed that CsbHLHs and CsMYBs belonging to the same evolutionary clade share common features at both gene and amino acidic level. Candidate regulatory genes for key metabolic pathways leading to flavonoid and cannabinoid synthesis in Cannabis were also retrieved. Furthermore, a candidate gene approach was used to identify structural enzyme-coding genes for flavonoid and cannabinoid synthesis. Taken as a whole, this work represents a valuable resource of candidate genes for further investigation of the C. sativa cannabinoid and flavonoid metabolic pathways for genomic studies and breeding programs.

Structure and expression analysis of seven salt-related ERF genes of Populus

Ethylene response factors (ERFs) are plant-specific transcription factors (TFs) that play important roles in plant growth and stress defense and have received a great amount of attention in recent years. In this study, seven ERF genes related to abiotic stress tolerance and response were identified in plants of the Populus genus. Systematic bioinformatics, including sequence phylogeny, genome organisation, gene structure, gene ontology (GO) annotation, etc. were detected. Expression-pattern of these seven ERF genes were analyzed using RT-qPCR and cross validated using RNA-Seq. Data from a phylogenetic tree and multiple alignment of protein sequences indicated that these seven ERF TFs belong to three subfamilies and contain AP2, YRG, and RAYD conserved domains, which may interact with downstream target genes to regulate the plant stress response. An analysis of the structure and promoter region of these seven ERF genes showed that they have multiple stress-related motifs and cis-elements, which may play roles in the plant stress-tolerance process through a transcriptional regulation mechanism; moreover, the cellular_component and molecular_function terms associated with these ERFs determined by GO annotation supported this hypothesis. In addition, the spatio-temporal expression pattern of these seven ERFs, as detected using RT-qPCR and RNA-seq, suggested that they play a critical role in mediating the salt response and tolerance in a dynamic and tissue-specific manner. The results of this study provide a solid basis to explore the functions of the stress-related ERF TFs in Populus abiotic stress tolerance and development process.

Light and ripening-regulated BBX protein-encoding genes in Solanum lycopersicum

Light controls several aspects of plant development through a complex signalling cascade. Several B-box domain containing proteins (BBX) were identified as regulators of Arabidopsis thaliana seedling photomorphogenesis. However, the knowledge about the role of this protein family in other physiological processes and species remains scarce. To fill this gap, here BBX protein encoding genes in tomato genome were characterised. The robust phylogeny obtained revealed how the domain diversity in this protein family evolved in Viridiplantae and allowed the precise identification of 31 tomato SlBBX proteins. The mRNA profiling in different organs revealed that SlBBX genes are regulated by light and their transcripts accumulation is directly affected by the chloroplast maturation status in both vegetative and fruit tissues. As tomato fruits develops, three SlBBXs were found to be upregulated in the early stages, controlled by the proper chloroplast differentiation and by the PHYTOCHROME (PHY)-dependent light perception. Upon ripening, other three SlBBXs were transcriptionally induced by RIPENING INHIBITOR master transcriptional factor, as well as by PHY-mediated signalling and proper plastid biogenesis. Altogether, the results obtained revealed a conserved role of SlBBX gene family in the light signalling cascade and identified putative members affecting tomato fruit development and ripening.

Transcriptome sequencing and metabolite profiling analyses provide comprehensive insight into molecular mechanisms of flower development in Dendrobium officinale (Orchidaceae)

This research provides comprehensive insight into the molecular networks and molecular mechanisms underlying D. officinale flower development. Flowers are complex reproductive organs and play a crucial role in plant propagation, while also providing sustenance for insects and natural bioactive metabolites for humans. However, knowledge about gene regulation and floral metabolomes in flowers is limited. In this study, we used an important orchid species (Dendrobium officinale), whose flowers can be used to make herbal tea, to perform transcriptome sequencing and metabolic profiling of early- and medium-stage flower buds, as well as opened flowers, to provide comprehensive insight into the molecular mechanisms underlying flower development. A total of 8019 differentially expressed genes (DEGs) and 239 differentiated metabolites were found. The transcription factors that were identified and analyzed belong exclusively to the MIKC-type MADS-box proteins and auxin responsive factors that are known to be involved in flower development. The expression of genes involved in chlorophyll and carotenoid biosynthesis strongly matched the metabolite accumulation patterns. The genes related to flavonoid and polysaccharide biosynthesis were active during flower development. Interestingly, indole-3-acetic acid and abscisic acid, whose trend of accumulation was inverse during flower development, may play an important role in this process. Collectively, the identification of DEGs and differentiated metabolites could help to illustrate the regulatory networks and molecular mechanisms important for flower development in this orchid.

Genome-wide analysis of PmTCP4 transcription factor binding sites by ChIP-Seq during pistil abortion in Japanese apricot

The TCP4 transcription factor plays an important role in plant growth and development, especially in flower development. PmTCP4 is involved in the process of pistil abortion in Japanese apricot, but its molecular mechanism, particularly the DNA binding sites and co-regulatory genes, are quite unknown. Therefore, to identify the genome-wide binding sites of PmTCP4 transcription factors and their co-regulatory genes, chromatin immunoprecipitation sequencing (ChIP-Seq) was carried out. ChIP-Seq data produced the maximum enriched peaks in two Japanese apricot cultivars 'Daqiandi' (DQD) and 'Longyan' (LY), which showed that the majority of DNA-protein interactions are relevant and have a significant function in binding sites. Moreover, 720 and 251 peak-associated genes regulated by PmTCP4 were identified in DQD and LY, respectively, and most of them were involved in the flower and pistil development process. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that photosynthesis and oxidative phosphorylation were the most enriched pathways in both cultivars and all identified genes related to these pathways were down-regulated. This study will provide a reference for a better understanding of the PmTCP4 regulatory mechanism during pistil abortion in Japanese apricot.

Genome Wide Identification and Comparative Analysis of the Serpin Gene Family in Brachypodium and Barley

Serpins (serine protease inhibitors) constitute one of the largest and most widely distributed superfamilies of protease inhibitors and have been identified in nearly all organisms. To gain significant insights, a comprehensive in silico analysis of the serpin gene family was carried out in the model plant for temperate grasses Brachypodium distachyon and barley Hordeum vulgare using bioinformatic tools at the genome level for the first time. We identified a total of 27 BdSRPs and 25 HvSRP genes in Brachypodium and barley, respectively, showing an unexpectedly high gene number in these model plants. Gene structure, conserved motifs and phylogenetic comparisons of serpin genes supported the role of duplication events in the expansion and evolution of serpin gene family. Further, purifying selection pressure was found to be a main driving force in the evolution of serpin genes. Genome synteny analysis indicated that BdSRP genes were present in syntenic regions of barley, rice, sorghum and maize, suggesting that they evolved before the divergence of these species from common ancestor. The distinct expression pattern in specific tissues further suggested a specialization of functions during development and in plant defense. These results suggest that the LR serpins (serpins with Leu-Arg residues at P2-P1') identified here can be utilized as candidates for exploitation in disease resistance, pest control and preventing stress-induced cell death. Additionally, serpins were identified that could lead to further research aimed at validating and functionally characterizing the role of potential serpin genes from other plants.

Comparative phenotypic and transcriptomic analysis of Victoria and flame seedless grape cultivars during berry ripening

Grape berry development is a highly coordinated and intricate process. Herein, we analyzed the phenotypic and transcriptomic patterns of Victoria (VT) and Flame Seedless (FS) grape varieties during berry development. Physiological analysis and transcriptomic sequencing were performed at four berry developmental phases. VT berry size was comparatively larger to the FS variety. At maturity, 80 days postanthesis (DPA), the FS soluble solids were 61.8% higher than VT. Further, 4889 and 2802 differentially expressed genes were identified from VT and FS 40 DPA to 80 DPA development stages, respectively. VvSWEET15, VvHXK, and MYB44 genes were up‐regulated during the postanthesis period, while bHLH14, linked to glucose metabolism, was gradually down‐regulated during berry development. These genes may have significant roles in berry development, ripening, and sugar accumulation.

HvHOX9, a novel homeobox leucine zipper transcription factor, positively regulates aluminum tolerance in Tibetan wild barley

Aluminum (Al) toxicity is the primary limiting factor of crop production on acid soils. Tibetan wild barley germplasm is a valuable source of potential genes for breeding barley with acid and Al tolerance. We performed microRNA and RNA sequencing using wild (XZ16, Al-tolerant; XZ61, Al-sensitive) and cultivated (Dayton, Al-tolerant) barley. A novel homeobox-leucine zipper transcription factor, HvHOX9, was identified as a target gene of miR166b and functionally characterized. HvHOX9 was up-regulated by Al stress in XZ16 (but unchanged in XZ61 and Dayton) and was significantly induced only in root tip. Phylogenetic analysis showed that HvHOX9 is most closely related to wheat TaHOX9 and orthologues of HvHOX9 are present in the closest algal relatives of Zygnematophyceae. Barley stripe mosaic virus-induced gene silencing of HvHOX9 in XZ16 led to significantly increased Al sensitivity but did not affect its sensitivity to other metals and low pH. Disruption of HvHOX9 did not change Al concentration in the root cell sap, but led to more Al accumulation in root cell wall after Al exposure. Silencing of HvHOX9 decreased H+ influx after Al exposure. Our findings suggest that miR166b/HvHOX9 play a critical role in Al tolerance by decreasing root cell wall Al binding and increasing apoplastic pH for Al detoxification in the root.

The R2R3-MYB gene family in banana (Musa acuminata): Genome-wide identification, classification and expression patterns

The R2R3-MYB genes comprise one of the largest transcription factor gene families in plants, playing regulatory roles in plant-specific developmental processes, defense responses and metabolite accumulation. To date MYB family genes have not yet been comprehensively identified in the major staple fruit crop banana. In this study, we present a comprehensive, genome-wide analysis of the MYB genes from Musa acuminata DH-Pahang (A genome). A total of 285 R2R3-MYB genes as well as genes encoding three other classes of MYB proteins containing multiple MYB repeats were identified and characterised with respect to structure and chromosomal organisation. Organ- and development-specific expression patterns were determined from RNA-Seq data. For 280 M. acuminata MYB genes for which expression was found in at least one of the analysed samples, a variety of expression patterns were detected. The M. acuminata R2R3-MYB genes were functionally categorised, leading to the identification of seven clades containing only M. acuminata R2R3-MYBs. The encoded proteins may have specialised functions that were acquired or expanded in Musa during genome evolution. This functional classification and expression analysis of the MYB gene family in banana establishes a solid foundation for future comprehensive functional analysis of MaMYBs and can be utilized in banana improvement programmes.

Ethylene biosynthesis and response factors are differentially modulated during the interaction of peach petals with Monilinia laxa or Monilinia fructicola

Monilinia spp. may infect stone fruit at any growth stage, although susceptibility to brown rot depends on both host properties and climatological conditions. This said, no studies deciphering the host response in the interaction between peach blossoms and Monilinia spp. are yet available. This study presents an in-depth characterization of the role of ethylene in the interaction of 'Merrill O'Henry' peach petals (Prunus persica (L.) Batch) with Monilinia laxa and M. fructicola. We investigated the physiological responses of the host and the fungi to the application of ethylene and 1-methylcyclopropene (1-MCP) as well as the molecular patterns associated with the biosynthetic and ethylene-dependent responses during the interaction of both Monilinia species with the host. The incidence of both species was differentially affected by 1-MCP and ethylene; M. laxa was favoured by the enhanced host ethylene production associated with the treatments whereas M. fructicola reduced its infection capacity. Such differences were host-dependent as treatments did not affect growth or colony morphology of Monilinia spp. Besides, host ethylene production was altered in M. laxa inoculated petals, either by the fungus or the host itself. Molecular analysis revealed some important ERFs that could be involved in the different ability of both species to activate a cascade response of peach petals against these pathogens.

Large-Scale Cloning and Comparative Analysis of TaNAC Genes in Response to Stripe Rust and Powdery Mildew in Wheat (Triticum aestivum L.)

The NAM, ATAF1/2, and CUC2 (NAC) transcription factors (TFs) constitute the largest plant-specific TF superfamily, and play important roles in various physiological processes, including stress responses. Stripe rust and powdery mildew are the most damaging of the fungal diseases that afflict wheat (Triticum aestivum L.). However, studies on Triticum aestivum NAC (TaNAC)s’ role in resistance to the two diseases are still limited, especially in an overall comparative analysis of TaNACs responding or not to fungal stress. In the present study, 186 TaNAC transcripts were obtained from the resistant hexaploid wheat line N9134 under fungal stress, and 180 new transcripts were submitted to GenBank. Statistical results show that 35.1% (54/154) of TaNAC genes responded to stripe rust and powdery mildew in the seedling stage. “Abnormal” coding transcripts of differentially expressed (DE)-TaNAC genes in wheat responding to fungal stress were found in a significantly higher proportion (24/117 vs. 8/69, p = 0.0098) than in non-DE-NACs. This hinted that the alternative splicing of TaNAC genes was active in transcriptional or post-transcriptional regulation during plant-pathogen interactions. Full-length NAC proteins were classified into nine groups via phylogenetic analysis. Multiple-sequence alignment revealed diversity in the C-terminal structural organization, but the differentially expressed gene (DEG)-encoding proteins enriched in Subgroups VI and VII were conserved, with WV[L/V]CR amino acid residues in Motif 7 following the NAM domain. Our data that showed TaNAC TFs responded to fungal disease, which was affected by expression levels and by the regulation of multifarious transcript variants. These data for TaNAC responses to stripe rust and/or powdery mildew and their numerous structural variants provide a good resource for NAC function–mechanism analysis in the context of biotic-stress tolerance in wheat.

Characterization and Comparative Analysis of RWP-RK Proteins from Arachis duranensis, Arachis ipaensis, and Arachis hypogaea

RWP-RK proteins are important factors involved in nitrate response and gametophyte development in plants, and the functions of RWP-RK proteins have been analyzed in many species. However, the characterization of peanut RWP-RK proteins is limited. In this study, we identified 16, 19, and 32 RWP-RK members from Arachis duranensis, Arachis ipaensis, and Arachis hypogaea, respectively, and investigated their evolution relationships. The RWP-RK proteins were classified into two groups, RWP-RK domain proteins and NODULE-INCEPTION-like proteins. Chromosomal distributions, gene structures, and conserved motifs of RWP-RK genes were compared among wild and cultivated peanuts. In addition, we identified 12 orthologous gene pairs from the two wild peanut species, 13 from A. duranensis and A. hypogaea, and 13 from A. ipaensis and A. hypogaea. One, one, and seventeen duplicated gene pairs were identified within the A. duranensis, A. ipaensis, and A. hypogaea genomes, respectively. Moreover, different numbers of cis-acting elements in the RWP-RK promoters were found in wild and cultivated species (87 in A. duranensis, 89 in A. ipaensis, and 92 in A. hypogaea), and as a result, many RWP-RK genes showed distinct expression patterns in different tissues. Our study will provide useful information for further functional and evolutionary analysis of the RWP-RK genes.

Overexpression of Maize ZmMYB59 Gene Plays a Negative Regulatory Role in Seed Germination in Nicotiana tabacum and Oryza sativa

MYB transcription factors are involved in many biological processes, including metabolism, stress response and plant development. In our previous work, ZmMYB59 was down-regulated by deep sowing during maize seed germination. However, there are few reports on seed germination regulated by MYB proteins. In this study, to examine its functions during seed germination, Agrobacterium-mediated transformation was exploited to generate ZmMYB59 overexpression (OE) tobacco and rice. In T2 generation transgenic tobacco, germination rate, germination index, vigor index and hypocotyl length were significantly decreased by 25.0-50.9, 34.5-54.4, 57.5-88.3, and 21.9-31.3% compared to wild-type (WT) lines. In T2 generation transgenic rice, above corresponding parameters were notably reduced by 39.1-53.8, 51.4-71.4, 52.5-74.0, and 28.3-41.5%, respectively. On this basis, antioxidant capacity and endogenous hormones were determined. The activities of catalase, peroxidase, superoxide dismutase, ascorbate peroxidase of OE lines were significantly lower than those of WT, suggesting that ZmMYB59 reduced their oxidation resistance. As well, ZmMYB59 overexpression extremely inhibited the synthesis of gibberellin A1 (GA1) and cytokinin (CTK), and promoted the synthesis of abscisic acid (ABA) concurrently. Taken together, it proposed that ZmMYB59 was a negative regulator during seed germination in tobacco and rice, which also contributes to illuminate the molecular mechanisms regulated by MYB transcription factors.

Genome-wide identification and expression analyses of R2R3-MYB transcription factor genes from two Orchid species

MYB transcription factors play important roles in different plant biological processes during plant growth, development and stress response. In this study, 101 (DoMYB1-101) and 99 (PaMYB1-99) R2R3-MYB genes were identified in the genomes of Dendrobium officinale and Phalaenopsis aphrodite, respectively. To classify the isolated candidate genes, the R2R3-MYB genes from A. thaliana were selected as references. As a result, all identified DoMYB and PaMYB genes were classified into 22 subfamilies. Phylogenetic analysis revealed that S21 had the largest number of members of all the subfamilies. The numbers of introns, exons and conserved sequences in all of the identified genes are different. In addition, 20 DoMYB genes from six subfamilies were selected for further analysis of tissue-specific expression and responses to various abiotic stresses treatments. The results showed that all of the DoMYB genes in S4 and S19 subfamilies exhibited the highest relative expression levels in flowers. And five DoMYB genes including DoMYB31, DoMYB40, DoMYB49, DoMYB52 and DoMYB54, responded to the stress response. These results may provide useful information for further studies of the R2R3-MYB gene family.

De novo assembly of transcriptome and genome-wide identification reveal GA3 stress-responsive WRKY transcription factors involved in fiber formation in jute (Corchorus capsularis)

BACKGROUND

WRKY is a group of transcription factors (TFs) that play a vital role in plant growth, development, and stress tolerance. To date, none of jute WRKY (CcWRKY) genes have been identified, even if jute (Corchorus capsularis) is one of the most important natural fiber crops in the world. Little information about the WRKY genes in jute is far from sufficient to understand the molecular mechanism of bast fiber biosynthesis.

RESULTS

A total of 244,489,479 clean reads were generated using Illumina paired-end sequencing. De novo assembly yielded 90,982 unigenes with an average length of 714 bp. By sequence similarity searching for known proteins, 48,896 (53.74%) unigenes were annotated. To mine the CcWRKY TFs and identify their potential function, the search for CcWRKYs against the transcriptome data of jute was performed, and a total of 43 CcWRKYs were identified in this study. The gene structure, phylogeny, conserved domain and three-dimensional structure of protein were analyzed by bioinformatics tools of GSDS2.0, MEGA7.0, DNAMAN5.0, WebLogo 3 and SWISS-MODEL respectively. Phylogenetic analysis showed that 43 CcWRKYs were divided into three groups: I, II and III, containing 9, 28, and 6 members respectively, according to the WRKY conserved domain features and the evolution analysis with Arabidopsis thaliana. Gene structure analysis indicated that the number of exons of these CcWRKYs varied from 3 to 11. Among the 43 CcWRKYs, 10, 2, 2, and 14 genes showed higher expression in leaves, stem sticks, stem barks, and roots at the vigorous vegetative growth stage, respectively. Moreover, the expression of 21 of 43 CcWRKYs was regulated significantly with secondary cell wall biosynthesis genes using FPKM and RT-qPCR by GA₃ stress to a typical GA₃ sensitive dwarf germplasm in comparison to an elite cultivar in jute. The Cis-element analysis showed that promoters of these 21 CcWRKYs had 1 to 4 cis-elements involved in gibberellin-responsiveness, suggesting that they might regulate the development of bast fiber in response to GA₃ stress.

CONCLUSIONS

A total of 43 CcWRKYs were identified in jute for the first time. Analysis of phylogenetic relationship and gene structure revealed that these CcWRKYs might have a functional diversity. Expression analysis showed 21 TFs as GA₃ stress responsive genes. The identification of these CcWRKYs and the characterization of their expression pattern will provide a basis for future clarification of their functions in bast fiber development in jute.

Genome-wide identification and characterization of SRO gene family in wheat: Molecular evolution and expression profiles during different stresses

SRO (SIMILAR TO RCD ONE), a type of plant-specific small protein family, play important roles in plant growth and development, as well as in response to biotic/abiotic stresses. Although characterization of SROs have been performed in model plants, little is known about their function in wheat, especially under stress conditions. In this study, 30 SRO genes were identified from the wheat genome (TaSRO). They were phylogenetically separated into two groups with distinct structures. The cis-regulatory elements in the promoter region of TaSROs were analyzed and numerous elements functionally associated with stress responding and hormones were interpreted, implying the reason for induction expression patterns of TaSROs during abiotic and biotic stresses in wheat. Whole-genome replication events in the SRO gene family of wheat and seven other species (Arabidopsis thaliana, rice, maize, barley, soybean, upland cotton, and cucumber) were analyzed, resulting in 1, 12, 9, 23, 6, 5, and 3 of gene pairs, respectively. The tissue-specific expression pattern profiling revealed that most TaSROs are highly expressed in one or more tissues and may play an important role in wheat growth and development. In addition, qRT-PCR results further confirmed that these TaSRO genes are involved in wheat stress response. In summary, our study laid a theoretical basis for molecular function deciphering of TaSROs, especially in plant hormones and biotic/abiotic stress responses.

Genomic Characterization and Expression Analysis of Basic Helix-Loop-Helix (bHLH) Family Genes in Traditional Chinese Herb Dendrobium officinale

Dendrobium officinale Kimura et Migo is of great importance as a traditional Chinese herb due to its abundant metabolites. The family of basic helix-loop-helix (bHLH) transcription factors widely exists in plants and plays an essential role in plant growth and development, secondary metabolism as well as responses to environmental changes. However, there is limited information on bHLH genes in D. officinale. In the present study, a total of 98 putative DobHLH genes were identified at the genomic level, which could be classified into 18 clades. Gene structures and conserved motifs in DobHLH genes showed high conservation during their evolution. The conserved amino acids and DNA bindings of DobHLH proteins were predicted, both of which are pivotal for their function. Furthermore, gene expression from eight tissues showed that some DobHLH genes were ubiquitously expressed while other DobHLH genes were expressed in the specific tissues. Expressional changes of DobHLH genes under MeJA and ABA treatments were detected by qRT-PCR. The protein-protein interactions between DobHLHs were predicted and several interactions were confirmed by yeast two hybrid. Therefore, our results here contribute to the understanding of bHLH genes in D. officinale and lay a foundation for the further functional study of its biological processes.

Genome-Wide Analysis of the Role of NAC Family in Flower Development and Abiotic Stress Responses in Cleistogenes songorica

Plant-specific NAC (NAM, ATAF, CUC) transcription factor (TF) family plays important roles in biological processes such as plant growth and response to stress. Nevertheless, no information is known about NAC TFs in Cleistogenes songorica, a prominent xerophyte desert grass in northwestern China. In this study, 162 NAC genes were found from the Cleistogenes songorica genome, among which 156 C. songoricaNAC (CsNAC) genes (96.3%) were mapped onto 20 chromosomes. The phylogenetic tree constructed by CsNAC and rice NAC TFs can be separated into 14 subfamilies. Syntenic and Ka/Ks analyses showed that CsNACs were primarily expanded by genomewide replication events, and purifying selection was the primary force driving the evolution of CsNAC family genes. The CsNAC gene expression profiles showed that 36 CsNAC genes showed differential expression between cleistogamous (CL) and chasmogamous (CH) flowers. One hundred and two CsNAC genes showed differential expression under heat, cold, drought, salt and ABA treatment. Twenty-three CsNAC genes were commonly differentially expressed both under stress responses and during dimorphic floret development. Gene Ontology (GO) annotation, coexpression network and qRT-PCR tests revealed that these CsNAC genes may simultaneously regulate dimorphic floret development and the response to stress. Our results may help to characterize the NAC transcription factors in C. songorica and provide new insights into the functional research and application of the NAC family in crop improvement, especially in dimorphic floret plants.

Genome-Wide Characterization and Comparative Analysis of MYB Transcription Factors in Ganoderma Species

Numerous studies in plants have shown the vital roles of MYB transcription factors in signal transduction, developmental regulation, biotic/abiotic stress responses and secondary metabolism regulation. However, less is known about the functions of MYBs in Ganoderma. In this study, five medicinal macrofungi of genus Ganoderma were subjected to a genome-wide comparative analysis of MYB genes. A total of 75 MYB genes were identified and classified into four types: 1R-MYBs (52), 2R-MYBs (19), 3R-MYBs (2) and 4R-MYBs (2). Gene structure analysis revealed varying exon numbers (3-14) and intron lengths (7-1058 bp), and noncanonical GC-AG introns were detected in G. lucidum and G. sinense. In a phylogenetic analysis, 69 out of 75 MYB genes were clustered into 15 subgroups, and both single-copy orthologous genes and duplicated genes were identified. The promoters of the MYB genes harbored multiple cis-elements, and specific genes were co-expressed with the G. lucidum MYB genes, indicating the potential roles of these MYB genes in stress response, development and metabolism. This comprehensive and systematic study of MYB family members provides a reference and solid foundation for further functional analysis of MYB genes in Ganoderma species.

The genome of kenaf (Hibiscus cannabinus L.) provides insights into bast fibre and leaf shape biogenesis

Kenaf is an annual crop that is widely cultivated as a source of bast (phloem) fibres, the phytoremediation of heavy metal-contaminated farmlands and textile-relevant compounds. Leaf shape played a unique role in kenaf improvement, due to the inheritance as a single locus and the association with fibre development in typical lobed-leaf varieties. Here we report a high-quality genome assembly and annotation for var. 'Fuhong 952' with 1078 Mbp genome and 66 004 protein-coding genes integrating single-molecule real-time sequencing, a high-density genetic map and high-throughput chromosome conformation capture techniques. Gene mapping assists the identification of a homeobox transcription factor LATE MERISTEM IDENTITY 1 (HcLMI1) gene controlling lobed-leaf. Virus-induced gene silencing (VIGS) of HcLMI1 in a lobed-leaf variety was critical to induce round (entire)-like leaf formation. Candidate genes involved in cell wall formation were found in quantitative trait loci (QTL) for fibre yield and quality-related traits. Comparative genomic and transcriptome analyses revealed key genes involved in bast fibre formation, among which there are twice as many cellulose synthase A (CesA) genes due to a recent whole-genome duplication after divergence from Gossypium. Population genomic analysis showed two recent population bottlenecks in kenaf, suggesting domestication and improvement process have led to an increase in fibre biogenesis and yield. This chromosome-scale genome provides an important framework and toolkit for sequence-directed genetic improvement of fibre crops.

Genome-wide identification and evolution of Dof transcription factor family in cultivated and ancestral cotton species

The DNA-binding with one finger (Dof) proteins are transcription factors involved in many biological processes in plants. To predict the evolutionary pattern, a genome-wide in-silico analysis of Dof TFs family in diploid (Gossypium arboreum and Gossypium raimondii) and allotetraploid (Gossypium hirsutum and Gossypium barbadense) cotton species were carried out. In G. arboreum, we have identified 58 non-redundant genes encoding Dof proteins renamed as GaDof (G. arboreum Dof), 55 Dof genes were identified in G. raimondii (GrDof), 89 were predicted ffrom G. hirsutum (GhDof) and the highest, 110 Dof genes were identified in G. barbadense (GbDof). The phylogenetic analysis, physical location, gene structure, conserved domain analyses were also investigated for G. arboreum, G. raimondii, and G. hirsutum. The gene expression pattern in G. hirsutum, at different growth stages, revealing the probable involvement of some GhDof genes in growth and development. These genes may improve seed germination and growth in cotton.

Transcriptional regulation and expression network responding to cadmium stress in a Cd-tolerant perennial grass Poa Pratensis

Kentucky bluegrass has good capability to absorb and accumulate cadmium (Cd) through developed root system, thus having potential phytoremediation function in Cd contaminated soils. Understanding the molecular mechanisms of Cd tolerance and accumulation in this species will be crucial to generating novel Cd-tolerance cultivars through genetic improvement, while it has not well documented yet. In the present study, comparative transcriptome analysis was performed for the seedlings of high Cd-tolerant genotype (M) and low Cd-tolerant genotype (R) under Cd stress. A total of 7022 up-regulated and 1033 down-regulated transcripts were identified in M genotype, whereas, only 850 up-regulated and 846 down-regulated transcripts were detected in R. Further transcriptional regulation analysis in M genotype showed that Dof, MADS25, BBR-BPC, B3, bZIP23 and MYB30 might be the hub transcription factors in response to Cd stress due to the orchestrated multiple functional genes associated with carbohydrate, lipid and secondary metabolism, as well as signal transduction. Differential expressed genes involved in auxin, ethylene, brassinosteroid and ABA signalling formed signal transduction cascades, which interacted with hub transcription factors, thereby finally orchestrated the expression of multiple genes associated with cell wall and membrane stability, cell elongation and Cd tolerance, including IAAs, ARFs, SnRK2, PP2C, PIFs, BES1/BZR1, CCR, CAD, FATB, fabF and HACD. Additionally, post-transcriptional modification of CIPKs, MAPKs, WAXs, UBCs, and E3 ubiquitin ligases were identified and also involved in plant signalling pathways and abiotic resistance. The study could contribute to our understanding the transcriptional regulation and complex internal network associated with Cd tolerance in Kentucky bluegrass.