Genetic interactions in the control of flowering time and reproductive structure development in tomato (Solanum lycopersicum)

Strigolactones promote flowering by inducing the miR319-LA-SFT module in tomato

Strigolactones are a class of phytohormones with various functions in plant development, stress responses, and in the interaction with (micro)organisms in the rhizosphere. While their effects on vegetative development are well studied, little is known about their role in reproduction. We investigated the effects of genetic and chemical modification of strigolactone levels on the timing and intensity of flowering in tomato (Solanum lycopersicum L.) and the molecular mechanisms underlying such effects. Results showed that strigolactone levels in the shoot, whether endogenous or exogenous, correlate inversely with the time of anthesis and directly with the number of flowers and the transcript levels of the florigen-encoding gene SINGLE FLOWER TRUSS (SFT) in the leaves. Transcript quantifications coupled with metabolite analyses demonstrated that strigolactones promote flowering in tomato by inducing the activation of the microRNA319-LANCEOLATE module in leaves. This, in turn, decreases gibberellin content and increases the transcription of SFT. Several other floral markers and morpho-anatomical features of developmental progression are induced in the apical meristems upon treatment with strigolactones, affecting floral transition and, more markedly, flower development. Thus, strigolactones promote meristem maturation and flower development via the induction of SFT both before and after floral transition, and their effects are blocked in plants expressing a miR319-resistant version of LANCEOLATE. Our study positions strigolactones in the context of the flowering regulation network in a model crop species.

EARLY FLOWERING is a dominant gain-of-function allele of FANTASTIC FOUR 1/2c that promotes early flowering in tomato

Flowering time, an important factor in plant adaptability and genetic improvement, is regulated by various genes in tomato (Solanum lycopersicum). In this study, we characterized a tomato mutant, EARLY FLOWERING (EF), that developed flowers much earlier than its parental control. EF is a dominant gain-of-function allele with a T-DNA inserted 139 bp downstream of the stop codon of FANTASTIC FOUR 1/2c (FAF1/2c). The transcript of SlFAF1/2c was at elevated levels in the EF mutant. Overexpressing SlFAF1/2c in tomato plants phenocopied the early flowering trait of the EF mutant. Knocking out SlFAF1/2c in the EF mutant reverted the early flowering phenotype of the mutant to the normal flowering time of the wild-type tomato plants. SlFAF1/2c promoted the floral transition by shortening the vegetative phase rather than by reducing the number of leaves produced before the emergence of the first inflorescence. The COP9 signalosome subunit 5B (CSN5B) was shown to interact with FAF1/2c, and knocking out CSN5B led to an early flowering phenotype in tomato. Interestingly, FAF1/2c was found to reduce the accumulation of the CSN5B protein by reducing its protein stability. These findings imply that FAF1/2c regulates flowering time in tomato by reducing the accumulation and stability of CSN5B, which influences the expression of SINGLE FLOWER TRUSS (SFT), JOINTLESS (J) and UNIFLORA (UF). Thus, a new allele of SlFAF1/2c was discovered and found to regulate flowering time in tomato.

The auxin-responsive transcription factor SlDOF9 regulates inflorescence and flower development in tomato

Understanding the mechanisms underlying differentiation of inflorescence and flower meristems is essential towards enlarging our knowledge of reproductive organ formation and to open new prospects for improving yield traits. Here, we show that SlDOF9 is a new modulator of floral differentiation in tomato. CRISPR/Cas9 knockout strategy uncovered the role of SlDOF9 in controlling inflorescence meristem and floral meristem differentiation via the regulation of cell division genes and inflorescence architecture regulator LIN. Tomato dof9-KO lines have more flowers in both determinate and indeterminate cultivars and produce more fruit upon vibration-assisted fertilization. SlDOF9 regulates inflorescence development through an auxin-dependent ARF5-DOF9 module that seems to operate, at least in part, differently in Arabidopsis and tomato. Our findings add a new actor to the complex mechanisms underlying reproductive organ differentiation in flowering plants and provide leads towards addressing the diversity of factors controlling the transition to reproductive organs.

Germplasm Resources and Strategy for Genetic Breeding of Lycium Species: A Review

Lycium species (goji), belonging to Solanaceae, are widely spread in the arid to semiarid environments of Eurasia, Africa, North and South America, among which most species have affinal drug and diet functions, resulting in their potential to be a superior healthy food. However, compared with other crop species, scientific research on breeding Lycium species lags behind. This review systematically introduces the present germplasm resources, cytological examination and molecular-assisted breeding progress in Lycium species. Introduction of the distribution of Lycium species around the world could facilitate germplasm collection for breeding. Karyotypes of different species could provide a feasibility analysis of fertility between species. The introduction of mapping technology has discussed strategies for quantitative trait locus (QTL) mapping in Lycium species according to different kinds of traits. Moreover, to extend the number of traits and standardize the protocols of trait detection, we also provide 1,145 potential traits (275 agronomic and 870 metabolic) in different organs based on different reference studies on Lycium, tomato and other Solanaceae species. Finally, perspectives on goji breeding research are discussed and concluded. This review will provide breeders with new insights into breeding Lycium species.

CGFS-type glutaredoxin mutations reduce tolerance to multiple abiotic stresses in tomato

Sessile organisms such as plants have adopted diverse reactive oxygen species (ROS) scavenging mechanisms to mitigate damage under abiotic stress conditions. Though CGFS-type glutaredoxin (GRX) genes are important regulators of ROS homeostasis, each of their functions in crop plants have not yet been well understood. We performed a targeted mutagenesis analysis of four CGFS-type GRXs (SlGRXS14, SlGRXS15, SlGRXS16, and SlGRXS17) in tomato plants (Solanum lycopersicum) using a multiplex clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 system and found that Slgrxs mutants were more sensitive to various abiotic stresses compared with the wild-type tomatoes. Slgrxs15 mutants were embryonic lethal. Single, double, and triple combinations of Slgrxs14, 16, and 17 mutants were examined under heat, chilling, drought, heavy metal toxicity, nutrient deficiency, and short photoperiod stresses. Slgrxs14 and 17 mutants showed hypersensitivity to almost all stresses while Slgrxs16 mutants were affected by chilling stress and showed milder sensitivity to other stresses. Additionally, Slgrxs14 and 17 mutants showed delayed flowering time. Our results indicate that the CGFS-type SlGRXs have specific roles against abiotic stresses, providing valuable resources to develop tomato and, possibly, other crop species that are tolerant to multiple abiotic stresses by genetic engineering.

Transcriptome Analysis of Flower Development and Mining of Genes Related to Flowering Time in Tomato (Solanum lycopersicum)

Flowering is a morphogenetic process in which angiosperms shift from vegetative growth to reproductive growth. Flowering time has a strong influence on fruit growth, which is closely related to productivity. Therefore, research on crop flowering time is particularly important. To better understand the flowering period of the tomato, we performed transcriptome sequencing of early flower buds and flowers during the extension period in the later-flowering “Moneymaker” material and the earlier-flowering “20965” homozygous inbred line, and we analyzed the obtained data. At least 43.92 million clean reads were obtained from 12 datasets, and the similarity with the tomato internal reference genome was 92.86–94.57%. Based on gene expression and background annotations, 49 candidate genes related to flowering time and flower development were initially screened, among which the greatest number belong to the photoperiod pathway. According to the expression pattern of candidate genes, the cause of early flowering of “20965” is predicted. The modes of action of the differentially expressed genes were classified, and the results show that they are closely related to hormone regulation and participated in a variety of life activities in crops. The candidate genes we screened and the analysis of their expression patterns provide a basis for future functional verification, helping to explore the molecular mechanism of tomato flowering time more comprehensively.

Characterization of two TERMINAL FLOWER1 homologs PgTFL1 and PgCENa from pomegranate (Punica granatum L.)

FLOWERING LOCUS T (FT) and TERMINAL FLOWER1/CENTRORADIALIS (TFL1/CEN) are the key regulators of flowering time in plants with FT promoting flowering and TFL1 repressing flowering. TFL1 also controls floral meristem identity and its maintenance. In this study we have characterized two pomegranate (Punica granatum L.) TFL1/CEN-like genes designated as PgTFL1 and PgCENa. The expression of PgTFL1 and PgCENa fluctuated through alternate pruning and flowering cycles, being highly expressed during the vegetative phase (immediately after pruning) and decreasing gradually in the months thereafter such that their lowest levels, especially for PgCENa coincided with the flowering phase. Both the genes are able to functionally suppress the Arabidopsis tfl1-14 mutant flowering defect. Their expression in Arabidopsis resulted in delayed flowering time, increased plant height and leaf number, branches and shoot buds as compared with wild type, suggesting that PgTFL1 and PgCENa are bonafide homologs of TFL1. However, both the genes show distinct expression patterns, being expressed differentially in vegetative shoot apex and floral bud samples. While PgTFL1 expression was low in vegetative shoot apex and high in flower bud, PgCENa expression showed the opposite trend. These results suggest that the two TFL1s in pomegranate may be utilized to control distinct developmental processes, namely repression of flowering by PgCENa and development and growth of the reproductive tissues by PgTFL1 via distinct temporal and developmental regulation of their expression.

Natural and induced loss of function mutations in SlMBP21 MADS-box gene led to jointless-2 phenotype in tomato

Abscission is the mechanism by which plants disconnect unfertilized flowers, ripe fruits, senescent or diseased organs from the plant. In tomato, pedicel abscission is an important agronomic factor that controls yield and post-harvest fruit quality. Two non-allelic mutations, jointless (j) and jointless-2 (j-2), controlling pedicel abscission zone formation have been documented but only j-2 has been extensively used in breeding. J was shown to encode a MADS-box protein. Using a combination of physical mapping and gene expression analysis we identified a positional candidate, Solyc12g038510, associated with j-2 phenotype. Targeted knockout of Solyc12g038510, using CRISPR/Cas9 system, validated our hypothesis. Solyc12g038510 encodes the MADS-box protein SlMBP21. Molecular analysis of j-2 natural variation revealed two independent loss-of-function mutants. The first results of an insertion of a Rider retrotransposable element. The second results of a stop codon mutation that leads to a truncated protein form. To bring new insights into the role of J and J-2 in abscission zone formation, we phenotyped the single and the double mutants and the engineered alleles. We showed that J is epistatic to J-2 and that the branched inflorescences and the leafy sepals observed in accessions harboring j-2 alleles are likely the consequences of linkage drags.

Manipulation of flowering time and branching by overexpression of the tomato transcription factor SlZFP2

Flowering of higher plants is orchestrated by complex regulatory networks through integration of various environmental signals such as photoperiod, temperature, light quality and developmental cues. In Arabidopsis, transcription of the flowering integrator gene FLOWERING LOCUS T (FT) that several flowering pathways converge to is directly regulated by more than ten transcription factors. However, very little is known about the transcriptional regulation of the FT homolog SINGLE FLOWER TRUESS (SFT) in the day-neutral plant tomato (Solanum lycopersicum). Previously, we showed that the zinc finger transcription factor SlZFP2 plays important roles in regulation of seed germination and fruit ripening in tomato and also found that overexpression of SlZFP2 impacted flowering and branching. Here, we characterized in detail the early flowering and high branching phenotypes by overexpression of this transcription factor. Our data showed that overexpression of SlZFP2 accelerated flowering in an SFT-dependent manner as demonstrated by elevated SFT expression in the leaves and the transcription factor's binding ability to SFT promoter in vitro and in vivo. Furthermore, overexpression of the SlZFP2 gene in the sft plants failed to rescue the mutant's late flowering. Through analysis of grafting phenotype, growth response of branches to auxin application and transcriptome profiling by RNA sequencing, we also showed that overexpression of SlZFP2 affected shoot apical dominance through multiple regulatory pathways. Our results suggest that the transcription factor SlZFP2 has potential applications in genetic modification of plant architecture and flowering time for tomato production and other crops as well.

Genetic interactions of the unfinished flower development (ufd) mutant support a significant role of the tomato UFD gene in regulating floral organogenesis

Genetic interactions of UFD gene support its specific function during reproductive development of tomato; in this process, UFD could play a pivotal role between inflorescence architecture and flower initiation genes. Tomato (Solanum lycopersicum L.) is a major vegetable crop that also constitutes a model species for the study of plant developmental processes. To gain insight into the control of flowering and floral development, a novel tomato mutant, unfinished flower development (ufd), whose inflorescence and flowers were unable to complete their normal development was characterized using double mutant and gene expression analyses. Genetic interactions of ufd with mutations affecting inflorescence fate (uniflora, jointless and single flower truss) were additive and resulted in double mutants displaying the inflorescence structure of the non-ufd parental mutant and the flower phenotype of the ufd mutant. In addition, ufd mutation promotes an earlier inflorescence meristem termination. Taken together, both results indicated that UFD is not involved in the maintenance of inflorescence meristem identity, although it could participate in the regulatory system that modulates the rate of meristem maturation. Regarding the floral meristem identity, the falsiflora mutation was epistatic to the ufd mutation even though FALSIFLORA was upregulated in ufd inflorescences. In terms of floral organ identity, the ufd mutation was epistatic to macrocalyx, and MACROCALYX expression was differently regulated depending on the inflorescence developmental stage. These results suggest that the UFD gene may play a pivotal role between the genes required for flowering initiation and inflorescence development (such as UNIFLORA, FALSIFLORA, JOINTLESS and SINGLE FLOWER TRUSS) and those required for further floral organ development such as the floral organ identity genes.

Characterization of vegetative inflorescence (mc-vin) mutant provides new insight into the role of MACROCALYX in regulating inflorescence development of tomato

Inflorescence development is a key factor of plant productivity, as it determines flower number. Therefore, understanding the mechanisms that regulate inflorescence architecture is critical for reproductive success and crop yield. In this study, a new mutant, vegetative inflorescence (mc-vin), was isolated from the screening of a tomato (Solanum lycopersicum L.) T-DNA mutant collection. The mc-vin mutant developed inflorescences that reverted to vegetative growth after forming two to three flowers, indicating that the mutated gene is essential for the maintenance of inflorescence meristem identity. The T-DNA was inserted into the promoter region of the MACROCALYX (MC) gene; this result together with complementation test and expression analyses proved that mc-vin is a new knock-out allele of MC. Double combinations between mc-vin and jointless (j) and single flower truss (sft) inflorescence mutants showed that MC has pleiotropic effects on the reproductive phase, and that it interacts with SFT and J to control floral transition and inflorescence fate in tomato. In addition, MC expression was mis-regulated in j and sft mutants whereas J and SFT were significantly up-regulated in the mc-vin mutant. Together, these results provide new evidences about MC function as part of the genetic network regulating the development of tomato inflorescence meristem.

Capsicum annuum S (CaS) promotes reproductive transition and is required for flower formation in pepper (Capsicum annuum)

The genetic control of the transition to flowering has mainly been studied in model species, while few data are available in crop species such as pepper (Capsicum spp.). To elucidate the genetic control of the transition to flowering in pepper, mutants that lack flowers were isolated and characterized. Genetic mapping and sequencing allowed the identification of the gene disrupted in the mutants. Double mutants and expression analyses were used to characterize the relationships between the mutated gene and other genes controlling the transition to flowering and flower differentiation. The mutants were characterized by a delay in the initiation of sympodial growth, a delay in the termination of sympodial meristems and complete inhibition of flower formation. Capsicum annuum S (CaS), the pepper (Capsicum annuum) ortholog of tomato (Solanum lycopersicum) COMPOUND INFLORESCENCE and petunia (Petunia hybrida) EVERGREEN, was found to govern the mutant phenotype. CaS is required for the activity of the flower meristem identity gene Ca-ANANTHA and does not affect the expression of CaLEAFY. CaS is epistatic over other genes controlling the transition to flowering with respect to flower formation. Comparative homologous mutants in the Solanaceae indicate that CaS has uniquely evolved to have a critical role in flower formation, while its role in meristem maturation is conserved in pepper, tomato and petunia.

Inferring the gene network underlying the branching of tomato inflorescence

The architecture of tomato inflorescence strongly affects flower production and subsequent crop yield. To understand the genetic activities involved, insight into the underlying network of genes that initiate and control the sympodial growth in the tomato is essential. In this paper, we show how the structure of this network can be derived from available data of the expressions of the involved genes. Our approach starts from employing biological expert knowledge to select the most probable gene candidates behind branching behavior. To find how these genes interact, we develop a stepwise procedure for computational inference of the network structure. Our data consists of expression levels from primary shoot meristems, measured at different developmental stages on three different genotypes of tomato. With the network inferred by our algorithm, we can explain the dynamics corresponding to all three genotypes simultaneously, despite their apparent dissimilarities. We also correctly predict the chronological order of expression peaks for the main hubs in the network. Based on the inferred network, using optimal experimental design criteria, we are able to suggest an informative set of experiments for further investigation of the mechanisms underlying branching behavior.

Inflorescence development in tomato: gene functions within a zigzag model

Tomato is a major crop plant and several mutants have been selected for breeding but also for isolating important genes that regulate flowering and sympodial growth. Besides, current research in developmental biology aims at revealing mechanisms that account for diversity in inflorescence architectures. We therefore found timely to review the current knowledge of the genetic control of flowering in tomato and to integrate the emerging network into modeling attempts. We developed a kinetic model of the tomato inflorescence development where each meristem was represented by its "vegetativeness" (V), reflecting its maturation state toward flower initiation. The model followed simple rules: maturation proceeded continuously at the same rate in every meristem (dV); floral transition and floral commitment occurred at threshold levels of V; lateral meristems were initiated with a gain of V (ΔV) relative to the V level of the meristem from which they derived. This last rule created a link between successive meristems and gave to the model its zigzag shape. We next exploited the model to explore the diversity of morphotypes that could be generated by varying dV and ΔV and matched them with existing mutant phenotypes. This approach, focused on the development of the primary inflorescence, allowed us to elaborate on the genetic regulation of the kinetic model of inflorescence development. We propose that the lateral inflorescence meristem fate in tomato is more similar to an immature flower meristem than to the inflorescence meristem of Arabidopsis. In the last part of our paper, we extend our thought to spatial regulators that should be integrated in a next step for unraveling the relationships between the different meristems that participate to sympodial growth.

CaBLIND regulates axillary meristem initiation and transition to flowering in pepper

Plant architecture is a major motif in plant diversity. The shape of the plant is regulated by genes that have been found to have similar or related functions in different species. However, changes in gene regulation or their recruitment to additional developmental pathways contribute to the wide range of plant patterns. Our aim was to unravel the genetic mechanisms governing the unique architecture of pepper (Capsicum annuum) and to determine whether these genetic factors have conserved functions in other plant species. We describe the pepper CaBLIND (CaBL) gene that is orthologous to the tomato (Solanum lycopersicum) BLIND (BL) and to the Arabidopsis thaliana REGULATOR OF AXILLARY MERISTEMS (RAX). We identified two allelic Cabl mutants that show dramatic reduction in axillary meristem initiation. In addition, Cabl exhibits late flowering and ectopic vegetative growth during the reproductive phase. Double-mutant and expression analyses suggest that CaBL functions independently of FASCICULATE, the pepper ortholog of SELF PRUNING in regulating sympodial growth, but is epistatic to FASCICULATE in controlling axillary meristem formation. Furthermore, CaBL operates independently of CaREVOLUTA and CaLATERAL SUPPRESSOR in regulating axillary branching. Our results provide evidence of CaBL's conserved function with BL and RAX genes in regulating axillary meristem initiation early in development. In addition, similar to BL but opposite to RAX, CaBL acts to promote the transition from vegetative to reproductive phase. However, in contrast to BL and RAX, CaBL is co-opted to play a role in suppressing vegetative growth during the reproductive phase in pepper.

Flowering response of the uniflora:blind:self-pruning and jointless:uniflora:self-pruning tomato (Solanum lycopersicum) triple mutants

Tomato (Solanum lycopersicum) is a day-neutral plant with a sympodial growth habit. Triple mutants were produced in order to investigate the genetic interactions in the flowering regulation of this species in the initial and sympodial segments. The jointless:uniflora:self-pruning and uniflora:blind:self-pruning triple mutants were produced by crossing the jointless:uniflora and uniflora:blind double mutants with the uniflora:self-pruning double mutant. The phenotype of the triple mutants was characterized and the expression of the affected genes was studied in the uniflora (uf) mutant through semi-quantitative reverse transcriptase polymerase chain reaction (PCR). The triple mutants produced solitary flowers, as their uf parent, instead of inflorescences. They were both late flowering in the initial segment with a flowering time intermediate between their parental double mutants. The flowering time of the sympodial segments was delayed in the jointless:uniflora:self-pruning triple mutant while the uniflora:blind:self-pruning triple mutant did not initiate sympodial segments. The expression of the studied genes was not markedly affected by the uf mutation. These results suggest that floral transition of the primary shoot and of sympodial segments is regulated differently in tomato. The UNIFLORA (UF) gene acts upstream of the other investigated genes in controlling reproductive structure and flowering time of the initial segment although their expression does not seem to be affected by the uf mutation. In the sympodial segments, the self-pruning determinate phenotype is strengthened by the blind mutation and suppressed by the jointless mutation.

Genetic effects on the biomass partitioning and growth of Pisum and Lycopersicon

We examined a series of eight pea genotypes differing in three naturally occurring allelic mutations, i.e., af (afila), st (stipules reduced), and tl (tendril-less) and three species, five cultivars, and one interspecific hybrid of tomato differing in SP (SELF-PRUNING) allele composition to determine whether different phenotypes ontogenetically express different biomass partitioning patterns compared to the isometric partitioning pattern and an interspecific 3/4 scaling "rule" governing annual growth with respect to body mass. The slopes and "elevations" (i.e., α and log β, respectively) of log-log linear regression curves of bivariate plots of leaf, stem, and root dry mass and of annual growth vs. total body mass were used to assess pattern homogeneity. The annual growth of all pea and tomato phenotypes complied with the 3/4 growth rule. The biomass partitioning patterns of all tomato phenotypes were statistically indistinguishable from the isometric pattern as were those of the pea wild type and three single-mutant genotypes. However, significant departures from the isometric (and pea wild type) biomass allocation pattern were observed for three genotypes, all of which were homozygous for the af allele. These results open the door to explore the heritability and genetics underlying the allometry of biomass partitioning patterns and growth.

The rice StMADS11-like genes OsMADS22 and OsMADS47 cause floral reversions in Arabidopsis without complementing the svp and agl24 mutants

During floral induction and flower development plants undergo delicate phase changes which are under tight molecular control. MADS-box transcription factors have been shown to play pivotal roles during these transition phases. SHORT VEGETATIVE PHASE (SVP) and AGAMOUS LIKE 24 (AGL24) are important regulators both during the transition to flowering and during flower development. During vegetative growth they exert opposite roles on floral transition, acting as repressor and promoter of flowering, respectively. Later during flower development they act redundantly as negative regulators of AG expression. In rice, the orthologues of SVP and AGL24 are OsMADS22, OsMADS47, and OsMADS55 and these three genes are involved in the negative regulation of brassinosteroid responses. In order to understand whether these rice genes have maintained the ability to function as regulators of flowering time in Arabidopsis, complementation tests were performed by expressing OsMADS22 and OsMADS47 in the svp and agl24 mutants. The results show that the rice genes are not able to complement the flowering-time phenotype of the Arabidopsis mutants, indicating that they are biologically inactive in Arabidopsis. Nevertheless, they cause floral reversions, which mimic the SVP and AGL24 floral overexpressor phenotypes. Yeast two-hybrid analysis suggests that these floral phenotypes are probably the consequence of protein interactions between OsMADS22 and OsMADS47 and other MADS-box proteins which interfere with the formation of complexes required for normal flower development.

Interaction study of MADS-domain proteins in tomato

MADS-domain proteins are important transcription factors involved in many biological processes of plants. Interactions between MADS-domain proteins are essential for their functions. In tomato (Solanum lycopersicum), the number of MIKC(c)-type MADS-domain proteins identified has totalled 36, but a large-scale interaction assay is lacking. In this study, 22 tomato MADS-domain proteins were selected from six functionally important subfamilies of the MADS-box gene family, to create the first large-scale tomato protein interaction network. Compared with Arabidopsis and petunia (Petunia hybrida), protein interaction patterns in tomato displayed both conservation and divergence. The majority of proteins that can be identified as putative orthologues exhibited conserved interaction patterns, and modifications were mostly found in genes underlining traits unique to tomato. JOINTLESS and RIN, characterized for their roles in abscission zone development and fruit ripening, respectively, showed enlarged interaction networks in comparison with their Arabidopsis and petunia counterparts. Novel interactions were also found for members of the expanded subfamilies, such as those represented by AP1/FUL and AP3/PI MADS-domain proteins. In search for higher order complexes, TM5 was found to be the preferred bridge among the five SEP-like proteins. Additionally, 16 proteins with the MADS-domain removed were used to assess the role of the MADS-domain in protein-protein interactions. The current work provides important knowledge for further functional and evolutionary study of the MADS-box genes in tomato.