First RNA-seq approach to study fruit set and parthenocarpy in zucchini (Cucurbita pepo L.)

[EN] Background: Zucchini fruit set can be limited due to unfavourable environmental conditions in off-seasons crops that caused ineffective pollination/fertilization. Parthenocarpy, the natural or artificial fruit development without fertilization, has been recognized as an important trait to avoid this problem, and is related to auxin signalling. Nevertheless, differences found in transcriptome analysis during early fruit development of zucchini suggest that other complementary pathways could regulate fruit formation in parthenocarpic cultivars of this species. The development of next-generation sequencing technologies (NGS) as RNA-sequencing (RNA-seq) opens a new horizon for mapping and quantifying transcriptome to understand the molecular basis of pathways that could regulate parthenocarpy in this species. The aim of the current study was to analyze fruit transcriptome of two cultivars of zucchini, a non-parthenocarpic cultivar and a parthenocarpic cultivar, in an attempt to identify key genes involved in parthenocarpy. Results: RNA-seq analysis of six libraries (unpollinated, pollinated and auxin treated fruit in a non-parthenocarpic and parthenocarpic cultivar) was performed mapping to a new version of C. pepo transcriptome, with a mean of 92% success rate of mapping. In the non-parthenocarpic cultivar, 6479 and 2186 genes were differentially expressed (DEGs) in pollinated fruit and auxin treated fruit, respectively. In the parthenocarpic cultivar, 10,497 in pollinated fruit and 5718 in auxin treated fruit. A comparison between transcriptome of the unpollinated fruit for each cultivar has been performed determining that 6120 genes were differentially expressed. Annotation analysis of these DEGs revealed that cell cycle, regulation of transcription, carbohydrate metabolism and coordination between auxin, ethylene and gibberellin were enriched biological processes during pollinated and parthenocarpic fruit set. Conclusion: This analysis revealed the important role of hormones during fruit set, establishing the activating role of auxins and gibberellins against the inhibitory role of ethylene and different candidate genes that could be useful as markers for parthenocarpic selection in the current breeding programs of zucchini.Research worked is supported by the project RTA2014-00078 from the Spanish Institute of Agronomy Research INIA (Instituto Nacional de Investigacion y Tecnologia Agraria y Alimentaria) and also PP.AVA.AVA201601.7, FEDER y FSE (Programa Operativo FSE de Andalucia 2007-2013 "Andalucia se mueve con Europa"). TPV is supported by a FPI scholarship from RTA2011-00044-C02-01/02 project of INIA. The funding agencies were not involved in the design of the study, collection, analysis, and interpretation of data and in writing the manuscript.Pomares-Viciana, T.; Del Rio-Celestino, M.; Roman, B.; Die, J.; Picó Sirvent, MB.; Gómez, P. (2019). First RNA-seq approach to study fruit set and parthenocarpy in zucchini (Cucurbita pepo L.). BMC Plant Biology. 19:1-20. https://doi.org/10.1186/s12870-019-1632-2S12019Varga A, Bruinsma J. Tomato. In: Monselise SP, editor. CRC Handbook of Fruit Set and Development. Boca Raton: CRC Press; 1986. p. 461–80.Nepi M, Cresti L, Guarnieri M, Pacini E. Effect of relative humidity on water content, viability and carbohydrate profile of Petunia hybrid and Cucurbita pepo pollen. Plant Syst Evol. 2010;284:57–64.Gustafson FG. Parthenocarpy: natural and artificial. Bot Rev. 1942;8:599–654.Robinson RW, Reiners S. Parthenocarpy in summer squash. Hortscience. 1999;34:715–7.Pomares-Viciana T, Die J, Del Río-Celestino M, Román B, Gómez P. Auxin signalling regulation during induced and parthenocarpic fruit set in zucchini. Mol Breeding. 2017;37:56.Ozga JA, Reinecke DM. Hormonal interactions in fruit development. J Plant Growth Regul. 2003;22:73–81.Kim IS, Okubo H, Fujieda K. Endogenous levels of IAA in relation to parthenocarpy in cucumber (Cucumis sativus L). Sci Hortic. 1992;52:1–8.Olimpieri I, Siligato F, Caccia R, Mariotti L, Ceccarelli N, Soressi GP, et al. Tomato fruit set driven by pollination or by the parthenocarpic fruit allele are mediated by transcriptionally regulated gibberellin biosynthesis. Planta. 2007;226:877–88.Cui L, Zhang T, Li J, Lou Q, Chen J. Cloning and expression analysis of Cs-TIR1/AFB2: the fruit development-related genes of cucumber (Cucumis sativus L.). Acta Physiol Plant. 2014;36:139–49.De Jong M, Wolters-Arts J, Feron R, Mariani C, Vriezen WH. The Solanum lycopersicum auxin response factor 7 (SlARF7) regulates auxin signalling during tomato fruit set and development. Plant J. 2009;57:160–70.Wang H, Jones B, Li Z, Frasse P, Delalande C, Regad F, Chaabouni S, Latché A, Pech JC, Bouzayen M. The tomato aux/IAA transcription factor IAA9 is involved in fruit development and leaf morphogenesis. Plant Cell. 2005;17(10):2676–92.Goetz M, Vivian-Smith A, Johnson SD, Koltunow AM. AUXIN RESPONSE FACTOR 8 is a negative regulator of fruit initiation in Arabidopsis. Plant Cell. 2006;18(8):1873–86.Mazzucato A, Cellini F, Bouzayen M, Zouine M, Mila I, Minoia S et al. A TILLING allele of the tomato aux/IAA9 gene offers new insights into fruit set mechanisms and perspectives for breeding seedless tomatoes. Mol Breeding. 2015; 35(22):1-15.Blanca J, Cañizares J, Roig C, Ziarsolo P, Nuez F, Picó B. Transcriptome characterization and high throughput SSRs and SNPs discovery in Cucurbita pepo (Cucurbitaceae). BMC Genomics. 2011;12:104.Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63.Da Fonseca RR, Albrechtsen A, Themudo GE, Ramos-Madrigal J, Sibbesen JA, Maretty L, et al. Next-generation biology: sequencing and data analysis approaches for non-model organisms. Mar Genomics. 2016;30:3–13.Conesa A, Madrigal P, Tarazona S, Gomez-Cabrero D, Cervera A, McPherson A, et al. A survey of best practices for RNA-seq data analysis. Genome Biol. 2016;17:13.Li J, Cui ZWJ, Zhang T, Guo Q, Xu J, Li J, et al. Transcriptome comparison of global distinctive features between pollination and parthenocarpic fruit set reveals transcriptional phytohormone cross-talk in cucumber (Cucumis sativus L). Plant Cell Physiol. 2014;55(7):1325–42.Fu L, Niu B, Zhu Z, Wu S, Li W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics. 2012;28(23):3150–2.Montero-Pau J, Blanca J, Bombarely A, Ziarsolo P, Esteras C, Martí-Gómez C, et al. De novo assembly of the zucchini genome reveals a whole genome duplication associated with the origin of the Cucurbita genus. Plant Biotechnol J. 2017. https://doi.org/10.1111/pbi.12860 .Vriezen WH, Feron R, Maretto F, Keijman J, Mariani C. Changes in tomato ovary transcriptome demonstrate complex hormonal regulation of fruit set. New Phytol. 2008;177:60–76.Tang N, Deng W, Hu G, Hu N, Li Z. Transcriptome profiling reveals the regulatory mechanism underlying pollination dependent and parthenocarpic fruit set mainly mediated by auxin and gibberellin. PLoS One. 2015;10(4):e0125355.Li J, Yan S, Yang W, Li Y, Xia M, Chen Z, et al. Transcriptomic analysis reveals the roles of microtubule-related genes and transcription factors in fruit length regulation in cucumber (Cucumis sativus L.). Sci Rep. 2015;26(5):8031.Mironov V, De Veylder L, Van Montagu M, Inze D. Cyclin-dependent kinases and cell division in plants- the nexus. Plant Cell. 1999;11(4):509–22.Perrot-Rechenmann C. Cellular responses to auxin: division versus expansion. Cold Spring Harb Perspect Biol. 2010;2(5):a001446.De Veylder L, Beeckman T, Beemster GT, Krols L, Terras F, Landrieu I, et al. Functional analysis of cyclin-dependent kinase inhibitors of Arabidopsis. Plant Cell. 2001;13:1653–68.Nieuwland J, Menges M, Murray JAH. The plant cyclins. In: Inze D, editor. Cell cycle control and plant development, vol. 2007. Oxford: Wiley-Blackwell Publishing; 2007. p. 33–61.Menges M, Samland AK, Planchais S, Murray JA. The D-type cyclin CYCD3;1 is limiting for the G1-to-S-phasetransition in Arabidopsis. Plant Cell. 2006;18:893–906.Boruc J, Mylle E, Duda M, De Clercq R, Rombauts S, Geelen D, et al. Systematic localization of the Arabidopsis core cell cycle proteins reveals novel cell division complexes. Plant Physiol. 2010;152(2):553–65.Sampedro J, Cosgrove DJ. The expansin superfamily. Genome Biol. 2005;6:242.Esmon CA, Tinsley AG, Ljung K, Sandberg G, Hearne LB, Liscum E. A gradient of auxin and auxin-dependent transcription precedes tropic growth responses. Proc Natl Acad Sci. 2006;103:236–41.De Folter S, Busscher J, Colombo L, Losa A, Angenent GC. Transcript profiling of transcription factors genes during siliques development in Arabidopsis. Plant Mol Bio. 2004;56:351–3662004.Son O, Cho HY, Kim MR, Lee H, Lee MS, Song E, et al. Induction of a homeodomain-leucine zipper gene by auxin is inhibited by cytokinin in Arabidopsis roots. Biochem Biophys Res Commun. 2005;326(1):203–9.Olsson ASB, Engstroem P, Seoderman E. The homeobox genes ATHB12 and ATHB7 encode potential regulators of growth in response to water deficit in Arabidopsis. Plant Mol Biol. 2004;55:663–77.Merrow SB, Hopp RJ. Storage effects on winter squashes. Associations between the sugar and starch content of and the degree of preference for winter squashes. J Agric Food Chem. 1961;9:321–6.Berg JM, Tymoczko JL, Stryer L. Carbohydrates. In: Freeman WH, editor. Biochemistry. 5th ed. New York: W H Freeman; 2002.Prabhakar V, Löttgert T, Gigolashvili T, Bell K, Flügge UI, Häusler RE. Molecular and functional characterization of the plastid-localized phosphoenolpyruvate enolase (ENO1) from Arabidopsis thaliana. FEBS Lett. 2009;583(6):983–91.Rius SP, Casati P, Iglesias AA, Gomez-Casati DF. Characterization of Arabidopsis lines deficient in GAPC-1, a cytosolic NAD-dependent glyceraldehyde-3-phosphate dehydrogenase. Plant Physiol. 2008;148(3):1655–67.Van der Linde K, Gutsche N, Leffers HM, Lindermayr C, Müller B, Holtgrefe S, et al. Regulation of plant cytosolic aldolase functions by redox-modifications. Plant Physiol Biochem. 2011;49(9):946–57.Lim H, Cho MH, Jeon JS, Bhoo SH, Kwon YK, Hahn TR. Altered expression of pyrophosphate: fructose-6-phosphate 1-phosphotransferase affects the growth of transgenic Arabidopsis plants. Mol Cells. 2009;27(6):641–9.Baud S, Wuillème S, Dubreucq B, De Almeida A, Vuagnat C, Lepiniec L, et al. Function of plastidial pyruvate kinases in seeds of Arabidopsis thaliana. Plant J. 2007;52:405–19.De Jong M, Mariani C, Vriezen WH. The role of auxin and gibberellin in tomato fruit set. J Exp Bot. 2009;60(5):1523–32.Martínez C, Manzano S, Megías Z, Garrido D, Picó B, Jamilena M. Involvement of ethylene biosynthesis and signalling in fruit set and early fruit development in zucchini squash (Cucurbita pepo L.). BMC Plant Biol. 2013;13:139.Serrani JC, Fos M, Atarés A, Garcia-martinez JL. Effect of gibberellin and auxin on parthenocarpic fruit growth induction in the cv. micro-tom of tomato. J Plant Growth Regul. 2007;26:211–21.Mapelli S. Changes in cytokinin in the fruits of parthenocarpic and normal tomatoes. Plant Sci Lett. 1981;22:227–33.Ulmasov T, Hagen G, Guilfoyle TJ. Activation and repression of transcription by auxin-response factors. Proc Natl Acad Sci U S A. 1999;96:5844–9.Ulmasov T, Hagen G, Guilfoyle TJ. Dimerization and DNA binding of auxin response factors. Plant J. 1999;19:309–19.Tiwari SB, Hagen G, Guilfoyle TJ. Aux/IAA proteins contain a potent transcriptional repression domain. Plant Cell. 2004;16:533–43.Switzenberg JA, Beaudry RM, Grumet R. Effect of CRC:: etr1-1 transgene expression on ethylene production, sex expression, fruit set and fruit ripening in transgenic melon (Cucumis melo L.). Transgenic Res. 2015;24(3):497-507.Nitsch LM, Oplaat C, Feron R, Ma Q, Wolters-Arts M, Hedden P, et al. Abscisic acid levels in tomato ovaries are regulated by LeNCED1 and SlCYP707A1. Planta. 2009;229(6):1335–46.Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-seq. Nat Methods. 2008;5(7):621–8.Robinson MD, McCarthy DJ, Smyth GK. Edger: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2008;26(1):139–40.Raza K, Mishra A. A novel anticlustering filtering algorithm for the prediction of genes as a drug target. Am J Bio Engi. 2012;2(5):206–11.Van Iterson M, Boer JM, Menezes RX. Filtering, FDR and power. BMCBioinformatics. 2010;11:450.Conesa A, Götz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.Berardini TZ, Reiser L, Li D, Mezheritsky Y, Muller R, Strait E, Huala E. The Arabidopsis information resource: making and mining the “gold standard” annotated reference plant genome. Genesis. 2015. https://doi.org/10.1002/dvg.22877 .Bairoch A, Apweiler R. The SWISS-PROT protein sequence database and its supplement TrEMBL. Nucleic Acids Res. 2000;28(1):45–8.Johnson M, Zaretskaya I, Raytselis Y, Merezhuk Y, McGinnis S, Madden TL. NCBI BLAST: a better web interface. Nucleic Acids Res. 2008;36:W5–9.Wyatt LE, Strickler SR, Mueller LA, Mazourek M. An acorn squash (Cucurbita pepo ssp. ovifera) fruit and seed transcriptome as a resource for the study of fruit traits in Cucurbita. Hortic Res. 2015;2:14070. https://doi.org/10.1038/hortres.2014.70 .Zhang A, Ren GA, Sun YA, Guo H, Zhang SA, Zhang FA, et al. A high-density genetic map for anchoring genome sequences and identifying QTLs associated with dwarf vine in pumpkin (Cucurbita maxima Duch.). BMC Genomics. 2015;16:1101.Finn RD, Attwood TK, Babbit AB, Bork P, Bridge AJ, Chang HY. InterPro in 2017-beyond protein family and domain annotations. Nucleic Acids Res. https://doi.org/10.1093/nar/gkw1107 .Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Sherlock G. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25–9.Kanehisa M, Araki M, Goto S, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:480–4

Genes involved in ethylene and gibberellins metabolism are required for endosperm-limited germiantion of Sisymbrium officinales L. Seeds

The rupture of the seed coat and that of the endosperm were found to be two sequential events in the germination of Sisymbrium officinale L. seeds, and radicle protrusion did not occur exactly in the micropylar area but in the neighboring zone. The germination patterns were similar both in the presence of gibberellins (GA4+7) and in presence of ethrel. The analysis of genes involved in GAs synthesis and breakdown demonstrated that (1) SoGA2ox6 expression peaked just prior to radicle protrusion (20–22 h), while SoGA3ox2 and SoGA20ox2 expression was high at early imbibition (6 h) diminishing sharply thereafter; (2) the accumulation of SoGA20ox2 transcript was strongly inhibited by paclobutrazol (PB) as well as by inhibitors of ET synthesis and signaling (IESS) early after imbibition (6 h), while SoGA3ox2 and SoGA2ox6 expression was slowly depressed as germination progressed; (3) ethrel and GA4+7 positively or negatively affected expression of SoGA3ox2, SoGA20ox2, and SoGA2ox6, depending on the germination period studied. Regarding genes involved in ET synthesis, our results showed that SoACS7 was expressed, just prior to radicle emergence while SoACO2 expression slowly increased as germination progressed. Both genes were strongly inhibited by PB but were almost unaffected by externally added ethrel or GA4+7. These results suggest that GAs are more important than ET during the early stages of imbibition, while ET is more important at the late phases of germination of S. officinale L. seed

Characterization of the cork oak transcriptome dynamics during acorn development

Background: Cork oak (Quercus suber L.) has a natural distribution across western Mediterranean regions and is a keystone forest tree species in these ecosystems. The fruiting phase is especially critical for its regeneration but the molecular mechanisms underlying the biochemical and physiological changes during cork oak acorn development are poorly understood. In this study, the transcriptome of the cork oak acorn, including the seed, was characterized in five stages of development, from early development to acorn maturation, to identify the dominant processes in each stage and reveal transcripts with important functions in gene expression regulation and response to water. Results: A total of 80,357 expressed sequence tags (ESTs) were de novo assembled from RNA-Seq libraries representative of the several acorn developmental stages. Approximately 7.6 % of the total number of transcripts present in Q. suber transcriptome was identified as acorn specific. The analysis of expression profiles during development returned 2,285 differentially expressed (DE) transcripts, which were clustered into six groups. The stage of development corresponding to the mature acorn exhibited an expression profile markedly different from other stages. Approximately 22 % of the DE transcripts putatively code for transcription factors (TF) or transcriptional regulators, and were found almost equally distributed among the several expression profile clusters, highlighting their major roles in controlling the whole developmental process. On the other hand, carbohydrate metabolism, the biological pathway most represented during acorn development, was especially prevalent in mid to late stages as evidenced by enrichment analysis. We further show that genes related to response to water, water deprivation and transport were mostly represented during the early (S2) and the last stage (S8) of acorn development, when tolerance to water desiccation is possibly critical for acorn viability. Conclusions: To our knowledge this work represents the first report of acorn development transcriptomics in oaks. The obtained results provide novel insights into the developmental biology of cork oak acorns, highlighting transcripts putatively involved in the regulation of the gene expression program and in specific processes likely essential for adaptation. It is expected that this knowledge can be transferred to other oak species of great ecological value.Fundação para a Ciência e a Tecnologi

Regulation of the Fruit-Specific PEP Carboxylase SlPPC2 Promoter at Early Stages of Tomato Fruit Development

The SlPPC2 phosphoenolpyruvate carboxylase (PEPC; EC 4.1.1.31) gene from tomato (Solanum lycopersicum) is differentially and specifically expressed in expanding tissues of developing tomato fruit. We recently showed that a 1966 bp DNA fragment located upstream of the ATG codon of the SlPPC2 gene (GenBank AJ313434) confers appropriate fruit-specificity in transgenic tomato. In this study, we further investigated the regulation of the SlPPC2 promoter gene by analysing the SlPPC2 cis-regulating region fused to either the firefly luciferase (LUC) or the β-glucuronidase (GUS) reporter gene, using stable genetic transformation and biolistic transient expression assays in the fruit. Biolistic analyses of 5′ SlPPC2 promoter deletions fused to LUC in fruits at the 8th day after anthesis revealed that positive regulatory regions are mostly located in the distal region of the promoter. In addition, a 5′ UTR leader intron present in the 1966 bp fragment contributes to the proper temporal regulation of LUC activity during fruit development. Interestingly, the SlPPC2 promoter responds to hormones (ethylene) and metabolites (sugars) regulating fruit growth and metabolism. When tested by transient expression assays, the chimeric promoter:LUC fusion constructs allowed gene expression in both fruit and leaf, suggesting that integration into the chromatin is required for fruit-specificity. These results clearly demonstrate that SlPPC2 gene is under tight transcriptional regulation in the developing fruit and that its promoter can be employed to drive transgene expression specifically during the cell expansion stage of tomato fruit. Taken together, the SlPPC2 promoter offers great potential as a candidate for driving transgene expression specifically in developing tomato fruit from various tomato cultivars

Petiole hyponasty: an ethylene-driven, adaptive response to changes in the environment

Hyponastic (upwardly bending) growth by leaves is a response of numerous plant species to adverse environmental conditions. This review summarises current knowledge on hyponasty with a particular focus on the role of ethylene in regulating this phenomenon and its possible adaptive significance