1. Introduction
Transcription factors (TFs) are key regulators of gene expression that control various biological processes in plants, including growth, development, metabolism, and responses to environmental stimuli [1,2]. Among the numerous plant TF families, the APETALA2/Ethylene-Responsive Factor (AP2/ERF) superfamily represents one of the largest and most functionally diverse groups [3,4]. This superfamily is characterized by the presence of the highly conserved AP2/ERF DNA-binding domain, typically consisting of 60–70 amino acids [5,6]. Based on the number and structural features of the AP2/ERF domains, it is conventionally classified into four major families: AP2 (containing two AP2/ERF domains), ERF (containing a single AP2/ERF domain), RAV (containing one AP2/ERF domain and one B3 domain), and Soloist (a divergent group). The ERF family is further divided into the ERF and DREB (Dehydration-Responsive Element-Binding protein) subfamilies [4,7,8]. This classification system provides a framework for understanding the functional diversity of AP2/ERF proteins.
The AP2/ERF superfamily has been extensively characterized in model plants and several crops, revealing their crucial roles in various biological processes. The AP2 family members primarily regulate floral meristem determinacy, floral organ identity, seed development, and fruit ripening [9,10]. ERF family members are predominantly involved in responses to biotic and abiotic stresses, including pathogen infection, drought, salinity, and extreme temperatures [11,12]. RAV family members function in leaf senescence, flowering time control, and responses to ethylene and brassinosteroid hormones [13]. The Soloist family, though small, has been implicated in various stress responses, with members from different species shown to enhance tolerance to salt stress and pathogen resistance [14,15]. The functional significance of these genes makes them attractive targets for crop improvement through molecular breeding and genetic engineering approaches.
The Asteraceae (Compositae) family is one of the largest and most successful plant families, comprising approximately 10% of all flowering plants with an estimated 32,000–34,000 species distributed worldwide [16]. The remarkable evolutionary success of this family is largely attributed to its unique inflorescence structure, the capitulum, which represents one of the most significant innovations in plant reproductive morphology [17,18]. The capitulum, often perceived as a single flower, is actually a complex inflorescence consisting of numerous individual florets arranged on a common receptacle and surrounded by protective bracts [19]. A distinctive characteristic of many Asteraceae species is the presence of dimorphic florets within a single capitulum: peripheral ray florets with enlarged, zygomorphic corollas that often function primarily for attraction, and central disc florets with actinomorphic corollas that are typically hermaphroditic and ensure reproduction [20]. In Arabidopsis and most model plant species, the outermost floral organs are referred to as sepals and petals, following the classical floral architecture. However, in Asteraceae, these structures have undergone significant evolutionary modifications resulting in distinct terminology. The first whorl organs are referred to as pappus in Asteraceae and often appear as scales, bristles, or hair-like structures that aid in seed dispersal [20,21]. Similarly, the second whorl organs are termed corolla in Asteraceae, which, in ray florets, often fuse to form the characteristic ligule structure [22]. This specialized terminology reflects the unique evolutionary adaptations in Asteraceae floral architecture and highlights the morphological innovations that contribute to the success of this diverse plant family [19].
The molecular mechanisms underlying capitulum development and floret differentiation in Asteraceae have garnered substantial research interest in recent years. Investigations into model Asteraceae species, such as sunflower (Helianthus annuus), gerbera (Gerbera hybrida), and chrysanthemum (Chrysanthemum morifolium), have revealed complex regulatory networks involving multiple transcription factor families [23,24]. Notably, CYCLOIDEA2 (CYC2)-like TCP transcription factors are crucial for determining ray and disc floret identity and establishing flower symmetry in Asteraceae [25,26]. Based on the classic ABCDE model of floral organ development, MADS-box genes, particularly A-, B-, C-, and E-class genes, have been implicated in the differentiation of ray and disc florets in Asteraceae [27,28]. The regulatory networks of TCP proteins, in conjunction with MADS-box transcription factors, are vital for establishing floral dorsoventral asymmetry and regulating cell proliferation and differentiation, thereby collectively shaping floral architecture [29,30].
Within the classic ABCDE model of floral development, AP2 occupies a unique position as the only non-MADS-box transcription factor, fulfilling the A-class function that establishes the foundation for proper floral organ patterning and identity determination [31,32]. In Arabidopsis, mutations in the AP2 gene result in homeotic transformations where sepals are converted to carpel-like structures and petals to stamen-like structures [33]. Despite the well-established role of AP2 genes in flower development in model species, their function in Asteraceae capitulum development and floret differentiation remains largely unexplored. The limited studies available suggest that AP2 genes in Asteraceae may have evolved specialized functions related to the unique architecture of the capitulum and dimorphic florets. In chrysanthemum, the AP2/ERF member WRINKLED1 (WRI1) shows differential expression between contrasting ray floret types, suggesting its involvement in determining ray floret morphology. These AP2/ERF proteins likely mediate auxin and ethylene signaling to regulate the distinctive characteristics of ray florets [28]. This functional divergence of AP2 genes in Asteraceae compared to non-Asteraceae species may represent an evolutionary adaptation contributing to the diversity and success of the capitulum inflorescence across the family.
Marigold (Tagetes erecta), native to South America, is an economically important ornamental plant in the Asteraceae family, valued for its distinctive flowers, medicinal properties, and as a source of carotenoids, particularly lutein [34,35]. The inflorescence of marigold is a capitulum, characterized by peripheral ray florets and central disc florets [36]. Specifically, ray florets possess enhanced corollas but lack stamens, while disc florets contain complete floral structures [37]. This distinctive floral structure makes marigold an excellent model system for studying floral development and organ differentiation in Asteraceae. The recent publication of the chromosome-scale genome assembly of marigold [38] has provided an unprecedented opportunity to investigate the molecular mechanisms underlying its development and adaptive traits. Previous studies in marigold have identified and characterized several MADS-box genes and CYC2 genes involved in floral development [39,40,41,42], providing valuable insights into the genetic control of flower organ formation. However, a systematic analysis of the AP2/ERF superfamily—particularly AP2 genes that function as key regulators of flower organ development alongside MADS-box genes—has not yet been conducted in marigold.
In this study, we conducted a genome-wide identification and comprehensive characterization of AP2/ERF superfamily genes in marigold. We analyzed their phylogenetic relationships, gene structures, protein motifs, chromosomal distribution, and duplication patterns. Additionally, we investigated the expression profiles of AP2 family genes across different tissues and developmental stages, with a particular focus on floral organs in ray and disc florets. To better understand the regulatory networks involving AP2 family genes, we also performed protein interaction network analysis based on homology with Arabidopsis proteins, revealing potential interactions between AP2 family, MADS-box, and TCP family proteins. Our findings provide valuable insights into the evolutionary history and potential functions of AP2/ERF superfamily genes in marigold, establishing a foundation for future functional studies and genetic improvement of this economically important plant.
2. Materials and Methods
2.1. Plant Materials
The marigold inbred line V5, exhibiting typical Asteraceae capitulum structure with peripheral ray florets and central disc florets, was developed through 5 generations of self-crossing of cultivar ‘Vanilla’. Plants were grown in the Floriculture Teaching Experimental Base at Huazhong Agricultural University, Wuhan, Hubei Province, China (30°28′36.5″ N, 114°21′59.4″ E) under natural field conditions during the fall of 2024, with standard management practices.
2.2. Identification, Physicochemical Property Analysis, and Subcellular Localization of AP2/ERF Genes in Marigold
Marigold genome and protein sequences were obtained from the NCBI database (
Physicochemical properties of the identified AP2/ERF proteins, including amino acid length, molecular weight, theoretical isoelectric point (pI), instability index, and aliphatic index, were predicted using the ExPASy ProtParam tool (
2.3. Phylogenetic Analysis of AP2/ERF Genes
Arabidopsis AP2/ERF protein sequences were obtained from The Arabidopsis Information Resource (TAIR) database (
2.4. Gene Structure Analysis and Motif Identification
The exon–intron structures of marigold AP2/ERF genes were analyzed and visualized using TBtools V2.225 software based on the genome annotation file [51]. Conserved motifs in AP2/ERF proteins were identified using the MEME Suite v5.5.7 online tool (
2.5. Chromosomal Mapping and Synteny Analysis
Based on the marigold genome annotation file, a chromosome localization map of AP2/ERF genes was constructed via TBtools V2.225 software. Core parameters were input, including chromosome assignment, start/end positions of 177 genes, and corresponding chromosome physical lengths. Gene density for each chromosome was calculated using 500 kb genetic intervals. A gradient color scale was employed to characterize gene density distribution. Finally, a standardized chromosome distribution map was generated [53]. According to the linear distribution order of genes on chromosomes, members of this family were systematically named from TeAP2/ERF001 to TeAP2/ERF177.
For comparative genomic analysis, genome annotation files for Arabidopsis (Arabidopsis thaliana) and lettuce (Lactuca sativa) were obtained from The Arabidopsis Information Resource (
2.6. RNA Extraction from Various Tissues and Organs
During the flowering period, we collected samples for tissue-specific expression analysis, including roots, stems, leaves, receptacles, bracts, seeds, floral organs from disc florets (pappus, corollas, stamens, pistils, and ovaries) and ray florets (pappus, corollas, pistils, and ovaries). Additionally, we also collected flowers buds at different developmental stages: flower buds with diameter of 0–1 mm (ray florets and outermost disc florets primordia formation), flower buds with diameter of 4–5 mm (stamen development in outermost disc florets), flower buds with diameter of 9–10 mm (floral organ morphological completion) [40]. All plant materials were rapidly frozen in liquid nitrogen and stored at −80 °C for RNA extraction.
Total RNA was isolated using the SteadyPure RNA Extraction Kit (Accurate Biotechnology, Changsha, China). RNA concentration was determined using a NanoDrop 2000c spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA), and RNA integrity was assessed by 1% agarose gel electrophoresis. First-strand cDNA was synthesized using the Evo M-MLV RT Mix Kit (Accurate Biotechnology, Changsha, China) according to the manufacturer’s instructions.
2.7. Expression Analysis of AP2 Family Genes
Quantitative real-time PCR (qRT-PCR) was performed to analyze the expression patterns of AP2 family genes across different tissues and developmental stages. PCR reactions were conducted using the SYBR Green Premix Pro Taq HS qPCR Kit III (Low Rox Plus) (Accurate Biotechnology, Changsha, China) on the QuantStudio™ 6 Flex Real-Time PCR System (Thermo Fisher Scientific, Waltham, MA, USA). The marigold β-actin gene was used as an internal reference for expression normalization [36]. Gene-specific primers were designed using Primer Premier 5.0 software (
2.8. Protein–Protein Interaction Network Analysis
To gain insights into the potential regulatory relationships of AP2 family proteins in marigold, we analyzed their putative protein–protein interactions using the STRING database v12.0 (
3. Results
3.1. Identification and Classification of AP2/ERF Superfamily Genes in Marigold
To identify AP2/ERF superfamily genes in the marigold genome, we performed a comprehensive analysis using the Hidden Markov Model (HMM) profile of the AP2/ERF domain (PF00847). After validation through conserved domain analysis, we identified a total of 177 AP2/ERF genes in the marigold genome, which were subsequently named TeAP2/ERF001 to TeAP2/ERF177 based on their chromosomal positions (Supplementary Table S2).
The identified AP2/ERF proteins varied considerably in their physicochemical properties. The length of these proteins ranged from 52 to 624 amino acids, with molecular weights ranging from 6.00 to 69.10 kDa. The theoretical isoelectric points (pI) varied from 4.23 to 9.87, with 107 proteins exhibiting acidic properties (pI < 7) and 70 displaying basic properties (pI > 7). The instability index analysis revealed that approximately 52% of the AP2/ERF proteins were unstable (instability index > 50), while the remaining 48% were stable. Subcellular localization prediction indicated that most of TeAP2/ERF proteins were localized in the nucleus, consistent with their function as transcription factors (Supplementary Table S3).
Based on the number and configuration of AP2/ERF domains, the 177 marigold TeAP2/ERF genes were classified into four families: AP2 (28 genes), ERF (143 genes), RAV (4 genes), and Soloist (2 genes) (Table 1). Members of the AP2 family contain two AP2 domains, while ERF family members possess a single AP2 domain, and RAV family members have one AP2 domain and one B3 domain.
3.2. Phylogenetic Analysis of AP2/ERF Superfamily Members in Marigold
To understand the evolutionary relationships of AP2/ERF genes in marigold, we constructed a maximum-likelihood phylogenetic tree using the full-length protein sequences of AP2/ERF superfamily members from marigold and Arabidopsis (Figure 1). The resulting tree clearly separated the AP2/ERF proteins into distinct clades corresponding to the AP2, ERF, RAV, and Soloist families, confirming our initial classification.
Within the AP2 family, 28 marigold AP2 proteins formed three distinct groups, consistent with the classification in Arabidopsis. Based on their sequence similarities and phylogenetic positions, these groups likely represent functional equivalents of the AP2, ANT (AINTEGUMENTA), and AIL (AINTEGUMENTA-LIKE) groups described in other plant species [55]. In the ERF family, containing the ERF and DREB subfamilies, were further divided into 12 groups (ERF-B1 to B6 and DREB-A1 to A6), following the nomenclature proposed by Nakano et al. [3] (Table 1). The clustering pattern of marigold AP2/ERF proteins with their Arabidopsis counterparts suggests functional conservation across these evolutionarily distant species. Additionally, the uneven distribution of marigold AP2/ERF proteins across different groups, with notable expansions in the ERF-B3 and DREB-A4 groups, suggests differential selection pressures on specific groups during the evolutionary history of marigold.
3.3. Gene Structure and Conserved Motif Analysis
To gain insights into the structural diversity and functional evolution of AP2/ERF genes in marigold, we analyzed their exon–intron structures and conserved motifs (Supplementary Figure S1). The analysis revealed distinct structural patterns among different families and groups, which generally corresponded with their phylogenetic relationships.
The AP2 family genes exhibited the most complex gene structures, containing four to eight exons. In contrast, most ERF subfamily and DREB subfamily genes generally had one to two exons, RAV family genes consistently contained one exon, and Soloist family genes had six exons. These structural patterns are consistent with those observed in other plant species [56], suggesting evolutionary conservation of gene structures within each family.
The distribution of introns within the AP2/ERF domain was also family-specific. AP2 family genes typically contained introns within their AP2/ERF domains, while the majority of ERF subfamily and DREB subfamily genes lacked introns in this region. This differential intron pattern likely reflects the ancient divergence between the AP2 family and ERF family and may have functional implications for the regulation of gene expression and alternative splicing.
Conserved motif analysis using the MEME program identified 10 distinct motifs (designated as Motifs 1–10) across the marigold AP2/ERF proteins (Supplementary Figures S1 and S4). Motifs 1, 2, 3, and 4 were identified as the core components of the AP2/ERF domain and were present in all AP2/ERF proteins, confirming their conserved nature across the superfamily. The distribution of the remaining motifs showed family- and group-specific patterns. For instance, Motifs 5 and 8 were specifically present in the AP2 family. Motifs 9 and 10 were exclusively found in B3 of ERF subfamily proteins. Motif 6 was only present in the DREB subfamily and Motif 7 was only present in A1 and A4 of the DREB subfamily.
3.4. Chromosomal Distribution and Duplication Analysis
To investigate the genomic organization and evolutionary expansion of the AP2/ERF superfamily genes in marigold, we mapped their chromosomal positions and analyzed duplication events (Figure 2). The 177 AP2/ERF genes were unevenly distributed across 11 of the 12 marigold chromosomes. Chromosome 1 harbored the highest number of AP2/ERF genes (38 genes), while chromosome 12 contained no AP2/ERF genes.
Analysis of gene duplication events identified 56 pairs of segmentally duplicated TeAP2/ERF genes (Figure 3), predominantly involving ERF subfamily and DREB subfamily members. The prevalence of segmental duplications suggests that whole-genome or segmental duplication events have contributed significantly to the expansion of the AP2/ERF superfamily members in marigold. Calculation of non-synonymous (Ka) to synonymous (Ks) substitution ratios for the duplicated gene pairs revealed that all pairs had Ka/Ks ratios less than 1 (Supplementary Table S5), such as TeAP2/ERF011-TeAP2/ERF043 (Ka/Ks = 0.101) and TeAP2/ERF052-TeAP2/ERF107 (Ka/Ks = 0.093). This indicates that the AP2/ERF superfamily genes are subject to strong purifying selection within the species, with highly conserved functions where mutations likely disrupt their biological functions. Some gene pairs show Ka/Ks ratios approaching 0.5, such as TeAP2/ERF059-TeAP2/ERF119 (Ka/Ks = 0.463) and TeAP2/ERF033-TeAP2/ERF064 (Ka/Ks = 0.366). While still less than 1, these values suggest that certain families may be experiencing weaker purifying selection or contain local sites with positive selection signals that require further validation.
To understand the evolutionary relationships of AP2/ERF genes between marigold and other plant species, we performed synteny analysis with Arabidopsis (a model dicot) and lettuce (a related Asteraceae species) (Figure 4). We identified 49 syntenic AP2/ERF gene pairs between marigold and Arabidopsis and 93 pairs between marigold and lettuce. The higher number of syntenic pairs between marigold and lettuce reflects their closer evolutionary relationship as members of the Asteraceae family. The syntenic gene pairs were predominantly distributed in the AP2 and ERF family members, suggesting that these families may have experienced stronger selection pressure for functional conservation during the divergence of these species.
3.5. Expression Patterns of AP2 Family Genes Reveal Their Roles in Tissue-Specific and Floral Organ Development
AP2 family members function as A-class genes in the classical ABCDE model of floral organ development, playing essential roles in the specification of sepals and petals [32,57,58]. Additionally, given our research focus on understanding the molecular mechanisms controlling the distinctive capitulum architecture in marigold, we prioritized the 28 AP2 family genes and analyzed their expression patterns across different tissues and developmental stages using qRT-PCR.
We analyzed the temporal expression patterns of AP2 family genes during flower bud development, examining flower buds at different developmental stages (0–1 mm, 4–5 mm, and 9–10 mm diameter) (Figure 5 and Figure S2). Several genes showed stage-specific expression patterns. For instance, TeAP2/ERF014, TeAP2/ERF030, TeAP2/ERF038, TeAP2/ERF050, and TeAP2/ERF145 were highly expressed in early-stage flower buds, suggesting their involvement in floral organ initiation, while TeAP2/ERF009, TeAP2/ERF026, TeAP2/ERF036, and TeAP2/ERF042 showed increasing expression as flower development progressed, indicating their potential roles in later stages of floral development.
Based on their predominant expression patterns of the 28 AP2 family genes across different tissues (Figure 6 and Figure S3), these genes could be classified into four groups: (1) those with high expression in vegetative tissues (roots, stems, and leaves), (2) those predominantly expressed in reproductive tissues (receptacle and flowers), (3) those specifically expressed in seeds, and (4) those with broad expression across multiple tissues.
Several AP2 family genes, including TeAP2/ERF030, TeAP2/ERF036, TeAP2/ERF048, TeAP2/ERF138, and TeAP2/ERF154, showed high expression in roots. TeAP2/ERF028, TeAP2/ERF050, and TeAP2/ERF137, were highly expressed in the stems, while TeAP2/ERF015 and TeAP2/ERF114 exhibited high expression in leaves. Twelve genes, including TeAP2/ERF002, TeAP2/ERF009, TeAP2/ERF014, and TeAP2/ERF035 etc., were highly expressed in receptacles. In bracts, TeAP2/ERF026, TeAP2/ERF042, and TeAP2/ERF111 showed notably high expression levels. Interestingly, TeAP2/ERF038, TeAP2/ERF103, and TeAP2/ERF137 showed particularly high expression in seeds, suggesting potential roles in seed development or maturation.
Given the distinctive floral architecture of marigold, characterized by peripheral ray florets and central disc florets, we further investigated the expression patterns of AP2 family genes in different floral organs of these two floret types (Figure 6, Figure 7 and Figure S3). In the corollas, several genes showed differential expression between the two floret types. Notably, TeAP2/ERF114 and TeAP2/ERF145 exhibited significantly higher expression in ray floret corollas compared to disc floret corollas. This elevated expression suggests their potential role in determining the enlarged, showy corolla morphology characteristic of ray florets. In disc florets, which uniquely possess stamens, several genes showed high stamen-specific expression. TeAP2/ERF036 showed the highest expression in stamens among all floral organs, while TeAP2/ERF103 also exhibited strong expression in this male reproductive organ. The stamen-specific expression of these genes suggests their involvement in stamen development or function, which may contribute to the fertility of disc florets.
When comparing the same organs between ray and disc florets, we observed both consistent and differential expression patterns. In pappus, TeAP2/ERF015, TeAP2/ERF137, and TeAP2/ERF151 were highly expressed in both floret types, while TeAP2/ERF042 and TeAP2/ERF111 showed notably higher expression in ray floret pappus. In pistils, TeAP2/ERF035 and TeAP2/ERF154 exhibited high expression in both floret types, indicating conserved functions in pistil development. For ovaries, TeAP2/ERF103 and TeAP2/ERF151 showed consistently high expression in both floret types, suggesting their roles in ovule or seed development. Several genes showed broader expression patterns across multiple floral organs. TeAP2/ERF114 exhibited high expression in both pappus and corollas of ray florets, consistent with the A-class function in the floral organ identity model. Similarly, TeAP2/ERF100, TeAP2/ERF103, and TeAP2/ERF111 were expressed across multiple organs in both floret types, indicating their potential involvement in general floral development processes.
3.6. Protein Interaction Network Analysis of AP2 Family, MADS-Box, and TCP Family Transcription Factors
Based on our expression profile analysis results, which revealed distinct expression patterns and significant differential expression of the 28 AP2 family genes across tissues and between ray and disc florets, we specifically focused on these genes for protein–protein interaction network based on homology with Arabidopsis proteins (Figure 8). The resulting network revealed extensive predicted interactions between AP2 family proteins and members of the MADS-box and TCP transcription factors. AP2 emerged as a central hub in this network, showing strong interactions with multiple MADS-box proteins. The TeAP2/ERF026, TeAP2/ERF035, TeAP2/ERF114, and TeAP2/ERF154 identified in our study correspond to AP2 in Arabidopsis, highlighting their potential conservation of function.
These interactions suggest potential cooperative functions between AP2 family and MADS-box proteins in regulating flower development in marigold. Notably, TCP1, a homolog of CYC2-like genes which are known to regulate floral symmetry in Asteraceae, was also integrated into this network. The predicted interactions between AP2 family, MADS-box, and TCP family proteins provide molecular evidence for the potential cross-talk between these three transcription factor families in controlling floral organ specification and floret identity in marigold. The strength of these predicted interactions, indicated by the thickness of connecting lines in the network, suggests varying degrees of functional association, with particularly strong connections observed between AP2 protein and several MADS-box proteins (including AP3, AGL8, and SOC1).
4. Discussion
4.1. Evolutionary Relationships and Structural Features of AP2/ERF Genes in Marigold
In this study, we conducted a comprehensive genome-wide analysis of the AP2/ERF superfamily members in marigold, an important ornamental plant with distinctive heterogamous capitulum inflorescence, and investigated their potential roles in floral development, particularly in the differentiation of ray and disc florets.
We identified 177 AP2/ERF genes in the marigold genome, which is comparable to the number found in other plant species such as Arabidopsis (145) [7], rice (170) [59], sunflower (288) [60], lettuce (228), and tomato (134) [61]. The classification of TeAP2/ERF genes into different families based on domain architecture and phylogenetic analysis is consistent with previous studies in other plants [7,8]. Notably, the total number of ERF subfamily and DREB subfamily genes (143) far exceeds the number of AP2 family genes (28), and the predominance of ERF family genes (including both ERF and DREB subfamilies) aligns with findings in other species [3,62,63]. This pattern highlights the evolutionary conservation and functional significance of the ERF family members across plant lineages and reflects the central role of ERF and DREB genes in diverse stress response pathways [11]. Comparison between marigold and other Asteraceae species like lettuce reveals interesting patterns in the composition of AP2/ERF superfamily genes (Table 1). While the total number of AP2/ERF superfamily genes in marigold (177) is somewhat lower than in lettuce (228), the proportional distribution across families shows similarities, particularly in the expansion of the ERF family relative to the AP2 family. These observations suggest conserved evolutionary patterns within the Asteraceae family, possibly reflecting shared adaptive pressures.
The gene structure and motif composition analysis revealed distinct patterns among different subfamilies, reflecting their evolutionary divergence and functional specialization. The AP2 family genes typically contained four to eight exons, while most ERF and DREB subfamily genes had only one to two exons. This structural difference has been observed across multiple plant species [56,60] and likely represents an ancient divergence in the evolution of the AP2/ERF superfamily members. The presence of multiple introns in AP2 family genes could provide additional opportunities for alternative splicing and post-transcriptional regulation, potentially contributing to their functional complexity in developmental processes [64]. These family-specific motifs likely contribute to the functional specificity of different AP2/ERF proteins and may serve as transcriptional activation domains, nuclear localization signals, or protein–protein interaction sites [65,66,67]. A comparative analysis of gene structure between marigold and lettuce revealed striking conservation patterns within the AP2/ERF superfamily members [43]. In both Asteraceae species, AP2 family members consistently displayed more complex gene architecture with multiple introns (ranging from four to eight in marigold and four to thirteen in lettuce). By contrast, the vast majority of ERF family members exhibited a markedly simpler structure, predominantly lacking introns entirely. This pronounced structural dichotomy, preserved across these evolutionarily related species, provides compelling evidence for an ancient divergence event that separated the AP2 family members from the ERF family members long before the diversification of the Asteraceae family.
Our analysis of chromosomal distribution and duplication events revealed that segmental duplications have played a major role in the expansion of the AP2/ERF superfamily genes in marigold [68]. The prevalence of segmental duplications suggests that whole-genome or large-scale chromosomal duplication events have been particularly important in the evolutionary history of marigold’s AP2/ERF superfamily genes. The collinearity relationships in the marigold AP2/ERF superfamily genes are primarily driven by purifying selection, while elevated Ka/Ks ratios in specific gene pairs may reflect functional specialization or adaptive evolutionary events that require further verification through functional experiments.
The synteny analysis between marigold and other plant species provided evolutionary evidence for the conservation and divergence of AP2/ERF genes. The higher number of syntenic pairs between marigold and lettuce (93 pairs), compared to marigold and Arabidopsis (49 pairs), reflects their closer evolutionary relationship as members of the Asteraceae family. This observation aligns with previous studies showing higher synteny among closely related species [69]. The preferential conservation of syntenic relationships in certain families, particularly the AP2 and ERF families, suggests stronger functional constraints on these families during evolution, possibly due to their roles in essential developmental and adaptive processes.
4.2. Diverse Expression Patterns Reveal Functional Specialization of AP2 Family Genes in Marigold
The expression analysis of AP2 family genes revealed diverse spatial and temporal expression patterns, suggesting their involvement in various developmental processes in marigold. The tissue-specific expression patterns observed for many AP2 family genes provide valuable clues about their potential functions. The high expression of TeAP2/ERF030, TeAP2/ERF036, TeAP2/ERF048, TeAP2/ERF138, and TeAP2/ERF154 in roots suggests potential roles in root development or root-specific responses, similar to functions reported for certain AP2 genes in Arabidopsis [70] and rice [71]. Similarly, the seed-specific expression of TeAP2/ERF038, TeAP2/ERF130, and TeAP2/ERF137 indicates their potential roles in seed development or maturation, consistent with the functions of AP2 genes in seed development in other plants [72,73,74].
Several genes showed specific expression patterns in Asteraceae-specific structures such as the receptacle and bracts, which are crucial components of the capitulum [75]. TeAP2/ERF002, TeAP2/ERF009, TeAP2/ERF014, TeAP2/ERF035, etc., were highly expressed in receptacles, suggesting their potential roles in the development of this specialized structure that supports the florets. Similarly, TeAP2/ERF026, TeAP2/ERF042, and TeAP2/ERF111 showed high expression in bracts, indicating their possible involvement in the development of these protective structures surrounding the capitulum.
One of the most intriguing findings of our study is the differential expression of several AP2 family genes between ray and disc florets, suggesting their potential roles in determining floret identity or organ development in the complex inflorescence of marigold. The expression of TeAP2/ERF036 and TeAP2/ERF103 specifically in stamens of disc florets, while they are absent in sterile ray florets, suggests their potential roles in stamen development or function. This pattern could be relevant to understanding the evolutionary origin of ray florets, which are thought to have evolved from disc florets through the suppression of stamen development and corolla fusion [19]. Further functional studies of these genes could provide valuable insights into the molecular mechanisms underlying the evolution of the complex capitulum structure in Asteraceae. The preferential expression of TeAP2/ERF111, TeAP2/ERF114, and TeAP2/ERF145 in both pappus and corollas of disc and ray florets aligns with the canonical A-class function in the ABCDE model of floral organ development [76]. However, the additional expression of TeAP2/ERF111 in pistils of ray florets suggests an expanded functional role beyond the classical A-function typically limited to pappus and corollas [77]. The significantly higher expression of TeAP2/ERF145 in ray floret corollas compared to disc floret corollas is particularly interesting, as it may contribute to the development of the enlarged, showy corollas characteristic of ray florets in marigold. Similarly, the differential expression of TeAP2/ERF103, TeAP2/ERF114, and TeAP2/ERF13 between ray and disc florets suggests their potential involvement in establishing the distinct morphological and functional identities of these floret types [78].
4.3. Integration of AP2 Family, MADS-Box, and CYC2 Protein Regulatory Networks in Marigold Floral Development
The complex floral architecture of the Asteraceae family, characterized by the capitulum inflorescence with distinct ray and disc florets, is likely regulated by intricate genetic networks involving multiple transcription factor families. Our findings on the expression patterns of AP2 family genes, together with previous studies on MADS-box [39] and CYC2 [37] genes in marigold, provide an opportunity to explore the potential integration of these transcription factor networks in controlling floral development (Figure 9). Previous studies in marigold have identified and characterized several MADS-box genes involved in floral organ development, including A-class genes (TeAP1/FUL-like), B-class genes (TeAP3/PI-like), C/D-class genes (TeAG-like), and E-class genes (TeSEP-like). These genes showed distinct expression patterns in ray and disc florets, with some genes exhibiting floret type-specific expression. For instance, TeMADS49, TeMADS45, and TeMADS54 (A-class genes) were highly expressed in pappus of ray and disc florets, similar to the expression patterns we observed for TeAP2/ERF151. This suggests potential co-regulation or interaction between A-class MADS-box genes and A-class AP2 genes in specifying pappus identities in both floret types.
Interestingly, some MADS-box genes, such as TeMADS10 and TeMADS33 (A-class), were specifically expressed in the stamens of disc florets, similar to the expression of TeAP2/ERF103 in our study. This parallel expression pattern suggests a potential cooperative role of these genes in regulating stamen development specifically in disc florets.
Similarly, CYC2 genes in marigold have been shown to play crucial roles in regulating floral symmetry and organ development, particularly in the differentiation of ray and disc florets. Seven CYC2 genes were identified in marigold (TeCYC2a, TeCYC2b, TeCYC2c, TeCYC2d, TeCYC2e1, TeCYC2e2, and TeCYC2g). Among these, TeCYC2b, TeCYC2c, and TeCYC2d were specifically expressed in ray floret corollas, exhibiting an expression pattern similar to the ray floret-specific expression of TeAP2/ERF114 and TeAP2/ERF145 observed in our study. Notably, functional analysis revealed that TeCYC2c inhibits stamen development when overexpressed in Arabidopsis, providing a molecular explanation for the characteristic absence of stamens in marigold ray florets. The expression patterns of TeAP2 genes and TeCYC2 genes in marigold florets suggest potential interactions or regulatory relationships between these transcription factor families. For instance, the ray floret-specific expression of TeCYC2b, TeCYC2c, and TeCYC2d, together with the higher expression of TeAP2/ERF145 in ray floret corollas, suggests a possible genetic interaction in specifying ray floret identity and corolla development. Given the distinctive expression pattern of TeAP2/ERF145 and its potential role in determining ray floret morphology, this gene deserves particular attention in future functional studies to elucidate its precise regulatory mechanisms. Similarly, the expression of TeCYC2a, TeCYC2e1, and TeCYC2e2 in disc floret corollas, contrasted with the notably low expression of all TeAP2 family genes in this tissue, suggests a distinctive regulatory mechanism governing floret differentiation. The absence of TeAP2 family genes expression in disc floret corollas likely contributes to the development of their characteristic tubular morphology. This expression pattern indicates that while CYC family members are active in disc florets, the lack of AP2 family member activity may be a critical factor in determining disc floret identity. Conversely, the potential expression of TeAP2 family genes in ray florets could be instrumental in establishing their distinctive asymmetric structure and elongated ligule formation, highlighting how differential transcription factor activity may drive the morphological diversification of floret types within the composite inflorescence of Asteraceae.
Our protein–protein interaction network analysis further supports this integrated regulatory model, revealing extensive predicted interactions between AP2 family, MADS-box, and TCP family proteins (Figure 8). The emergence of AP2 as a central hub in this network, with strong connections to multiple MADS-box proteins (including AP3, AGL8, and SOC1) and TCP1 (a homolog of CYC2-like genes), provides molecular evidence for the coordinated function of these transcription factor families in floral development. The varying strengths of these predicted interactions suggest differential regulatory relationships, which may contribute to the distinctive development of ray and disc florets in marigold.
In other Asteraceae species, interactions between MADS-box, CYC2, and AP2 family members have been reported to regulate floral development. In chrysanthemum, the A-class gene ClAP2 was found to particularly accumulate in ray floret corollas and form heterodimers with CYC2-like proteins, suggesting their cooperative role in ray floret identity and development [29,79]. Similarly, in gerbera, a regulatory network involving CINCINNATA-like TCP proteins functions upstream of GhCYC to specify ray flower identity. The E-class MADS-box transcription factor GRCD5 was shown to activate GhCYC3 expression during ligule development, while the C-class MADS-box transcription factor GAGA1 appears to contribute to stamen development upstream of GhCYC3 [30]. These findings, together with our results on differential expression of AP2 family genes in marigold florets, demonstrate how intricate genetic networks involving AP2 family, MADS-box, and CYC2 family members collectively regulate the complex floral architecture of the Asteraceae family, with particular importance in controlling heteromorphic flower type identity and organ development.
In the context of the ABCDE model of floral organ development, our results suggest both conservation and divergence of AP2 family member functions in marigold compared to model plants. The expression patterns of several AP2 family genes in sepals and petals are consistent with the A-class function in the classical model. However, their expression in additional organs suggests potential neofunctionalization or subfunctionalization in the complex floral architecture of marigold. This functional divergence might be related to the evolution of the unique capitulum structure in Asteraceae, which represents one of the most successful adaptations in flowering plants [80].
5. Conclusions
In this study, we conducted a comprehensive genome-wide analysis of the AP2/ERF superfamily members in marigold, identifying 177 members classified into four families: AP2 (28), ERF (143), RAV (4), and Soloist (2). These genes exhibited family-specific structural features, conserved motifs, and chromosomal distributions, with significant expansion through segmental duplications under strong purifying selection. Expression profiling revealed tissue-specific and developmental stage-specific patterns, with several TeAP2 family genes showing differential expression between ray and disc florets that may contribute to the distinctive floral architecture of marigold. Protein interaction network analysis suggested cross-talk among AP2 family, MADS-box, and TCP family proteins, providing insights into how these regulatory networks collectively shape the complex capitulum structure. Our findings establish a foundation for future functional studies of TeAP2/ERF superfamily members in marigold and contribute to understanding the molecular mechanisms underlying floral development in Asteraceae.
Conceptualization, H.L., G.C., M.B. and Y.H.; methodology, H.L.; software, G.C. and C.L.; validation, H.L. and G.C.; formal analysis, H.L., G.C. and S.H.; investigation, H.L., G.C. and S.H.; resources, H.L. and G.C.; data curation, H.L., G.C. and C.L.; writing—original draft preparation, H.L. and G.C.; writing—review and editing, M.B. and Y.H.; visualization, G.C.; supervision, M.B. and Y.H.; project administration, M.B. and Y.H.; funding acquisition, Y.H. All authors have read and agreed to the published version of the manuscript.
The original data underlying this study are available within the article and its
The authors declare no conflicts of interest.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Figure 1 Phylogenetic analysis and classification of AP2/ERF superfamily members in marigold and Arabidopsis. The maximum-likelihood phylogenetic tree was constructed using full-length protein sequences of AP2/ERF superfamily members from marigold (Tagetes erecta) and Arabidopsis (Arabidopsis thaliana). The different colored arcs represent different families: AP2 family (AP2, ANT, and AIL, purple shades), RAV family (brown), Soloist family (orange), DREB subfamily (A1–A6, blue shades), and ERF subfamily (B1–B6, green shades). Marigold proteins are labeled as TeAP2/ERF followed by numbers, while Arabidopsis proteins are shown with their original gene names. Red and green spots indicate AP2/ERF superfamily members from marigold and Arabidopsis, respectively.
Figure 2 Chromosomal distribution of AP2/ERF superfamily genes in marigold. The 177 AP2/ERF genes identified in this study are distributed across 11 of the 12 marigold chromosomes (TeORChr1-TeORChr11). Chromosome numbers are labeled on the left side of each chromosome, with their physical lengths presented in megabases (Mb) using the scale bar on the left. Gene IDs are labeled at their corresponding chromosome positions according to their actual locations. Gene density for each chromosome was calculated using 500 kb genetic intervals, and a gradient color scale was used to characterize gene density distribution, where the red to blue gradient corresponds to high-density to low-density regions (red: gene-rich regions; blue: gene-sparse regions). Blank areas represent genomic regions lacking gene distribution information.
Figure 3 Segmental duplications of TeAP2/ERF genes in marigold. Circular representation of the 56 segmental duplication pairs of TeAP2/ERF genes mapped across the 11 marigold chromosomes (Chr1–Chr11). Chromosomes are depicted as colored arcs with gene density represented by the red histogram track along each chromosome. TeAP2/ERF genes are labeled at their respective chromosome positions. Red curved lines connect segmentally duplicated AP2/ERF gene pairs, revealing extensive inter-chromosomal duplication events, while gray lines represent other duplication gene pairs of non-AP2/ERF genes across the genome. The inner heat map bands on each chromosome illustrate gene density with a blue to red gradient (representing low to high density).
Figure 4 Synteny analysis of AP2/ERF genes between marigold and related species. Syntenic relationships of AP2/ERF genes between (A) marigold and Arabidopsis and (B) marigold and lettuce. Chromosomes of Arabidopsis (Chr1–Chr5) and lettuce (Chr1–Chr9) are represented by orange bars, while green bars represent the corresponding marigold (TeChr1–TeChr12) AP2/ERF genes. Red curves connect syntenic gene pairs between the species, with the intensity of red indicating the strength of syntenic relationships. The grey lines represent other collinear gene pairs of non-AP2/ERF genes across genomes.
Figure 5 Expression profiles of AP2 family genes during flower bud development in marigold. Expression intensity is represented by circle size and color depth, with larger, darker circles indicating higher expression levels. The genes are clustered based on their expression patterns.
Figure 6 Expression profiles of AP2 family genes in different marigold tissues. The color scale represents normalized expression levels from low (light pink) to high (dark burgundy). Hierarchical clustering of genes (left) groups those with similar expression patterns.
Figure 7 Differential expression patterns of AP2 family genes in floral organs of ray and disc florets in marigold. To clarify the expression characteristics of AP2 family genes in marigold inflorescences, this study selected 12 candidate genes with normalized expression values > 0.3 in at least one tissue region of either ray florets or disc florets for visualization analysis according to
Figure 8 Protein–protein interaction network of AP2 family, MADS-box, and TCP family transcription factors based on homology with Arabidopsis. Nodes represent transcription factors, with node size and color intensity (from yellow to deep red) indicating the degree of connectivity. Line thickness represents the confidence level of predicted interactions, with thicker lines indicating stronger interaction scores. The TeAP2 labels (TeAP2XXX) represent marigold AP2 family transcription factors corresponding to their Arabidopsis homologs in the network.
Figure 9 Expression patterns of AP2 family members alongside MADS-box and CYC2 transcription factor family members in different organs of marigold. The schematic shows tissue-specific expression of TeAP2 family genes (labeled as AP2XXX) in comparison with MADS-box (AP1, FUL, AP3/PI, AG, SEP, AGL6, AGL15, AGL17, BS, MIKC*) and CYC2 (CYC2a, CYC2b, CYC2c, CYC2d, CYC2e1, CYC2e2, CYC2g) family genes across different marigold organs. Genes highly expressed in each tissue are grouped in colored ovals, with the size of the ovals representing the number of genes (1–6) as indicated in the legend.
Comparative analysis of AP2/ERF superfamily member distribution in marigold, lettuce, and Arabidopsis.
Family | Subfamily | Group | Tagetes erecta | Lactuca sativa | Arabidopsis |
---|---|---|---|---|---|
RAV | 4 | 4 | 6 | ||
Soloist | 2 | 2 | 1 | ||
AP2 | AP2 | 18 | 20 | 10 | |
ANT | 4 | 3 | 1 | ||
AIL | 6 | 6 | 7 | ||
ERF | DREB | A1 | 8 | 14 | 6 |
A2 | 3 | 12 | 8 | ||
A3 | 1 | 1 | 1 | ||
A4 | 24 | 29 | 16 | ||
A5 | 16 | 14 | 15 | ||
A6 | 6 | 11 | 10 | ||
ERF | B1 | 19 | 19 | 15 | |
B2 | 6 | 6 | 5 | ||
B3 | 37 | 45 | 17 | ||
B4 | 3 | 9 | 7 | ||
B5 | 2 | 10 | 8 | ||
B6 | 18 | 23 | 12 | ||
Total | 177 | 228 * | 145 |
* The set of 228 genes analyzed in the comparative section was identified from the lettuce (Lactuca sativa) genome, specifically from the “Genome assembly Lsat_Salinas_v8” (GCA_002870075.3), which is available at the National Center for Biotechnology Information (NCBI) database (
Supplementary Materials
The following supporting information can be downloaded at:
1. Riechmann, J.L.; Heard, J.; Martin, G.; Reuber, L.; Jiang, C.Z.; Keddie, J.; Adam, L.; Pineda, O.; Ratcliffe, O.J.; Samaha, R.R.
2. Liu, L.; White, M.J.; Macrae, T.H. Transcription factors and their genes in higher plants. Eur. J. Biochem.; 1999; 262, pp. 247-257. [DOI: https://dx.doi.org/10.1046/j.1432-1327.1999.00349.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10336605]
3. Nakano, T.; Suzuki, K.; Fujimura, T.; Shinshi, H. Genome-wide analysis of the ERF gene family in Arabidopsis and rice. Plant Physiol.; 2006; 140, pp. 411-432. [DOI: https://dx.doi.org/10.1104/pp.105.073783] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16407444]
4. Licausi, F.; Ohme-Takagi, M.; Perata, P. APETALA2/Ethylene Responsive Factor (AP2/ERF) transcription factors: Mediators of stress responses and developmental programs. New Phytol.; 2013; 199, pp. 639-649. [DOI: https://dx.doi.org/10.1111/nph.12291]
5. Buttner, M.; Singh, K.B. Arabidopsis thaliana ethylene-responsive element binding protein (AtEBP), an ethylene-inducible, GCC box DNA-binding protein interacts with an ocs element binding protein. Proc. Natl. Acad. Sci. USA; 1997; 94, pp. 5961-5966. [DOI: https://dx.doi.org/10.1073/pnas.94.11.5961]
6. Ohme-Takagi, M.; Shinshi, H. Ethylene-inducible DNA binding proteins that interact with an ethylene-responsive element. Plant Cell; 1995; 7, pp. 173-182. [DOI: https://dx.doi.org/10.1105/tpc.7.2.173]
7. Sakuma, Y.; Liu, Q.; Dubouzet, J.G.; 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, pp. 998-1009. [DOI: https://dx.doi.org/10.1006/bbrc.2001.6299]
8. Zhuang, J.; Cai, B.; Peng, R.H.; Zhu, B.; Jin, X.F.; Xue, Y.; Gao, F.; Fu, X.Y.; Tian, Y.S.; Zhao, W.
9. Krizek, B. AINTEGUMENTA and AINTEGUMENTA-LIKE6 act redundantly to regulate Arabidopsis floral growth and patterning. Plant Physiol.; 2009; 150, pp. 1916-1929. [DOI: https://dx.doi.org/10.1104/pp.109.141119]
10. Jose, R.J.; Bailey, L.J.; Mai, Q.A.; Wu, S.L.; Hon, C.T.; Chapman, E.J.; Ditta, G.S.; Estelle, M.; Yanofsky, M.F. microRNA regulation of fruit growth. Nat. Plants; 2015; 1, 15036. [DOI: https://dx.doi.org/10.1038/nplants.2015.36]
11. Mizoi, J.; Shinozaki, K.; Yamaguchi-Shinozaki, K. AP2/ERF family transcription factors in plant abiotic stress responses. Biochim. Biophys. Acta; 2012; 1819, pp. 86-96. [DOI: https://dx.doi.org/10.1016/j.bbagrm.2011.08.004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21867785]
12. Huang, P.Y.; Catinot, J.; Zimmerli, L. Ethylene response factors in Arabidopsis immunity. J. Exp. Bot.; 2016; 67, pp. 1231-1241. [DOI: https://dx.doi.org/10.1093/jxb/erv518] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26663391]
13. Matias-Hernandez, L.; Aguilar-Jaramillo, A.E.; Marin-Gonzalez, E.; Suarez-Lopez, P.; Pelaz, S. RAV genes: Regulation of floral induction and beyond. Ann. Bot.; 2014; 114, pp. 1459-1470. [DOI: https://dx.doi.org/10.1093/aob/mcu069]
14. Li, X.; Yang, R.; Liang, Y.; Gao, B.; Li, S.; Bai, W.; Oliver, M.J.; Zhang, D. The ScAPD1-like gene from the desert moss Syntrichia caninervis enhances resistance to Verticillium dahliae via phenylpropanoid gene regulation. Plant J.; 2023; 113, pp. 75-91. [DOI: https://dx.doi.org/10.1111/tpj.16035] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36416176]
15. Sun, Z.M.; Zhou, M.L.; Dan-Wang,; Tang, Y.X.; Lin, M.; Wu, Y.M. Overexpression of the Lotus corniculatus Soloist Gene LcAP2/ERF107 Enhances Tolerance to Salt Stress. Protein Pept. Lett.; 2016; 23, pp. 442-449. [DOI: https://dx.doi.org/10.2174/0929866523666160322152914]
16. Roeble, L.; van Benthem, K.J.; Weigelt, P.; Kreft, H.; Knope, M.L.; Mandel, J.R.; Vargas, P.; Etienne, R.S.; Valente, L. Island biogeography of the megadiverse plant family Asteraceae. Nat. Commun.; 2024; 15, 7276. [DOI: https://dx.doi.org/10.1038/s41467-024-51556-7]
17. Panero, J.L.; Crozier, B.S. Macroevolutionary dynamics in the early diversification of Asteraceae. Mol. Phylogenet. Evol.; 2016; 99, pp. 116-132. [DOI: https://dx.doi.org/10.1016/j.ympev.2016.03.007]
18. Mandel, J.R.; Dikow, R.B.; Siniscalchi, C.M.; Thapa, R.; Watson, L.E.; Funk, V.A. A fully resolved backbone phylogeny reveals numerous dispersals and explosive diversifications throughout the history of Asteraceae. Proc. Natl. Acad. Sci. USA; 2019; 116, pp. 14083-14088. [DOI: https://dx.doi.org/10.1073/pnas.1903871116]
19. Zhang, T.; Elomaa, P. Development and evolution of the Asteraceae capitulum. New Phytol.; 2024; 242, pp. 33-48. [DOI: https://dx.doi.org/10.1111/nph.19590]
20. Elomaa, P.; Zhao, Y.; Zhang, T. Flower heads in Asteraceae-recruitment of conserved developmental regulators to control the flower-like inflorescence architecture. Hortic. Res.; 2018; 5, 36. [DOI: https://dx.doi.org/10.1038/s41438-018-0056-8]
21. Funk, V.; Susanna, A.; Stuessy, T.; Bayer, R. Systematics, Evolution, and Biogeography of Compositae; International Association for Plant Taxonomy: Vienna, Austria, 2009; ISBN 978-3-9501754-3-1
22. Chen, J.; Shen, C.Z.; Guo, Y.P.; Rao, G.Y. Patterning the Asteraceae Capitulum: Duplications and Differential Expression of the Flower Symmetry CYC2-Like Genes. Front. Plant Sci.; 2018; 9, 551. [DOI: https://dx.doi.org/10.3389/fpls.2018.00551] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29922305]
23. Tahtiharju, S.; Rijpkema, A.S.; Vetterli, A.; Albert, V.A.; Teeri, T.H.; Elomaa, P. Evolution and diversification of the CYC/TB1 gene family in Asteraceae—A comparative study in Gerbera (Mutisieae) and sunflower (Heliantheae). Mol. Biol. Evol.; 2012; 29, pp. 1155-1166. [DOI: https://dx.doi.org/10.1093/molbev/msr283]
24. Cheng, H.; Yu, Y.; Zhai, Y.; Wang, L.; Wang, L.; Chen, S.; Chen, F.; Jiang, J. An ethylene-responsive transcription factor and a B-box protein coordinate vegetative growth and photoperiodic flowering in chrysanthemum. Plant Cell Environ.; 2023; 46, pp. 440-450. [DOI: https://dx.doi.org/10.1111/pce.14488] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36367211]
25. Broholm, S.K.; Tahtiharju, S.; Laitinen, R.A.; Albert, V.A.; Teeri, T.H.; Elomaa, P. A TCP domain transcription factor controls flower type specification along the radial axis of the Gerbera (Asteraceae) inflorescence. Proc. Natl. Acad. Sci. USA; 2008; 105, pp. 9117-9122. [DOI: https://dx.doi.org/10.1073/pnas.0801359105]
26. Fambrini, M.; Pugliesi, C. CYCLOIDEA 2 Clade Genes: Key Players in the Control of Floral Symmetry, Inflorescence Architecture, and Reproductive Organ Development. Plant Mol. Biol. Rep.; 2017; 35, pp. 20-36. [DOI: https://dx.doi.org/10.1007/s11105-016-1005-z]
27. Laitinen, R.A.; Broholm, S.; Albert, V.A.; Teeri, T.H.; Elomaa, P. Patterns of MADS-box gene expression mark flower-type development in Gerbera hybrida (Asteraceae). BMC Plant Biol.; 2006; 6, 11. [DOI: https://dx.doi.org/10.1186/1471-2229-6-11]
28. Pu, Y.; Liao, M.; Li, J.; Tian, Y.; Wang, Z.; Song, X.; Dai, S. Floral Development Stage-Specific Transcriptomic Analysis Reveals the Formation Mechanism of Different Shapes of Ray Florets in Chrysanthemum. Genes; 2023; 14, 766. [DOI: https://dx.doi.org/10.3390/genes14030766] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36981036]
29. Wen, X.; Qi, S.; Huang, H.; Wu, X.; Zhang, B.; Fan, G.; Yang, L.; Hong, Y.; Dai, S. The expression and interactions of ABCE-class and CYC2-like genes in the capitulum development of Chrysanthemum lavandulifolium and C. × morifolium. Plant Growth Regul.; 2019; 88, pp. 205-214. [DOI: https://dx.doi.org/10.1007/s10725-019-00491-5]
30. Zhao, Y.; Broholm, S.K.; Wang, F.; Rijpkema, A.S.; Lan, T.; Albert, V.A.; Teeri, T.H.; Elomaa, P. TCP and MADS-Box Transcription Factor Networks Regulate Heteromorphic Flower Type Identity in Gerbera hybrida. Plant Physiol.; 2020; 184, pp. 1455-1468. [DOI: https://dx.doi.org/10.1104/pp.20.00702]
31. Sang, Q.; Vayssieres, A.; O’Maoileidigh, D.S.; Yang, X.; Vincent, C.; Bertran, G.D.O.E.; Cerise, M.; Franzen, R.; Coupland, G. MicroRNA172 controls inflorescence meristem size through regulation of APETALA2 in Arabidopsis. New Phytol.; 2022; 235, pp. 356-371. [DOI: https://dx.doi.org/10.1111/nph.18111]
32. Lu, R.; Hu, S.; Feng, J.; Liu, Z.; Kang, C. The AP2 transcription factor BARE RECEPTACLE regulates floral organogenesis via auxin pathways in woodland strawberry. Plant Cell; 2024; 36, pp. 4970-4987. [DOI: https://dx.doi.org/10.1093/plcell/koae270] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/39367407]
33. Jofuku, K.D.; den Boer, B.G.; Van Montagu, M.; Okamuro, J.K. Control of Arabidopsis flower and seed development by the homeotic gene APETALA2. Plant Cell; 1994; 6, pp. 1211-1225. [DOI: https://dx.doi.org/10.1105/tpc.6.9.1211] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/7919989]
34. Walia, S.; Kumar, R. Wild marigold (Tagetes minuta L.) an important industrial aromatic crop: Liquid gold from the Himalaya. J. Essent. Oil Res.; 2020; 32, pp. 373-393. [DOI: https://dx.doi.org/10.1080/10412905.2020.1813211]
35. Vasudevan, P.; Kashyap, S.; Sharma, S. Tagetes: A multipurpose plant. Bioresour. Technol.; 1997; 62, pp. 29-35. [DOI: https://dx.doi.org/10.1016/S0960-8524(97)00101-6]
36. Ai, Y.; Zhang, Q.; Wang, W.; Zhang, C.; Cao, Z.; Bao, M.; He, Y. Transcriptomic Analysis of Differentially Expressed Genes during Flower Organ Development in Genetic Male Sterile and Male Fertile Tagetes erecta by Digital Gene-Expression Profiling. PLoS ONE; 2016; 11, e0150892. [DOI: https://dx.doi.org/10.1371/journal.pone.0150892]
37. Zhu, Y.; Liu, Y.; Wang, W.; Li, H.; Liu, C.; Dou, L.; Wei, L.; Cheng, W.; Bao, M.; Yi, Q.
38. Xin, H.; Ji, F.; Wu, J.; Zhang, S.; Yi, C.; Zhao, S.; Cong, R.; Zhao, L.; Zhang, H.; Zhang, Z. Chromosome-scale genome assembly of marigold (Tagetes erecta L.): An ornamental plant and feedstock for industrial lutein production. Hortic. Plant J.; 2023; 9, pp. 1119-1130. [DOI: https://dx.doi.org/10.1016/j.hpj.2023.04.001]
39. Liu, C.; Wang, F.; Li, R.; Zhu, Y.; Zhang, C.; He, Y. Marigold (Tagetes erecta) MADS-Box Genes: A Systematic Analysis and Their Implications for Floral Organ Development. Agronomy; 2024; 14, 1889. [DOI: https://dx.doi.org/10.3390/agronomy14091889]
40. Ai, Y.; Zhang, C.; Sun, Y.; Wang, W.; He, Y.; Bao, M. Characterization and Functional Analysis of Five MADS-Box B Class Genes Related to Floral Organ Identification in Tagetes erecta. PLoS ONE; 2017; 12, e0169777. [DOI: https://dx.doi.org/10.1371/journal.pone.0169777]
41. Zhang, C.; Sun, Y.; Yu, X.; Li, H.; Bao, M.; He, Y. Functional Conservation and Divergence of Five AP1/FUL-like Genes in Marigold (Tagetes erecta L.). Genes; 2021; 12, 2011. [DOI: https://dx.doi.org/10.3390/genes12122011]
42. Zhang, C.; Wei, L.; Wang, W.; Qi, W.; Cao, Z.; Li, H.; Bao, M.; He, Y. Identification, characterization and functional analysis of AGAMOUS subfamily genes associated with floral organs and seed development in Marigold (Tagetes erecta). BMC Plant Biol.; 2020; 20, 439. [DOI: https://dx.doi.org/10.1186/s12870-020-02644-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32967618]
43. Park, S.; Shi, A.; Meinhardt, L.W.; Mou, B. Genome-wide characterization and evolutionary analysis of the AP2/ERF gene family in lettuce (Lactuca sativa). Sci. Rep.; 2023; 13, 21990. [DOI: https://dx.doi.org/10.1038/s41598-023-49245-4] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/38081919]
44. Wang, Y.; Zhang, J.; Hu, Z.; Guo, X.; Tian, S.; Chen, G. Genome-Wide Analysis of the MADS-Box Transcription Factor Family in Solanum lycopersicum. Int. J. Mol. Sci.; 2019; 20, 2961. [DOI: https://dx.doi.org/10.3390/ijms20122961]
45. Li, X.; Li, W.; Huo, J.; Li, L.; Chen, B.; Guo, Z.; Ma, Z. Identification and expression analysis of citrate synthase 3 gene family members in apple. Chin. J. Biotechnol.; 2024; 40, pp. 137-149. [DOI: https://dx.doi.org/10.13345/j.cjb.230166]
46. Reiser, L.; Subramaniam, S.; Zhang, P.; Berardini, T. Using the Arabidopsis Information Resource (TAIR) to Find Information About Arabidopsis Genes. Curr. Protoc.; 2022; 2, e574. [DOI: https://dx.doi.org/10.1002/cpz1.574]
47. Minh, B.Q.; Schmidt, H.A.; Chernomor, O.; Schrempf, D.; Woodhams, M.D.; von Haeseler, A.; Lanfear, R. IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Mol. Biol. Evol.; 2020; 37, pp. 1530-1534. [DOI: https://dx.doi.org/10.1093/molbev/msaa015]
48. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.; von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat. Methods; 2017; 14, pp. 587-589. [DOI: https://dx.doi.org/10.1038/nmeth.4285]
49. Hoang, D.T.; Chernomor, O.; von Haeseler, A.; Minh, B.Q.; Vinh, L.S. UFBoot2: Improving the Ultrafast Bootstrap Approximation. Mol. Biol. Evol.; 2018; 35, pp. 518-522. [DOI: https://dx.doi.org/10.1093/molbev/msx281]
50. Letunic, I.; Bork, P. Interactive Tree of Life (iTOL) v6: Recent updates to the phylogenetic tree display and annotation tool. Nucleic Acids Res.; 2024; 52, pp. W78-W82. [DOI: https://dx.doi.org/10.1093/nar/gkae268]
51. Zhao, Y.; Wang, T.; Wan, S.; Tong, Y.; Wei, Y.; Li, P.; Hu, N.; Liu, Y.; Chen, H.; Pan, X.
52. Hwa, K.; Subramani, B.; Lenka, G.; Chang, S. Sequence Motif Analysis of Phycobiliproteins in Cyanobacterial Genome. Proceedings of the 2011 International Conference on Bioscience, Biochemistry and Bioinformatics; Singapore, 26–28 February 2011.
53. Zhang, D.; Zeng, B.; He, Y.; Li, J.; Yu, Z. Genome-wide identification and comparative analysis of the AP2/ERF gene family in Prunus dulcis and Prunus tenella: Expression of PdAP2/ERF genes under freezing stress during dormancy. BMC Genom.; 2025; 26, 95. [DOI: https://dx.doi.org/10.1186/s12864-025-11275-9]
54. Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods; 2001; 25, pp. 402-408. [DOI: https://dx.doi.org/10.1006/meth.2001.1262] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11846609]
55. Ahmed, S.; Rashid, M.A.R.; Zafar, S.A.; Azhar, M.T.; Waqas, M.; Uzair, M.; Rana, I.A.; Azeem, F.; Chung, G.; Ali, Z.
56. 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.; 2015; 6, 1247. [DOI: https://dx.doi.org/10.3389/fpls.2015.01247] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26834762]
57. Krizek, B.A. Aintegumenta and Aintegumenta-Like6 regulate auxin-mediated flower development in Arabidopsis. BMC Res. Notes; 2011; 4, 176. [DOI: https://dx.doi.org/10.1186/1756-0500-4-176]
58. Yant, L.; Mathieu, J.; Dinh, T.T.; Ott, F.; Lanz, C.; Wollmann, H.; Chen, X.; Schmid, M. Orchestration of the floral transition and floral development in Arabidopsis by the bifunctional transcription factor APETALA2. Plant Cell; 2010; 22, pp. 2156-2170. [DOI: https://dx.doi.org/10.1105/tpc.110.075606]
59. 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. Bioinform. Online; 2012; 8, pp. 321-355. [DOI: https://dx.doi.org/10.4137/EBO.S9369]
60. Najafi, S.; Sorkheh, K.; Nasernakhaei, F. Characterization of the APETALA2/Ethylene-responsive factor (AP2/ERF) transcription factor family in sunflower. Sci. Rep.; 2018; 8, 11576. [DOI: https://dx.doi.org/10.1038/s41598-018-29526-z]
61. Yang, H.; Sun, Y.; Wang, H.; Zhao, T.; Xu, X.; Jiang, J.; Li, J. Genome-wide identification and functional analysis of the ERF2 gene family in response to disease resistance against Stemphylium lycopersici in tomato. BMC Plant Biol.; 2021; 21, 72. [DOI: https://dx.doi.org/10.1186/s12870-021-02848-3]
62. Song, X.; Li, Y.; Hou, X. Genome-wide analysis of the AP2/ERF transcription factor superfamily in Chinese cabbage (Brassica rapa ssp. pekinensis). BMC Genom.; 2013; 14, 573. [DOI: https://dx.doi.org/10.1186/1471-2164-14-573]
63. Girardi, C.L.; Rombaldi, C.V.; Dal Cero, J.; Nobile, P.M.; Laurens, F.; Bouzayen, M.; Quecini, V. Genome-wide analysis of the AP2/ERF superfamily in apple and transcriptional evidence of ERF involvement in scab pathogenesis. Sci. Hortic.; 2013; 151, pp. 112-121. [DOI: https://dx.doi.org/10.1016/j.scienta.2012.12.017]
64. Li, Z.; Wang, S.; Cheng, J.; Su, C.; Zhong, S.; Liu, Q.; Fang, Y.; Yu, Y.; Lv, H.; Zheng, Y.
65. Fujimoto, S.Y.; Ohta, M.; Usui, A.; Shinshi, H.; Ohme-Takagi, M. Arabidopsis ethylene-responsive element binding factors act as transcriptional activators or repressors of GCC box-mediated gene expression. Plant Cell; 2000; 12, pp. 393-404. [DOI: https://dx.doi.org/10.1105/tpc.12.3.393] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10715325]
66. 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, pp. 470-478. [DOI: https://dx.doi.org/10.1093/nar/27.2.470]
67. Nole-Wilson, S.; Krizek, B.A. DNA binding properties of the Arabidopsis floral development protein AINTEGUMENTA. Nucleic Acids Res.; 2000; 28, pp. 4076-4082. [DOI: https://dx.doi.org/10.1093/nar/28.21.4076]
68. Cannon, S.B.; Mitra, A.; Baumgarten, A.; Young, N.D.; May, G. The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol.; 2004; 4, 10. [DOI: https://dx.doi.org/10.1186/1471-2229-4-10]
69. Cheng, F.; Wu, J.; Fang, L.; Wang, X. Syntenic gene analysis between Brassica rapa and other Brassicaceae species. Front. Plant Sci.; 2012; 3, 198. [DOI: https://dx.doi.org/10.3389/fpls.2012.00198] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22969786]
70. Ye, B.B.; Shang, G.D.; Pan, Y.; Xu, Z.G.; Zhou, C.M.; Mao, Y.B.; Bao, N.; Sun, L.; Xu, T.; Wang, J.W. AP2/ERF Transcription Factors Integrate Age and Wound Signals for Root Regeneration. Plant Cell; 2020; 32, pp. 226-241. [DOI: https://dx.doi.org/10.1105/tpc.19.00378]
71. Kitomi, Y.; Ito, H.; Hobo, T.; Aya, K.; Kitano, H.; Inukai, Y. The auxin responsive AP2/ERF transcription factor CROWN ROOTLESS5 is involved in crown root initiation in rice through the induction of OsRR1, a type-A response regulator of cytokinin signaling. Plant J.; 2011; 67, pp. 472-484. [DOI: https://dx.doi.org/10.1111/j.1365-313X.2011.04610.x]
72. Jofuku, K.D.; Omidyar, P.K.; Gee, Z.; Okamuro, J.K. Control of seed mass and seed yield by the floral homeotic gene APETALA2. Proc. Natl. Acad. Sci. USA; 2005; 102, pp. 3117-3122. [DOI: https://dx.doi.org/10.1073/pnas.0409893102]
73. Ohto, M.A.; Fischer, R.L.; Goldberg, R.B.; Nakamura, K.; Harada, J.J. Control of seed mass by APETALA2. Proc. Natl. Acad. Sci. USA; 2005; 102, pp. 3123-3128. [DOI: https://dx.doi.org/10.1073/pnas.0409858102] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15708976]
74. Yu, B.; He, X.; Tang, Y.; Chen, Z.; Zhou, L.; Li, X.; Zhang, C.; Huang, X.; Yang, Y.; Zhang, W.
75. Bello, M.A.; Alvarez, I.; Torices, R.; Fuertes-Aguilar, J. Floral development and evolution of capitulum structure in Anacyclus (Anthemideae, Asteraceae). Ann. Bot.; 2013; 112, pp. 1597-1612. [DOI: https://dx.doi.org/10.1093/aob/mcs301] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23287557]
76. Wollmann, H.; Mica, E.; Todesco, M.; Long, J.A.; Weigel, D. On reconciling the interactions between APETALA2, miR172 and AGAMOUS with the ABC model of flower development. Development; 2010; 137, pp. 3633-3642. [DOI: https://dx.doi.org/10.1242/dev.036673]
77. Krogan, N.T.; Hogan, K.; Long, J.A. APETALA2 negatively regulates multiple floral organ identity genes in Arabidopsis by recruiting the co-repressor TOPLESS and the histone deacetylase HDA19. Development; 2012; 139, pp. 4180-4190. [DOI: https://dx.doi.org/10.1242/dev.085407]
78. Kerstens, M.; Schranz, M.E.; Bouwmeester, K. Phylogenomic analysis of the APETALA2 transcription factor subfamily across angiosperms reveals both deep conservation and lineage-specific patterns. Plant J.; 2020; 103, pp. 1516-1524. [DOI: https://dx.doi.org/10.1111/tpj.14843]
79. Ding, L.; Song, A.; Zhang, X.; Li, S.; Su, J.; Xia, W.; Zhao, K.; Zhao, W.; Guan, Y.; Fang, W.
80. Zhang, G.; Yang, J.; Zhang, C.; Jiao, B.; Panero, J.L.; Cai, J.; Zhang, Z.R.; Gao, L.M.; Gao, T.; Ma, H. Nuclear phylogenomics of Asteraceae with increased sampling provides new insights into convergent morphological and molecular evolution. Plant Commun.; 2024; 5, 100851. [DOI: https://dx.doi.org/10.1016/j.xplc.2024.100851]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
The APETALA2/Ethylene-Responsive Factor (AP2/ERF) superfamily is one of the largest transcription factor families in plants, playing diverse roles in development, stress response, and metabolic regulation. Despite their ecological and economic importance, AP2/ERF genes remain uncharacterized in marigold (Tagetes erecta), a valuable ornamental and medicinal plant in the Asteraceae family known for its unique capitulum-type inflorescence with distinct ray and disc florets. Here, we conducted a comprehensive genome-wide analysis of the AP2/ERF superfamily in marigold and identified 177 AP2/ERF genes distributed across 11 of the 12 chromosomes. Phylogenetic analysis revealed their classification into the AP2 (28 genes), ERF (143 genes), RAV (4 genes), and Soloist (2 genes) families based on domain architecture. Gene structure and motif composition analyses demonstrated group-specific patterns that correlated with their evolutionary relationships. Chromosome mapping and synteny analyses revealed that segmental duplications significantly contributed to AP2/ERF superfamily gene expansion in marigold, with extensive collinearity observed between marigold and other species. Expression profiling across different tissues and developmental stages indicated distinct spatio-temporal expression patterns, with several genes exhibiting tissue-specific expression in Asteraceae-specific structures. In floral organs, TeAP2/ERF145 exhibited significantly higher expression in ray floret corollas compared to disc florets, while TeAP2/ERF103 showed stamen-specific expression in disc florets. Protein interaction network analysis revealed AP2 as a central hub with extensive predicted interactions with MADS-box and TCP family proteins. These findings suggest that AP2 family genes may collaborate with MADS-box and CYC2 genes in regulating the characteristic floral architecture of marigold, establishing a foundation for future functional studies and molecular breeding efforts to enhance ornamental and agricultural traits in this economically important plant.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer