Genome-wide analysis of the distribution of AP2/ERF transcription factors reveals duplication and CBFs genes elucidate their potential function in Brassica oleracea

The cold-responsive C-repeat binding factors in Betula platyphylla Suk. positively regulate cold tolerance

Cold stress is one of the most destructive abiotic stresses limiting plant growth and development. CBF (C-repeat binding factor) transcription factors and their roles in cold response have been identified in Arabidopsis as well as several other plant species. However, the biological functions and related molecular mechanisms of CBFs in birch (Betula platyphylla Suk.) remain undetermined. In this study, five cold-responsive BpCBF genes, BpCBF1, BpCBF2, BpCBF7, BpCBF10 and BpCBF12 were cloned. Via protoplast transformation, BpCBF7 was found to be localized in nucleus. The result of yeast one hybrid assay validated the binding of BpCBF7 to the CRT/DRE (C-repeat/dehydration responsive element) elements in the promoter of BpERF1.1 gene. By overexpressing and repressing BpCBFs in birch plants, it was proven that BpCBFs play positive roles in the cold tolerance. At the metabolic level, BpCBFs OE lines had lower ROS accumulation, as well as higher activities of antioxidant enzymes (SOD, POD and CAT) and higher accumulation of protective substances (soluble sugar, soluble protein and proline). Via yeast one hybrid and co-transformation of effector and reporter vectors assay, it was proven that BpCBF7 can regulate the expression of BpERF5 and BpZAT10 genes by directly binding to their promoters. An interacting protein of BpCBF7, BpWRKY17, was identified by yeast two hybrid library sequencing and the interaction was validated with in vivo methods. These results indicates that BpCBFs can increase the cold tolerance of birch plants, partly by gene regulation and protein interaction. This study provides a reference for the research on CBF transcription factors and genetic improvement of forest trees upon abiotic stresses.

Genome-wide association analysis identifies a candidate gene controlling seed size and yield in Xanthoceras sorbifolium Bunge

Yellow horn (Xanthoceras sorbifolium Bunge) is a woody oilseed tree species whose seed oil is rich in unsaturated fatty acids and rare neuronic acids, and can be used as a high-grade edible oil or as a feedstock for biodiesel production. However, the genetic mechanisms related to seed yield in yellow horn are not well elucidated. This study identified 2 164 863 SNP loci based on 222 genome-wide resequencing data of yellow horn germplasm. We conducted genome-wide association study (GWAS) analysis on three core traits (hundred-grain weight, single-fruit seed mass, and single-fruit seed number) that influence seed yield for the years 2022 and 2020, and identified 399 significant SNP loci. Among these loci, the Chr10_24013014 and Chr10_24012613 loci caught our attention due to their consistent associations across multiple analyses. Through Sanger sequencing, we validated the genotypes of these two loci across 16 germplasms, confirming their consistency with the GWAS analysis results. Downstream of these two significant loci, we identified a candidate gene encoding an AP2 transcription factor protein, which we named XsAP2. RT-qPCR analysis revealed high expression of the XsAP2 gene in seeds, and a significant negative correlation between its expression levels and seed hundred-grain weight, as well as single-fruit seed mass, suggesting its potential role in the normal seed development process. Transgenic Arabidopsis lines with the overexpressed XsAP2 gene exhibited varying degrees of reduction in seed size, number of seeds per silique, and number of siliques per plant compared with wild-type Arabidopsis. Combining these results, we hypothesize that the XsAP2 gene may have a negative regulatory effect on seed yield of yellow horn. These results provide a reference for the molecular breeding of high-yielding yellow horn.

Genome-Wide Identification of the CBF Gene Family and ICE Transcription Factors in Walnuts and Expression Profiles under Cold Conditions

Cold stress impacts woody tree growth and perennial production, especially when the temperature rapidly changes in late spring. To address this issue, we conducted the genome-wide identification of two important transcription factors (TFs), CBF (C-repeat binding factors) and ICE (inducers of CBF expression), in three walnut (Juglans) genomes. Although the CBF and ICE gene families have been identified in many crops, very little systematic analysis of these genes has been carried out in J. regia and J. sigillata. In this study, we identified a total of 16 CBF and 12 ICE genes in three Juglans genomes using bioinformatics analysis. Both CBF and ICE had conserved domains, motifs, and gene structures, which suggests that these two TFs were evolutionarily conserved. Most ICE genes are located at both ends of the chromosomes. The promoter cis-regulatory elements of CBF and ICE genes are largely involved in light and phytohormone responses. Based on 36 RNA sequencing of leaves from four walnut cultivars ('Zijing', 'Lvling', 'Hongren', and 'Liao1') under three temperature conditions (8 °C, 22 °C, and 5 °C) conditions in late spring, we found that the ICE genes were expressed more highly than CBFs. Both CBF and ICE proteins interacted with cold-related proteins, and many putative miRNAs had interactions with these two TFs. These results determined that CBF1 and ICE1 play important roles in the tolerance of walnut leaves to rapid temperature changes. Our results provide a useful resource on the function of the CBF and ICE genes related to cold tolerance in walnuts.

Genome-wide analysis of the DREB family genes and functional identification of the involvement of BrDREB2B in abiotic stress in wucai (Brassica campestris L.)

Dehydration responsive element binding protein (DREB) is a significant transcription factor class known to be implicated in abiotic stresses. In this study, we systematically conducted a genome-wide identification and expression analysis of the DREB gene family, including gene structures, evolutionary relationships, chromosome distribution, conserved domains, and expression patterns. A total of 65 DREB family gene members were identified in Chinese cabbage (Brassica rapa L.) and were classified into five subgroups based on phylogenetic analysis. Through analysis of the conserved domains of BrDREB family genes, only one exon existed in the gene structure. Through the analysis of cis-acting elements, these genes were mainly involved in hormone regulation and adversity stress. In order to identify the function of BrDREB2B, overexpressed transgenic Arabidopsis was constructed. After different stress treatments, the germination rate, root growth, survival rate, and various plant physiological indicators were measured. The results showed that transgenic Arabidopsis thaliana plants overexpressing BrDREB2B exhibited enhanced tolerance to salt, heat and drought stresses. Taken together, our results are the first to report the BrDREB2B gene response to drought and heat stresses in Chinese cabbage and provide a basis for further studies to determine the function of BrDREBs in response to abiotic stresses.

Genome-wide identification of CBF genes and their responses to cold acclimation in Taraxacum kok-saghyz

C-repeat binding factors (CBFs) are transcription factors that are known to play important roles in plant cold acclimation. They are highly conserved in most higher plants. Taraxacum kok-saghyz (TKS) is an herb native to China and Kazakhstan and is well-known for its production of rubber silk with industrial and economic value. To understand cold acclimation mechanisms, we conducted a genome-wide discovery of the CBF family genes in TKS and revealed ten CBF genes. A bioinformatic analysis of the CBF genes was carried out to analyze the phylogenetic relationship, protein conservative motifs, protein physicochemical properties, gene structure, promoter cis-acting elements, and the gene expression patterns under cold acclimation and control conditions. It was found that most of these genes were highly responsive at the late stage of cold acclimation, indicating that they play important roles in the cold acclimation processes of TKS. This study provides a theoretical basis for the study of the molecular functions of the CBF gene family in TKS, and a useful guidance for the genetic improvement of the cold tolerance traits of TKS and other plants, including crops.

Expansion and Molecular Characterization of AP2/ERF Gene Family in Wheat (Triticum aestivum L.)

The AP2/ERF is a large protein family of transcription factors, playing an important role in signal transduction, plant growth, development, and response to various stresses. AP2/ERF super-family is identified and functionalized in a different plant but no comprehensive and systematic analysis in wheat (Triticum aestivum L.) has been reported. However, a genome-wide and functional analysis was performed and identified 322 TaAP2/ERF putative genes from the wheat genome. According to the phylogenetic and structural analysis, TaAP2/ERF genes were divided into 12 subfamilies (Ia, Ib, Ic, IIa, IIb, IIc, IIIa, IIIb, IIIc, IVa, IVb, and IVc). Furthermore, conserved motifs and introns/exons analysis revealed may lead to functional divergence within clades. Cis-Acting analysis indicated that many elements were involved in stress-related and plant development. Chromosomal location showed that 320 AP2/ERF genes were distributed among 21 chromosomes and 2 genes were present in a scaffold. Interspecies microsynteny analysis revealed that maximum orthologous between Arabidopsis, rice followed by wheat. Segment duplication events have contributed to the expansion of the AP2/ERF family and made this family larger than rice and Arabidopsis. Additionally, AP2/ERF genes were differentially expressed in wheat seedlings under the stress treatments of heat, salt, and drought, and expression profiles were verified by qRT-PCR. Remarkably, the RNA-seq data exposed that AP2/ERF gene family might play a vital role in stress-related. Taken together, our findings provided useful and helpful information to understand the molecular mechanism and evolution of the AP2/ERF gene family in wheat.

Genome-wide investigation and expression analysis of APETALA-2 transcription factor subfamily reveals its evolution, expansion and regulatory role in abiotic stress responses in Indica Rice (Oryza sativa L. ssp. indica)

Rice is an important cereal crop that serves as staple food for more than half of the world population. Abiotic stresses resulting from changing climatic conditions are continuously threating its yield and production. Genes in APETALA-2 (AP2) family encode transcriptional regulators implicated during regulation of developmental processes and abiotic stress responses but their identification and characterization in indica rice was still missing. In this context, twenty-six genes distributed among eleven chromosomes in Indica rice encoding AP2 transcription-factor subfamily were identified and their diverse haplotypes were studied. Phylogenetic analysis of OsAP2 TF family-members grouped them into three clades indicating conservation of clades among cereals. Segmental duplications were observed to be principal route of evolution, supporting the higher positive selection-pressure, which were estimated to be originated about 10.57 to 56.72 million years ago (MYA). Conserved domain analysis and intron-exon distribution pattern of identified OsAP2s revealed their exclusive distribution among the specific clades of the phylogenetic tree. Moreover, the members of osa-miR172 family were also identified potentially targeting four OsAP2 genes. The real-time quantitative expression profiling of OsAP2s under heat stress conditions in contrasting indica rice genotypes revealed the differential expression pattern of OsAP2s (6 genes up-regulated and 4 genes down-regulated) in stress- and genotype-dependent manner. These findings unveiled the evolutionary pathways of AP2-TF in rice, and can help the functional characterization under developmental and stress responses.

Evolutionary history of the C-repeat binding factor/dehydration-responsive element-binding 1 (CBF/DREB1) protein family in 43 plant species and characterization of CBF/DREB1 proteins in Solanum tuberosum

Background

Plants are easily affected by temperature variations, and high temperature (heat stress) and low temperature (cold stress) will lead to poor plant development and reduce crop yields. Therefore, it is very important to identify resistance genes for improving the ability of plants to resist heat stress or cold stress by using modern biotechnology. Members of the C-repeat binding factor/Dehydration responsive element-binding 1 (CBF/DREB1) protein family are related to the stress resistance of many plant species. These proteins affect the growth and development of plants and play vital roles during environmental stress (cold, heat, drought, salt, etc.). In this study, we identified CBF/DREB1 genes from 43 plant species (including algae, moss, ferns, gymnosperms, angiosperms) by using bioinformatic methods to clarify the characteristics of the CBF/DREB1 protein family members and their functions in potato under heat and cold stresses.

Results

In this study, we identified 292 CBF/DREB1 proteins from 43 plant species. However, no CBF/DREB1 protein was found in algae, moss, ferns, or gymnosperms; members of this protein family exist only in angiosperms. Phylogenetic analysis of all the CBF/DREB1 proteins revealed five independent groups. Among them, the genes of group I do not exist in eudicots and are found only in monocots, indicating that these genes have a special effect on monocots. The analysis of motifs, gene duplication events, and the expression data from the PGSC website revealed the gene structures, evolutionary relationships, and expression patterns of the CBF/DREB1 proteins. In addition, analysis of the transcript levels of the 8 CBF/DREB1 genes in potato (Solanum tuberosum) under low-temperature and high-temperature stresses showed that these genes were related to temperature stresses. In particular, the expression levels of StCBF3 and StCBF4 in the leaves, stems, and roots significantly increased under high-temperature conditions, which suggested that StCBF3 and StCBF4 may be closely related to heat tolerance in potato.

Conclusion

Overall, members of the CBF/DREB1 protein family exist only in angiosperms and plays an important role in the growth and development of plants. In addition, the CBF/DREB1 protein family is related to the heat and cold resistance of potato. Our research revealed the evolution of the CBF/DREB1 family, and is useful for studying the precise functions of the CBF/DREB1 proteins when the plants are developing and are under temperature stress.

Genome-wide AP2/ERF gene family analysis reveals the classification, structure, expression profiles and potential function in orchardgrass (Dactylis glomerata)

The AP2/ERF transcription factor (TF) family is of great importance in developmental regulation and responses to stress and pathogenic stimuli. Orchardgrass (Dactylis glomerata), a perennial cold-season forage of high quality in the world's temperate zones, contributes to grazing land through mixed sowing with alfalfa (Medicago sativa) and white clover (Trifolium repens). However, little is known about AP2/ERF TFs in orchardgrass. In this study, 193 AP2/ERF genes were classified into five subfamilies and 13 subgroups through phylogenetic analysis. Chromosome structure analysis showed that AP2/ERF family genes in orchardgrass were distributed on seven chromosomes and specific conservative sequences were found in each subgroup. Exon-intron structure and motifs in the same subgroup were almost identical, and the unique motifs contributed to the classification and functional annotation of DgERFs. Expression analysis showed tissue-specific expression of DgERFs in roots and flowers, with most DgERFs widely expressed in roots. The expression levels of each subgroup (subgroups Vc, VIIa, VIIIb, IXa, and XIa) were high at the before-heading and heading stages (BH_DON and H_DON). In addition, 12 DgERFs in various tissues and five DgERFs associated with abiotic stresses were selected for qRT-PCR analysis showed that four dehydration-responsive element binding (DREB) genes and one ERF subfamily gene in orchardgrass were regulated with PEG, heat and salt stresses. DgERF056 belonged to ERF subfamily was involved in the processes of flowering and development stage. This study systematic explored the DgERFs at the genome level for the first time, which lays a foundation for a better understanding of AP2/ERF gene function in Dactylis glomerata and other types of forage.

DREB Genes from Common Bean (Phaseolus vulgaris L.) Show Broad to Specific Abiotic Stress Responses and Distinct Levels of Nucleotide Diversity

We analyzed the nucleotide variability and the expression profile of DREB genes from common bean, a crop of high economic and nutritional value throughout the world but constantly affected by abiotic stresses in cultivation areas. As DREB genes have been constantly associated with abiotic stress tolerance, we systematically categorized 54 putative PvDREB genes distributed in the common bean genome. It involved from AP2 domain location and amino acid conservation analysis (valine at the 14^th position) to the identification of conserved motifs within peptide sequences representing six subgroups (A-1 to A-6) of PvDREB proteins. Four genes (PvDREB1F, PvDREB2A, PvDREB5A, and PvDREB6B) were cloned and analyzed for their expression profiles under abiotic stresses and their nucleotide and amino acid diversity in genotypes of Andean and Mesoamerican origin, showing distinct patterns of expression and nucleotide variability. PvDREB1F and PvDREB5A showed high relative inducibilities when genotypes of common bean were submitted to stresses by drought, salt, cold, and ABA. PvDREB2A inducibility was predominantly localized to the stem under drought. PvDREB6B was previously described as an A-2 (DREB2) gene, but a detailed phylogenetic analysis and its expression profile clearly indicated it belongs to group A-6. PvDREB6B was found as a cold- and dehydration-responsive gene, mainly in leaves. Interestingly, PvDREB6B also showed a high nucleotide and amino acid diversity within its coding region, in comparison to the others, implicating in several nonsynonymous amino acid substitutions between Andean and Mesoamerican genotypes. The expression patterns and nucleotide diversity of each DREB found in this study revealed fundamental characteristics for further research aimed at understanding the molecular mechanisms associated with drought, salt, and cold tolerance in common bean, which could be performed based on association mapping and functional analyses.

A tuber mustard AP2/ERF transcription factor gene, BjABR1, functioning in abscisic acid and abiotic stress responses, and evolutionary trajectory of the ABR1 homologous genes in Brassica species

The AP2/ERF superfamily of transcription factors is one of the largest transcription factor families in plants and plays an important role in plant development processes and stress responses. In this study, BjABR1, an AP2/ERF superfamily gene, from tuber mustard (Brassica juncea var. tumida Tsen et Lee), sharing high amino acid sequence similarity with the AtABR1 (Arabidopsis thaliana AP2-like abscisic acid repressor 1) gene, were performed functional research, and the ABR1 homologous genes in Brassica species were identified and performed phylogenetic analysis. The promoter sequence of BjABR1 contained many phytohormone- and stress-related cis-elements; ABA (abscisic acid) and abiotic stresses can induce BjABR1 expression in tuber mustard; overexpression of BjABR1 in Arabidopsis can alleviate plant sensitivity to ABA and salt and osmotic stresses, and the alleviation may be due to changes in stress/ABA-induced gene expression. These results indicated that BjABR1 functions in ABA and abiotic stress responses. By BLAST searches against the genome database of five Brassica species (three diploids, B. rapa, B. nigra, and B. oleracea, and two allotetraploid, B. juncea and B. napus) using the protein sequence of AtABR1, 3, 3, 3, 6, and 5 ABR1 homologous genes in B. nigra, B. rapa, B. oleracea, B. juncea, and B. napus were identified, respectively, and they shared high sequence similarity. By sequence analysis, annotation mistakes of the protein-coding regions of two ABR1 homologous genes, GSBRNA2T00134741001 and BjuB007684, were found and corrected. Then, the evolution analysis of these ABR1 homologous genes showed that the ancestor of the three diploid species had three ABR1 homologous genes and each diploid inherited all the three genes from their ancestor; then, allotetraploid B. juncea inherited all the six genes from B. rapa and B. nigra with no gene lost, while allotetraploid B. napus inherited all the three genes from B. oleracea and two genes from B. rapa with one gene lost, indicating that ABR1 homologous genes possessed greater hereditary conservation in Brassica species. The ABR1 homologous genes between B. rapa and B. oleracea shared much higher sequence similarity compared to that of B. nigra in diploid species, indicating that ABR1 homologous genes in B. nigra had experienced more rapid evolution, and B. rapa and B. oleracea may share closer relationship compared to B. nigra. Moreover, the spatial and temporal expression analysis of six ABR1 homologous genes of tuber mustard showed that they possessed different expression models. These results imply that ABR1 homologous genes are important to Brassica plants, and they may possess similar function in ABA and abiotic stress responses but play a role in different tissues and growing stages of plant. This study will provide the foundation to the functional research of ABR1 homologous genes in the Brassica species and help to reveal and understand the evolution mechanisms of Brassica species.

Characterization of the APETALA2/Ethylene-responsive factor (AP2/ERF) transcription factor family in sunflower

One of the most prominent families of genes in plants is the AP2/ERF which play an important role in regulating plant growth and responses to various stresses. In this research, a genome-wide survey was conducted to recognize the AP2/ERF genes in sunflower (Helianthus annuus L.), and a total of 288 HaAP2/ERF was obtained. Phylogenetic analysis divided them into four sub-families, including 248 ERF, 4 RAV and 35 AP2, and one subgroup of the Soloist family. Localization of chromosome, gene structure, the conserved motif, gene ontology, interaction networks, homology modeling, the modeling of cis-regulatory elements and the analysis of events in the duplication of genes were carried out for HaAP2/ERF genes. Finally, 9AP2/ERF genes were chosen to confirm the gene expression of the selected genes in leaf and root tissues in various abiotic stress conditions by qPCR. The results confirmed that AP2/ERFs genes could effectively resist abiotic stress. Also, proline content was studied under drought, salinity, cold and heat stress. The results indicated that proline was increased under abiotic stress. This research has been done for the first time to determine the HaAP2/ERF family, which prepared valuable data for the evolutionary and practical research regarding AP2/ERF in sunflower.

Mitochondrial Biogenesis in Diverse Cauliflower Cultivars under Mild and Severe Drought. Impaired Coordination of Selected Transcript and Proteomic Responses, and Regulation of Various Multifunctional Proteins

Mitochondrial responses under drought within Brassica genus are poorly understood. The main goal of this study was to investigate mitochondrial biogenesis of three cauliflower (Brassica oleracea var. botrytis) cultivars with varying drought tolerance. Diverse quantitative changes (decreases in abundance mostly) in the mitochondrial proteome were assessed by two-dimensional gel electrophoresis (2D PAGE) coupled with liquid chromatography-tandem mass spectrometry (LC-MS/MS). Respiratory (e.g., complex II, IV (CII, CIV) and ATP synthase subunits), transporter (including diverse porin isoforms) and matrix multifunctional proteins (e.g., components of RNA editing machinery) were diversely affected in their abundance under two drought levels. Western immunoassays showed additional cultivar-specific responses of selected mitochondrial proteins. Dehydrin-related tryptic peptides (found in several 2D spots) immunopositive with dehydrin-specific antisera highlighted the relevance of mitochondrial dehydrin-like proteins for the drought response. The abundance of selected mRNAs participating in drought response was also determined. We conclude that mitochondrial biogenesis was strongly, but diversely affected in various cauliflower cultivars, and associated with drought tolerance at the proteomic and functional levels. However, discussed alternative oxidase (AOX) regulation at the RNA and protein level were largely uncoordinated due to the altered availability of transcripts for translation, mRNA/ribosome interactions, and/or miRNA impact on transcript abundance and translation.

Genome-Wide Identification, Characterization, and Expression Profiling of Glutathione S-Transferase (GST) Family in Pumpkin Reveals Likely Role in Cold-Stress Tolerance

Plant growth and development can be adversely affected by cold stress, limiting productivity. The glutathione S-transferase (GST) family comprises important detoxifying enzymes, which play major roles in biotic and abiotic stress responses by reducing the oxidative damage caused by reactive oxygen species. Pumpkins (Cucurbitamaxima) are widely grown, economically important, and nutritious; however, their yield can be severely affected by cold stress. The identification of putative candidate genes responsible for cold-stress tolerance, including the GST family genes, is therefore vital. For the first time, we identified 32 C. maxima GST (CmaGST) genes using a combination of bioinformatics approaches and characterized them by expression profiling. These CmaGST genes represent seven of the 14 known classes of plant GSTs, with 18 CmaGSTs categorized into the tau class. The CmaGSTs were distributed across 13 of pumpkin's 20 chromosomes, with the highest numbers found on chromosomes 4 and 6. The large number of CmaGST genes resulted from gene duplication; 11 and 5 pairs of CmaGST genes were segmental- and tandem-duplicated, respectively. In addition, all CmaGST genes showed organ-specific expression. The expression of the putative GST genes in pumpkin was examined under cold stress in two lines with contrasting cold tolerance: cold-tolerant CP-1 (C. maxima) and cold-susceptible EP-1 (Cucurbita moschata). Seven genes (CmaGSTU3, CmaGSTU7, CmaGSTU8, CmaGSTU9, CmaGSTU11, CmaGSTU12, and CmaGSTU14) were highly expressed in the cold-tolerant line and are putative candidates for use in breeding cold-tolerant crop varieties. These results increase our understanding of the cold-stress-related functions of the GST family, as well as potentially enhancing pumpkin breeding programs.

Characterization of OglDREB2A gene from African rice (Oryza glaberrima), comparative analysis and its transcriptional regulation under salinity stress

In this study, AP2 DNA-binding domain-containing transcription factor, OglDREB2A, was cloned from the African rice (Oryza glaberrima) and compared with 3000 rice genotypes. Further, the phylogenetic and various structural analysis was performed using in silico approaches. Further, to understand its allelic variation in rice, SNPs and indels were detected among the 3000 rice genotypes which indicated that while coding region is highly conserved, yet noncoding regions such as UTR and intron contained most of the variation. Phylogenetic analysis of the OglDREB2A sequence in different Oryza as well as in diverse eudicot species revealed that DREB from various Oryza species were diversed much earlier than other genes. Further, structural features and in silico analyses provided insights into different properties of OglDREB2A protein. The neutrality test on the coding region of OglDREB2A from different genotypes of O. glaberrima showed the lack of selection in this gene. Among the different developmental stages, it was upregulated at tillering and flag leaf under salinity treatment indicating its positive role in seedling and reproductive stage tolerance. Real-time PCR analysis also indicated the conserve expression pattern of this gene under salinity stress across the three different Oryza species having different degree of salinity tolerance.

RAV transcription factors are essential for disease resistance against cassava bacterial blight via activation of melatonin biosynthesis genes

With 1 AP2 domain and 1 B3 domain, 7 MeRAVs in apetala2/ethylene response factor (AP2/ERF) gene family have been identified in cassava. However, the in vivo roles of these remain unknown. Gene expression assays showed that the transcripts of MeRAVs were commonly regulated after Xanthomonas axonopodis pv manihotis (Xam) and MeRAVs were specifically located in plant cell nuclei. Through virus-induced gene silencing (VIGS) in cassava, we found that MeRAV1 and MeRAV2 are essential for plant disease resistance against cassava bacterial blight, as shown by the bacterial propagation of Xam in plant leaves. Through VIGS in cassava leaves and overexpression in cassava leave protoplasts, we found that MeRAV1 and MeRAV2 positively regulated melatonin biosynthesis genes and the endogenous melatonin level. Further investigation showed that MeRAV1 and MeRAV2 are direct transcriptional activators of 3 melatonin biosynthesis genes in cassava, as evidenced by chromatin immunoprecipitation-PCR in cassava leaf protoplasts and electrophoretic mobility shift assay. Moreover, cassava melatonin biosynthesis genes also positively regulated plant disease resistance. Taken together, this study identified MeRAV1 and MeRAV2 as common and upstream transcription factors of melatonin synthesis genes in cassava and revealed a model of MeRAV1 and MeRAV2-melatonin biosynthesis genes-melatonin level in plant disease resistance against cassava bacterial blight.

Transcription Factors and Their Roles in Signal Transduction in Plants under Abiotic Stresses

In agricultural production, abiotic stresses are known as the main disturbance leading to negative impacts on crop performance. Research on elucidating plant defense mechanisms against the stresses at molecular level has been addressed for years in order to identify the major contributors in boosting the plant tolerance ability. From literature, numerous genes from different species, and from both functional and regulatory gene categories, have been suggested to be on the list of potential candidates for genetic engineering. Noticeably, enhancement of plant stress tolerance by manipulating expression of Transcription Factors (TFs) encoding genes has emerged as a popular approach since most of them are early stress-responsive genes and control the expression of a set of downstream target genes. Consequently, there is a higher chance to generate novel cultivars with better tolerance to either single or multiple stresses. Perhaps, the difficult task when deploying this approach is selecting appropriate gene(s) for manipulation. In this review, on the basis of the current findings from molecular and post-genomic studies, our interest is to highlight the current understanding of the roles of TFs in signal transduction and mediating plant responses towards abiotic stressors. Furthermore, interactions among TFs within the stress-responsive network will be discussed. The last section will be reserved for discussing the potential applications of TFs for stress tolerance improvement in plants.

Transcription Factors and Their Roles in Signal Transduction in Plants under Abiotic Stresses

In agricultural production, abiotic stresses are known as the main disturbance leading to negative impacts on crop performance. Research on elucidating plant defense mechanisms against the stresses at molecular level has been addressed for years in order to identify the major contributors in boosting the plant tolerance ability. From literature, numerous genes from different species, and from both functional and regulatory gene categories, have been suggested to be on the list of potential candidates for genetic engineering. Noticeably, enhancement of plant stress tolerance by manipulating expression of Transcription Factors (TFs) encoding genes has emerged as a popular approach since most of them are early stress-responsive genes and control the expression of a set of downstream target genes. Consequently, there is a higher chance to generate novel cultivars with better tolerance to either single or multiple stresses. Perhaps, the difficult task when deploying this approach is selecting appropriate gene(s) for manipulation. In this review, on the basis of the current findings from molecular and post-genomic studies, our interest is to highlight the current understanding of the roles of TFs in signal transduction and mediating plant responses towards abiotic stressors. Furthermore, interactions among TFs within the stress-responsive network will be discussed. The last section will be reserved for discussing the potential applications of TFs for stress tolerance improvement in plants.

Overexpression of ERF1-V from Haynaldia villosa Can Enhance the Resistance of Wheat to Powdery Mildew and Increase the Tolerance to Salt and Drought Stresses

The APETALA 2/Ethylene-responsive element binding factor (AP2/ERF) transcription factor gene family is widely involved in the biotic and abiotic stress regulation. Haynaldia villosa (VV, 2n = 14), a wild species of wheat, is a potential gene pool for wheat improvement. H. villosa confers high resistance to several wheat diseases and high tolerance to some abiotic stress. In this study, ERF1-V, an ethylene-responsive element-binding factor gene of the AP2/ERF transcription factor gene family from wild H. villosa, was cloned and characterized. Sequence and phylogenetic analysis showed that ERF1-V is a deduced B2 type ERF gene. ERF1-V was first identified as a Blumeria graminis f. sp. tritici (Bgt) up-regulated gene, and later found to be induced by drought, salt and cold stresses. In responses to hormones, ERF1-V was up-regulated by ethylene and abscisic acid, but down-regulated by salicylic acid and jasmonic acid. Over expression of ERF1-V in wheat could improve resistance to powdery mildew, salt and drought stress. Chlorophyll content, malondialdehyde content, superoxide dismutase and peroxidase activity were significantly differences between the recipient Yangmai158 and the transgenic plants following salt treatment. Furthermore, the expression levels of some stress responsive genes were differences after drought or salt treatments. Although ERF1-V was activated by the constitutive promoter, the agronomic traits, including flowering time, plant height, effective tiller number, spikelet number per spike and grain size, did not changed significantly. ERF1-V is a valuable gene for wheat improvement by genetic engineering.

Transcriptomic changes reveal gene networks responding to the overexpression of a blueberry DWARF AND DELAYED FLOWERING 1 gene in transgenic blueberry plants

BACKGROUND

Constitutive expression of the CBF/DREB1 for increasing freezing tolerance in woody plants is often associated with other phenotypic changes including dwarf plant and delayed flowering. These phenotypic changes have been observed when Arabidopsis DWARF AND DELAYED FLOWERING 1 (DDF1) was overexpressed in A. thaliana plants. To date, the DDF1 orthologues have not been studied in woody plants. The aim of this study is to investigate transcriptomic responses to the overexpression of blueberry (Vaccinium corymbosum) DDF1 (herein, VcDDF1-OX).

RESULTS

The VcDDF1-OX resulted in enhanced freezing tolerance in tetraploid blueberry plants and did not result in significant changes in plant size, chilling requirement, and flowering time. Comparative transcriptome analysis of transgenic 'Legacy-VcDDF1-OX' plants containing an overexpressed VcDDF1 with non-transgenic highbush blueberry 'Legacy' plants revealed the VcDDF1-OX derived differentially expressed (DE) genes and transcripts in the pathways of cold-response, plant flowering, DELLA proteins, and plant phytohormones. The increase in freezing tolerance was associated to the expression of cold-regulated genes (CORs) and the ethylene pathway genes. The unchanged plant size, dormancy and flowering were due to the minimal effect of the VcDDF1-OX on the expression of DELLA proteins, flowering pathway genes, and the other phytohormone genes related to plant growth and development. The DE genes in auxin and cytokinin pathways suggest that the VcDDF1-OX has also altered plant tolerance to drought and high salinity.

CONCLUSION

A DDF1 orthologue in blueberry functioned differently from the DDF1 reported in Arabidopsis. The overexpression of VcDDF1 or its orthologues is a new approach to increase freezing tolerance of deciduous woody plant species with no obvious effect on plant size and plant flowering time.

Special trends in CBF and DREB2 groups in Eucalyptus gunnii vs Eucalyptus grandis suggest that CBF are master players in the trade-off between growth and stress resistance

Annotation of the Eucalyptus grandis genome showed a large amplification of the dehydration-responsive element binding 1/C-repeat binding factor (DREB1/CBF) group without recent DREB2 gene duplication compared with other plant species. The present annotation of the CBF and DREB2 genes from a draft of the Eucalyptus gunnii genome sequence reveals at least one additional CBF copy in the E. gunnii genome compared with E. grandis, suggesting that this group is still evolving, unlike the DREB2 group. This study aims to investigate the redundancy/neo- or sub-functionalization of the duplicates and the relative involvement of the two groups in abiotic stress responses in both E. grandis and E. gunnii (lower growth but higher cold resistance). A comprehensive transcriptional analysis using high-throughput quantitative real-time polymerase chain reaction (qRT-PCR) was performed on leaves, stems and roots from the two Eucalyptus species after cold, heat or drought treatment. A large CBF cluster accounted for most of the cold response in all the organs, whereas heat and drought responses mainly involved a small CBF cluster and the DREB2 genes. In addition, CBF putative target genes, known to be involved in plant tolerance and development, were found to be cold-regulated. The higher transcript amounts of both the CBF and target genes in the cold tolerant E. gunnii contrasted with the higher CBF induction rates in the fast growing E. grandis. Altogether, the present results, in agreement with previous data about Eucalyptus transgenic lines over-expressing CBF, suggest that these factors, which promote both stress protection and growth limitation, participate in the trade-off between growth and resistance in this woody species.

Comparative genome-wide phylogenetic and expression analysis of SBP genes from potato (Solanum tuberosum)

Flowering time is a very important phase in transition to reproductive stage of life in higher plants. SQUAMOSA promoter-binding protein (SBP) gene family encodes plant-specific transcription factors that are involved in regulation of several developmental processes, especially flowering. Although SBP-box genes have been identified in different plants, there have been no study indicating the regulatory effect of SBP box in potato flowering. Here, we report for the first time the identification and characterization of SBP-box transcription factors as well as determination of expression level of SBP-box genes in Solanum tuberosum L. an important crop worldwide. Fifteen different SBP-box transcription factor genes were identified in potato genome. They were found to be distributed in nine chromosomes and eight of them included miR156 and miRNA157 target sites. Characterization of amino acid sequences were performed and protein interactions were predicted. In addition, expression levels of five S. tuberosum SBP-box genes were analysed by both in silico and qRT-PCR. All these results provide a better understanding of functional role of SBP-box gene family members in flowering time in potato.

Genome-Wide Identification of AP2/ERF Transcription Factors in Cauliflower and Expression Profiling of the ERF Family under Salt and Drought Stresses

The AP2/ERF transcription factors (TFs) comprise one of the largest gene superfamilies in plants. These TFs perform vital roles in plant growth, development, and responses to biotic and abiotic stresses. In this study, 171 AP2/ERF TFs were identified in cauliflower (Brassica oleracea L. var. botrytis), one of the most important horticultural crops in Brassica. Among these TFs, 15, 9, and 1 TFs were classified into the AP2, RAV, and Soloist family, respectively. The other 146 TFs belong to ERF family, which were further divided into the ERF and DREB subfamilies. The ERF subfamily contained 91 TFs, while the DREB subfamily contained 55 TFs. Phylogenetic analysis results indicated that the AP2/ERF TFs can be classified into 13 groups, in which 25 conserved motifs were confirmed. Some motifs were group- or subgroup- specific, implying that they are significant to the functions of the AP2/ERF TFs of these clades. In addition, 35 AP2/ERF TFs from the 13 groups were selected randomly and then used for expression pattern analysis under salt and drought stresses. The majority of these AP2/ERF TFs exhibited positive responses to these stress conditions. In specific, Bra-botrytis-ERF054a, Bra-botrytis-ERF056, and Bra-botrytis-CRF2a demonstrated rapid responses. By contrast, six AP2/ERF TFs were showed to delay responses to both stresses. The AP2/ERF TFs exhibiting specific expression patterns under salt or drought stresses were also confirmed. Further functional analysis indicated that ectopic overexpression of Bra-botrytis-ERF056 could increase tolerance to both salt and drought treatments. These findings provide new insights into the AP2/ERF TFs present in cauliflower, and offer candidate AP2/ERF TFs for further studies on their roles in salt and drought stress tolerance.

Ethylene responsive transcription factor ERF109 retards PCD and improves salt tolerance in plant

BACKGROUND

The ultimate goal of this work was to detect the role of transcription factors (TFs) concordantly expressed with genes related to programmed cell death (PCD) during PCD and salt stress. This work was based on the hypothesis that TFs and their driven genes likely co-express under different stimuli. The conserved superfamily ethylene responsive factor (AP2/ERF) draw attention of the present study as it participates in the response to biotic and abiotic stimuli as well as to program cell death (PCD).

RESULTS

RNA-Seq analysis was done for tobacco (N. benthamiana) leaves exposed to oxalic acid (OA) at 20 mM for 0, 2, 6, 12 and 24 h to induce PCD. Genes up-regulated after 2 h of OA treatment with known function during PCD were utilized as landmarks to select TFs with concordant expression. Knockdown mutants of these TFs were generated in tobacco via virus induced gene silencing (VIGS) in order to detect their roles during PCD. Based on the results of PCD assay, knockout (KO) T-DNA insertion mutants of Arabidopsis as well as over-expression lines of two selected TFs, namely ERF109 and TFIID5, analogs to those in tobacco, were tested under salt stress (0, 100, 150 and 200 mM NaCl).

CONCLUSIONS

Results of knockdown mutant tobacco cells confirmed the influence of these two TFs during PCD. Knockout insertion mutants and over-expression lines indicated the role of ERF109 in conferring salt tolerance in Arabidopsis.

De Novo Assembly and Transcriptome Analysis of Bulb Onion (Allium cepa L.) during Cold Acclimation Using Contrasting Genotypes

Bulb onion (Allium cepa) is the second most widely cultivated and consumed vegetable crop in the world. During winter, cold injury can limit the production of bulb onion. Genomic resources available for bulb onion are still very limited. To date, no studies on heritably durable cold and freezing tolerance have been carried out in bulb onion genotypes. We applied high-throughput sequencing technology to cold (2°C), freezing (-5 and -15°C), and control (25°C)-treated samples of cold tolerant (CT) and cold susceptible (CS) genotypes of A. cepa lines. A total of 452 million paired-end reads were de novo assembled into 54,047 genes with an average length of 1,331 bp. Based on similarity searches, these genes were aligned with entries in the public non-redundant (nr) database, as well as KEGG and COG database. Differentially expressed genes (DEGs) were identified using log₁₀ values with the FPKM method. Among 5,167DEGs, 491 genes were differentially expressed at freezing temperature compared to the control temperature in both CT and CS libraries. The DEG results were validated with qRT-PCR. We performed GO and KEGG pathway enrichment analyses of all DEGs and iPath interactive analysis found 31 pathways including those related to metabolism of carbohydrate, nucleotide, energy, cofactors and vitamins, other amino acids and xenobiotics biodegradation. Furthermore, a large number of molecular markers were identified from the assembled genes, including simple sequence repeats (SSRs) 4,437 and SNP substitutions of transition and transversion types of CT and CS. Our study is the first to provide a transcriptome sequence resource for Allium spp. with regard to cold and freezing stress. We identified a large set of genes and determined their DEG profiles under cold and freezing conditions using two different genotypes. These data represent a valuable resource for genetic and genomic studies of Allium spp.

Genome-Wide Analysis of APETALA2/Ethylene-Responsive Factor (AP2/ERF) Gene Family in Barley (Hordeum vulgare L.)

APETALA2/Ethylene-Responsive Factor (AP2/ERF) gene family is plant specific transcription factor. It plays critical roles in development process, tolerance to biotic and abiotic stresses, and responses to plant hormones. However, limited data are available on the contributions of AP2/ERF gene family in barley (Hordeum vulgare L.). In the present study, 121 HvAP2/ERF genes in barley were identified by using bioinformatics methods. A total of 118 HvAP2/ERF (97.5%) genes were located on seven chromosomes. According to phylogenetic classification of AP2/ERF family in Arabidopsis, HvAP2/ERF proteins were divided into AP2 (APETALA2), RAV (Related to ABI3/VP), DREB (dehydration responsive element binding), ERF (ethylene responsive factors) and soloist sub families. The analysis of duplication events indicated that tandem repeat and segmental duplication contributed to the expansion of the AP2/ERF family in barley. HvDREB1s/2s genes displayed various expression patterns under abiotic stress and phytohormone. Taken together, the data generated in this study will be useful for genome-wide analysis to determine the precise role of the HvAP2/ERF gene during barley development, abiotic stress and phytohormone responses with the ultimate goal of improving crop production.

Glutathione Transferases Superfamily: Cold-Inducible Expression of Distinct GST Genes in Brassica oleracea

Plants, as sessile organisms, can suffer serious growth and developmental consequences under cold stress conditions. Glutathione transferases (GSTs, EC 2.5.1.18) are ubiquitous and multifunctional conjugating proteins, which play a major role in stress responses by preventing oxidative damage by reactive oxygen species (ROS). Currently, understanding of their function(s) during different biochemical and signaling pathways under cold stress condition remain unclear. In this study, using combined computational strategy, we identified 65 Brassica oleracea glutathione transferases (BoGST) and characterized them based on evolutionary analysis into 11 classes. Inter-species and intra-species duplication was evident between BoGSTs and Arabidopsis GSTs. Based on localization analyses, we propose possible pathways in which GST genes are involved during cold stress. Further, expression analysis of the predicted putative functions for GST genes were investigated in two cold contrasting genotypes (cold tolerance and susceptible) under cold condition, most of these genes were highly expressed at 6 h and 1 h in the cold tolerant (CT) and cold susceptible (CS) lines, respectively. Overall, BoGSTU19, BoGSTU24, BoGSTF10 are candidate genes highly expressed in B. oleracea. Further investigation of GST superfamily in B. oleracea will aid in understanding complex mechanism underlying cold tolerance in plants.

Characterization and abiotic stress-responsive expression analysis of SGT1 genes in Brassica oleracea

SGT1 genes are involved in enhancing plant responses to various biotic and abiotic stresses. Brassica oleracea is known to contain two types of SGT1 genes, namely suppressor of G2 allele of SKP1 and suppressor of GCR2. In this study, through systematic analysis, four putative SGT1 genes were identified and characterized in B. oleracea. In phylogenetic analysis, the genes clearly formed separate groups, namely BolSGT1a, BolSGT1b (both suppressor of G2 allele of SKP1 types), and BolSGT1 (suppressor of GCR2). Functional domain analysis and organ-specific expression patterns suggested possible roles for BolSGT1 genes during stress conditions. BolSGT1 genes showed significant changes in expression in response to heat, cold, drought, salt, or ABA treatment. Interaction network analysis supported the expression analysis, and showed that the BolSGT1a and BolSGT1b genes are strongly associated with co-regulators during stress conditions. However, the BolSGT1 gene did not show any strong association. Hence, BolSGT1 might be a stress resistance-related gene that functions without a co-regulator. Our results show that BolSGT1 genes are potential target genes to improve B. oleracea resistance to abiotic stresses such as heat, cold, and salt.

Dynamic transcriptome analysis reveals AP2/ERF transcription factors responsible for cold stress in rapeseed (Brassica napus L.)

The APETALA2/ethylene response factor (AP2/ERF) transcription factor (TF) superfamily plays an important regulatory role in signal transduction of the plant responses to various stresses including low temperature. Significant progress has been made in understanding the mechanism of cold resistance in Brassica napus, an important oilseed crop. However, comprehensive studies on the induction and activity of these TFs under low temperature have been lacking. In this study, 132 AP2/ERF genes were identified by transcriptome sequencing of rapeseed leaves exposed to 0, 2, 6, 12, and 24 h of low (4 °C) temperature stress. The genes were classified into 4 subfamilies (AP2, DREB, ERF, and RAV) and 13 subgroups, among which the DREB subfamily and ERF subfamily contained 114 genes, no genes were assigned to soloist or DREB A3 subgroups. One hundred and eighteen genes were located on chromosomes A1 to C9. GO functional analysis and promoter sequence analysis revealed that these genes are involved in many molecular pathways that may enhance cold resistance in plants, such as the low-temperature responsiveness, methyl jasmonate, abscisic acid, and ethylene-responsiveness pathways. Their expression patterns revealed dynamic control at different times following initiation of cold stress; the RAV and DREB subfamilies were expressed at the early stage of cold stress, whereas the AP2 subfamily was expressed later. Quantitative PCR analyses of 13 cold-induced AP2/ERF TFs confirmed the accuracy of above results. This study is the first dynamic analysis of the AP2/ERF TFs responsible for cold stress in rapeseed. These findings will serve as a reference for future functional research on transcription in rapeseed.

Genome-Wide Survey of Flavonoid Biosynthesis Genes and Gene Expression Analysis between Black- and Yellow-Seeded Brassica napus

Flavonoids, the compounds that impart color to fruits, flowers, and seeds, are the most widespread secondary metabolites in plants. However, a systematic analysis of these loci has not been performed in Brassicaceae. In this study, we isolated 649 nucleotide sequences related to flavonoid biosynthesis, i.e., the Transparent Testa (TT) genes, and their associated amino acid sequences in 17 Brassicaceae species, grouped into Arabidopsis or Brassicaceae subgroups. Moreover, 36 copies of 21 genes of the flavonoid biosynthesis pathway were identified in Arabidopsis thaliana, 53 were identified in Brassica rapa, 50 in Brassica oleracea, and 95 in B. napus, followed the genomic distribution, collinearity analysis and genes triplication of them among Brassicaceae species. The results showed that the extensive gene loss, whole genome triplication, and diploidization that occurred after divergence from the common ancestor. Using qRT-PCR methods, we analyzed the expression of 18 flavonoid biosynthesis genes in 6 yellow- and black-seeded B. napus inbred lines with different genetic background, found that 12 of which were preferentially expressed during seed development, whereas the remaining genes were expressed in all B. napus tissues examined. Moreover, 14 of these genes showed significant differences in expression level during seed development, and all but four of these (i.e., BnTT5, BnTT7, BnTT10, and BnTTG1) had similar expression patterns among the yellow- and black-seeded B. napus. Results showed that the structural genes (BnTT3, BnTT18, and BnBAN), regulatory genes (BnTTG2 and BnTT16) and three encoding transfer proteins (BnTT12, BnTT19, and BnAHA10) might play an crucial roles in the formation of different seed coat colors in B. napus. These data will be helpful for illustrating the molecular mechanisms of flavonoid biosynthesis in Brassicaceae species.

Origination, Expansion, Evolutionary Trajectory, and Expression Bias of AP2/ERF Superfamily in Brassica napus

The AP2/ERF superfamily, one of the most important transcription factor families, plays crucial roles in response to biotic and abiotic stresses. So far, a comprehensive evolutionary inference of its origination and expansion has not been available. Here, we identified 515 AP2/ERF genes in B. napus, a neo-tetraploid forming ~7500 years ago, and found that 82.14% of them were duplicated in the tetraploidization. A prominent subgenome bias was revealed in gene expression, tissue-specific, and gene conversion. Moreover, a large-scale analysis across plants and alga suggested that this superfamily could have been originated from AP2 family, expanding to form other families (ERF, and RAV). This process was accompanied by duplicating and/or alternative deleting AP2 domain, intragenic domain sequence conversion, and/or by acquiring other domains, resulting in copy number variations, alternatively contributing to functional innovation. We found that significant positive selection occurred at certain critical nodes during the evolution of land plants, possibly responding to changing environment. In conclusion, the present research revealed origination, functional innovation, and evolutionary trajectory of the AP2/ERF superfamily, contributing to understanding their roles in plant stress tolerance.

Genome-wide investigation and expression analysis of AP2-ERF gene family in salt tolerant common bean

Apetala2-ethylene-responsive element binding factor (AP2-ERF) superfamily with common AP2-DNA binding domain have developmentally and physiologically important roles in plants. Since common bean genome project has been completed recently, it is possible to identify all of the AP2-ERF genes in the common bean genome. In this study, a comprehensive genome-wide in silico analysis identified 180 AP2-ERF superfamily genes in common bean (Phaseolus vulgaris). Based on the amino acid alignment and phylogenetic analyses, superfamily members were classified into four subfamilies: DREB (54), ERF (95), AP2 (27) and RAV (3), as well as one soloist. The physical and chemical characteristics of amino acids, interaction between AP2-ERF proteins, cis elements of promoter region of AP2-ERF genes and phylogenetic trees were predicted and analyzed. Additionally, expression levels of AP2-ERF genes were evaluated by in silico and qRT-PCR analyses. In silico micro-RNA target transcript analyses identified nearly all PvAP2-ERF genes as targets of by 44 different plant species' miRNAs were identified in this study. The most abundant target genes were PvAP2/ERF-20-25-62-78-113-173. miR156, miR172 and miR838 were the most important miRNAs found in targeting and BLAST analyses. Interactome analysis revealed that the transcription factor PvAP2-ERF78, an ortholog of Arabidopsis At2G28550, was potentially interacted with at least 15 proteins, indicating that it was very important in transcriptional regulation. Here we present the first study to identify and characterize the AP2-ERF transcription factors in common bean using whole-genome analysis, and the findings may serve as a references for future functional research on the transcription factors in common bean.

Genome-wide characterization and expression analysis of common bean bHLH transcription factors in response to excess salt concentration

Members of basic helix-loop-helix (bHLH) gene family found in all eukaryotes play crucial roles in response to stress. Though, most eukaryotes carry the proteins of this family, biological functions of the most bHLH family members are not deeply evaluated in plants. In this study, we conducted a comprehensive genome-wide analysis of bHLH transcription factors in salt tolerant common bean. We identified 155 bHLH protein-encoding genes (PvbHLH) by using in silico comparative genomics tools. Based on the phylogenetic tree, PvbHLH genes were classified into 8 main groups with 21 subfamilies. Exon-intron analysis indicated that proteins belonging to same main groups exhibited a closely related gene structure. While, the PvbHLH gene family has been mainly expanded through segmental duplications, a total of 11 tandem duplication were detected. Genome-wide expression analysis of bHLH genes showed that 63 PvbHLH genes were differentially expressed in at least one tissue. Three of them displayed higher expression values in both leaf and root tissues. The in silico micro-RNA target transcript analyses revealed that totally 100 PvHLH genes targeted by 86 plant miRNAs. The most abundant transcripts, which were targeted by all 18 plant miRNA, were belonging to PvHLH-22 and PvHLH-44 genes. The expression of 16 PvbHLH genes in the root and leaf tissues of salt-stressed common bean was evaluated using qRT-PCR. Among them, two of PvbHLHs, PvbHLH-54, PvbHLH-148, were found to be up-regulated in both tissues in correlation with RNA-seq measurements. The results of this study could help improve understanding of biological functions of common bean bHLH family under salt stress. Additionally, it may provide basic resources for analyzing bHLH protein function for improving economic, agronomic and ecological benefit in common bean and other species.

Gene-to-metabolite network for biosynthesis of lignans in MeJA-elicited Isatis indigotica hairy root cultures

Root and leaf tissue of Isatis indigotica shows notable anti-viral efficacy, and are widely used as "Banlangen" and "Daqingye" in traditional Chinese medicine. The plants' pharmacological activity is attributed to phenylpropanoids, especially a group of lignan metabolites. However, the biosynthesis of lignans in I. indigotica remains opaque. This study describes the discovery and analysis of biosynthetic genes and AP2/ERF-type transcription factors involved in lignan biosynthesis in I. indigotica. MeJA treatment revealed differential expression of three genes involved in phenylpropanoid backbone biosynthesis (IiPAL, IiC4H, Ii4CL), five genes involved in lignan biosynthesis (IiCAD, IiC3H, IiCCR, IiDIR, and IiPLR), and 112 putative AP2/ERF transcription factors. In addition, four intermediates of lariciresinol biosynthesis were found to be induced. Based on these results, a canonical correlation analysis using Pearson's correlation coefficient was performed to construct gene-to-metabolite networks and identify putative key genes and rate-limiting reactions in lignan biosynthesis. Over-expression of IiC3H, identified as a key pathway gene, was used for metabolic engineering of I. indigotica hairy roots, and resulted in an increase in lariciresinol production. These findings illustrate the utility of canonical correlation analysis for the discovery and metabolic engineering of key metabolic genes in plants.