Genomic and transcriptomic analysis of the AP2/ERF superfamily in Vitis vinifera
© Licausi et al; licensee BioMed Central Ltd. 2010
Received: 5 August 2010
Accepted: 20 December 2010
Published: 20 December 2010
The AP2/ERF protein family contains transcription factors that play a crucial role in plant growth and development and in response to biotic and abiotic stress conditions in plants. Grapevine (Vitis vinifera) is the only woody crop whose genome has been fully sequenced. So far, no detailed expression profile of AP2/ERF-like genes is available for grapevine.
An exhaustive search for AP2/ERF genes was carried out on the Vitis vinifera genome and their expression profile was analyzed by Real-Time quantitative PCR (qRT-PCR) in different vegetative and reproductive tissues and under two different ripening stages.
One hundred and forty nine sequences, containing at least one ERF domain, were identified. Specific clusters within the AP2 and ERF families showed conserved expression patterns reminiscent of other species and grapevine specific trends related to berry ripening. Moreover, putative targets of group IX ERFs were identified by co-expression and protein similarity comparisons.
The grapevine genome contains an amount of AP2/ERF genes comparable to that of other dicot species analyzed so far. We observed an increase in the size of specific groups within the ERF family, probably due to recent duplication events. Expression analyses in different aerial tissues display common features previously described in other plant systems and introduce possible new roles for members of some ERF groups during fruit ripening. The presented analysis of AP2/ERF genes in grapevine provides the bases for studying the molecular regulation of berry development and the ripening process.
The AP2/ERF superfamily is one of the largest groups of transcription factors in plants . It includes all genes coding for at least one APETALA2 (AP2) domain and can be further separated into the ETHYLENE RESPONSE FACTOR (ERF), the AP2, and the RAV families. The AP2 domain, which is responsible for DNA binding, was named after the Arabidopsis thaliana APETALA2 protein identified by Jofuku et al. . Subsequently, tobacco (Nicotiana tabacum) transcription factors containing a motif related to the AP2 domain were also identified . Since these proteins are able to bind an ethylene responsive DNA element (AGCCGCC), they were classified as ERFs (Ethylene Responsive Factors) and their DNA binding domain was named after them .
The subdivision of the AP2/ERF group into families is based on the number of AP2 domains present in the proteins, together with the presence of other DNA binding domains. The AP2 family contains proteins with a double, tandem repeated, AP2 domain . The ERF family genes code for transcription factors with a single AP2 domain. Lastly, the RAV family encodes proteins possessing a single AP2 domain plus an additional B3 domain, which is also present in other, non-ERF transcription factors .
Two major efforts to define a nomenclature for ERF sequences have been published. Sakuma et al. divided the Arabidopsis ERF family into two subfamilies based on the amino acid sequence of the DNA binding domain: the CBF/DREB subfamily (group A) and the ERF-like subfamily (group B), both further subdivided in six groups. Nakano et al.  on the other hand, proposed an innovative classification for Arabidopsis ERFs, based on the amino acid sequence of the whole proteins, splitting the family into ten groups.
In the AP2 family, a binding element gCAC(A/G)N(A/T)TcCC(a/g)ANG(c/t) has been reported for the protein AINTEGUMENTA (ANT) , and the same sequence was shown to be bound by five other AP2 members . This DNA element is quite distinct from the consensus sequence CCGA/CC bound by the DREB/ERF group and, surprisingly, is not composed of two similar half sites, since the AP2 proteins contain a double AP2 domain .
After the release of the whole genomic sequences of several plant organisms, including Arabidopsis, rice (Oryza sativa) and poplar (Populus tricocarpa), the AP2/ERF transcription factor superfamily was analyzed, both to place each member in an organized nomenclature system, and to provide maps of their expression. To date, two full genome sequences of highly homozygous and heterozygous grapevine (Vitis vinifera) Pinot noir genotypes have been carried out [[20, 21] and ]. These two milestones provided a useful genomic platform to study this fruit crop. As AP2/ERF transcription factors are involved in flower development and tolerance to biotic and abiotic stresses, their superfamily represents one of the best pools to investigate, when searching for important grapevine traits. Since genome scale analyses of the transcriptional response to the development and environmental stimuli require a precise and complete annotation in order to provide reliable and exhaustive data, we decided to annotate the ERF family members and create a qRT-PCR platform that allows investigating their expression profile.
A recent study suggested that 132 genes encoding AP2/ERF proteins are present in the grapevine genome . However, a higher number of AP2/ERF genes is present in Arabidopsis and poplar genomes (147 and 202 sequences, respectively). We therefore re-screened the grapevine genome for AP2/ERF sequences, adopting two different strategies, in order to accurately identify AP2/ERF-like sequences in the Vitis vinifera genome. In this study we provide a characterization of the grapevine AP2/ERF transcription factor superfamily and demonstrate that it is in fact composed of almost 149 genes. Moreover, using qRT-PCR platform encompassing the whole ERF/AP2 superfamily, we show how AP2/ERF-like genes are expressed in both vegetative and reproductive tissues at different developmental stages, and we infer roles and putative targets for some of these genes. Overall, our analysis suggests that AP2/ERF proteins play a strong role in ripening-related processes.
Identification of the AP2/ERF family transcription factors in Vitis vinifera
As result of an extensive search for AP2-domain containing proteins, 149 distinct AP2/ERF putative TFs were identified (Additional file 1, Table S1). One hundred and twenty two genes encoding for proteins with a single AP2/ERF domain were assigned to the ERF superfamily. The AP2 family was grouped into 20 genes, that could be identified due to the tandem repeated double AP2/ERF motif. Six genes, containing a single AP2/ERF DNA binding domain together with a B3 type domain, were assigned to the RAV family. A single AP2 protein, GIDVvP00018355001, showed a low similarity to the other ERF sequences although it is homologous to the Arabidopsis ERF transcription factor At4g13040. A similar gene has also been identified in P. trichocarpa and named Soloist . Surprisingly, the Vitis gene annotated as Soloist in the grapevine genome by Zhuang et al.  is not homologous to the poplar Soloist. Zhuang et al.  named protein sequence GSVIVP00025602001 as Soloist, although its similarity to the sequences of Arabidopsis At4g13040 and PtSoloist is rather low. On the other hand, GIDVvP00018355001 is the closest homolog of both At4g13040 and PtSoloist in the grapevine genome (Additional file 2, Figure S2).
Previous annotations of AP2/ERF genes in poplar  and grapevine  followed the nomenclature proposed by Sakuma et al. , based on a homology of the DNA binding domain alone. However, Nakano et al.  proposed an alternative method, based on the presence of domains that were different from the DNA binding domain. Therefore we subdivided the grapevine ERF genes into 11 groups, according to their similarity to the Arabidopsis ERF sequences. Amino acid motifs located outside the DNA binding domain are conserved among the Arabidopsis and Vitis ERF proteins (Additional file 3, Table S3).
A comparison of AP2/ERF families and groups between monocot (Oryza sativa) and dicot (Arabidopsis thaliana, Populus trichocarpa, Vitis vinifera) species.
Vitis / Arabidopsis
Vitis / Poplar
Vitis / Rice
Accuracy of protein predictions
Proteins sequences from the Genoscope database have been deduced using an automated method . However, this methodology is prone to errors . We found similar issues in grapevine predictions. In fact, some gene annotations and intron junctions did not match those of the Arabidopsis and poplar homologs. We found unrealistic introns of 8.4 and 35 Kb predicted in GSVIVP00009519001 and in GSVIVT0009456001 gene models respectively. Moreover the gene models GSVIVT00021812001, GSVIVP00013482001, GSVIVT00019482001, and GSVIVP00007524001 were predicted to encode proteins that contained two or more ERF domains. Since both of these sequences were homologous to Arabidopsis proteins belonging to the ERF family and not to the AP2 family, we decided to split each gene model into sequences encoding a single ERF domain protein, using the closest Arabidopsis protein as a template (GSVIVT00021812001: VvERF014 and VvERF016; GSVIVP00013482001: VvERF029 and VvERF030; GSVIVT00019482001: VvERF115, VvERF116 and VvERF119; GSVIVP00007524001: VvERF101 and VvERF102). In order to test whether the prediction of these dubious genes was correct, we compared the production of amplicons corresponding to single exons with that of the predicted cDNA derived from exon-joining. The absence of an amplicon that corresponded to exon joining, together with the correct amplification of the single exon, was interpreted as an indication of misprediction. The gene model GSVIVT00034010001, which is supported by cDNA report evidence, was used as a positive control. No amplification was obtained using primers annealing to the N- and C- terminus of GSVIVT00009519001, GSVIVT0009456001 and GSVIVT00021812001, although an amplicon of the expected size was obtained using primers that anneal within the exons (Additional file 4, Figure S4). Since the present results revealed the possibility that other gene structure predictions were partly incorrect, all dubious sequences were checked and, where necessary, corrected, using Arabidopsis and poplar closest homologs as templates (Additional file 1, Table S1).
In order to investigate the evolutionary relatedness of the identified sequences, together with the ERF genes encoded by the other fully sequenced plant species, we performed a sequence-based phylogenetic analysis.
Chromosome position of the identified Vitis vinifera AP2/ERF genes
Expression analyses of the AP2/ERF genes in aerial tissues of V. vinifera cv. Corvina
A broad variability of expression patterns within families and groups was observed: for instance, in the ERF group V, VvERF042 and VvERF043 were highly expressed in leaf, inflorescence and stem tissues while VvERF046 and VvERF047 reached high levels in the flesh of ripe berries (Figure 3). Very few genes exhibited a specific single tissue expression: VvERF103 mRNA was barely detectable except in leaf tissues, VvERF114 was highly expressed in flesh tissues of the berry at the veraison stage and VvAP2-03 showed a strong induction in skin tissues during ripening (Figure 4).
Three ERF genes displayed a high expression in leaves, stems and inflorescences and low levels in the berry tissues considered: VvERF027 (group III of the ERF family), VvERF042 and VvERF043 (ERF group V) (Figure 3), together with a big group of AP2 genes (VvAP2-09, VvAP2-11, VvAP2-12, VvAP215, VvAP2-17, VvAP2-18 and VvAP2-19) (Figure 4). Surprisingly, none of the AP2-coding genes, which have been reported to play a role in flower development, were exclusively expressed in the inflorescence tissues (Figure 4).
Expression pattern similarity did not mimic the sequence similarity of ERF group IX. For instance, VvERF084 and VvERF085, which share more than 95% similarity, exhibited different expression patterns in fruit tissues and VvERF109 displayed stem specificity whereas its closest paralog, VvERF110, was poorly expressed in the same tissue.
Ripening regulated AP2/ERF genes in grape berry
AP2 and ERF related transcription factors are assumed to play a role in fruit development during the ripening process in several species [[28, 29] and ]. To identify AP2/ERF genes whose expression changes during the transition from veraison to full ripeness (stage III according to Deluc et al.), we analyzed the mRNA levels of 149 AP2/ERF genes identified in skin and flesh tissues.
Differentially regulated AP2/ERF genes during ripening in the skin tissues.
Genes up-regulated in the transition from veraison to ripe in the skin tissues
Genes down-regulated in the transition from veraison to ripe in the skin tissues
Fold Change (log2)
Fold Change (log2)
Differentially regulated AP2/ERF genes during ripening in the flesh tissues.
Genes up-regulated in the transition from veraison to ripe in the flesh tissues
Genes down-regulated in the transition from veraison to ripe in the flesh tissues
Fold Change (log2)
Fold Change (log2)
To test this hypothesis, we analyzed the expression of genes encoding for β-1,3 glucanases (VvGlucβ-1, VvGlucβ-2), class IV chitinase (VvChit4-1) and pathogenesis related proteins (VvPR10.1 and VvPR10L-b) using RT-qPCR. VvGlucβ-1 and VvGlucβ-2, VvChit4-1, VvPR10.1 and VvPR10L-b were induced in the transition from veraison to full ripeness in the skin of the berry (Additional file 8, Figure S8). The same trend was observed for several ERF-IX genes including VvERF086, VvERF093 and VvERF098 (Figure 3 and Additional file 8, Figure S8). Moreover, the expression profile of ERF-IX VvERF078 and VvERF085 is highly correlated with that of genes encoding chitinases, β 1,3-glucanases and thaumatin during different stages of berry development in the dataset published by Deluc et al.  (Additional file 9, Figure S9).
Only the VII group, which contains three members, was not represented in any of the DEG lists (Figure 3, Tables 2, 3). Four AP2 genes were induced and two were repressed in the skin during the ripening phase considered, whereas in the flesh only one gene was up-regulated and two were down-regulated (Tables 2, 3).
Comparison of the AP2/ERF superfamily across sequenced plant genomes
In this study, a comprehensive search for genes encoding for AP2/ERF transcription factors encoded in the Vitis vinifera genome was carried out, leading to the identification of 149 genes. Previously, 132 AP2/ERF genes had been annotated, as reported by Zhuang et al.. Eighteen newly identified genes encoded transcription factors belonging to groups II (VvERF013), III (VvERF028 and VvERF029), V (VvERF041), VIII (VvERF060, VvERF061, VvERF062, VvERF063, VvERF071, VvERF073 and VvERF074), IX (VvERF094) and X (VvERF118 and VvERF121) of the ERF family, to the AP2 family (VvAP2-12 and VvAP2-17) and the RAV family (VvRAV3 and VvRAV5) (Additional file 1, Table S1). Our analysis did not identify GSVIVP00019513001 as an ERF gene due to its short sequence and the presence of two long introns. It is likely that this locus constitutes a pseudogene. This number of genes is comparable with that of Arabidopsis but is only 70% of the number found in poplar and 80% of the rice AP2/ERF superfamily (Table 1).
The availability of the complete genome sequence of three dicot species and one monocot enabled a comparison of individual families and groups using Nakano et al.'s criteria . Several ERF groups (II, III, VI, VII, VIII, IX, X, VI-L and Xb-L) and the AP2 family contained more members in poplar than in grape and Arabidopsis (Table 1). A feature common to both woody species Vitis vinifera and Populus trichocarpa is the lower number of genes in group I in the ERF family and the higher number in groups V and IX. In Medicago truncatula, WXP1, a group I ERF, contributes to drought tolerance by regulating the accumulation of cuticular waxes  while, in Arabidopsis, this function is performed by the SHYNE clade of ERF genes (SHN1-3), all belonging to group V . Therefore, it is possible that the function of group I and V overlaps and the increase in group V in woody dicot species poises the lower number of group I genes. The expression of homologs of the Arabidopsis SHN clade, VvERF043 and VvERF044, in fruit skin tissues (Figure 3) is not unexpected, since waxy layers are produced to cover the grape berries .
Compared to all the species considered in this study, the grapevine genome has very small ERF-group II and VII. In Arabidopsis, the ERF-VII genes are highly expressed, suggesting to be regulated at the post-transcriptional level. Interestingly, also VvERF057 and VvERF059 are highly expressed in all the tissues analysed in this study.
ERF-IX is the largest ERF group in poplar and grape genomes, since it is almost twice as big as Arabidopsis and rice (Table 1). Of the 40 ERF-IX genes identified, 17 exhibited a high sequence similarity not limited to the coding sequence but also extended to the 5'and 3' surrounding regions (Additional file 6, Figure S6). However, since the expression profiles of these genes did not overlap completely, it is likely that their expression is driven by different promoter sequences. Interestingly, members of group IX displayed a particular clustered organization along chromosomes 7 and 16 (Figure 2), similar to the one observed in the poplar genome for the same group, although in this species the tandem repetition does not exceed four genes.
The tandem repetition of almost identical coding sequences suggests that these duplication events in the grapevine genome are quite recent. Moreover, it is tempting to speculate on the molecular mechanisms that led to the formation of such a long cluster of almost identical sequences, maybe involving the interaction of transposons or retroposons with the genomic sequence encoding ERF-IX genes. Magnani et al.  showed that a cyanobacterium (Trichodesmium erythraeum) possesses a HNH endonuclease protein with an AP2 DNA binding domain and suggested that the AP2 sequence approached the plant genome as a lateral transfer from a viral or bacterial donor. Gene duplication is more likely to be retained for gene families involved in signal transduction and transcriptional regulation . Nakano et al. suggested that an increased number of genes in groups III and IX might improve plant fitness since these genes are involved in the regulation of environmental stress responses. However, many other ERF groups, which do not have a high number of members, play a role in stress responses . Therefore, the functional role of the ERF-IX group alone does not seem enough to explain the high number of genes in this family and their repetition on chromosomes. Further analyses of the nucleotidic and amino acid sequences of the genes belonging to this family are likely to shed light on this topic.
Genes belonging to the RAV family are highly conserved among dicot species, which contain six members, and not very different in rice (5 members) .
The number of AP2 genes is very similar between Arabidopsis (18) and Vitis (19) and lower than that of poplar (29 members) (Table 1). The generally higher number of AP2/ERF genes in the poplar genome is likely to be a consequence of a recent whole-genome duplication in the poplar lineage . Curiously, a divergent member of the AP2 family, whose poplar homolog has been called Soloist, is maintained in a single copy in all plant genomes analyzed so far. The conservativeness of the sequence, plus the fact that it is only weakly related to other AP2/ERF genes, suggests that this gene diverged from the AP2 family early in the evolution of the plant species. Surprisingly, it also maintained its single copy feature during the separation of the monocot and dicot and the spread of the eurosid clade .
Differential expression of AP2/ERF genes in aerial tissues of Vitis vinifera
The development of an RT-qPCR platform (Additional file 7, Table S7) encompassing the complete number of AP2/ERF genes previously identified, together with knowledge about the molecular function of their homolog in other plant species, may help in associating a role with some of these transcription factors.
Traditionally, the expression of ERF genes has been associated with the molecular response to ethylene, from which their acronym derives. Although AP2/ERF transcription factors are regulated by a number of physical-chemical stimuli, many ERF genes are indeed ethylene responsive and therefore have been suggested as being involved in the ripening process in climacteric fruits. Wine grape is considered to be a non-climacteric fruit but ethylene still plays an important role during development and ripening processes such as anthocyanin accumulation  and alcohol dehydrogenase expression .
Although many ERF genes have been shown to play a role in the response to abiotic and biotic stresses, only a few genes have been described as regulators of organ development. Two homologues of the development-related ERFs AtLEP and AtBOL, VvERF067 and VvERF068, are constantly expressed in the berry skin tissues, whereas a vitis homolog of the cytokinin responsive factors exhibits leaf-specific expression.
Many members of the groups II, III and IX (VvERF085, VvERF087 and VvERF089) were highly expressed in all the tissues examined (Figure 3). In Arabidopsis, tobacco (Nicotiana tabacum) and tomato (Solanum lycopersicum), ERFs belonging to groups II and III play a role in response to abiotic stress [ and ] whereas group IX includes genes involved in the response to pathogens [ and ]. The relatively high expression of members of these groups in grapevine tissues, under non stress conditions, may represent a default defence system, which was already present before the pathogen attack or the environmental stress took place, activated only when a trigger signal is perceived.
AP2 and RAV genes are the main actors in the determination of cell fate and development from meristems to fully developed organs. However, with the exception of VvAP2-5, the expression of AP2 and RAV genes in the tissues considered was not restricted to a single organ. On the other hand, VvAP2-5 is mainly expressed in the inflorescence tissues (Figure 4).
Contrary to what was expected, the expression homolog of the AP2 gene, VvAP2-11, VvAP2-14 and VvAP2-15, was not restricted to floral organs but extended to vegetative and fruit tissues (Figure 4). VvAP2-9, VvAP2-11, VvAP2-12 and VvAP2-17, which are homologs of the Arabidopsis AIL-5, AP2, TOE1 and ANT genes respectively, had a similar expression pattern, with the highest expression in the leaf, inflorescence and stem, but not in the fruit tissues (Figure 4). No member of the AP2 or RAV family showed a fruit-specific expression pattern.
Several ERF-IX members (VvERF074, VvERF079, VvERF085, VvERF086, VvERF098, VvERF099, and VvERF100) reached high mRNA levels in the skin tissues of the berry, but not in the flesh tissues (Figure 3). This differential expression between skin and flesh among members of group IX may reflect the need for a constantly activated defence in the skin of the berry, which is continuously exposed to pathogens attack. This was not the case however in the flesh tissues, where bacteria and fungi can only be present once the physical barrier of the skin has been penetrated.
Expression changes of grapevine AP2/ERF genes in the veraison to full ripeness transition
Skin and flesh tissues from grape berries differ greatly during the transition from veraison to full ripeness with respect to ERF-gene expression. In the skin, 31 genes were up-regulated at least four-fold and 18 were down-regulated. In the flesh tissues, only 18 AP2/ERF genes were up-regulated and 30 down-regulated (Figure 5).
Group IX contained most of the genes induced or repressed by ripening in both skin and flesh tissues of the grapevine berry (Figure 5). Members of the same group are up-regulated in apple and plum fruit during ripening [28, 46]. Deluc et al.  identified a small number of ERF IX as being differentially regulated throughout the development of the grape berry. This was either due to the sensitivity limitation of the microarray technique or because the transcriptome of skin and flesh tissues was not analyzed separately. The ripening-related induction of many members of group IX suggests that the skin cells need to be prepared against pathogen attack in a phase during which the berry is enriched with sugars and therefore represents the ideal substrate for pathogen growth. This hypothesis is supported by the fact that VvGlucβ-1, VvGlucβ-2, VvChit4-1, VvPR10.1 and VvPR10L-b were up-regulated during ripening, preferentially in the skin tissues, together with ERFs belonging to group IX (Additional file 8, Figure S8). In addition, in a more extensive series of developmental stages in whole berries (Deluc et al. ), the expression of stress-related genes correlated with that of ERF-IX members (Additional file 9, Figure S9). Moreover, increased levels of PR proteins and chitinases during the transition from veraison to ripening have been reported for grape Barbera cultivar . In wine, PR proteins represent one of the most abundant classes of proteins . Since they can negatively affect wine clarity and stability, their regulation may represent an important trait for improving wine quality.
The ERF-V VvERF042 and VvERF046, homologs of the Arabidopsis SHN1 gene, were induced in both skin and flesh tissues, while VvERF044 was repressed. An up-regulation of group V members during berry development was also reported by Deluc et al. . In Arabidopsis, SHN1 regulates lipid biosynthetic pathways towards cutin synthesis , while a barley SHN1-homolog regulates lipid anabolism to generate organ adhesion . Of the direct targets of SHN1, Kannangara et al.  identified the Long Acyl-coA synthetase 2 (LACS2), which is required for cutin biosynthesis in wild-type plants. Interestingly, in Vitis vinifera cv. Cabernet, VvLACS7, a homolog of the Arabidopsis LACS2, was reported to be up-regulated during ripening  and its expression highly correlates with VvERF046 in the same subset (Spearman correlation coefficient 0.92). An equal number of ERF-III genes was induced and repressed in the skin, whereas only two were up-regulated in the flesh tissues (Tables 2, 3). Group III encodes transcription factors involved in the crosstalk between biotic and abiotic stress  and therefore may play a similar role as group IX during the last phases of ripening in the berry. Members of group IV, which were also involved in stress response regulation, were induced in the skin but not in the flesh (Table 2). The ripening related induction of group VIII-ERF has already been reported in plums [ and ]. Also in the grape berry, three members were up-regulated in the skin tissues (Table 2).
Previous studies on ripening regulated transcripts identified members of the group VII as being associated with ripening [ and ]. However, in our transcriptional analysis, none of the three ERF-VII genes showed a differential expression in the berry skin or flesh in the transition from veraison to full ripeness (Figure 3). It is possible that the ripening-related induction of these genes occurs at earlier stages of the ripening process.
Two AP2 genes, VvAP2-5 and VvAP2-16, homologs to the Arabidopsis PTL3 (AIL-6) and PTL4 (BBM), respectively were strongly up-regulated in skin tissues during ripening, whereas an AIL-9 homolog, VvAP2-9, was repressed. In Arabidopsis, The PTL proteins act in a dose-dependent way to regulate stem cell maintenance and meristem boundaries in Arabidopsis, probably by PIN-dependent auxin distribution [ and ]. Their function in fruit development has not yet been studied; however they represent interesting candidates as regulators of the transition from different ripening stages. The RAV gene VvRAV2 was up-regulated in the skin tissues (Table 2). Its closest Arabidopsis homolog, RAV2, has been proposed to act as repressor, through the RLFGV motif located at the C-terminus . Since VvRAV2 maintained the presence of this motif, it is possible to hypothesize a similar molecular activity.
In summary, in this study we identified 149 AP2/ERF genes in the grape genome and characterized their expression patterns in aerial vegetative, inflorescence and berry skin and flesh at two different ripening stages, veraison and full ripening. A comparison of homologs from other species, whose genome has been sequenced, together with their expression profiles, may help in an understanding of the role of these transcription factors in perennial plants. Vitis vinifera represent the only woody plant whose genome is fully sequenced and whose commercial value is due to fruit production. Unveiling the role of AP2/ERF transcription factors in the developmental and ripening processes in this species, may help molecular breeders in improving fruit quality.
Identification of Vitis AP2/ERF genomic sequences
A search of the Vitis vinifera genome database was performed in order to find all members of the AP2/ERF family. A double strategy to obtain every gene of the AP2/ERF family in the genome was used. The sequences of all members of the ERF family in the genome of Arabidopsis thaliana were downloaded from the DATF database  and the amino acid sequence of one or most representative members (i.e. the maximum number of different conserved motifs distinctive of each group) for each group defined by Nakano et al.  were used as queries to search the grapevine genome database (Genoscope, CEA - Institut de génomique, France) using the BLAT program . We also used the consensus sequence of the ERF domain to search the Vitis genome database on the NCBI web site. Every sequence identified was subsequently checked against the Arabidopsis and Poplar protein databases to confirm that it belonged to the AP2/ERF superfamily. As a final quality check, we confirmed the presence of the AP2 domain in every AP2/ERF Vitis vinifera gene candidate using SMART .
Phylogenetic and molecular evolutionary analyses were conducted using MEGA version 4 . To generate a phylogenetic tree, complete AP2/ERF predicted proteins Vitis vinifera were aligned using the ClustalW algorithm version 2.0 . The neighbour-joining method was used to construct different trees, using the pair-wise deletion option. The reliability of the obtained trees was tested using bootstrapping with 1000 replicates.
Gene prediction testing
Genomic sequences predicted to encode ERF genes that contained introns longer than 10 kb or that do not exist in Arabidopsis or poplar were considered as not reliable and four were (GSVIVT00021812001, GSVIVP00013482001, GSVIVT00019482001, and GSVIVP00007524001) chosen as examples. A primer pair annealing to the beginning and the end of the predicted coding sequence was used to check the reliability of the prediction. Primer pairs annealing within exons of the predicted coding sequence were used as controls to test that the gene considered was indeed expressed in the cDNA sample used. GSVIVT00034010001 was used as a positive control, since its coding sequence has also been predicted by EST analyses and its structure is strongly conserved by its homologues in different plant species such as potato, rice, Arabidopsis, sorrel (Licausi et al. unpublished).
Gene Expression Analyses
Total RNA, extracted using the RNeasy kit (Qiagen) according to the manufacturer's instructions, was subjected to DNase treatment using the TURBO DNA-free kit (Ambion). Five micrograms of each sample were reverse transcribed into cDNA using the Superscript III reverse transcriptase kit (Invitrogen). Real-time PCR amplification was carried out with the ABI Prism 6900 sequence detection system (Applied Biosystems), using a power sybr-green master mix (Applied Biosystems) according to the manufacturer's instructions. Gene-specific primers (Additional file 7, Table S7) were designed on regions not interrupted by introns using Quantprime software , tested for the specificity of each amplification reaction. This was done by analyzing the dissociation curves of each amplicon and the primer efficiency was calculated using LinReg . For all primer pairs and for all samples, the primer efficiency was between 1.7 and 2.0. Data were analyzed using 6900 SDS software 1.3 (Applied Biosystems). Dissociation curves for each amplicon were analyzed in order to verify the specificity of each amplification reaction. Transcript levels were normalized against the average of the grapevine housekeeping (HK) genes: EF1- gene (BQ799343), the ubiquitin gene (VvUB; CF406001), the actin gene (AB073011) and the glyceraldehydes dehydrogenase gene (EF192466). Relative gene expression in the transcriptomic analysis considering the different plant tissues was obtained by the formula CtHK-CtGene, where Ct is the cycle number at which a reaction reaches a specified fluorescence level. Hierarchical clustering of the gene expression data was performed using R software. For the relative gene expression level related to the ripening, the ΔΔCt  method was employed. Statistical analysis of ERF group overrepresentation among differentially expressed genes was performed applying the Fisher's Exact test .
List of abbreviations used
AIntegumenta Like 5
Cold responsive element binding factor
class IV chitinase
Differentiallly Expressed Genes
Drought Responsive Elemenet Binding protein
Ethylene Responsive Factor
β 1-3, Glucanase
Long AcylCoA synthase
Pathogen Related protein 10.1
Pathogen Related protein b like
Related To APETALA2
Related to ABI3/VP1
Real Time quantitative Polymerase Chain Reaction
Target of EAT1 1
Funding for the research described here was generously provided by Scuola Superiore Sant'Anna, by the BACCA Project funded by the ORVIT Consortium, by the Completamento del Centro di Genomica Funzionale Vegetale Project funded by the Cariverona Bank Foundation, and by the Structural and Functional Characterization of the Grapevine Genome (Vigna) Project funded by the Italian Ministry of Agricultural and Forestry Policies. We thank Pasqua Vini e Cantine (Verona, Italy) for allowing us to sample material from its vineyard.
- Wessler SR: Homing into the origin of the AP2 DNA binding domain. Trends in Plant Sci. 2005, 10: 54-56. 10.1016/j.tplants.2004.12.007.View Article
- Jofuku KD, den Boer BGW, van Montagu M, Okamuro JK: Control of Arabidopsis flower and seed development by the homeotic gene APETALA2. Plant Cell. 1994, 6: 1211-1225. 10.1105/tpc.6.9.1211.PubMed CentralPubMedView Article
- Ohme-Takagi M, Shinshi H: Ethylene-inducible DNA binding proteins that interact with an ethylene-responsive element. Plant Cell. 1995, 7: 173-182. 10.1105/tpc.7.2.173.PubMed CentralPubMedView Article
- Ecker JR: The ethylene signal transduction pathway in plants. Science. 1995, 268: 667-675. 10.1126/science.7732375.PubMedView Article
- Okamuro JK, Caster B, Villarroel R, Van Montagu M, Jofuku KD: The AP2 domain APETALA2 defines a large new family of DNA binding proteins in Arabidopsis. Proc Natl Acad Sci USA. 1997, 94: 7076-7081. 10.1073/pnas.94.13.7076.PubMed CentralPubMedView Article
- 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-478. 10.1093/nar/27.2.470.PubMed CentralPubMedView Article
- 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. Biochemical d Biophyisical Res Comm. 2002, 290: 998-1009. 10.1006/bbrc.2001.6299.View Article
- Nakano T, Suzuki K, Fujimura T, Shinshi H: Genome-wide analysis of the ERF gene family in Arabidopsis and rice. Plant Phys. 2006, 140: 411-432. 10.1104/pp.105.073783.View Article
- Masaki T, Mitsui N, Tsukagoshi H, Nishii T, Morikami A, Nakamura K: Activator of Spomin::LUC1/WRINKLED1 of Arabidopsis thaliana transactivates sugar-inducible promoters. Plant Cell Phys. 2005, 46: 547-556. 10.1093/pcp/pci072.View Article
- Welsch R, Maass D, Voegel T, Della Penna D, Beyer P: Transcription factor RAP2.2 and its interacting partner SINAT2: stable elements in the carotenogenesis of Arabidopsis leaves. Plant Phys. 2007, 145: 1073-1085. 10.1104/pp.107.104828.View Article
- Jiang C, Iu B, Singh J: Requirement of a CCGAC cis-acting element for cold induction of the BN115 gene from winter Brassica napus. Plant Mol Biol. 1996, 30: 679-684. 10.1007/BF00049344.PubMedView Article
- Nole-Wilson S, Krizek BA: DNA binding properties of the Arabidopsis floral development protein AINTEGUMENTA. Nucleic Acids Res. 2000, 21: 4076-4082. 10.1093/nar/28.21.4076.View Article
- Gong W, He K, Covington M, Dinesh-Kumarb SP, Snyder M, Harmerc SL, Zhua Y-X, Deng XW: The development of protein microarrays and their applications in DNA-protein and protein-protein interaction analyses of Arabidopsis transcription factors. Mol Plant. 2007, 1: 27-41. 10.1093/mp/ssm009.PubMedView Article
- Krizek BA: AINTEGUMENTA utilizes a mode of DNA recognition distinct from that used by proteins containing a single AP2 domain. Nucleic Acids Res. 2003, 31: 1859-1868. 10.1093/nar/gkg292.PubMed CentralPubMedView Article
- Elliott RC, Betzner AS, Huttner E, Oakes MP, Tucker WQJ, Gerentes D, Perez P, Smyth DR: AINTEGUMENTA, an APETALA2-like gene of Arabidopsis with pleiotropic roles in ovule development and floral organ growth. Plant Cell. 1996, 8: 155-168. 10.1105/tpc.8.2.155.PubMed CentralPubMedView Article
- Stockinger EJ, Gilmour SJ, Thomashow MF: Arabidopsis thaliana CBF1 encodes an AP2 domain-containing transcriptional activator that binds to the C-repeat/DRE a cis-acting DNA regulatory element that stimulates transcription in response to low temperature and water deficit. Proc Natl Acad Sci USA. 1997, 94: 1035-1040. 10.1073/pnas.94.3.1035.PubMed CentralPubMedView Article
- Liu Q, Kasuga M, Sakuma Y, Abe H, Miura S, Yamaguchi-Shinozaki K, Shinozaki K: Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought- and low-temperature-responsive gene expression, respectively, in Arabidopsis. Plant Cell. 1998, 10: 1391-1406. 10.1105/tpc.10.8.1391.PubMed CentralPubMedView Article
- Yamamoto S, Suzuki K, Shinshi H: Elicitor-responsive, ethylene-independent activation of GCC box-mediated transcription that is regulated by both protein phosphorylation and dephosphorylation in cultured tobacco cells. Plant J. 1999, 20: 571-579. 10.1046/j.1365-313X.1999.00634.x.PubMedView Article
- Boutilier K, Offringa R, Sharma VK, Kieft H, Ouellet T, Zhang LM, Hattori J, Liu CM, van Lammeren AAM, Miki BLA, Custers JBM, van Lookeren Campagne MM: Ectopic expression of BABY BOOM triggers a conversion from vegetative to embryonic growth. Plant Cell. 2002, 14: 1737-1749. 10.1105/tpc.001941.PubMed CentralPubMedView Article
- Jaillon O, Aury J-M, Noel B, Policrit A, Clepet C, Casagrande A, Choisne N, Aubourg C, Vitulo N, Jubin C, Vezzi A, Legeai F, Hugueney P, Dasilva C, Horner D, Mica E, Jublot D, Poulain J, Bruyère C, Billault A, Segurens B, Gouyvenoux M, Ugarte E, Cattonaro F, Anthouard V, Vico V, Del Fabbro C, Alaux M, Di Gaspero G, Dumas V, et al: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007, 449: 463-467. 10.1038/nature06148.PubMedView Article
- Troggio M, Malacarne G, Coppola G, Segala C, Cartwright DA, Pindo M, Stefanini M, Mank R, Moroldo M, Morgante M, Grando MS, Velasco R: A dense single-nucleotide polymorphism-based genetic linkage map of grapevine (Vitis vinifera L.) anchoring Pinot Noir bacterial artificial chromosome contigs. Genetics. 2007, 176: 2637-2650. 10.1534/genetics.106.067462.PubMed CentralPubMedView Article
- Velasco R, Zharkikh A, Troggio M, Cartwright DA, Cestaro C, Pruss D, Pindo P, FitzGerald LM, Vezzulli S, Reid J, Malacarne G, Iliev D, Coppola C, Wardell B, Micheletti D, Macalma T, Facci M, Mitchell JT, Perazzolli M, Eldredge G, Gatto G, Oyzerski R, Moretto M, Gutin N, Stefanini M, Chen Y, Segala C, Davenport C, Demattè L, Mraz A, et al: A high quality draft consensus sequence of the genome of a heterozygous grapevine variety. PLoS ONE. 2007, 2: e1326-10.1371/journal.pone.0001326.PubMed CentralPubMedView Article
- Zhuang J, Peng R-H, Cheng Z-M(M), Zhang J, Cai B, Zhang Z, Gao F, Zhu B, Fu X-Y, Jin X-F, Chen J-M, Qiao Y-S, Xiong A-S, Yao Q-H: Genome-wide analysis of the putative AP2/ERF family genes in Vitis vinifera. Scientia Horticolturae. 2009, 123: 73-81. 10.1016/j.scienta.2009.08.002.View Article
- Zhuang J, Cai B, Peng RH, Zhu B, Jin XF, Xue Y, Gao F, Fu XY, Tian YS, Zhao W, Qiao YS, Zhang Z, Xiong AS, Yao QH: Genome-wide analysis of the AP2/ERF gene family in Populus trichocarpa. Biochemical Biophysical Res Comm. 2008, 371: 468-474. 10.1016/j.bbrc.2008.04.087.View Article
- Howe KL, Chothia T, Durbin R: GAZE: a generic framework for the integration of gene-prediction data by dynamic programming. Genome Res. 2002, 12: 418-1427. 10.1101/gr.149502.View Article
- Díaz-Riquelme J, Lijavetzky D, Martínez-Zapater JM, Carmona MJ: Genome-wide analysis of MIKCC-type MADS box genes in grapevine. Plant Phys. 2009, 149: 354-369.View Article
- Czechowski T, Bari RP, Stitt M, Scheible W-R, Udvardi MK: Real-time RT-PCR profiling of over 1400 Arabidopsis transcription factors: unprecedented sensitivity reveals novel root- and shoot-specific genes. Plant J. 2002, 38: 366-379. 10.1111/j.1365-313X.2004.02051.x.View Article
- Wang A, Tan D, Takahashi A, Li TZ, Harada T: MdERFs, two ethylene-response factors involved in apple fruit ripening. J Exp Bot. 2007, 58: 3743-3748. 10.1093/jxb/erm224.PubMedView Article
- Chen G, Hu Z, Grierson D: Differential regulation of tomato ethylene responsive factor LeERF3b, a putative repressor, and the activator Pti4 in ripening mutants and in response to environmental stresses. J Plant Phys. 2007, 6: 662-670.
- El-Sharkawy I, Kim WS, El-Kereamy A, Jayasankar S, Svircev AMD, Brown CW: Isolation and characterization of four ethylene signal transduction elements in plums (Prunus salicina L.). J Exp Bot. 2007, 58: 3631-3643. 10.1093/jxb/erm213.PubMedView Article
- Deluc LG, Grimplet J, Wheatley MD, Tillett RL, Quilici DR, Osborne C, Schooley DA, Schlauch KA, Cushman JC, Cramer GR: Transcriptomic and metabolite analyses of Cabernet Sauvignon grape berry development. BMC Genomics. 2007, 8: 429-10.1186/1471-2164-8-429.PubMed CentralPubMedView Article
- Pré M, Atallah M, Champion C, De Vos M, Pieterse CMJ, Memelink J: The AP2/ERF domain transcription factor ORA59 integrates jasmonic acid and ethylene signals in plant defense. Plant Phys. 2008, 147: 1347-1357.View Article
- Zhang JY, Broeckling CD, Sumner LW, Wang Z-Y: Heterologous expression of two Medicago truncatula putative ERF transcription factor genes, WXP1 and WXP2, in Arabidopsis led to increased leaf wax accumulation and improved drought tolerance, but differential response in freezing tolerance. Plant Mol Biol. 2007, 64: 265-278. 10.1007/s11103-007-9150-2.PubMedView Article
- Aharoni A, Dixit S, Jetter R, Thoenes E, van Arkel G, Pereira A: The SHINE clade of AP2 domain transcription factors activates wax biosynthesis, alters cuticle properties, and confers drought tolerance when overexpressed in Arabidopsis. Plant Cell. 2004, 16: 2463-2480. 10.1105/tpc.104.022897.PubMed CentralPubMedView Article
- Rogiers SY, Hatfield JM, Jaudzems VG, White RG, Keller M: Grape berry cv. Shiraz epicuticular wax and transpiration during ripening and preharvest weight loss. American J Enology Viticulture. 2004, 55: 121-127.
- Magnani E, Sjölanderc K, Hakea S: From endonucleases to transcription factors: evolution of the AP2 DNA binding domain in plants. Plant Cell. 2004, 16: 2265-2277. 10.1105/tpc.104.023135.PubMed CentralPubMedView Article
- Blanc G, Wolfe KH: Functional divergence of duplicated genes formed by polyploidy during Arabidopsis evolution. Plant Cell. 2004, 16: 1679-1691. 10.1105/tpc.021410.PubMed CentralPubMedView Article
- Swaminathan K, Peterson K, Jack T: The plant B3 superfamily. Trends in Plant Sci. 2008, 13: 647-655. 10.1016/j.tplants.2008.09.006.View Article
- Tuskan GA, DiFazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A, Schein J, Sterck L, Aerts A, Bhalerao RR, Bhalerao RP, Blaudez D, Boerjan B, Brun A, Brunner A, Busov V, Campbell M, Carlson J, Chalot M, Chapman J, Chen G-L, Cooper D, Coutinho PM, Couturier J, Covert S, Cronk Q, Van de Peer Y: The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science. 2006, 313: 1596-1604. 10.1126/science.1128691.PubMedView Article
- El-Kereamy A, Chervin C, Roustan J, Cheynier V, Souquet JM, Moutoumet M, Raynal J, Ford CM, Latché A, Pech JC, Bouzayen M: Exogenous ethylene stimulates the long-term expression of genes related to anthocyanin biosynthesis in grape berries. Phys Plantarum. 2003, 119: 175-182. 10.1034/j.1399-3054.2003.00165.x.View Article
- Tesniere C, Pradal M, El-Kereamy A, Torregrosa L, Chatelet P, Roustan JP, Chervin C: Involvement of ethylene signalling in a non-climacteric fruit: new elements regarding the regulation of ADH expression in grapevine. J Exp Bot. 2004, 55: 2235-2240. 10.1093/jxb/erh244.PubMedView Article
- Gilmour SJ, Zarka DG, Stockinger EJ, Salazar MP, Houghton JM, Thomashow MF: Low temperature regulation of the Arabidopsis CBF family of AP2 transcriptional activators as an early step in cold-induced COR gene expression. Plant J. 1998, 16: 433-442. 10.1046/j.1365-313x.1998.00310.x.PubMedView Article
- Tsutsui T, Kato W, Asada Y, Sako K, Sato Takeo, Sonoda Y, Kidokoro S, Yamaguchi-Shinozaki K, Tamaoki M, Arakawa K, et al: DEAR1, a transcriptional repressor of DREB protein that mediates plant defense and freezing stress responses in Arabidopsis. J Plant Res. 2009, 122: 633-643. 10.1007/s10265-009-0252-6.PubMedView Article
- Berrocal-Lobo M, Molina A, Solano R: Constitutive expression of ETHYLENE-RESPONSE-FACTOR 1 in Arabidopsis confers resistance to several necrotrophic fungi. Plant J. 2002, 29: 23-32. 10.1046/j.1365-313x.2002.01191.x.PubMedView Article
- Gu YQ, Wildermuth MC, Chakravarthy S, Loh YT, Yang C, He X, Han Y, Martin GB: Tomato transcription factor Pti4, Pti5, and Pti6 activate defence responses when expressed in Arabidopsis. Plant Cell. 2002, 14: 817-831. 10.1105/tpc.000794.PubMed CentralPubMedView Article
- El-Sharkawy S, Sherif I, Mila M, Bouzayen M, Jayasankar S: Molecular characterization of seven genes encoding ethylene-responsive transcriptional factors during plum fruit development and ripening. J Exp Bot. 2009, 60: 907-922. 10.1093/jxb/ern354.PubMed CentralPubMedView Article
- Negri AS, Prinsi B, Rossoni M, Failla O, Scienza A, Cocucci M, Espen L: Proteome changes in the skin of the grape cultivar Barbera among different stages of ripening. BMC Genomics. 2009, 9: 378-10.1186/1471-2164-9-378.View Article
- Kannangara R, Branigan C, Liu Y, Penfielda T, Raob V, Mouillec G, Höftec H, Pauly M, Riechmann JL, Brouna P: The transcription factor WIN1/SHN1 regulates cutin biosynthesis in Arabidopsis thaliana. Plant Cell. 2007, 19: 1278-1294. 10.1105/tpc.106.047076.PubMed CentralPubMedView Article
- Taketa S, Amano S, Tsujino Y, Sato T, Saisho D, Kakeda K, Nomura M, Suzuki T, Matsumoto T, Sato K, Kanamori H, Kawasaki S, Takeda K: Barley grain with adhering hulls is controlled by an ERF family transcription factor gene regulating a lipid biosynthesis pathway. Proc Natl Acad Sci USA. 2008, 105: 4062-4067. 10.1073/pnas.0711034105.PubMed CentralPubMedView Article
- Sun S, Yu JP, Chen F, Zhao TJ, Fang XH, Li YQ, Sui SF: TINY, a dehydration-responsive element (DRE)-binding protein-like transcription factor connecting the DRE- and ethylene-responsive element-mediated signaling pathways in Arabidopsis. J Exp Bot. 2008, 283: 6261-6271.
- Ziliotto F, Begheldo M, Rasori A, Bonghi C, Tonutti P: Transcriptome profiling of ripening nectarine (Prunus persica L. Batsch) fruit treated with 1-MCP. J Exp Bot. 2008, 59: 2781-2791. 10.1093/jxb/ern136.PubMed CentralPubMedView Article
- Galinha C, Hofhuis H, Luijten M, Willemsen V, Blilou I, Heidstra R, Scheres B: PLETHORA proteins as dose-dependent master regulators of Arabidopsis root development. Nature. 2007, 449: 1053-1057. 10.1038/nature06206.PubMedView Article
- Krizek BA: AINTEGUMENTA and AINTEGUMENTA-like 6 act redundantly to regulate Arabidopsis floral growth and patterning. Plant Phys. 2009, 150: 1916-1929. 10.1104/pp.109.141119.View Article
- Ikeda M, Ohme-Takagi M: A novel group of transcriptional repressors in Arabidopsis. Plant Cell Phys. 2009, 50: 970-975. 10.1093/pcp/pcp048.View Article
- Guo A, He K, Liu D, Bai S, Gu X, Wei L, Luo J: DATF: a database of Arabidopsis transcription factors. Bioinformatics. 2005, 21: 2568-2569. 10.1093/bioinformatics/bti334.PubMedView Article
- Kent WJ: BLAT - The BLAST-Like Alignment Tool. Genome Res. 2002, 4: 656-664.View Article
- Letunic I, Tobias Doerks T, Bork P: SMART 6: recent updates and new developments. Nucleic Acids Res. 2008, 37: D229-D232. 10.1093/nar/gkn808.PubMed CentralPubMedView Article
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.PubMedView Article
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.PubMedView Article
- Arvidsson S, Kwasniewski M, Riano-Pachon DM, Mueller-Roeber B: QuantPrime - a flexible tool for reliable high-throughput primer design for quantitative PCR. BMC Bioinformatics. 2008, 9: 465-10.1186/1471-2105-9-465.PubMed CentralPubMedView Article
- Ramakers C, Ruijter JM, Lekanne Deprez LH, Moorman AFM: Assumption-free analysis of quantitative real-time PCR data. Neurosci Letts. 2003, 339: 62-66. 10.1016/S0304-3940(02)01423-4.View Article
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.PubMedView Article
- Agresti A: A Survey of Exact Inference for Contingency Tables. Statistical Sci. 1992, 7: 131-153. 10.1214/ss/1177011454.View Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<url>http://creativecommons.org/licenses/by/2.0</url>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.