Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes

Transcriptome profiles identify the common responsive genes to drought stress in two Elymus species

Elymus, the largest genus of the Triticeae Dumort, is a forage grass in the Qinghai-Tibetan Plateau, where the climate has gradually become increasingly dry in recent years. To understand the mechanisms of the response to drought stress in Elymus species, we first investigated physiological and biochemical responses to polyethylene glycol (PEG-6000) simulated drought stress in two Elymus species, Elymus nutans and Elymus sibiricus, and found that E. nutans was more tolerant to drought stress than E. sibiricus. De novo transcriptome analysis of these two Elymus species treated with or without 10 % PEG-6000 revealed that a total of 1695 unigenes were commonly regulated by drought treatment in these two Elymus species, with 1614 unigenes up-regulated and 81 unigenes down-regulated. The coexpressed differentially expressed genes (DEGs) were enriched in regulation of transcription and gene expression in the GO database. KEGG pathway analysis indicated plant hormone signaling transduction were mostly enriched in co-expressed DEGs. Furthermore, genes annotated in the plant hormone signaling transduction were screened from co-expressed DEGs, and found that abscisic acid plays the major role in the drought stress tolerance of Elymus. Meanwhile, transcription factors screened from co-expressed DEGs were mainly classified into the ERF subfamily and WRKY, DREB, and HSF family members. Our results provide further reference for studying the response mechanism and culturing highly tolerant grasses of the Elymus species under drought stress.

Understanding the early cold response mechanism in IR64 indica rice variety through comparative transcriptome analysis

Background

Cellular reprogramming in response to environmental stress involves alteration of gene expression, changes in the protein and metabolite profile for ensuring better stress management in plants. Similar to other plant species originating in tropical and sub-tropical areas, indica rice is highly sensitive to low temperature that adversely affects its growth and grain productivity. Substantial work has been done to understand cold induced changes in gene expression in rice plants. However, adequate information is not available for early gene expression, especially in indica variety. Therefore, a transcriptome profile was generated for cold shock treated seedlings of IR64 variety to identify early responsive genes.

Results

The functional annotation of early DEGs shows enrichment of genes involved in altered membrane rigidity and electrolytic leakage, the onset of calcium signaling, ROS generation and activation of stress responsive transcription factors in IR64. Gene regulatory network suggests that cold shock induced Ca2+ signaling activates DREB/CBF pathway and other groups of transcription factors such as MYB, NAC and ZFP; for activating various cold-responsive genes. The analysis also indicates that cold induced signaling proteins like RLKs, RLCKs, CDPKs and MAPKK and ROS signaling proteins. Further, several late-embryogenesis-abundant (LEA), dehydrins and low temperature-induced-genes were upregulated under early cold shock condition, indicating the onset of water-deficit conditions. Expression profiling in different high yielding cultivars shows high expression of cold-responsive genes in Heera and CB1 indica varieties. These varieties show low levels of cold induced ROS production, electrolytic leakage and high germination rate post-cold stress, compared to IR36 and IR64. Collectively, these results suggest that these varieties may have improved adaptability to cold stress.

Conclusions

The results of this study provide insights about early responsive events in Oryza sativa l.ssp. indica cv IR64 in response to cold stress. Our data shows the onset of cold response is associated with upregulation of stress responsive TFs, hydrophilic proteins and signaling molecules, whereas, the genes coding for cellular biosynthetic enzymes, cell cycle control and growth-related TFs are downregulated. This study reports that the generation of ROS is integral to the early response to trigger the ROS mediated signaling events during later stages.

Genome-wide identification and characterization of UDP-glucose dehydrogenase family genes in moso bamboo and functional analysis of PeUGDH4 in hemicellulose synthesis

Uridine diphosphate glucose dehydrogenases (UGDHs) are critical for synthesizing many nucleotide sugars and help promote the carbohydrate metabolism related to cell wall synthesis. In plants, UGDHs are encoded by a small gene family. Genome-wide analyses of these genes have been conducted in Glycine max and Arabidopsis thaliana, however, the UGDH gene family has not been comprehensively and systematically investigated in moso bamboo (Phyllostachys edulis), which is a special woody grass monocotyledonous species. In this study, we identified nine putative PeUGDH genes. Furthermore, analysis of gene duplication events and divergences revealed that the expansion of the PeUGDH family was mainly due to segmental and tandem duplications approximately 4.76-83.16 million years ago. An examination of tissue-specific PeUGDH expression indicated that more than 77% of the genes were predominantly expressed in the stem. Based on relative expression levels among PeUGDH members in different tissues in moso bamboo, PeUGDH4 was selected for detailed analysis. The results of subcellular localization indicated that PeUGDH4-GFP fusion proteins was observed to be localized in the cytoplasm. The ectopic overexpression of PeUGDH4 in Arabidopsis significantly increased the contents of hemicellulose and soluble sugar, suggesting that PeUGDH4 acts as a key enzyme involved in bamboo cell wall synthesis.

Genome-wide identification and characterization of cucumber bHLH family genes and the functional characterization of CsbHLH041 in NaCl and ABA tolerance in Arabidopsis and cucumber

BACKGROUND

The basic/helix-loop-helix (bHLH) transcription factor family exists in all three eukaryotic kingdoms as important participants in biological growth and development. To date, the comprehensive genomic and functional analyses of bHLH genes has not been reported in cucumber (Cucumis sativus L.).

RESULTS

Here, a total of 142 bHLH genes were identified and classified into 32 subfamilies according to the conserved motifs, phylogenetic analysis and gene structures in cucumber. The sequences of CsbHLH proteins were highly conserved based on the results of multiple sequence alignment analyses. The chromosomal distribution, synteny analysis, and gene duplications of these 142 CsbHLHs were further analysed. Many elements related to stress responsiveness and plant hormones were present in the promoter regions of CsbHLH genes based on a cis-element analysis. By comparing the phylogeny of cucumber and Arabidopsis bHLH proteins, we found that cucumber bHLH proteins were clustered into different functional clades of Arabidopsis bHLH proteins. The expression analysis of selected CsbHLHs under abiotic stresses (NaCl, ABA and low-temperature treatments) identified five CsbHLH genes that could simultaneously respond to the three abiotic stresses. Tissue-specific expression profiles of these five genes were also analysed. In addition, 35S:CsbHLH041 enhanced the tolerance to salt and ABA in transgenic Arabidopsis and in cucumber seedlings, suggesting CsbHLH041 is an important regulator in response to abiotic stresses. Lastly, the functional interoperability network among the CsbHLH proteins was analysed.

CONCLUSION

This study provided a good foundation for further research into the functions and regulatory mechanisms of CsbHLH proteins and identified candidate genes for stress resistance in cucumber.

Advances in AP2/ERF super-family transcription factors in plant

In the whole life process, many factors including external and internal factors affect plant growth and development. The morphogenesis, growth, and development of plants are controlled by genetic elements and are influenced by environmental stress. Transcription factors contain one or more specific DNA-binding domains, which are essential in the whole life cycle of higher plants. The AP2/ERF (APETALA2/ethylene-responsive element binding factors) transcription factors are a large group of factors that are mainly found in plants. The transcription factors of this family serve as important regulators in many biological and physiological processes, such as plant morphogenesis, responsive mechanisms to various stresses, hormone signal transduction, and metabolite regulation. In this review, we summarized the advances in identification, classification, function, regulatory mechanisms, and the evolution of AP2/ERF transcription factors in plants. AP2/ERF family factors are mainly classified into four major subfamilies: DREB (Dehydration Responsive Element-Binding), ERF (Ethylene-Responsive-Element-Binding protein), AP2 (APETALA2) and RAV (Related to ABI3/VP), and Soloists (few unclassified factors). The review summarized the reports about multiple regulatory functions of AP2/ERF transcription factors in plants. In addition to growth regulation and stress responses, the regulatory functions of AP2/ERF in plant metabolite biosynthesis have been described. We also discussed the roles of AP2/ERF transcription factors in different phytohormone-mediated signaling pathways in plants. Genomic-wide analysis indicated that AP2/ERF transcription factors were highly conserved during plant evolution. Some public databases containing the information of AP2/ERF have been introduced. The studies of AP2/ERF factors will provide important bases for plant regulatory mechanisms and molecular breeding.

Gene Regulation via the Combination of Transcription Factors in the INDETERMINATE DOMAIN and GRAS Families

INDETERMINATE DOMAIN (IDD) family proteins are plant-specific transcription factors. Some Arabidopsis IDD (AtIDD) proteins regulate the expression of SCARECROW (SCR) by interacting with GRAS family transcription factors SHORT-ROOT (SHR) and SCR, which are involved in root tissue formation. Some AtIDD proteins regulate genes involved in the synthesis (GA3ox1) or signaling (SCL3) of gibberellic acid (GA) by interacting with DELLA proteins, a subfamily of the GRAS family. We analyzed the DNA binding properties and protein–protein interactions of select AtIDD proteins. We also investigated the transcriptional activity of the combination of AtIDD and GRAS proteins (AtIDD proteins combined with SHR and SCR or with REPRESSOR of ga1-3 (RGA)) on the promoters of SCR,SCL3, and GA3ox1 by conducting a transient assay using Arabidopsis culture cells. Our results showed that the SCR promoter could be activated by the IDD and RGA complexes and that the SCL3 and GA3ox1 promoters could be activated by the IDD, SHR, and SCR complexes, indicating the possibility that these complexes regulate and consequently coordinate the expression of genes involved in GA synthesis (GA3ox1), GA signaling (SCL3), and root formation (SCR).

Characterisation of the ERF102 to ERF105 genes of Arabidopsis thaliana and their role in the response to cold stress

The four phylogenetically closely related ERF102 to ERF105 transcription factors of Arabidopsis thaliana are regulated by different stresses and are involved in the response to cold stress. The ETHYLENE RESPONSE FACTOR (ERF) genes of Arabidopsis thaliana form a large family encoding plant-specific transcription factors. Here, we characterise the four phylogenetically closely related ERF102/ERF5, ERF103/ERF6, ERF104 and ERF105 genes. Expression analyses revealed that these four genes are similarly regulated by different hormones and abiotic stresses. Analyses of tissue-specific expression using promoter:GUS reporter lines revealed their predominant expression in root tissues including the root meristem (ERF103), the quiescent center (ERF104) and the root vasculature (all). All GFP-ERF fusion proteins were nuclear-localised. The analysis of insertional mutants, amiRNA lines and 35S:ERF overexpressing transgenic lines indicated that ERF102 to ERF105 have only a limited impact on regulating shoot and root growth. Previous work had shown a role for ERF105 in the cold stress response. Here, measurement of electrolyte leakage to determine leaf freezing tolerance and expression analyses of cold-responsive genes revealed that the combined activity of ERF102 and ERF103 is also required for a full cold acclimation response likely involving the CBF regulon. These results suggest a common function of these ERF genes in the response to cold stress.

The PalERF109 transcription factor positively regulates salt tolerance via PalHKT1;2 in Populus alba var. pyramidalis

Salinity restricts the growth of trees to varying extents, but the regulatory mechanisms involved in their varying salt tolerance are largely unknown. In an effort to elucidate these mechanisms, we identified a total of 99 genes in the Ethylene Responsive Factor (ERF) family of transcription factors and examined their expression patterns under salt stress in Populus alba var. pyramidalis. We found that a B4 group gene, PalERF109, was rapidly induced by salt treatment and preferentially expressed in stems and petioles, where it is probably involved in transport of ions and water in xylem. Overexpression of PalERF109 enhanced the salt tolerance of the poplar, and further analysis showed that it directly upregulated a high-affinity K+transporter (HKT) gene, PalHKT1;2. The results clearly indicate that PalERF109 enhances salt tolerance at least partially through direct activation of PalHKT1;2 and extends understanding of the roles of ERF genes in tree stress responses.

Genome-Wide Characterization of GRAS Family and Their Potential Roles in Cold Tolerance of Cucumber (Cucumis sativus L.)

Cucumber (Cucumis sativus L.) is one of the most important cucurbit vegetables but is often subjected to stress during cultivation. GRAS (gibberellic acid insensitive, repressor of GAI, and scarecrow) genes encode a family of transcriptional factors that regulate plant growth and development. In the model plant Arabidopsis thaliana, GRAS family genes function in formation of axillary meristem and root radial structure, phytohormone (gibberellin) signal transduction, light signal transduction and abiotic/biological stress. In this study, a gene family was comprehensively analyzed from the aspects of evolutionary tree, gene structure, chromosome location, evolutionary and expression pattern by means of bioinformatics; 37 GRAS gene family members have been screened from cucumber. We reconstructed an evolutionary tree based on multiple sequence alignment of the typical GRAS domain and conserved motif sequences with those of other species (A. thaliana and Solanum lycopersicum). Cucumber GRAS family was divided into 10 groups according to the classification of Arabidopsis and tomato genes. We conclude that tandem and segmental duplication have played important roles in the expansion and evolution of the cucumber GRAS (CsaGRAS) family. Expression patterns of CsaGRAS genes in different tissues and under cold treatment, combined with gene ontology annotation and interaction network analysis, revealed potentially different functions for CsaGRAS genes in response to cold tolerance, with members of the SHR, SCR and DELLA subfamilies likely playing important roles. In conclusion, this study provides valuable information and candidate genes for improving cucumber tolerance to cold stress.

Genome-Wide Analysis of Basic Helix-Loop-Helix Transcription Factors to Elucidate Candidate Genes Related to Fruit Ripening and Stress in Banana (Musa acuminata L. AAA Group, cv. Cavendish)

The basic helix-loop-helix (bHLH) proteins are a superfamily of transcription factors (TFs) that can bind to specific DNA target sites, playing a central role in a wide range of metabolic, physiological, and developmental processes in higher organisms. However, no systemic analysis of bHLH TFs has been reported in banana, a typical climacteric fruit in tropical and subtropical regions. In our study, 259 MabHLH TF genes were identified in the genome of Musa acuminata (A genome), and phylogenetic analysis indicated that these MabHLHs could be classified into 23 subfamilies with the bHLHs from rice and Arabidopsis. The amino acid sequences of the bHLH domain in all MabHLH protein sequences were quite conserved, especially Arg-12, Arg-13, Leu-23, and Leu-79. Distribution mapping results showed that 258 MabHLHs were localized on the 11 chromosomes in the M. acuminata genome. The results indicated that 40.7% of gene duplication events were located in collinear fragments, and segmental duplications might have played a key role in the expansion of MabHLHs. Moreover, the expression profiles of MabHLHs in different fruit development and ripening stages and under various abiotic and biotic stresses were investigated using available RNA-sequencing data to obtain fruit development, ripening-specific, and stress-responsive candidate genes. Finally, a co-expression network of MabHLHs was constructed by weighted gene co-expression network analysis to elucidate the MabHLHs that might participate in important metabolic biosynthesis pathways in banana during development and the response to stress. A total of 259 MabHLHs were identified, and their sequence features, conserved domains, phylogenetic relationships, chromosomal distributions, gene duplications, expression profiles, and co-expression networks were investigated. This study systematically identified the MabHLHs in the M. acuminata genome at the genome-wide level, providing important candidate genes for further functional analysis. These findings improve our understanding of the molecular basis of developmental and stress tolerance in an important banana cultivar.

Genome-wide identification and expression of YABBY genes family during flower development in Punica granatum L

The plant-specific YABBY transcription factors have important biological roles in plant morphogenesis, growth and development. In this study, we identified six YABBY genes in pomegranate (Punica granatum) and characterized their expression pattern during flower development. Six PgYABBY genes were divided into five subfamilies (YAB1/3, YAB2, INO, CRC, and YAB5), based on protein sequence, motifs and similarity of exon-intron structure. Next, analysis of putative cis-acting element showed that PgYABBYs contained lots of hormone response and stress response elements. Subsequently, gene function prediction and protein-protein network analysis showed that PgYABBYs were associated with the development of apical meristem, flower, carpel, and ovule. Analysis of PgYABBY genes expression in various structures and organs suggested that PgYABBYs were highly activated in flower, leaf and seed coat. Analysis of expression during flower development in pomegranate showed that PgINO might play critical role in regulating the differentiation of flowers. This study provided a theoretical basis for function research and utilization of YABBY genes in pomegranate.

Mutations in orthologous PETALOSA TOE-type genes cause a dominant double-flower phenotype in phylogenetically distant eudicots

The double-flower phenotype has been selected by humans for its attractiveness in various plant species and it is of great commercial value for the ornamental market. In this study we investigated the genetic determinant of the dominant double-flower trait in carnation, petunia, and Rosa rugosa, and identified mutant alleles of TARGET OF EAT (TOE)-type genes characterized by a disruption of the miR172 target sequence and of the C-terminal portion of the encoded protein. Despite the phylogenetic distance between these eudicots, which diverged in the early Cretaceous, the orthologous genes carrying these mutations all belong to a single TOE-type subgroup, which we name as PETALOSA (PET). Homology searches allowed us to identify PET sequences in various other species. To confirm the results from naturally occurring mutations, we used CrispR-Cas9 to induce lesions within the miR172 target site of Nicotiana tabacum PET genes, and this resulted in the development of supernumerary petaloid structures. This study describes pet alleles in economically important ornamental species and provides evidence about the possibility of identifying and engineering PET genes to obtain the desirable double-flower trait in different plants.

Transcriptome profiling reveals insertional mutagenesis suppressed the expression of candidate pathogenicity genes in honeybee fungal pathogen, Ascosphaera apis

Chalkbrood disease is caused by Ascosphaera apis which severely affects honeybee brood. Spore inoculation experiments shown pathogenicity varies among different strains and mutants, however, the molecular mechanism of pathogenicity is unclear. We sequenced, assembled and annotated the transcriptomes of wild type (SPE1) and three mutants (SPE2, SPE3 and SPE4) with reduced pathogenicity that were constructed in our previous study. Illumina sequencing generated a total of 394,910,604 clean reads and de novo Trinity-based assembled into 12,989 unigenes, among these, 9,598 genes were successfully annotated to known proteins in UniProt database. A total of 172, 3,996, and 650 genes were up-regulated and 4,403, 2,845, and 3,016 genes were down-regulated between SPE2-SPE1, SPE3-SPE1, and SPE4-SPE1, respectively. Overall, several genes with a potential role in fungal pathogenicity were detected down-regulated in mutants including 100 hydrolytic enzymes, 117 transcriptional factors, and 47 cell wall related genes. KEGG pathway enrichment analysis reveals 216 genes involved in nine pathways were down-regulated in mutants compared to wild type. The down-regulation of more pathways involved in pathogenicity in SPE2 and SPE4 than SPE3 supports their lower pathogenicity during in-vitro bioassay experiment. Expression of 12 down-regulated genes in mutants was validated by quantitative real time PCR. This study provides valuable information on transcriptome variation caused by mutation for further functional validation of candidate pathogenicity genes in A. apis.

Envisioning the immune interactome in Arabidopsis

During plant-pathogen interaction, immune targets were regulated by protein-protein interaction events such as ligand-receptor/co-receptor, kinase-substrate, protein sequestration, activation or repression via post-translational modification and homo/oligo/hetro-dimerisation of proteins. A judicious use of molecular machinery requires coordinated protein interaction among defence components. Immune signalling in Arabidopsis can be broadly represented in successive or simultaneous steps; pathogen recognition at cell surface, Ca2+ and reactive oxygen species signalling, MAPK signalling, post-translational modification, transcriptional regulation and phyto-hormone signalling. Proteome wide interaction studies have shown the existence of interaction hubs associated with physiological function. So far, a number of protein interaction events regulating immune targets have been identified, but their understanding in an interactome view is lacking. We focussed specifically on the integration of protein interaction signalling in context to plant-pathogenesis and identified the key targets. The present review focuses towards a comprehensive view of the plant immune interactome including signal perception, progression, integration and physiological response during plant pathogen interaction.

Transcription factors as molecular switches to regulate drought adaptation in maize

Here, we reviewed major transcription factors of maize that confer drought stress tolerance, and their target genes and involved signaling pathway. Transcription factors in maize can be promising candidates for improving comprehensive resistance of multiple environmental stimuli. Adverse environmental stress is the main influencing factor affecting plant growth and reproduction, which poses tremendous threats to sustainable agriculture development and crops productivity worldwide. Among various abiotic stress factors, drought is the most vital adversity with the characteristics of frequent occurrences, long duration, and globality. Maize (Zea mays L.) is a major source of food supply for human being and livestock and recently for biofuel. Maize is the crop that is highly susceptible to drought stress. Drought stress tolerance in plants is quite complex, and it is not ideal to improve crop drought tolerance through a single resistant gene. Transcription factors participate in the regulation of plant growth and development, morphogenesis, and various environmental stress responses via regulating the expression level of their target stress-responsive genes independently or cross talk with other transcription factors, thereby the comprehensive resistance of multiple stresses in crops is improved. This review aims to summarize the major drought-tolerant transcription factors in maize and their regulatory network. With the continuous identification of maize transcription factors, more will be demonstrated to confer drought tolerance either in maize or other crops. It is expected that the transcription factors will greatly enrich the functional gene resources and will be a benefit to drought-tolerant maize cultivars breeding.

Analysis of differentially expressed genes and pathways associated with male sterility lines in watermelon via bulked segregant RNA-seq

Genic male sterility (GMS) is a common and important trait, which is widely used for the production of hybrid seeds. However, the molecular mechanism of GMS in watermelon remains poorly understood. In this study, we comparatively analyzed the transcriptome profiles of sterile and fertile floral buds using the bulked segregant analysis (BSA) and transcriptome sequencing (RNA-seq). A total of 2507 differentially expressed genes (DEGs) including 593 up-regulated and 1914 down-regulated, were identified to be related to male sterility in watermelon line Se18. Gene ontology (GO) analysis showed that 57 GO terms were significantly enriched, while Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis revealed plant hormone signal transduction, glycolysis/gluconeogenesis, starch and sucrose metabolism, plant-pathogen interaction, phenylpropanoid biosynthesis pathways were obviously enriched. Furthermore, the efficiency of the RNA-seq analysis was validated by quantitative real-time PCR (qRT-PCR). Among the DEGs, some valuable candidate genes involved in pollen development were identified, such as gene Cla000029, a bHLH transcription factor and homologous to MS1 in Arabidopsis. Moreover, other DEGs including MYB gene Cla012590 (MYB26), Cla017100 (MYB21), etc., also provide useful information for further understanding the function of key genes involved in pollen development. This study provides new insights into the global network of male sterility in watermelon.

Epigenetic memory and priming in plants

In nature, plants are regularly exposed to biotic and abiotic stress conditions. These conditions create potential risks for survival. Plants have evolved in order to compete with these stress conditions through physiological adjustments that are based on epigenetic background. Thus, the ecological signals create different levels of stress memory. Recent studies have shown that this stress-induced environmental memory is mediated by epigenetic mechanisms that have fundamental roles in the aspect of controlling gene expression via DNA methylation, histone modifications and, small RNAs and these modifications could be transmitted to the next generations. Thus, they provide alternative mechanisms to constitute stress memories in plants. In this review, we summarized the epigenetic memory mechanisms related with biotic and abiotic stress conditions, and relationship between priming and epigenetic memory in plants by believing that it can be useful for analyzing memory mechanisms and see what is missing out in order to develop plants more resistant and productive under diverse environmental cues.

Integrated Analysis of mRNA and microRNA Elucidates the Regulation of Glycyrrhizic Acid Biosynthesis in Glycyrrhiza uralensis Fisch

Licorice (Glycyrrhiza) is a staple Chinese herbal medicine in which the primary bioactive compound is glycyrrhizic acid (GA), which has important pharmacological functions. To date, the structural genes involved in GA biosynthesis have been identified. However, the regulation of these genes in G. uralensis has not been elucidated. In this study, we performed a comprehensive analysis based on the transcriptome and small RNAome by high-throughput sequencing. In total, we identified 18 structural GA genes and 3924 transporter genes. We identified genes encoding 2374 transporters, 1040 transcription factors (TFs), 262 transcriptional regulators (TRs) and 689 protein kinases (PKs), which were coexpressed with at least one structural gene. We also identified 50,970 alternative splicing (AS) events, in which 17 structural genes exhibited AS. Finally, we also determined that miRNAs potentially targeted 4 structural genes, and 318, 8, and 218 miRNAs potentially regulated 150 TFs, 34 TRs, and 88 PKs, respectively, related to GA. Overall, the results of this study helped to elucidate the gene expression and regulation of GA biosynthesis in G. uralensis, provided a theoretical basis for the synthesis of GA via synthetic biology, and laid a foundation for the cultivation of new varieties of licorice with high GA content.

Dynamics of Global Gene Expression and Regulatory Elements in Growing Brachypodium Root System

Root systems are dynamic and adaptable organs that play critical roles in plant development. However, how roots grow and accumulate biomass during plant life cycle and in relation to shoot growth phenology remains understudied. A comprehensive time-dependent root morphological analysis integrated with molecular signatures is then required to advance our understanding of root growth and development. Here we studied Brachypodium distachyon rooting process by monitoring root morphology, biomass production, and C/N ratios during developmental stages. To provide insight into gene regulation that accompanies root growth, we generated comprehensive transcript profiles of Brachypodium whole-root system at four developmental stages. Our data analysis revealed that multiple biological processes including trehalose metabolism and various families of transcription factors (TFs) were differentially expressed in root system during plant development. In particular, the AUX/IAA, ERFs, WRKY, NAC, and MADS TF family members were upregulated as plant entered the booting/heading stage, while ARFs and GRFs were downregulated suggesting these TF families as important factors involved in specific phases of rooting, and possibly in regulation of transition to plant reproductive stages. We identified several Brachypodium candidate root biomass-promoting genes and cis-regulatory elements for further functional validations and root growth improvements in grasses.

Identification and Characterization of microRNAs in the Developing Seed of Linseed Flax (Linum usitatissimum L.)

Seed development plays an important role during the life cycle of plants. Linseed flax is an oil crop and the seed is a key organ for fatty acids synthesis and storage. So it is important to understand the molecular mechanism of fatty acid biosynthesis during seed development. In this study, four small RNA libraries from early seeds at 5, 10, 20 and 30 days after flowering (DAF) were constructed and used for high-throughput sequencing to identify microRNAs (miRNAs). A total of 235 miRNAs including 114 known conserved miRNAs and 121 novel miRNAs were identified. The expression patterns of these miRNAs in the four libraries were investigated by bioinformatics and quantitative real-time polymerase chain reaction (qPCR) analysis. It was found that several miRNAs, including Lus-miRNA156a was significantly correlated with seed development process. In order to confirm the actual biological function of Lus-miRNA156a, over-expression vector was constructed and transformed to Arabidopsis. The phenotypes of homozygous transgenic lines showed decreasing of oil content and most of the fatty acid content in seeds as well as late flowering time. The results provided a clue that miRNA156a participating the fatty acid biosynthesis pathway and the detailed molecular mechanism of how it regulates the pathway needs to be further investigated.

Two promoter regions confer heat-induced activation of SlDREBA4 in Solanum lycopersicum

Dehydration-responsive element binding (DREB) transcription factors activate the expression of downstream functional genes in combination with a dehydration-responsive element (DRE), and thereby improve the resistance of plants to abiotic stresses such as heat. However, the upstream regulatory mechanism of DREB genes under heat is unclear. A DREBA4 subfamily transcription factor (SlDREBA4), which is heat-responsive and improves heat resistance, was isolated from Solanum lycopersicum 'Microtom'. In this study, promoter truncation experiments were performed to verify changes in β-glucuronidase (GUS) enzyme activity and GUS gene expression levels in transgenic plants with different lengths of promoter fragments under heat and to identify specific regions in the promoter that respond to heat. Our results showed that the GUS reporter gene was constitutively expressed in tissues of the full-length promoter transgenic 'Microtom' plants, with higher expression in conducting tissues of root, stem, and leaf, as well as sepals of flowers and fruits. Under heat treatment, GUS enzyme activity and GUS gene expression levels in tissues of the full-length promoter transgenic plants increased. Promoter deletion analysis identified two positive regulatory regions (-1095 to -730 bp and -162 to -38 bp) responsible for the promoter's response to heat. These results indicated that the heat shock element (HSE) and MYC recognition sequences may cooperate in heat-induced activation of SlDREBA4 promoter.

Radial or Bilateral? The Molecular Basis of Floral Symmetry

In the plant kingdom, the flower is one of the most relevant evolutionary novelties. Floral symmetry has evolved multiple times from the ancestral condition of radial to bilateral symmetry. During evolution, several transcription factors have been recruited by the different developmental pathways in relation to the increase of plant complexity. The MYB proteins are among the most ancient plant transcription factor families and are implicated in different metabolic and developmental processes. In the model plant Antirrhinum majus, three MYB transcription factors (DIVARICATA, DRIF, and RADIALIS) have a pivotal function in the establishment of floral dorsoventral asymmetry. Here, we present an updated report of the role of the DIV, DRIF, and RAD transcription factors in both eudicots and monocots, pointing out their functional changes during plant evolution. In addition, we discuss the molecular models of the establishment of flower symmetry in different flowering plants.

Overexpressing the NAC transcription factor LpNAC13 from Lilium pumilum in tobacco negatively regulates the drought response and positively regulates the salt response

NACs are one of the largest transcription factor families in plants and are involved in the response to abiotic stress. A new stress-responsive NAC transcription factor gene, LpNAC13, was isolated from Lilium pumilum bulbs. The expression of LpNAC13 was induced by drought, salt, cold and ABA treatments. LpNAC13 overexpressing plants were generated to explore the function of LpNAC13 in response to drought and salt stress. Overexpression of LpNAC13 in tobacco displayed a reduced drought tolerance but exhibited an enhanced salt tolerance. The LpNAC13 overexpression plants had decreased antioxidant enzyme activities, content of proline and chlorophyll, increased MDA content under drought condition, the results in the LpNAC13 plants under salt condition were opposite to those under drought condition. The seed germination and root length assays of overexpression of LpNAC13 showed decreased sensitivity to ABA. Functional analyses demonstrate that LpNAC13 plays opposite roles in drought and salt stress tolerance, acting as a negative regulator of drought response but as a positive regulator of salt response in tobacco.

Overexpression of ethylene response factor ERF96 gene enhances selenium tolerance in Arabidopsis

Ethylene response factors (ERFs) are involved in the regulation of plant responses to biotic and abiotic stresses. Here we provide evidence for a role of ERF96, a member of the ERF transcription factor group IX, in selenite tolerance in Arabidopsis. ERF96 gene was rapidly up-regulated in response to selenite stress. Overexpression of ERF96 enhanced Arabidopsis resistance to selenite stress, while ERF96-silenced plants demonstrated wild-type (WT) resistance to selenite. In addition, the overexpression plants had significantly lower selenium (Se) content in shoots when subjected to selenite stress. Further investigation indicated that overexpression of ERF96 reduced transcript levels of selenite/phosphate transporters PHT1;1 and PHT2;1, which influenced Arabidopsis Se uptake and allocation in the presence of selenite. Moreover, our experiments showed that overexpression of ERF96 enhanced Arabidopsis antioxidant activity. Under selenite stress, ERF96-overexpressing lines exhibited the significant increases in catalase (CAT) and glutathione peroxidase (GPX) activities as well as the glutathione (GSH) content, while had a decrease in reactive oxygen species (ROS) accumulation compared to WT. Taken together, our results demonstrate that ERF96 plays a positive role in the regulation of selenite tolerance in Arabidopsis.

The bHLH gene family and its response to saline stress in Jilin ginseng, Panax ginseng C.A. Meyer

Basic helix-loop-helix (bHLH) gene family is a gene family of transcription factors that plays essential roles in plant growth and development, secondary metabolism and response to biotic and abiotic stresses. Therefore, a comprehensive knowledge of the bHLH gene family is paramount to understand the molecular mechanisms underlying these processes and develop advanced technologies to manipulate the processes efficiently. Ginseng, Panax ginseng C.A. Meyer, is a well-known medicinal herb; however, little is known  about the bHLH genes (PgbHLH) in the species. Here, we identified 137 PgbHLH genes from Jilin ginseng cultivar, Damaya, widely cultivated in Jilin, China, of which 50 are newly identified by pan-genome analysis. These 137 PgbHLH genes were phylogenetically classified into 26 subfamilies, suggesting their sequence diversification. They are alternatively spliced into 366 transcripts in a 4-year-old plant and involved in 11 functional subcategories of the gene ontology, indicating their functional differentiation in ginseng. The expressions of the PgbHLH genes dramatically vary spatio-temporally and across 42 genotypes, but they are still somehow functionally correlated. Moreover, the PgbHLH gene family, at least some of its genes, is shown to have roles in plant response to the abiotic stress of saline. These results provide a new insight into the evolution and functional differentiation of the bHLH gene family in plants, new bHLH genes to the PgbHLH gene family, and saline stress-responsive genes for genetic improvement in ginseng and other plant species.

APETALA 2 transcription factor CBX1 is a regulator of mycorrhizal symbiosis and growth of Lotus japonicus

An AP2 family gene CBX1 is involved in mycorrhizal symbiosis and growth of Lotus japonicus. APETALA 2 (AP2) transcriptional regulator is highly conserved in plants. CBX1 from Lotus japonicus is a member of AP2 family. AMF (Arbuscular mycorrhizal fungi) inoculation experiment demonstrated that expression of CBX1 was significantly induced by AMF. Further promoter analysis showed that the - 764 to - 498 bp region of the CBX1 promoter containing CTTC motif is the AMF responsive region. Functional analysis of cbx1 mutant suggested CBX1 is critical for mycorrhizal symbiosis, especially for arbuscule formation. Moreover, under noncolonized condition, overexpression of CBX1 reduced the root length of L. japonicus but increased the size of root system and shoot length, whereas cbx1 mutant reduced the root size and shoot length, but not effect on root length. In addition, cbx1 altered activity of monolignol biosynthetic gene and increased lignin levels. Collectively, these data indicated that CBX1 is a positive regulator of symbiotic activity and plays roles in the growth of L. japonicus.

Expression of a NGATHA1 Gene from Medicago truncatula Delays Flowering Time and Enhances Stress Tolerance

Shoot branching is one of the most variable determinants of crop yield, and the signaling pathways of plant branches have become a hot research topic. As an important transcription factor in the B3 family, NGATHA1 (NGA1), plays an important role in regulating plant lateral organ development and hormone synthesis and transport, but few studies of the role of this gene in the regulation of plant growth and stress tolerance have been reported. In this study, the NGA1 gene was isolated from Medicago truncatula (Mt) and its function was characterized. The cis-acting elements upstream of the 5′ end of MtNGA1 and the expression pattern of MtNGA1 were analyzed, and the results indicated that the gene may act as a regulator of stress resistance. A plant expression vector was constructed and transgenic Arabidopsis plants were obtained. Transgenic Arabidopsis showed delayed flowering time and reduced branching phenotypes. Genes involved in the regulation of branching and flowering were differentially expressed in transgenic plants compared with wild-type plants. Furthermore, transgenic plants demonstrated strong tolerances to salt- and mannitol-induced stresses, which may be due to the upregulated expression of NCED3 (NINE-CIS-EPOXYCAROTENOID DIOXYGENASE 3) by the MtNGA1 gene. These results provide useful information for the exploration and genetic modification use of MtNGA1 in the future.

Genome-Wide Identification and Comparative Analysis of MYB Transcription Factor Family in Musa acuminata and Musa balbisiana

MYB transcription factors (TFs) make up one of the most important TF families in plants. These proteins play crucial roles in processes related to development, metabolism, and stimulus-response; however, very few studies have been reported for the characterization of MYB TFs from banana. The current study identified 305 and 251 MYB genes from Musa acuminata and Musa balbisiana, respectively. Comprehensive details of MYBs are reported in terms of gene structure, protein domain, chromosomal localization, phylogeny, and expression patterns. Based on the exon-intron arrangement, these genes were classified into 12 gene models. Phylogenetic analysis of MYBs involving both species of banana, Oryza sativa, and Arabidopsis thaliana distributed these genes into 27 subfamilies. This highlighted not only the conservation, but also the gain/loss of MYBs in banana. Such genes are important candidates for future functional investigations. The MYB genes in both species exhibited a random distribution on chromosomes with variable densities. Estimation of gene duplication events revealed that segmental duplications represented the major factor behind MYB gene family expansion in banana. Expression profiles of MYB genes were also explored for their potential involvement in acetylene response or development. Collectively, the current comprehensive analysis of MYB genes in both species of banana will facilitate future functional studies.

Transcriptomic and metabolomic analysis provides insights into anthocyanin and procyanidin accumulation in pear

BACKGROUND

Pear is one of the most important fruit crops worldwide. Anthocyanins and procyanidins (PAs) are important secondary metabolites that affect the appearance and nutritive quality of pear. However, few studies have focused on the molecular mechanism underlying anthocyanin and PA accumulation in pear.

RESULTS

We conducted metabolome and transcriptome analyses to identify candidate genes involved in anthocyanin and PA accumulation in young fruits of the pear cultivar 'Clapp Favorite' (CF) and its red mutation cultivar 'Red Clapp Favorite' (RCF). Gene-metabolite correlation analyses revealed a 'core set' of 20 genes that were strongly correlated with 10 anthocyanin and seven PA metabolites. Of these, PcGSTF12 was confirmed to be involved in anthocyanin and PA accumulation by complementation of the tt19-7 Arabidopsis mutant. Interestingly, PcGSTF12 was found to be responsible for the accumulation of procyanidin A3, but not petunidin 3, 5-diglucoside, opposite to the function of AtGSTs in Arabidopsis. Transformation with PcGSTF12 greatly promoted or repressed genes involved in anthocyanin and PA biosynthesis, regulation, and transport. Electrophoretic mobility shift and luciferase reporter assays confirmed positive regulation of PcGSTF12 by PcMYB114.

CONCLUSION

These findings identify a core set of genes for anthocyanin and PA accumulation in pear. Of these, PcGSTF12, was confirmed to be involved in anthocyanin and PA accumulation. Our results also identified an important anthocyanin and PA regulation node comprising two core genes, PcGSTF12 and PcMYB114. These results provide novel insights into anthocyanin and PA accumulation in pear and represent a valuable data set to guide future functional studies and pear breeding.

Genome-Wide Identification, Classification, and Expression Profiling Reveals R2R3-MYB Transcription Factors Related to Monoterpenoid Biosynthesis in Osmanthus fragrans

Osmanthus fragrans is widely grown for the purpose of urban greening and the pleasant aroma emitted from its flowers. The floral scent is determined by several monoterpenoid volatiles, such as linalool and its oxides, which are a few of the most common volatiles and the main components of the essential oils in most sweet osmanthus cultivars. In addition, the relative contents of cis- and trans-linalool oxide (furan) may affect the aromas and quality of the essential oils. MYB proteins represent the largest family of transcription factors in plants and participate in regulating secondary metabolites. Several cis-elements, especially AC-rich regions, are known to be bound by 2R-MYBs and could be found in the promoter of the enzyme genes in the terpenoid metabolic pathway. However, there has to date been no investigation into the 2R-MYB family genes involved in regulating terpenoid biosynthesis in O. fragrans. Here, 243 non-redundant 2R-MYB proteins were grouped into 33 clusters based on the phylogeny and exon-intron distribution. These genes were unevenly distributed on 23 chromosomes. Ka/Ks analysis showed that the major mode of 2R-MYB gene evolution was purifying selection. Expression analysis indicated that 2R-MYB genes in O. fragrans exhibited varied expression patterns. A total of 35 OfMYBs representing the highest per kilobase per million mapped reads in the flower were selected for quantitative real-time PCR analysis. The correlation analysis between the expression level and the contents of fragrant compounds at different flowering stages suggested that OfMYB19/20 exhibited remarkably positive correlation with the accumulation of cis-linalool oxides. OfMYB51/65/88/121/137/144 showed significantly negative correlations with one or more linalool oxides. Characterization of these proteins revealed that OfMYB19 and OfMYB137 were localized in the nuclei, but did not show transcriptional activation in the yeast system, which suggested that they may be bound to other transcription factors to exert regulatory functions. These findings provide useful information for further functional investigation of the 2R-MYBs and offer a foundation for clarifying the 2R-MYB transcription factors involved in the molecular mechanism of the regulation of monoterpenoid biosynthesis in Osmanthus fragrans.

Transcriptome analysis uncovers the gene expression profile of salt-stressed potato (Solanum tuberosum L.)

Potato (Solanum tuberosum L.) is an important staple food worldwide. However, its growth has been heavily suppressed by salt stress. The molecular mechanisms of salt tolerance in potato remain unclear. It has been shown that the tetraploid potato Longshu No. 5 is a salt-tolerant genotype. Therefore, in this study we conducted research to identify salt stress response genes in Longshu No. 5 using a NaCl treatment and time-course RNA sequencing. The total number of differentially expressed genes (DEGs) in response to salt stress was 5508. Based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, it was found that DEGs were significantly enriched in the categories of nucleic acid binding, transporter activity, ion or molecule transport, ion binding, kinase activity and oxidative phosphorylation. Particularly, the significant differential expression of encoding ion transport signaling genes suggests that this signaling pathway plays a vital role in salt stress response in potato. Finally, the DEGs in the salt response pathway were verified by Quantitative real-time PCR (qRT-PCR). These results provide valuable information on the salt tolerance of molecular mechanisms in potatoes, and establish a basis for breeding salt-tolerant cultivars.

Ectopic Expression of AhGLK1b (GOLDEN2-like Transcription Factor) in Arabidopsis Confers Dual Resistance to Fungal and Bacterial Pathogens

GOLDEN2-LIKE (GLK) is a member of the myeloblastosis (MYB) family transcription factor and it plays an important role in the regulation of plastid development and stress tolerance. In this study, a gene named AhGLK1b was identified from a cultivated peanut showing down-regulation in response to low calcium with a complete open reading frame (ORF) of 1212 bp. The AhGLK1b has 99.26% and 96.28% sequence similarities with its orthologs in Arachis ipaensis and A. duranensis, respectively. In the peanut, the AhGLK1b was localized in the nucleus and demonstrated the highest expression in the leaf, followed by the embryo. Furthermore, the expression of AhGLK1b was induced significantly in response to a bacterial pathogen, Ralstonia solanacearum infection. Ectopic expression of AhGLK1b in Arabidopsis showed stronger resistance against important phytopathogenic fungi S. sclerotiorum. It also exhibited high resistance to infection of the bacterial pathogen Pst DC3000. AhGLK1b-expressing Arabidopsis induced defense-related genes including PR10 and Phox/Bem 1 (PBI), which are involved in multiple disease resistance. Taken together, the results suggest that AhGLK1b might be useful in providing dual resistance to fungal and bacterial pathogens as well as tolerance to abiotic stresses.

MYB Transcription Factors as Regulators of Secondary Metabolism in Plants

MYB transcription factors (TFs), as one of the largest gene families in plants, play important roles in multiple biological processes, such as plant growth and development, cell morphology and pattern building, physiological activity metabolism, primary and secondary metabolic reactions, and responses to environmental stresses. The function of MYB TFs in crops has been widely studied, but few studies have been done on medicinal plants. In this review, we summarized the MYB TFs that play important roles in secondary metabolism and emphasized the possible mechanisms underlying how MYB TFs are regulated at the protein, posttranscriptional, and transcriptional levels, as well as how they regulate the downstream target gene networks related to secondary metabolism in plants, especially in medicinal plants.

Specific functions for Mediator complex subunits from different modules in the transcriptional response of Arabidopsis thaliana to abiotic stress

Adverse environmental conditions are detrimental to plant growth and development. Acclimation to abiotic stress conditions involves activation of signaling pathways which often results in changes in gene expression via networks of transcription factors (TFs). Mediator is a highly conserved co-regulator complex and an essential component of the transcriptional machinery in eukaryotes. Some Mediator subunits have been implicated in stress-responsive signaling pathways; however, much remains unknown regarding the role of plant Mediator in abiotic stress responses. Here, we use RNA-seq to analyze the transcriptional response of Arabidopsis thaliana to heat, cold and salt stress conditions. We identify a set of common abiotic stress regulons and describe the sequential and combinatorial nature of TFs involved in their transcriptional regulation. Furthermore, we identify stress-specific roles for the Mediator subunits MED9, MED16, MED18 and CDK8, and putative TFs connecting them to different stress signaling pathways. Our data also indicate different modes of action for subunits or modules of Mediator at the same gene loci, including a co-repressor function for MED16 prior to stress. These results illuminate a poorly understood but important player in the transcriptional response of plants to abiotic stress and identify target genes and mechanisms as a prelude to further biochemical characterization.

NAC Family Transcription Factors in Carrot: Genomic and Transcriptomic Analysis and Responses to Abiotic Stresses

Carrot is an annual or biennial herbaceous plant of the Apiaceae family. Carrot is an important vegetable, and its fresh taproot, which contains rich nutrients, is the main edible part. In the life cycle of carrot, NAC family transcription factors (TFs) are involved in almost all physiological processes. The function of NAC TFs in carrot remains unclear. In this study, 73 NAC family TF members in carrot were identified and characterized using transcriptome and genome databases. These members were divided into 14 subfamilies. Multiple sequence alignment was performed, and the conserved domains, common motifs, phylogenetic tree, and interaction network of DcNAC proteins were predicted and analyzed. Results showed that the same group of NAC proteins of carrot had high similarity. Eight DcNAC genes were selected to detect their expression profiles under abiotic stress treatments. The expression levels of the selected DcNAC genes significantly increased under treatments with low temperature, high temperature, drought, and salt stress. Results provide potentially useful information for further analysis of the roles of DcNAC transcription factors in carrot.

Structural insights into target DNA recognition by R2R3-MYB transcription factors

As the largest group of MYB family transcription factors, R2R3-MYB proteins play essential roles during plant growth and development. However, the structural basis underlying how R2R3-MYBs recognize the target DNA remains elusive. Here, we report the crystal structure of Arabidopsis WEREWOLF (WER), an R2R3-MYB protein, in complex with its target DNA. Structural analysis showed that the third α-helices in both the R2 and R3 repeats of WER fit in the major groove of the DNA, specifically recognizing the DNA motif 5'-AACNGC-3'. In combination with mutagenesis, in vitro binding and in vivo luciferase assays, we showed that K55, N106, K109 and N110 are critical for the function of WER. Although L59 of WER is not involved in DNA binding in the structure, ITC analysis suggested that L59 plays an important role in sensing DNA methylation at the fifth position of cytosine (5mC). Like 5mC, methylation at the sixth position of adenine (6mA) in the AAC element also inhibits the interaction between WER and its target DNA. Our study not only unravels the molecular basis of how WER recognizes its target DNA, but also suggests that 5mC and 6mA modifications may block the interaction between R2R3-MYB transcription factors and their target genes.

Genome-wide survey of the bHLH super gene family in Brassica napus

BACKGROUND

The basic helix-loop-helix (bHLH) gene family is one of the largest transcription factor families in plants and is functionally characterized in diverse species. However, less is known about its functions in the economically important allopolyploid oil crop, Brassica napus.

RESULTS

We identified 602 potential bHLHs in the B. napus genome (BnabHLHs) and categorized them into 35 subfamilies, including seven newly separated subfamilies, based on phylogeny, protein structure, and exon-intron organization analysis. The intron insertion patterns of this gene family were analyzed and a total of eight types were identified in the bHLH regions of BnabHLHs. Chromosome distribution and synteny analyses revealed that hybridization between Brassica rapa and Brassica oleracea was the main expansion mechanism for BnabHLHs. Expression analyses showed that BnabHLHs were widely in different plant tissues and formed seven main patterns, suggesting they may participate in various aspects of B. napus development. Furthermore, when roots were treated with five different hormones (IAA, auxin; GA3, gibberellin; 6-BA, cytokinin; ABA, abscisic acid and ACC, ethylene), the expression profiles of BnabHLHs changed significantly, with many showing increased expression. The induction of five candidate BnabHLHs was confirmed following the five hormone treatments via qRT-PCR. Up to 246 BnabHLHs from nine subfamilies were predicted to have potential roles relating to root development through the joint analysis of their expression profiles and homolog function.

CONCLUSION

The 602 BnabHLHs identified from B. napus were classified into 35 subfamilies, and those members from the same subfamily generally had similar sequence motifs. Overall, we found that BnabHLHs may be widely involved in root development in B. napus. Moreover, this study provides important insights into the potential functions of the BnabHLHs super gene family and thus will be useful in future gene function research.

Phenotypic and Transcriptomic Analysis of Two Pinellia ternata Varieties T2 line and T2Plus line

Pinellia (Pinellia ternata (Thunb.) Breit.), as important medicinal plant, has been used to treat various ailments for a long time. The sixteen ploid plant (2n = 16 * 13 = 208) Pinellia T2Plus line was obtained from an octoploid (2n = 8 * 13 = 104) T2 line by chromosome-doubling technique. Compared with T2 line, the content of various medicinal components (polysaccharide, guanosine, adenosine and ephedrine) was increased in T2Plus line. In this study, the transcriptome of T2 line and T2Plus line were characterized by RNA sequencing (RNA-seq) technology. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enrichment analysis on differential expressed unigenes (DEGs) revealed that multiple metabolic pathway were enriched significantly, such as 'Starch and sucrose metabolism', 'Purine metabolism', 'Photosynthesis' and six transcription factors (MYB, WRKY, bHLH, lateral organ boundaries domain (LBD), homeodomain-zipper (HD-ZIP) and Ethylene-responsive factor (ERF)) play a key role in difference of transcriptome between T2 line and T2Plus line. These metabolic pathways and transcription factors may play an important role in the difference of medicinal components and epigenetic features between these two Pinellia cultivars. This conclusion provides a robust theoretical basis for the mechanism of the formation of medicinal ingredients in Pinellia cultivars.

Phylogenetic analyses of key developmental genes provide insight into the complex evolution of seeds

Gene duplication plays a decisive role in organismal diversification and in the appearance of novel structures. In plants the megagametophyte covered by the integuments, which after fertilization becomes the seed constitutes a novel structure: the ovule. In Arabidopsis thaliana, genetic mechanisms regulating ovule development, including the genetics underlying ovule initiation, ovule patterning and integument development, have been identified. Among seed plants, integuments are not only a novelty in evolution, but integuments also present an enormous morphological variation. This study is focused on the evolution of gene families that play a role in the proper morphological development of the integuments, BELL1 (BEL1), KANADIs (KAN1, KAN2, and KAN4/ATS), UNICORN (UCN) and SHORT INTEGUMENTS1 (SIN1). In Arabidopsis, BEL1 establishes the initiation of integument development. KAN1 and 2 act in the proper development of the outer integument. While ABERRANT TESTA SHAPE (ATS), is involved in the correct separation of both integuments. UCN acts in planar growth of the outer integument repressing ATS. SIN1 is involved in cell elongation in the integuments. The results of our analyses show that each of these genes has a different evolutionary history and that while gymnosperms appear to have a simpler ovule morphology, they have more homologues of these candidate genes than angiosperms. In addition, we present the conserved and novel motifs for each of these genes among seed plants and their selection constraints, which may be related to functional changes and to the diversity of ovule morphologies.

Pigment variation and transcriptional response of the pigment synthesis pathway in the S2309 triple-color ornamental kale (Brassica oleracea L. var. acephala) line

Ornamental kale is popular because of its colorful leaves and few studies have investigated the mechanism of color changes. In this study, an ornamental kale line (S2309) with three leaf colors was developed. Analysis of the anthocyanin, chlorophyll, and carotenoid contents and RNA-seq were performed on the three leaf color types. There was less chlorophyll in the white leaves and purple leaves than in the green leaves, and the anthocyanin content was greatest in the purple leaves. All the downregulated DEGs related to chlorophyll metabolism were detected only in the S2309_G vs. S2309_W comparison, which indicated that the decrease in chlorophyll content was caused mainly by the inhibition of chlorophyll biosynthesis during the leaf color change from green to white. Moreover, the expression of 19 DEGs involved in the anthocyanin biosynthesis pathway was upregulated. These results provide new insight into the mechanisms underlying the three-color formation.

The small secreted effector protein MiSSP7.6 of Laccaria bicolor is required for the establishment of ectomycorrhizal symbiosis

To establish and maintain a symbiotic relationship, the ectomycorrhizal fungus Laccaria bicolor releases mycorrhiza-induced small secreted proteins (MiSSPs) into host roots. Here, we have functionally characterized the MYCORRHIZA-iNDUCED SMALL SECRETED PROTEIN OF 7.6 kDa (MiSSP7.6) from L. bicolor by assessing its induced expression in ectomycorrhizae, silencing its expression by RNAi, and tracking in planta subcellular localization of its protein product. We also carried out yeast two-hybrid assays and bimolecular fluorescence complementation analysis to identify possible protein targets of the MiSSP7.6 effector in Populus roots. We showed that MiSSP7.6 expression is upregulated in ectomycorrhizal rootlets and associated extramatrical mycelium during the late stage of symbiosis development. RNAi mutants with a decreased MiSSP7.6 expression have a lower mycorrhization rate, suggesting a key role in the establishment of the symbiosis with plants. MiSSP7.6 is secreted, and it localizes both to the nuclei and cytoplasm in plant cells. MiSSP7.6 protein was shown to interact with two Populus Trihelix transcription factors. Furthermore, when coexpressed with one of the Trihelix transcription factors, MiSSP7.6 is localized to plant nuclei only. Our data suggest that MiSSP7.6 is a novel secreted symbiotic effector and is a potential determinant for ectomycorrhiza formation.

A basic helix-loop-helix protein (GhFP1) promotes fibre elongation of cotton (Gossypium hirsutum) by modulating brassinosteroid biosynthesis and signalling

Basic helix-loop-helix (bHLH) proteins are involved in transcriptional networks controlling a number of biological processes in plants. However, little information is known on the roles of bHLH proteins in cotton fibre development so far. Here, we show that a cotton bHLH protein (GhFP1) positively regulates fibre elongation. GhFP1 transgenic cotton and Arabidopsis plants were generated to study how GhFP1 regulates fibre cell elongation. Fibre length of the transgenic cotton overexpressing GhFP1 was significantly longer than that of wild-type, whereas suppression of GhFP1 expression hindered fibre elongation. Furthermore, overexpression of GhFP1 in Arabidopsis promoted trichome development. Expression of the brassinosteroid (BR)-related genes was markedly upregulated in fibres of GhFP1 overexpression cotton, but downregulated in GhFP1-silenced fibres. BR content in the transgenic fibres was significantly altered, relative to that in wild-type. Moreover, GhFP1 protein could directly bind to the promoters of GhDWF4 and GhCPD to activate expression of these BR-related genes. Therefore, our data suggest that GhFP1 as a positive regulator participates in controlling fibre elongation by activating BR biosynthesis and signalling. Additionally, homodimerisation of GhFP1 may be essential for its function, and interaction between GhFP1 and other cotton bHLH proteins may interfere with its DNA-binding activity.

In silico characterization, and expression analysis of rice golden 2-like (OsGLK) members in response to low phosphorous

The availability of phosphorus (P) affects productivity of rice. Under acidic soil conditions (pH < 5.5), P is rapidly immobilized in the soil. Several transcription factors play an important role in low Pi tolerance response, including MYB family members but their role in acidic soil is yet unknown. In this study, genome wide identification and characterization of golden 2-like (GLK) members belonging to GARP superfamily from rice (OsGLK) led to identification of 46 members distributed over 12 chromosomes. We assigned gene nomenclature, analyzed gene structure and identified mutant orthologs and phenotypes in maize and rice, respectively. On the basis of biological functions three categories viz., (a) two-component response regulator (five members), (b) putative transcription factor (21 members) and (c) phosphate starvation response (8 members) were identified. Phylogenetic analysis revealed a total of nine subgroups with MYB homeodomain-like and MYB CC-type domains conserved across members. Expression profiling of OsGLKs in response to 24 and 48 h of low Pi in four contrasting rice genotypes, revealed significantly higher expression of OsGLK10, OsGLK15, OsGLK22 and OsGLK30 in tolerant genotypes as compared to susceptible genotypes, suggesting their role in Pi starvation tolerance. Meta analyses and cis-regulatory elements (CREs) profiling of OsGLK showed diverse expression pattern in various tissues and organs and also modulation in response to various abiotic and biotic stresses. Our results highlight the versatile role of this diverse and complex GLK family, in particular to abiotic stress. These genes will form the basis of future studies on low Pi tolerance in acidic soils.

Genome-Wide Identification, Structure Characterization, and Expression Profiling of Dof Transcription Factor Gene Family in Wheat (Triticum aestivum L.)

DNA binding with one finger (Dof) proteins are plant-specific transcription factors with crucial roles in plant growth and stress response. Even so, little is known about them in wheat. In this study, 108 wheat Dof (TaDof) genes across 21 chromosomes were detected. Although variable in sequence length, molecular weight, and isoelectric point, all TaDof proteins contained conserved zinc-finger structures and were phylogenetically divided into 7 sub-groups. Exon/intron and motif analyses suggested that TaDof structures and conserved motifs were similar within sub-groups but diverse among sub-groups. Many segmental duplications were identified and Ka/Ks and inter-species synthetic analyses indicated that polyploidization was main reason for increased number of TaDofs. Prediction and experimental confirmation revealed that TaDofs functioned as transcription factors in the nucleus. Expression pattern profiling showed that TaDofs specifically affected growth and development, and biotic and abiotic stress responses. Wheat miRNAs and cis-regulator were predicted as essential players in molding TaDofs expression patterns. qRT-PCR analysis revealed that TaDofs were induced by salt and drought stresses. Customized annotation revealed that TaDofs were widely involved in phytohormone response, defense, growth and development, and metabolism. Our study provided a comprehensive understanding to wheat TaDofs.

The R2R3-MYB transcription factor family in Taxus chinensis: identification, characterization, expression profiling and posttranscriptional regulation analysis

The MYB transcription factor family is one of the largest gene families playing regulatory roles in plant growth and development. The MYB family has been studied in a variety of plant species but has not been reported in Taxus chinensis. Here we identified 72 putative R2R3-MYB genes in T. chinensis using a comprehensive analysis. Sequence features, conversed domains and motifs were characterized. The phylogenetic analysis showed TcMYBs and AtMYBs were clustered into 36 subgroups, of which 24 subgroups included members from T. chinensis and Arabidopsis thaliana, while 12 subgroups were specific to one species. This suggests the conservation and specificity in structure and function of plant R2R3-MYBs. The expression of TcMYBs in various tissues and different ages of xylem were investigated. Additionally, miRNA-mediated posttranscriptional regulation analysis revealed that TcMYBs were the targets of miR858, miR159 and miR828, suggesting the posttranscriptional regulation of MYBs is highly conserved in plants. The results provide a basis for further study the role of TcMYBs in the regulation of secondary metabolites of T. chinensis.

Genome-wide identification and analysis of B-BOX gene family in grapevine reveal its potential functions in berry development

BACKGROUND

The B-BOX (BBX) proteins are the class of zinc-finger transcription factors and can regulate plant growth, development, and endure stress response. In plants, the BBX gene family has been identified in Arabidopsis, rice, and tomato. However, no systematic analysis of BBX genes has been undertaken in grapevine.

RESULTS

In this study, 24 grapevine BBX (VvBBX) genes were identified by comprehensive bioinformatics analysis. Subsequently, the chromosomal localizations, gene structure, conserved domains, phylogenetic relationship, gene duplication, and cis-acting elements were analyzed. Phylogenetic analysis divided VvBBX genes into five subgroups. Numerous cis-acting elements related to plant development, hormone and/or stress responses were identified in the promoter of the VvBBX genes. The tissue-specific expressional dynamics of VvBBX genes demonstrated that VvBBXs might play important role in plant growth and development. The transcript analysis from transcriptome data and qRT-PCR inferred that 11 VvBBX genes were down-regulated in different fruit developmental stages, while three VvBBX genes were up-regulated. It is also speculated that VvBBX genes might be involved in multiple hormone signaling (ABA, ethylene, GA3, and CPPU) as transcriptional regulators to modulate berry development and ripening. VvBBX22 seems to be responsive to multiple hormone signaling, including ABA, ethylene GA3, and CPPU. Some VvBBX genes were strongly induced by Cu, salt, waterlogging, and drought stress treatment. Furthermore, the expression of VvBBX22 proposed its involvement in multiple functions, including leaf senescence, abiotic stress responses, fruit development, and hormone response.

CONCLUSIONS

Our results will provide the reference for functional studies of BBX gene family, and highlight its functions in grapevine berry development and ripening. The results will help us to better understand the complexity of the BBX gene family in abiotic stress tolerance and provide valuable information for future functional characterization of specific genes in grapevine.

QTL-Seq and Sequence Assembly Rapidly Mapped the Gene BrMYBL2.1 for the Purple Trait in Brassica rapa

Anthocyanins have strong antioxidant activity and are believed to be healthy for human beings. The Brassica rapa L. ssp. chinensis var. purpurea “Zicaitai” is rich in anthocyanins. We constructed an F₂ population of Zicaitai and “Caixin” (Brassica rapa ssp. parachinensis) and it shows clear segregation of the purple phenotype (i.e., variation in anthocyanin enrichment). Here, quantitative trait locus (QTL)-Seq was performed with two sample groups from the F₂ population: one exhibiting an intense purple phenotype and the other showed a completely green phenotype. The results showed that the QTL-Seq and linkage analysis located different major loci. This indicates that there are two major genetic factors that plays different roles in regulating anthocyanin enrichment in Zicaitai. This was further supported by the data simulation of an in silico F₂ population that QTL-Seq and linkage analysis can locate different major loci. Furthermore, the draft genomes of the two parents (Zicaitai and Caixin) were assembled and utilized to search for mutations in candidate genes. A ~100-bp insertion was found in the third exon of gene BrMYBL2.1 in Zicaitai. BrMYBL2.1 is a negative regulator of anthocyanin biosynthesis, while BrEGL3.2—previously located by linkage mapping—is a positive regulator. For these populations with multiple genes contributing large effects to a trait, a strategy of low depth re-sequencing of F₂ individuals followed by QTL-Seq analysis with the free combination of sample groups is proposed. Furthermore, draft-sequence assembly of parental genomes together with QTL mapping is suggested as an efficient means for fine-mapping genes rapidly in segregating populations.

Salinity tolerance in barley during germination- homologs and potential genes

Salinity affects more than 6% of the world's total land area, causing massive losses in crop yield. Salinity inhibits plant growth and development through osmotic and ionic stresses; however, some plants exhibit adaptations through osmotic regulation, exclusion, and translocation of accumulated Na+ or Cl-. Currently, there are no practical, economically viable methods for managing salinity, so the best practice is to grow crops with improved tolerance. Germination is the stage in a plant's life cycle most adversely affected by salinity. Barley, the fourth most important cereal crop in the world, has outstanding salinity tolerance, relative to other cereal crops. Here, we review the genetics of salinity tolerance in barley during germination by summarizing reported quantitative trait loci (QTLs) and functional genes. The homologs of candidate genes for salinity tolerance in Arabidopsis, soybean, maize, wheat, and rice have been blasted and mapped on the barley reference genome. The genetic diversity of three reported functional gene families for salt tolerance during barley germination, namely dehydration-responsive element-binding (DREB) protein, somatic embryogenesis receptor-like kinase and aquaporin genes, is discussed. While all three gene families show great diversity in most plant species, the DREB gene family is more diverse in barley than in wheat and rice. Further to this review, a convenient method for screening for salinity tolerance at germination is needed, and the mechanisms of action of the genes involved in salt tolerance need to be identified, validated, and transferred to commercial cultivars for field production in saline soil.

Genome-Wide Analysis of the MYB Transcription Factor Superfamily in Physcomitrella patens

MYB transcription factors (TFs) are one of the largest TF families in plants to regulate numerous biological processes. However, our knowledge of the MYB family in Physcomitrella patens is limited. We identified 116 MYB genes in the P. patens genome, which were classified into the R2R3-MYB, R1R2R3-MYB, 4R-MYB, and MYB-related subfamilies. Most R2R3 genes contain 3 exons and 2 introns, whereas R1R2R3 MYB genes contain 10 exons and 9 introns. N3R-MYB (novel 3RMYB) and NR-MYBs (novel RMYBs) with complicated gene structures appear to be novel MYB proteins. In addition, we found that the diversity of the MYB domain was mainly contributed by domain shuffling and gene duplication. RNA-seq analysis suggested that MYBs exhibited differential expression to heat and might play important roles in heat stress responses, whereas CCA1-like MYB genes might confer greater flexibility to the circadian clock. Some R2R3-MYB and CCA1-like MYB genes are preferentially expressed in the archegonium and during the transition from the chloronema to caulonema stage, suggesting their roles in development. Compared with that of algae, the numbers of MYBs have significantly increased, thus our study lays the foundation for further exploring the potential roles of MYBs in the transition from aquatic to terrestrial environments.

A population genomics appraisal suggests independent dispersals for bitter and sweet manioc in Brazilian Amazonia

Amazonia is a major world centre of plant domestication, but the genetics of domestication remains unclear for most Amazonian crops. Manioc (Manihot esculenta) is the most important staple food crop that originated in this region. Although manioc is relatively well-studied, little is known about the diversification of bitter and sweet landraces and how they were dispersed across Amazonia. We evaluated single nucleotide polymorphisms (SNPs) in wild and cultivated manioc to identify outlier SNPs putatively under selection and to assess the neutral genetic structure of landraces to make inferences about the evolution of the crop in Amazonia. Some outlier SNPs were in putative manioc genes possibly related to plant architecture, transcriptional regulation and responses to stress. The neutral SNPs revealed contrasting genetic structuring for bitter and sweet landraces. The outlier SNPs may be signatures of the genomic changes resulting from domestication, while the neutral genetic structure suggests independent dispersals for sweet and bitter manioc, possibly related to the earlier domestication and diversification of the former. Our results highlight the role of ancient peoples and current smallholders in the management and conservation of manioc genetic diversity, including putative genes and specific genetic resources with adaptive potential in the context of climate change.