- Research article
- Open Access
Genome-wide identification, phylogeny and expression analysis of AP2/ERF transcription factors family in Brachypodium distachyon
- Licao Cui†1,
- Kewei Feng†1,
- Mengxing Wang1,
- Meng Wang1,
- Pingchuan Deng1,
- Weining Song1, 2Email author and
- Xiaojun Nie1Email author
© The Author(s). 2016
- Received: 20 January 2016
- Accepted: 26 July 2016
- Published: 15 August 2016
The AP2/ERF transcription factor is one of the most important gene families in plants, which plays the vital role in regulating plant growth and development as well as in response to diverse stresses. Although AP2/ERFs have been thoroughly characterized in many plant species, little is known about this family in the model plant Brachypodium distachyon, especially those involved in the regulatory network of stress processes.
In this study, a comprehensive genome-wide search was performed to identify AP2/ERF gene family in Brachypodium and a total of 141 BdAP2/ERFs were obtained. Phylogenetic analysis classified them into four subfamilies, of which 112 belonged to ERF, four to RAV and 24 to AP2 as well as one to soloist subfamily respectively, which was in accordance with the number of AP2 domains and gene structure analysis. Chromosomal localization, gene structure, conserved protein motif and cis-regulatory elements as well as gene duplication events analysis were further performed to systematically investigate the evolutionary features of these BdAP2/ERF genes. Furthermore, the regulatory network between BdAP2/ERF and other genes were constructed using the orthology-based method, and 39 BdAP2/ERFs were found to be involved in the regulatory network and 517 network branches were identified. The expression profiles of BdAP2/ERF during development and under diverse stresses were investigated using the available RNA-seq and microarray data and ten tissue-specific and several stress-responsive BdAP2/ERF genes were identified. Finally, 11 AP2/ERF genes were selected to validate their expressions in different tissues and under different stress treatments using RT-PCR method and results verified that these AP2/ERFs were involved in various developmental and physiological processes.
This study for the first time reported the characteristics of the BdAP2/ERF family, which will provide the invaluable information for further evolutionary and functional studies of AP2/ERF in Brachypodium, and also contribute to better understanding the molecular basis for development and stresses tolerance in this model species and beyond.
- Abiotic stress
- Expression profiles
- Gene family
- Transcription factor
Plant growth, development and productivity are adversely affected by numerous abiotic stresses, such as drought, salt and heat. To survive and flourish under these environmental stresses, plants have developed a complicated response mechanism by repressing or inducing the expression of a series of genes with diverse functions. Transcription factors (TF), as an important group of regulatory proteins, play the central roles in regulation network and signaling pathways of plant development and in response to abiotic stresses. Among them, AP2/ERF (APETALA2/Ethylene Responsive Factor) superfamily is one of the biggest plant TFs, which distinguished by one or two highly conserved ethylene-responsive element-binding factor domains that consisted of 50–60 amino acids [1, 2]. Based on sequence similarities and repetitions of AP2 DNA-binding domains, it can be classified into AP2, ERF and RAV families . The members of AP2 family proteins contain two AP2/ERF domains and are further divided into AP2 and AINTEGUMENTA (ANT) monophyletic groups [4, 5], while the members of ERF subfamily possesses a AP2/ERF domain with the specific WLG motif and are subdivided into ten group , of which Group I to IV belong to the DREB subfamily and group V to X belong to the ERF subfamily. The ERF subfamily is characterized by an additional cis-acting element AGCCCGCC of the GCC-box in the promoter regions , whereas the DREB subfamily typically binds to dehydration-responsive element-binding factor, which has a core motif of CCGAC . The RAV family members containing the single AP2/ERF domain and a specific B3 DNA-binding motif . In addition, other members with an AP2-like domain but lacking additional motifs are often defined as Soloist.
Extensive studies have revealed the crucial role of the AP2/ERF genes playing in plant growth, development and stress responses [4, 9–11]. Generally, the AP2 subfamily members were the main factors involving in regulating organ architecture and development, such as leaf epidermal cell determinacy, spikelet meristem differentiation and floral organ patterning  as well as seed mass and seed yield [13, 14], while the RAV subfamily showed the important functions in plant hormone signal transduction, such as ethylene , Brassinosteroid , and also involved in response to biotic and aboitic stresses [17, 18]. Additionally, the DREB, together with other members in ERF subfamily mainly involved in response to biotic and abiotic stresses, such as water deficit , low and high temperature [20, 21] and high salinity .
B. distachyon, belong to Brachypodium tribe Poaceae family which has a close phylogentic relationships with the major cereal crops, including wheat, barley and rye. It has many favorable features, such as small genome (~300 Mb), diploid accessions, self-fertility, a short lifecycle and easy transformation, which make it an ideal model organism for functional genomic studies of temperate grasses, cereals and biofuel crops [23, 24] and now its genome has been completely sequenced . The available genome data facilitated the studies to reveal the gene function and regulation network in this species, and the study of B. distachyon will provide the vital clue for better understand the molecular mechanism of stress response and subsequently improve the abiotic stress tolerance of other cereal crops. So far, the AP2/ERF family has been identified in Arabidopsis , Bamboo , grapevine , maize , peach  and rice . However, to the best of our knowledge, the systematic identification of AP2/ERF family has not been performed in B. distachyon, limiting the further function analysis of this important gene family.
In this study, a genome-wide bioinformatics analysis was conducted to investigate the genomic organization, phylogenetic relationship and expression profiles of AP2/ERF genes in B. distachyon. The chromosomal localization, gene structures, cis-elements in the promoter region as well as gene duplication and evolutionary mechanisms were subsequently analyzed. By using RNA-seq and microarray expression data, the expression profiling of these identified AP2/ERF genes in different tissue as well as under cold and drought stresses was further investigated. Our study provided a basis for further study on the regulation roles of the AP2/ERF family playing in B. distachyon development and in response to biotic and abiotic stresses, which will not only provide the helpful information on the evolutionary mechanism of this TFs family in plant, but also contribute to revealing the molecular mechanism of development and stresses response in B. distachyon and other cereal crops.
Identification of AP2/ERF gene family in Brachpodium genome
The whole genome data of B. distachyon was available at Ensemble plants database (http://plants.ensembl.org/index.html). The predicted protein sequences were downloaded as the dataset for downstream analysis (v1.0.29). The AP2/ERF domain (PF00847) obtained from PFAM database (http://pfam.xfam.org/) was used as the query for Hidden Markov Model (HMM) search using HMMER 3.0 program with a pre-defined threshold of E <1e−5. Furhtermore, the AP2/ERF protein sequences ofArabidopsis and rice were obtained from the plant transcription factor database (http://plntfdb.bio.uni-potsdam.de/v3.0/) and then used as query to search against the Brachpodium protein dataset using the BLASTP program with an e-value of 1e-5 and identity of 50 % as the threshold. Furthermore, HMMER and BLAST hits were compared and parsed and then a self-blast of these sequences was performed to remove the redundancy and no any alternative splice variants were considered. After manual correcting, the putative BdAR2/ERF proteins were obtained. Then, the NCBI-CDD web server (http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) and SMART database (http://smart.embl-heidelberg.de/webcite) were used to further confirm the predicted BdAR2/ERF genes. The theoretical isoelectric point (PI) and molecular weight (MW) of the obtained proteins were conducted by the compute pI/Mw tool in the ExPASy server (http://www.expasy.org/). The subcellular localization prediction of each gene was predicted using the cello web server (http://cello.life.nctu.edu.tw/).
Multiple sequence alignment and phylogenetic analysis
Multiple sequence alignment was performed using Clustal X v2.0  with the default parameters. An un-rooted neighbor joining (NJ) tree with 1000 bootstrap replications was constructed using MEGA 6.0  based on the full-length protein alignment.
Chromosome distribution, gene structure and conserved motif analysis
The chromosome distribution of these genes were obtained from the genome annotation information, and then validated by BLASTN search. The exon-intron organizations and splicing phase of these predicted AR2/ERF genes were also investigated based on the annotation file of B. distachyon genome, and then graphically displayed by the Gene Structure Display Server (http://gsds.cbi.pku.edu.cn/). Conserved motifs or domains were predicted using the MEME Suite web server (http://meme-suite.org/), with the following parameters: maximum number of motifs set at 25 and optimum with of motifs set from 5 to 200 amino acids.
Promoter analysis and identification of miRNAs targets
The upstream 2 kb genomic DNA sequences of each predicted AR2/ERF genes were extracted from the B. distachyon genome, and then submitted to PLACE database (http://www.dna.affrc.go. jp/PLACE/) to identify the putative cis-regulatory elements in the promoter regions. Furthermore, all the identified AP2/ERF transcripts were searched against the published B. distachyon miRNAs in the miRBase using psRNATarget tool (http://plantgrn.noble.org/psRNATarget/) to predict the AR2/ERF targeted by miRNA.
Gene duplication and synteny analysis
Gene duplication events were identified manually using the method as described by Chen et al. . The segmental duplication events were characterized as copying the whole blocks of genes from one chromosome to another, while contiguous homologous genes with the original duplication on a single chromosome were defined as tandem duplications . For synteny analysis, duplications between B. distachyon AP2/ERF genes, as well as the synteny block of this family among B. distachyon and other 5 grass species (rice, maize, sorghum, foxtail millet and switchgrass) were obtained from the Plant Genome Duplication Database (http://chibba.pgml.uga.edu/duplication/) and the diagrams were visualized using the program Circos v0.67 .
Gene expression and network interaction analysis
Microarray data of B. distachyon were obtained from Gene Expression Omnibus (GEO) (http://www.ncbi.nlm.nih.gov/geo/) and EBI ArrayExpress (https://www.ebi.ac.uk/arrayexpress/) databases, and then used to detect the expression of the AR2/ERFs in different tissue and in response to abiotic stresses. Additionally, high throughput RNA sequencing data were also retrieved and downloaded from the SRA database (http://www.ncbi.nlm.nih.gov/sra) and then used to detect the differential expression of the AR2/ERF genes by FPKM analysis. A total of 9 RNA data of different tissues at different development stages were used, including anther, pistil, leaves (20 days), seed (5 and 10 days after pollination), endosperm(25 days after pollination), embryo(25 days after pollination), and inflorescence (early and emerging time). Finally, the interaction network which these putative AR2/EFR genes involved in were investigated based on the orthogous genes between B. distachyon and Arabidopsis using the AraNet V2 tool (http://www.inetbio.org/aranet/) .
Plant growth, stress treatment and RT-PCR analysis
Roots, stems, leaves and spikes were collected from two-months-old Bd21 genotype for RNA extraction and then used for organ-specific expression analysis. The 3 weeks old seedling were subjected to 4 °C, 20 % PEG, 150 mM NaCl conditions as cold,drought and salt treatments. After 24 hours treatment, the leaves of plant under these 3 stresses were collected for RNA isolation, respectively. Total RNA was isolated using RNAiso Reagent (TaKaRa, Dalian, China) according to the manufacturer’s instructions. Semi-quantitative RT-PCR was employed to determine the transcript levels of 11 randomly selected BdAP2/ERF genes following the method as described by Chen et al. . The primers are listed in Additional file 1: Table S1.
Identification of AP2/ERF family in Brachypodium
Characteristic features of AP2/ERF Transcription factor gene family identified in B. distachyon
Properties of AP2/ERF proteins
Start position (bp)
End Position (bp)
Protein length (aa)
Molecular weight (Da)
Summary of the abundance of each group of the AP2/ERF superfamily in B. distachyon, Arabidopsis and rice
a single group
Total AP2/ERF genes
genome size (Mbp)
The average number of AP2/ERF family genes per Mb (gene/MB)
The percentage of AP2/ERF family genes (%)
Chromosome distribution analysis found that the BdAP2/ERF genes were unevenly distributed on all of the five chromosomes of Brachypodium. In detail, 40 AP2/ERF genes located on the chromosome 3, representing the most abundant regions, followed by the chromosome 1, 2 and 5, with the number of 36, 32 and 18 respectively, while there were only 15 genes on the chromosome 4, which have the minimum number of AP2/ERFs. Interestingly, all the 4 RAV genes located on the chromosome 2, which may be a Brachypodium-specific feature. The putative proteins of BdAP2/ERFs ranged from 92 to 1338 amino acids in length, with molecular weights (Mw) ranging from 9.8 to 148.6 kDa and theoretical isoelectric points (PI) ranging from 4.33 to 11.63. Subcellular localization analysis indicated that majority of BdAP2/ERFs (108 out of 141, 76.5 %) localized in the nucleus, while 25 genes were predicted to be located in the chloroplast and the remaining 7 genes located in cytoplasmic, mitochondrial, plasma membrane and extra-celluar (Table 1). To further assess the actual existence of these genes identified in this study, all the available Brachypodium expressed sequence tags (EST) were used to search against these genes using the BlastN program. Results showed that most of the AP2/ERFs were supported by EST hits, only 36 genes (25.5 %, 36/111) showed no EST hits. In light of the limit of available ESTs, the not-supported BdAP2/ERF gene might not express under any the used conditions or express with very low level that cannot be detected experimentally.
Phylogenetic relationship, conserved motif and gene structure analysis
Cis-elements and miRNA targets analysis
In order to understand the possible biological functions and regulation network of these AP2/BdERFs involved in, 2 kb genomic sequences upstream of the 5′-UTR of BdAP2/ERF genes were extracted and used to identify cis-regulatory elements. A total of 276 putative cis-elements were found to be presented in at least one BdAP2/ERF gene and only 7 (GT1CONSENSUS, DOFCOREZM, EBOXBNNAPA, MYCCONSENSUSAT, CAATBOX1, CACTFTPPCA1, WRKY71OS) out of them were presented in the promoter region of all BdAP2/ERF genes (Additional file 1: Table S4). In addition, 32 cis-elements were detected as gene-specific, such as S2FSORPL21, ABREDISTBBNNAPA and ABREDISTBBNNAPA were unique to Bradi5g24360, Bradi3g58980 and Bradi5g17620, respectively. The different numbers and types of cis-elements presenting in BdAP2/ERF genes indicated the differential regulatory network which the BdAP2/ERF genes may involve in. Further analysis found that hormones-response (e.g. abscisic acid, gibberellins, auxin, jasmonic acid and ethylene), abiotic stress-related (e.g., drought, extreme temperatures, high salinity, wounding, and disease) and organogenesis-related cis-elements were abundantly presented in the promoter regions of BdAP2/ERF (Additional file 1: Table S5), which indicated that these AP2/ERF genes might have potential functions involving in regulating a variety of stresses response and hormone signaling transduction.
Furthermore, the putative microRNAs (miRNAs) targeted BdAP2/ERF genes were also detected in this study and a total of 8 BdAP2/ERFs were predicted to be targeted by seven miRNAs (Additional file 1: Table S6). Although miRNA inhibition mostly involved the transcript cleavage, the BdAP2/ERF006 was predicted to be inhibited to translation. Most predicted microRNA target sites located into CDS region but outside the AP2 domain, whereas for gene BdAP2/ERF051 the cleavage site located in the 3′UTR region. The miRNAs-AP2/ERF complex identified in this study would be useful in interpreting the post-transcriptional control of gene expression during various stress-induced physiological and cellular processes in B. distachyon as well as other cereal crops.
Gene duplication and synteny analyses of AP2/ERFs between B. distachyon and other three grass species
The substitution rate of non-synonymous (Ka) versus synonymous (Ks) was an effective measure to examine the positive selection pressure after duplication, wherein Ka/Ks =1 means neutral selection, Ka/Ks <1 stands for purifying selection, and Ka/Ks >1 signifies accelerated evolution with positive selection . Furthermore, the divergence rate of the tandem and segmental duplicated BdAP2/ERF genes was calculated to detect selection influence (Additional file 1: Table S7 and S8). The Ka/Ks ratio for tandem duplicated gene-pairs in B. distachyon AP2/ERF genes ranged from 0.23 to 0.51 with an average of 0.31, whereas Ka/Ks for segmental duplicated gene-pairs ranged from 0.19 to 0.85 with an average of 0.53. These results indicated that the duplicated BdAP2/ERF genes were under strong purifying selection pressure and had gone through substitution elimination and enormous selective constraint by natural selection during the process of evolution since their Ka/Ks ratios were estimated to be lower than one. In addition, the duplication event of these BdAP2/ERF tandem and segmental duplicated genes was estimated to have occurred around ~54 and ~61 Mya, respectively. Although the BdAP2/ERF gene-pairs of segmental (Ka/Ks = 0.53) and tandem duplication (Ka/Ks = 0.31) events are not under similar evolutionary positive selection pressure, both set of gene pairs revealed that these duplication events may take place simultaneously. Additionally, the Ka/Ks ratios of the orthologous gene-pairs between B. distachyon and other three grass species were also calculated (Additional file 1: Table S9, S10, S11). The average Ka/Ks value was maximum between B. distachyon and maize (0.47), followed by rice (0.44) and sorghum(0.43), suggesting the genes pairs between B. distachyon and those three grass species appeared to have undergone extensive intense purifying selection. The divergence time was about 47, 49 and 51Mya for rice, sorghum and maize, respectively. Therefore, it can be concluded that the segmental and tandem duplication events played a major role in evolution and functional divergent of AP2/ERF genes family in B. distachyon as well as other grass species.
Co-expression network between AP2/ERFs and other genes in B. distachyon
Expression profiles of BdAP2/ERF genes at different developmental stages and under stresses
Expression patterns of BdAP2/ERFs in various tissues and under stress treatment by semi-quantitative RT-PCR analysis
AP2/ERF superfamily is one of the largest groups of plant-specific transcription factors, which has been widely studied in diverse plant species, such as Arabidopsis, soybean, rice, maize, foxtail millet and switchgrass [1, 28, 30, 43, 44]. In this study, we performed a comprehensive search for AP2/ERF genes throughout Brachypodium genome, and 141 BdAP2/ERF genes were found, accounting for 0.45 % of all the Brachypodium genes, which was similar with the result in rice (0.43 %), maize (0.44 %) and foxtail millet (0.44 %) . While compared to other plants, the number of AP2/ERF in Brachpodium is much lower than that of rice (174), maize (184) as well as foxtail millet (171), and also slight lower than that of Arabidopsis (148) and grape (149). It has been revealed that the number of AP2/ERF gene family was mainly depending on the number of ERF family members . It’s found that there are 112 members in ERF family in Brachypodium, while 122, 132 and 158 in Arabidopsis, rice and maize, respectively. In contrast, the number of AP2 and RAV family members showed no significantly difference among them, with the value of 28, 24, 34 and 25 in Brachypodium, Arabidopsis, rice and maize respectively. Thus, the lower AP2/ERF gene abundance in Brachpodium may also due to the lower number of ERF and DREB subfamily. Furthermore, the gene density is 0.3972 AP2/ERF genes per Mb in B. distachyon, while the value for rice and Arabidopsis is 0.4047 and 1.1760 respectively. B. distachyon shows closer AP2/ERF density with rice than Arabidopsis, suggesting the specific evolutionary events might occur to regulate the retention and disposition of this gene family between Monocots and Eudicots.
It has been widely revealed that AP2/ERF transcription factors played crucial roles in regulating plant growth, development and response to diverse stresses as well as signal transduction pathway in plants . However, the function of BdAP2/ERFs is not well understood at present. In this study, the expression patterns of these genes in different tissues and under different stresses were systematically investigated to understand their potential function during development and stress response. Results found that a total of 138 BdAP2/ERFs were expressed in at least one tested tissue, indicating they widely involved in growth and development. Compared to AP2 family, the members in ERF family showed higher expression levels in these tissues. We found that the ERF family genes had less intron than AP2 family in Brachypodium, which may cause the quicker response and higher expression of ERF genes during development . At the same time, BdAP2/ERF genes also showed obvious spatial and temporal expression profiles. For example, BdAP2/ERF064 is specifically expressed in leaf, and BdAP2/ERF005 showed preferential expression in the emerging inflorescence. In addition, six genes having significantly higher expressed during pollination were also identified, which may play the vital role in embryo and endosperm development. AP2/ERF proteins could bind to GCC-box or DRE motifs through the ERF domain, and then regulated the target gene expression under stress conditions [47, 48]. Compare to control, 69 BdAP2/ERFs showed differential expression under cold stress and 16 showed differential expression under drought stress, respectively. Among them, BdAP2/ERF120 (Bradi4g35630) which is a member of DREB subfamily, showed significantly up-regulated under both cold and drought stresses. Previous study have reported Bradi4g35630 encoding a C-repeat binding factor 3-like protein, is a cold-responsive gene and over-expression of this gene could improve the drought, salt and cold tolerance in Brachypodium . Moreover, a total of 5 ACGTATERD1 (element of early responsive to dehydration), 1 DRE1COREZMRAB17 (element of responsive to drought) and 9 MYCCONSENSUSAT (element of responsive to dehydration and cold) cis-elements were identified in the promoter region of BdAP2/ERF120. In addition, BdAP2/ERF053 (Bradi2g27920) is found to be highly expressed in all three stress treatments, which contained 5 LTRECOREATCOR15 (core element of low temperature responsive), 6 EMBP1TAEM (element involving in ABA-mediated stress-signaling pathway) and 1 GT1GMSCAM4 (element required for salt-induced gene expression) cis-elements. We speculated that cis-elements were the vital regulators to control the spatial and temporal expression of the BdAP2/ERFs, which integrated other functional proteins with the AP2/ERF transcription factor to form the complex regulatory metabolic network during development and stress response processes . These identified tissue-specific and stress-induced BdAP2/ERF provided the valuable candidates for further functional studies of AP2/ERF genes in B. distachyon as well as in other cereal crops.
Our current study identified and characterized the AP2/ERF transcription factors in the model grass B. distachyon. By performing a genome-wide search, a total of 141 BdAP2/ERF genes were obtained. EST hits or full-length cDNA sequences confirmed their actual existence. The chromosome location, exon-intron structure and conserved motif composition as well as phylogenetic relationship of these BdAP2/ERFs were systematically analyzed and compared. BdAP2/ERFs could be classified into four subgroups in accordance with the number of AP2 domains and putative functions. Co-expression network analysis found that 39 BdAP2/ERFs were involved in regulating other B. distachyon genes, and 517 network branches were found. The expression profiles of BdAP2/ERF genes in various tissues as well as under cold and drought stresses were investigated, and several tissue-specific or stress-induced BdAP2/ERF genes were identified, which could considered as the candidates for further study of their function in plant development and stress response. Our study for the first time reported the organization, structure, evolutionary and expression features of the BdAP2/ERF family, which will facilitate the future functional analysis of BdAP2/ERF genes, and lay the foundation for better understanding the molecular mechanism of plant development and stress physiological processes in B. distachyon and beyond.
AP2/ERF, APETALA2/Ethylene responsive factor; DREB, dehydration responsive element binding protein; EST, expressed sequence tag; FPKM, fragments kilobase of exon model per millon mapped reads; Ka, substitution rate of non-synonymous; Ks, substitution rate of synonymous; MP, maximum parsimony; MW, molecular weight; NJ, neighbor joining; PI, isoelectric point; TF, transcription factors
This work was mainly funded by the National Natural Science Foundation of China (Grant No. 31401373) and partially supported by the Open Project Program of State Key Laboratory of Crop Stress Biology in Arid Areas, China (CSBAA2014002). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. We are grateful to Dr. Haifeng Li for providing the plant materials, and also grateful to Jianxin Bian and Le Wang for their constructive suggestion on data analysis as well as we would like to thank the anonymous reviewers for their constructive comments.
Availability of data and material
All of the datasets supporting the results of this article are included within the article and its Additional files. The deposited the phylogenetic data in our manuscript to Treebase to make it available publically with the accession NO. S19605. The access URL is http://purl.org/phylo/treebase/phylows/study/TB2:S19605.
CLC and FKW collected the public dataset, perform bioinformatics analysis and also drafted the manuscript. WMX contributed to data collection. WM and DPC contributed to data analysis and manuscript preparation. SWN provided the experimental coordination and reviewed the manuscript. NXJ conceived this study and prepared the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Sakuma Y, Liu Q, Dubouzet JG, Abe H, Shinozaki K, Yamaguchi-Shinozaki K. DNA-binding specificity of the ERF/AP2 domain of Arabidopsis DREBs, transcription factors involved in dehydration-and cold-inducible gene expression. Biochem Biophys Res Commun. 2002;290:998–1009.View ArticlePubMedGoogle Scholar
- Wessler SR. Homing into the origin of the AP2 DNA binding domain. Trends Plant Sci. 2005;10:54–6.View ArticlePubMedGoogle Scholar
- Nakano T, Suzuki K, Fujimura T, Shinshi H. Genome-wide analysis of the ERF gene family in Arabidopsis and rice. Plant Physiol. 2006;140:411–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Jofuku KD, Den Boer B, Van Montagu M, Okamuro JK. Control of Arabidopsis flower and seed development by the homeotic gene APETALA2. Plant Cell. 1994;6:1211–25.View ArticlePubMedPubMed CentralGoogle Scholar
- Elliott RC, Betzner AS, Huttner E. AINTEGUMENTA, an APETALA2-like gene of Arabidopsis with pleiotropic roles in ovule development and floral organ growth. Plant Cell. 1996;8:155–68.View ArticlePubMedPubMed CentralGoogle Scholar
- Ohme-Takagi M, Shinshi H. Ethylene-inducible DNA binding proteins that interact with an ethylene-responsive element. Plant Cell. 1995;7:173–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Hao D, Ohme-Takagi M, Sarai A. Unique mode of GCC box recognition by the DNA-binding domain of ethylene-responsive element-binding factor (ERF domain) in plant. J Biol Chem. 1998;273:26857–61.View ArticlePubMedGoogle Scholar
- Kagaya Y, Ohmiya K, Hattori T. RAV1, a novel DNA-binding protein, binds to bipartite recognition sequence through two distinct DNA-binding domains uniquely found in higher plants. Nucleic Acids Res. 1999;27:470–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Krizek B. AINTEGUMENTA and AINTEGUMENTA-LIKE6 act redundantly to regulate Arabidopsis floral growth and patterning. Plant Physiol. 2009;150:1916–29.View ArticlePubMedPubMed CentralGoogle Scholar
- Xu ZS, Chen M, Li LC, Ma YZ. Functions and application of the AP2/ERF transcription factor family in crop improvement. J Integr Plant Biol. 2011;53(7):570–85.View ArticlePubMedGoogle Scholar
- Gutterson N, Reuber TL. Regulation of disease resistance pathways by AP2/ERF transcription factors. Curr Opin Plant Biol. 2004;7(4):465–71.View ArticlePubMedGoogle Scholar
- Aukerman MJ, Sakai H. Regulation of flowering time and floral organ identity by a microRNA and its APETALA2-like target genes. Plant Cell. 2003;15:2730–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Taketa S, Amano S, Tsujino Y, et al. Barley grain with adhering hulls is controlled by an ERF family transcription factor gene regulating a lipid biosynthesis pathway. Proc Natl Acad Sci U S A. 2008;105:4062–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Jofuku KD, Omidyar PK, Gee Z, Okamuro JK. Control of seed mass and seed yield by the floral homeotic gene APETALA2. Proc Natl Acad Sci U S A. 2005;102:3117–22.View ArticlePubMedPubMed CentralGoogle Scholar
- Alonso JM, Stepanova AN, Leisse TJ, et al. Genome-wide insertional mutagenesis of Arabidopsis thaliana. Science. 2003;301:653–57.View ArticlePubMedGoogle Scholar
- Hu YX, Wang YH, Liu XF, Li JY. Arabidopsis RAV1 is down-regulated by brassinosteroid and may act as a negative regulator during plant development. Cell Res. 2004;14:8–15.View ArticlePubMedGoogle Scholar
- Sohn KH, Lee SC, Jung HW, Hong JK, Hwang BK. Expression and functional roles of the pepper pathogen-induced transcription factor RAV1 in bacterial disease resistance, and drought and salt stress tolerance. Plant Mol Biol. 2006;61:897–915.View ArticlePubMedGoogle Scholar
- Li CW, Su RC, Cheng CP, You SJ, Hsieh TH, Chao TC, Chan MT. Tomato RAV transcription factor is a pivotal modulator involved in the AP2/EREBP-mediated defense pathway. Plant Physiol. 2011;156:213–27.View ArticlePubMedPubMed CentralGoogle Scholar
- Ito Y, Katsura K, Maruyama K, Taji T, Kobayashi M, Seki M, et al. Functional analysis of rice DREB1/CBF-type transcription factors involved in cold-responsive gene expression in transgenic rice. Plant Cell Physiol. 2006;47:141–53.View ArticlePubMedGoogle Scholar
- Jaglo-Ottosen KR, Gilmour SJ, Zarka DG, Schabenberger O, Thomashow MF. Arabidopsis CBF1 overexpression induces COR genes and enhances freezing tolerance. Science. 1998;280:104–6.View ArticlePubMedGoogle Scholar
- Qin F, Kakimoto M, Sakuma Y, et al. Regulation and functional analysis of ZmDREB2A in response to drought and heat stresses in Zea mays L. Plant J. 2007;50:54–69.View ArticlePubMedGoogle Scholar
- Hong JP, Kim WT. Isolation and functional characterization of the Ca-DREBLP1 gene encoding a dehydration-responsive element binding-factor-like protein 1 in hot pepper (Capsicum annuum L. cv. Pukang). Planta. 2005;220:875–88.View ArticlePubMedGoogle Scholar
- Draper J, Mur LA, Jenkins G, Ghosh-Biswas GC, Bablak P, Hasterok R, Routledge AP. Brachypodium distachyon. A new model system for functional genomics in grasses. Plant Physiol. 2001;127:1539–55.View ArticlePubMedPubMed CentralGoogle Scholar
- Brkljacic J, Grotewold E, Scholl R, Mockler T, Garvin DF, Vain P, et al. Brachypodium as a model for the grasses: today and the future. Plant Physiol. 2011;57:3–13.View ArticleGoogle Scholar
- Vogel JP, Garvin DF, Mockler TC, et al. Genome sequencing and analysis of the model grass Brachypodium distachyon. Nature. 2010;463:763–8.View ArticleGoogle Scholar
- Wu H, Lv H, Li L, Liu J, Mu S, Li X, Gao J. Genome-wide analysis of the AP2/ERF transcription factors family and the expression patterns of DREB genes in Moso Bamboo (Phyllostachys edulis). PLoS One. 2015;5:e0126657.View ArticleGoogle Scholar
- Licausi F, Giorgi FM, Zenoni S, Osti F, Pezzotti M, Perata P. Genomic and transcriptomic analysis of the AP2/ERF superfamily in Vitis vinifera. BMC Genomics. 2010;11:719.View ArticlePubMedPubMed CentralGoogle Scholar
- Du H, Huang M, Zhang Z, Cheng S. Genome-wide analysis of the AP2/ERF gene family in maize waterlogging stress response. Euphytica. 2014;198:115–26.View ArticleGoogle Scholar
- Zhang C, Shangguan LF, Ma RJ, et al. Genome-wide analysis of the AP2/ERF superfamily in peach (Prunus persica). Genet Mol Res. 2012;11:4789–809.PubMedGoogle Scholar
- Rashid M, Guangyuan H, Guangxiao Y, Hussain J, Xu Y. AP2/ERF transcription factor in rice: genome-wide canvas and syntenic relationships between monocots and eudicots. Evol Bioinformatics Online. 2012;8:321.View ArticleGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, et al. Clustal W and clustal X version 2.0. Bioinformatics. 2007;23:2947–48.View ArticlePubMedGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–29.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen L, Hu W, Tan S, et al. Genome-wide identification and analysis of MAPK and MAPKK gene families in B. distachyon. PLoS One. 2012;7:e46744.View ArticlePubMedPubMed CentralGoogle Scholar
- Leister D. Tandem and segmental gene duplication and recombination in the evolution of plant disease resistance genes. Trends Genet. 2014;20:116–22.View ArticleGoogle Scholar
- Krzywinski M, Schein J, Birol I, et al. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19:1639–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee T, Yang S, Kim E, et al. AraNet v2: an improved database of co-functional gene networks for the study of Arabidopsis thaliana and 27 other nonmodel plant species. Nucleic Acids Res. 2015;43:D996–1002.View ArticlePubMedGoogle Scholar
- Mochida K, Yoshida T, Sakurai T, Yamaguchi-Shinozaki K, Shinozaki K, Tran LP. In Silico analysis of transcription factor repertoires and prediction of stress-responsive transcription factors from six major Gramineae plants. DNA Res. 2011;18(5):321–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu ZJ, Li XH, Liu ZW, Li H, Wang YX, Zhuang J. Transcriptome-based discovery of AP2/ERF transcription factors related to temperature stress in tea plant (Camellia sinensis). Funct Integr Genomics. 2015;15:741–52.View ArticlePubMedGoogle Scholar
- Liu L, White MJ, MacRae TH. Transcription factors and their genes in higher plants. Eur J Biochem. 1999;262:247–57.View ArticlePubMedGoogle Scholar
- Lynch M, Conery JS. The evolutionary fate and consequences of duplicate genes. Science. 2000;290:1151–5.View ArticlePubMedGoogle Scholar
- Huan Q, Mao Z, Zhang J, Xu Y. Chong K. Transcriptome-wide analysis of vernalization reveals conserved and species-specific mechanisms in Brachypodium distachyon. J Integr Plant Biol. 2013;55:696–709.Google Scholar
- Verelst W, Bertolini E, De Bodt S, Vandepoele K, Demeulenaere M, Pè ME, Inzé D. Molecular and physiological analysis of growth-limiting drought stress in Brachypodium distachyon leaves. Mol Plant. 2013;6:311–22.View ArticlePubMedGoogle Scholar
- Lata C, Mishra AK, Muthamilarasan M, Bonthala VS, Khan Y, Prasad M. Genome-wide investigation and expression profiling of AP2/ERF transcription factor superfamily in foxtail millet (Setaria italica L.). PLoS One. 2014;11:e113092.View ArticleGoogle Scholar
- Wuddineh WA, Mazarei M, Turner GB, et al. Identification and molecular characterization of the switchgrass AP2/ERF transcription factor superfamily, and overexpression of PvERF001 for improvement of biomass characteristics for biofuel. Front Bioeng Biotechnol. 2015;3:101.View ArticlePubMedPubMed CentralGoogle Scholar
- Shu Y, Liu Y, Zhang J, Song L, Guo C. Genome-wide analysis of the AP2/ERF superfamily genes and their responses to abiotic stress in Medicago truncatula. Front Plant Sci. 2016;6:1247.View ArticlePubMedPubMed CentralGoogle Scholar
- Mizoi J, Shinozaki K, Yamaguchi-Shinozaki K. AP2/ERF family transcription factors in plant abiotic stress responses. Biochim Biophys Acta. 2012;1819(2):86–96.View ArticlePubMedGoogle Scholar
- Shigyo M, Hasebe M, Ito M. Molecular evolution of the AP subfamily. Gene. 2006;366(2):256–65.View ArticlePubMedGoogle Scholar
- Lata C, Prasad M. Role of DREBs in regulation of abiotic stress responses in plants. J Exp Bot. 2011. doi:10.1093/jxb/err210.PubMed CentralGoogle Scholar
- Ryu JY, Hong SY, Jo SH, et al. Molecular and functional characterization of cold-responsive C-repeat binding factors from Brachypodium distachyon. BMC Plant Biol. 2014;14:15.View ArticlePubMedPubMed CentralGoogle Scholar
- Dietz KJ, Vogel MO, Viehhauser A. AP2/EREBP transcription factors are part of gene regulatory networks and integrate metabolic, hormonal and environmental signals in stress acclimation and retrograde signalling. Protoplasma. 2010;245:3–14.View ArticlePubMedGoogle Scholar