Expression profile analysis of early fruit development in iaaM-parthenocarpic tomato plants
BMC Research Notes volume 2, Article number: 143 (2009)
Fruit normally develops from the ovary after pollination and fertilization. However, the ovary can also generate seedless fruit without fertilization by parthenocarpy. Parthenocarpic fruit development has been obtained in tomato (Solanum lycopersicum) by genetic modification using auxin-synthesising gene(s) (DefH9-iaaM; DefH9-RI-iaaM) expressed specifically in the placenta and ovules.
We have performed a cDNA Amplified Fragment Length Polymorphism (cDNA-AFLP) analysis on pre-anthesis tomato flower buds (0.5 cm long) collected from DefH9-iaaM and DefH9-RI-iaaM parthenocarpic and wild-type plants, with the aim to identify genes involved in very early phases of tomato fruit development. We detected 212 transcripts differentially expressed in auxin-ipersynthesising pre-anthesis flower buds, 65 of them (31%) have unknown function. Several differentially expressed genes show homology to genes involved in protein trafficking and protein degradation via proteasome. These processes are crucial for auxin cellular transport and signaling, respectively.
The data presented might contribute to elucidate the molecular basis of the fruiting process and to develop new methods to confer parthenocarpy to species of agronomic interest. In a recently published work, we have demonstrated that one of the genes identified in this screening, corresponding to #109 cDNA clone, regulates auxin-dependent fruit initiation and its suppression causes parthenocarpic fruit development in tomato.
Tomato (Solanum lycopersicum) fruit represents an important component of the human diet due to high content in fibres, vitamins and antioxidants. The identification of genes that control growth and maturation of tomato fruit will allow their manipulation, by breeding and/or genetic engineering, to improve fruit quality. Biochemical and genetic aspects of late stages of fruit development, in particular ripening, have been widely investigated [1–4], while initial phases of fruit development have received less attention, despite their importance for both basic and applied research. Until now, few studies have been performed to investigate tomato fruit set and early development (1 to 15 days post anthesis) [5–7].
In Angiosperms, once a flower is pollinated and fertilization successfully takes place, ovary starts to grow and this is the first visible sign of fruit development . The earliest phase of fruit growth is referred to as fruit set or fruit initiation.
In parthenocarpic plants, fruit set and development occurs without fertilization leading to the production of seedless fruits. The first phytohormone shown to trigger parthenocarpic fruit development was auxin applied exogenously to tomato flowers . It is extensively demonstrated that one of the possible methods for achieving parthenocarpic fruit development, employing genetic engineering, is based on the ovary-specific expression (driven by DefH9 promoter) of iaaM and RI-iaaM genes, which code for an enzyme of the auxin biosynthetic pathway [10–13]. The two chimeric genes, DefH9-iaaM and DefH9-RI-iaaM, differ in vitro in their translational potential, and in planta in the level of auxin (total IAA) content in flower buds . RI-iaaM is an iaaM derivative, modified in its 5'ULR in order to down-regulate the level of expression. Both DefH9-iaaM and DefH9-RI-iaaM parthenocarpic tomato flower buds are characterised by an increased content of auxin in female gametophyte as compared with wild-type [10–12], but DefH9-iaaM-expressing flower buds contain 5 times higher IAA than DefH9-RI-iaaM-expressing flower buds . At pre-anthesis the iaaM-parthenocarpic flower buds appear morphologically identical to those of wild-type, except for an early enlargement of the ovary , indicating that high levels of auxin can derepress ovary growth before fertilization. Therefore, iaaM-parthenocarpic plants show precocious fruit growth. It is reasonable to assume that in parthenocarpic flower buds, the genetic program for fruit development has been already switched on during early stages of flower development, well before anthesis. On this premises, the iaaM-parthenocarpic flower buds represent a suitable experimental model to study molecular events taking place during the early phases of fruit growth, monitoring expression changes that occur mainly within the ovary.
Results and discussion
A cDNA-AFLP approach  was used to generate expression profiles of young flower buds isolated from four independent parthenocarpic UC82 tomato lines, two transgenic for DefH9-iaaM (#3 and #2; see  for a description of the lines) and two for DefH9-RI-iaaM gene (#s5 and #s6; see  for a description of the lines), and from wild-type plants. Flower buds under analysis (0.5 cm in length) were at a very early stage of development (6–7 days before anthesis).
Using 32 different primer combinations (BstT/C+n – Mse+n, where n represents selective nucleotide) more than 3000 cDNA fragments were generated. The AFLP fragments ranged in length from 50 to 500 bp; for each primer combination, approximately 100 bands were observed on polyacrylamide denaturing gel (Figure 1). The differentially expressed fragments were excised from the gels, re-amplified by PCR and sequenced. Good quality and unique sequences were obtained for 212 cDNA fragments (Additional files 1 and 2).
The expression profiles of the isolated fragments were reproducible in the independent transgenic lines chosen for the analysis (Figure 1). One hundred and thirty cDNAs out of the 212 differentially expressed were induced in iaaM-parthenocarpic lines as compared to wild-type, whilst the remaining 82 cDNAs were down-regulated.
We performed quantitative RT-PCR analysis to validate some of the differentially expressed gene (Additional file 3) representative of the different functional categories (Additional file 2), using RNA extracted from independent parthenocarpic lines different from those tested by cDNA-AFLP. The sequences of the 212 differentially expressed genes were compared with those in GenBank database http://www.ncbi.nlm.nih.gov, DNA Data Bank of Japan database http://www.ddbj.nig.ac.jp/searches-e.html and DFCI Tomato Unique Gene Indices http://www.tigr.org/tdb/tgi, using BlastN or BlastX homology search tool . The differentially expressed cDNA-AFLP clones are listed in additional file 2; for each clone the accession number, if available, and the e-value of the sequence that produced the highest identity score, together with the best S. lycopersicum Tentative Consensus (TC) sequence (DFCI Tomato Unique Gene Indices), are reported. A tentative annotation was assigned to about 70% of the cDNA-AFLP clones following GenBank and DFCI Tomato Unique Gene Indices databases. The cDNA clones were grouped in 10 functional categories (Figure 2) according to Universal Protein Resource (UniProt) http://www.uniprot.org and TAIR database http://www.arabidopsis.org.
The two most represented categories "metabolism" and "unknown function" contain 47 genes (22%) and 65 genes (31%), respectively (Figure 2 and additional file 2). The majority of genes of the first category include enzymes with primary metabolic roles (e.g. phospholipid, carbohydrate and amino-acid metabolism). Other genes belonging to this category code for enzymes involved in secondary metabolism, such as spermidine synthase (#376) and ornithine decarboxylase (#216) that are up-regulated in iaaM-parthenocarpic flower buds and participate in the polyamine biosynthetic pathway. This finding is consistent with previous data showing that polyamines play a role in early fruit development  and it is in accordance with the experimental evidence demonstrating that polyamine metabolism is affected by exogenous auxin treatment of tomato ovaries .
In the functional category "Cell rescue, Defense and Cell fate" most of the genes involved in biotic stress showed a decreased expression in parthenocarpic flower buds. This is in accordance with the observation that tomato unpollinated flowers are actively protected against pathogens whereas after induction of fruit growth the expression of defence genes is reduced .
iaaM-parthenocarpic tomato flower buds have higher auxin levels compared with wild-type flower buds . The increased indole-3-acetic acid (IAA) content in ovaries and placenta is considered to mimic the burst of auxin, which occurs in the carpel after pollination . Interestingly, one of the genes in the category "metabolism", cDNA clone #436 induced in iaaM-parthenocarpic flower buds, codes for a molybdenum cofactor sulfurase protein-like, an enzyme taking part in the molybdenum cofactor biosynthetic pathway . The molybdenum cofactor sulfurase ABA3 from Arabidopsis thaliana is necessary for post-translational activation of aldehyde oxidase (AOs) and xanthine dehydrogenase (XDH) . In plants, AOs form a multigene family whose members have broad substrate specificity for several aldehydes including indole-3-acetaldehyde and abscisic aldehyde suggesting for AOs a role in the biosynthesis of IAA and abscisic acid (ABA). Other experimental evidences linked AO activity to the biosynthesis of IAA. AO1 activity in IAA-overproducing Arabidopsis plants is 5 times higher as compared to wild-type plants . The up-regulation of #436 might contribute to increase IAA content observed in iaaM-parthenocarpic flower buds. However, we cannot exclude that the modulation of this gene might be related to ABA biosynthetic activity. In this regard, genes involved in ABA biosynthesis were shown to be strongly expressed in tomato ovary collected from flowers at anthesis .
The molecular events underlying the auxin-mediated derepression of ovary growth have not been fully elucidated. The augmented synthesis of auxin could affect IAA signaling pathway and/or modify the expression and localization of auxin transporters .
Auxin is perceived by intracellular located TIR1/AFB1–3 receptors that mediate the ubiquitination of AUX/IAA proteins that are then degraded by 26S-proteasome [23–25]. Most of the differentially expressed genes present in the "protein fate" functional category show homology to genes involved in proteolytic degradation via ubiquitination.
The establishment of auxin gradients in tissue/organs is crucial for auxin action [26–28] and requires the polar distribution of transporters in the plasma membrane. There is experimental evidence proving that the flexible asymmetric localization of IAA transporters is assured by cycling of these proteins between plasma membranes and endosome compartments [26–28], however many aspects of this regulatory process are still elusive.
In this regard, in the category "Plasmamembrane and cellular transport" a number of differentially expressed transcripts are homologous to genes implicated in intracellular protein transport. cDNA-AFLP clone #163 is similar to Rab8 of Nicotiana tabacum, putatively involved in protein movement between Golgi apparatus and the plasma membrane; #438 is homologous to SEC61 a protein translocator of endoplasmic reticulum (RE); #877 is similar to synthaxin-71, a component of the SNARE complex which functions in vesicles trafficking; #408 is annotated as μ-subunit of the clathrin adaptor complex. Furthermore, #303 clone encodes for a protein with not yet identified function, characterised by the clathrin-box motif (L [IVLMF]X [IVLMF] [DE]) (Eukaryotic Linear Motif (ELM) http://elm.eu.org). It would be interesting to investigate whether the differentially expressed genes found to be implicated in intracellular protein trafficking, play a role in the distribution of auxin transporters during early stage of fruit development. In fact, the treatment of tomato flower buds with inhibitors of polar auxin transport induces parthenocarpic fruit development .
Genes involved in fruit set and/or early fruit development are often transcriptionally regulated during the first phases of ovary growth. We have investigated in wild-type tomato plants the expression of some genes identified in this screening (cDNA clones #70, #855, #805 reported in Figure 3; #109, see ; #216 and #904, see ) at different stages of flower and fruit development. These genes show a drastic increase (#855, and #805) or decrease (#70, #109, #216, and #904) in transcript level at the stage of open flower. Considering their expression profiles, it would be reasonable to test these genes for their putative role in fruit development. For instance, cDNA clone #109 (named by us SlAucsia-1 gene; ) and the homologous gene SlAucsia-2 are highly expressed in flower buds and decrease dramatically (97%) after pollination/fertilization. Their suppression by RNA silencing causes parthenocarpic fruit development in tomato . Furthermore, Aucsia- silenced tomato plants exhibit other alterations such as reduced polar auxin transport in roots and increased sensitivity to 1-naphthylphthalamic acid, an inhibitor of polar auxin transport. In Arabidopsis thaliana Aucsia genes [GenBank:AK224828, GenBank:AK224647] are annotated as components of the endomembrane system. This finding argues in favour of the hypothesis that also other genes here presented can be candidate genes playing a role in fruit set.
In conclusion, this dataset represents a starting tool for the study of genes involved in fruit set in tomato and in other crops bearing fleshy fruits and could contribute to develop new methods to confer parthenocarpy.
The two transgenes DefH9-iaaM and DefH9-RI-iaaM were introduced in UC82 tomato plants. The UC82 is a typical cultivar used by the processing industry. Two independent DefH9-iaaM parthenocarpic UC82 lines (#2 and # 3) and two independent DefH9-RI-iaaM parthenocarpic UC82 lines (#s5 and #s6) were used in this work (for a description of the lines see ).
For a detailed description of the method used, see .
qRT-PCR Quantitative PCR Analysis
100 mg of pooled flower buds (collected from 4–5 plants) were ground in liquid nitrogen and total RNA was isolated. Starting from DNase-treated total RNA, first-strand cDNA was synthesized with oligo-dT primer and Superscript II (Invitrogen). The cDNA clones were amplified with gene-specific primers designed to give amplification products ranging from 100 to 150 bp, according to Applied Biosystem guidelines. Experiments were carried out using Platinum SYBR Green QPCR Supermix-UDG (Invitrogen) in ABI Prism 7000 Sequence Detection System (Applied Biosystems). The following cycling conditions were used: 2 min at 50°C, 2 min at 95°C, 40 cycles of 95°C for 30 sec, 56°C for 30 sec, 72°C for 30 sec and finally 72°C for 3 min. All quantitations were normalized to actin gene as endogenous control gene. Forward (F) and reverse (R) primers used for actin amplification are the following: F 5'-CCCGTTCAGCAGTGGTGGT-3' and R 5'-TACGAGGGTTATGCTTTGCC-3'. For each amplification reaction, the analysis of the product dissociation curve was performed to exclude the presence of non-specific amplification. Data from qRT-PCR experiments were analysed according to .
Northern Blot Analysis
Total RNAs were isolated with Trizol reagent (Invitrogen) and then 20 μg were separated on 1% agarose-formaldehyde denaturing gels. The gels were blotted overnight on Hybond N+ membrane (Amersham Biosciences) in 10× SSC. The DNA probes were labeled with [32P]CTP using "Ready to go DNA labeling beads (-dCTP)" (Amersham Biosciences). The membranes were hybridised overnight at 42°C in ULTRAhyb buffer (Ambion) in presence of 106 cpm mL-1 of labelled probe. The membranes were washed 2 times in 2× SSC containing 0.1% SDS for 5 min and 2 times in 0.1× SSC containing 0.1% SDS for 15 min at 42°C. Autoradiography was then performed using Kodak X-AR5 film. For each c-DNA AFLP clone tested, the probe was obtained by amplification of cDNA with the same primers adopted for qRT-PCR.
Moore S, Vrebalov J, Payton P, Giovannoni J: Use of genomics tools to isolate key ripening genes and analyse fruit maturation in tomato. Journal of Experimental Botany. 2002, 53: 2023-2030. 10.1093/jxb/erf057.
Giovannoni JJ: Genetic Regulation of Fruit Development and Ripening. Plant Cell. 2004, 16: S170-S180. 10.1105/tpc.019158.
Fei Z, Tang X, Alba RM, White JA, Ronning CM, Martin GB, Tanksley SD, Giovannoni JJ: Comprehensive EST analysis of tomato and comparative genomics of fruit ripening. Plant Journal. 2004, 40: 47-59. 10.1111/j.1365-313X.2004.02188.x.
Giovannoni JJ: Fruit ripening mutants yield insights into ripening control. Current Opinion in Plant Biology. 2007, 10: 283-289. 10.1016/j.pbi.2007.04.008.
Testa G, Caccia R, Tilesi F, Soressi G, Mazzucato A: Sequencing and characterization of tomato genes putatively involved in fruit set and early development. Sexual Plant Reproduction. 2002, 14: 269-277. 10.1007/s00497-001-0122-4.
Lemaire-Chamley M, Petit J, Garcia V, Just D, Baldet P, Germain V, Fagard M, Mouassite M, Cheniclet C, Rothan C: Changes in transcriptional profiles are associated with early fruit tissue specialization in tomato. Plant Physiology. 2005, 139: 750-769. 10.1104/pp.105.063719.
Vriezen WH, Feron R, Maretto F, Keijman J, Mariani C: Changes in tomato ovary transcriptome demonstrate complex hormonal regulation of fruit set. New Phytologist. 2008, 177: 60-76.
Gillaspy G, Ben-David H, Gruissem W: Fruits: A developmental perspective. Plant Cell. 1993, 5: 1439-1451. 10.1105/tpc.5.10.1439.
Gustafson F: Auxin distribution in fruits and its significance in fruit development. American Journal of Botany. 1939, 26: 189-194. 10.2307/2436487.
Rotino GL, Perri E, Zottini M, Sommer H, Spena A: Genetic engineering of parthenocarpic plants. Nature Biotechnology. 1997, 15: 1398-1401. 10.1038/nbt1297-1398.
Ficcadenti N, Sestili S, Pandolfini T, Cirillo C, Rotino GL, Spena A: Genetic engineering of parthenocarpic fruit development in tomato. Molecular Breeding. 1999, 5: 463-470. 10.1023/A:1009665409959.
Pandolfini T, Rotino GL, Camerini S, Defez R, Spena A: Optimization of transgene action at the post-transcriptional level: high quality parthenocarpic fruits in industrial tomatoes. BMC Biotechnology. 2002, 2: 1-10.1186/1472-6750-2-1.
Pandolfini T, Molesini B, Spena A: Parthenocarpy in crop plants. Fruit development and seed dispersal. Edited by: Ostergaard L. 2009, Oxford, Wiley Blackwell, 38: 326-345.
Molesini B, Pandolfini T, Rotino GL, Dani V, Spena A: Aucsia gene silencing causes parthenocarpic fruit development in tomato. Plant Physiology. 2009, 149: 534-548. 10.1104/pp.108.131367.
Breyne P, Zabeau M: Genome-wide expression analysis of plant cell cycle modulated genes. Current Opinion in Plant Biology. 2001, 4: 136-142. 10.1016/S1369-5266(00)00149-7.
Altschul SF, Madden TL, Schaffer AA, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Research. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
Egea-Cortines M, Mizrahi Y: Polyamines in cells division, fruit set and development, and seed germination. Biochemistry and Physiology of Polyamines in Plants. Edited by: Slocum RD, Flores HE. 1991, Boca Raton, FL, CRC Press, 143-158.
Alabadí D, Agüero MS, Pérez-Amador MA, Carbonell J: Arginase, arginine decarboxylase, ornithine decarboxylase, and polyamines in tomato ovaries: changes in unpollinated ovaries and parthenocarpic fruits induced by auxin or gibberellin. Plant Physiology. 1996, 112: 1237-1244.
Pandolfini T, Molesini B, Spena A: Molecular dissection of the role of auxin in fruit initiation. Trends in Plant Science. 2007, 12: 327-329. 10.1016/j.tplants.2007.06.011.
Mendel RR, Hänsch R: Molybdoenzymes and molybdenum cofactor in plants. Journal of Experimantal Botany. 2002, 53: 1689-1698. 10.1093/jxb/erf038.
Wollers S, Heidenreich T, Zarepour M, Zachmann D, Kraft C, Zhao Y, Mendel RR, Bittner F: Binding of sulfurated molybdenum cofactor to the C-terminal domain of ABA3 from Arabidopsis thaliana provides insight into the mechanism of molybdenum cofactor sulfuration. Journal of Biological Chemistry. 2008, 283: 9642-9650. 10.1074/jbc.M708549200.
Seo M, Akaba S, Oritani T, Delarue M, Bellini C, Caboche M, Koshiba T: Higher activity of an aldehyde oxidase in the auxin-overproducing superroot1 mutant of Arabidopsis thaliana. Plant Physiology. 1998, 116: 687-693. 10.1104/pp.116.2.687.
Leyser O: Dynamic integration of auxin transport and signalling. Current Biology. 2006, 16: 424-433. 10.1016/j.cub.2006.05.014.
Woodward AW, Bartel B: Auxin: regulation, action and interaction. Annals of Botany. 2005, 95: 707-735. 10.1093/aob/mci083.
Tan X, Calderon-Villalobos LIA, Sharon M, Zheng C, Robinson CV, Estelle M, Zheng N: Mechanism of auxin perception by the TIR1 ubiquitin ligase. Nature. 2007, 446: 640-645. 10.1038/nature05731.
Vieten A, Sauer M, Brewer PB, Friml J: Molecular and cellular aspects of auxin-transport-mediated development. Trends in Plant Science. 2007, 12: 160-168. 10.1016/j.tplants.2007.03.006.
Kleine-Vehn J, Dhonukshe R, Swarup P, Bennet M, Friml J: Subcellular trafficking of the Arabidopsis auxin influx carrier AUX1 uses a novel pathway distinct from PIN-1. Plant Cell. 2006, 18: 3171-3181. 10.1105/tpc.106.042770.
Dhonukshe P, Aniento F, Hwang I, Robinson DG, Mravec J, Stierhof YD, Friml J: Clathrin-mediated constitutive endocytosis of PIN auxin efflux carriers in Arabidopsis. Current Biology. 2007, 17: 520-527. 10.1016/j.cub.2007.01.052.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.
This work is part of the FIRB project RBAU01JTHS 2001–2004 of the MIUR (Italian Ministry of University and Research).
The authors declare that they have no competing interests.
BM performed cDNA-AFLP, qRT-PCR and Northern blot analysis and drafted the manuscript; GLR is responsible for the production of the transgenic lines; TP carried out bioinformatic analysis and drafted the manuscript; AS conceived the study.
Electronic supplementary material
Additional file 1: Sequences of the 212 differentially expressed cDNA-AFLP clones. A list of sequences corresponding to cDNA-AFLP clones differentially expressed in iaaM-parthenocarpic flower buds. (PDF 32 KB)
Additional file 2: Differentially expressed genes. The differentially expressed cDNA-AFLP clones are listed with their (if available) accession number, e-value of the sequence that produced the highest identity score, and with the best S. lycopersicum Tentative Consensus (TC) sequence (DFCI Tomato Unique Gene Indices). (XLS 64 KB)
Additional file 3: Validation of cDNA-AFLP expression pattern by qRT-PCR in iaaM-parthenocarpic flower buds. The data reported represent quantitative RT-PCR analysis of differentially expressed genes representative of the different functional categories performed using independent parthenocarpic lines different from those tested by cDNA-AFLP. (XLS 21 KB)
About this article
Cite this article
Molesini, B., Rotino, G.L., Spena, A. et al. Expression profile analysis of early fruit development in iaaM-parthenocarpic tomato plants. BMC Res Notes 2, 143 (2009). https://doi.org/10.1186/1756-0500-2-143