The complexity of the Fragaria x ananassa (octoploid) transcriptome by single-molecule long-read sequencing

Abstract

Strawberry (Fragaria x ananassa) is an allopolyploid species with diverse and complex transcripts. The regulatory mechanisms of fruit development and maturation have been extensively studied; however, little is known about the signaling mechanisms that direct this process in octoploid strawberry (Fragaria x ananassa). Here, we used long-read sequencing (LRS) technology and RNA-seq analysis to investigate the diversity and complexity of the polyploid transcriptome and differentially expressed transcripts along four successive fruit developmental stages of cultivated strawberry. We obtained a reference transcriptome with 119,897 unique full-length isoforms, including 2017 new isoforms and 2510 long noncoding RNAs. Based on the genome of the plausible progenitor (Fragaria vesca), 20,229 alternative splicing (AS) events were identified. Using this transcriptome, we found 17,485 differentially expressed transcripts during strawberry fruit development, including 527 transcription factors (TFs) belonging to 41 families. The expression profiles of all members of the auxin, ABA pathway, and anthocyanin biosynthesis gene families were also examined, and many of them were highly expressed at the ripe fruit stage, strongly indicating that the role of those genes is in the regulation of fruit ripening. We produce a high-quality reference transcriptome for octoploid strawberry, including much of the full-length transcript diversity, to help understand the regulatory mechanisms of fruit development and maturation of polyploid species, particularly via elucidation of the biochemical pathways involved in auxin, ABA, and anthocyanin biosynthesis.

References

[1] Nellist, C. F. et al. Quantitative trait loci controlling phytophthora cactorum resistance in the cultivated octoploid strawberry (Fragaria x ananassa). bioRxiv. https://doi.org/10.1101/249573 (2018).

[2] Giovannoni, J. J. Genetic regulation of fruit development and ripening. Plant Cell 16(Suppl), 170-180 (2004).

[3] Barry, C. S. & Giovannoni, J. J. Ethylene and fruit ripening. Plant Growth Regul. 26, 143-159 (2007).

[4] Ren, J. et al. Cloning and expression analysis of cDNAs for ABA 8’-hydroxylase during sweet cherry fruit maturation and under stress conditions. Plant Physiol. 167, 1486-1493 (2010).

[5] Aharoni, A. et al. Novel insight into vascular, stress, and auxin-dependent and -independent gene expression programs in strawberry, a non-climacteric fruit. Plant Physiol. 129, 1019-1031 (2002).

[6] Nitsch, J. P. Growth and morphogenesis of the strawberry as related to auxin. Am. J. Bot. 37, 211-215 (1950).

[7] Jia, H. F. et al. Abscisic acid plays an important role in the regulation of strawberry fruit ripening. Plant Physiol. 157, 188-199 (2011).

[8] Li, C. et al. Abscisic acid perception and signaling transduction in strawberry: a model for non-climacteric fruit ripening. Plant Signal. Behav. 6, 1950-1953 (2011).

[9] Hancock, J. F. (eds). Strawberries. (CAB International, Wallingford, Cambridge, New York. 1999).

[10] Hirakawa, H. et al. Dissection of the octoploid strawberry genome by deep sequencing of the genomes of Fragaria species. DNA Res. 21, 169-181 (2014).

[11] Yoo, M. J., Liu, X., Pires, J. C., Soltis, P. S. & Soltis, D. E. Nonadditive gene expression in polyploids. Annu. Rev. Genet. 48, 485-517 (2014).

[12] Wang, B. et al. Unveiling the complexity of the maize transcriptome by single-molecule long-read sequencing. Nat. Commun. https://doi.org/10.1038/ncomms11708 (2016).

[13] Abdel-Ghany, S. E. et al. A survey of the sorghum transcriptome using single-molecule long reads. Nat. Commun. https://doi.org/10.1038/ncomms11706 (2016).

[14] Lin, Y. et al. Comparative transcriptome profiling analysis of red- and white-fleshed strawberry (Fragaria×ananassa) provides new insight into the regulation of anthocyanins pathway. Plant Cell Physiol. https://doi.org/10.1093/pcp/pcy098 (2018).

[15] Hu, P. P. et al. Transcriptome profiling by RNA-Seq reveals differentially expressed genes related to fruit development and ripening characteristics in strawberries (Fragaria×ananassa). PeerJ. https://doi.org/10.7717/peerj.4976 (2018).

[16] Michael, T. P. & Vanburen, R. Progress, challenges and the future of crop genomes. Curr. Opin. Plant Biol. 24, 71-81 (2015).

[17] Li, Y. P., Dai, C., Hu, C. G., Liu, Z. C. & Kang, C. Y. Global identification of alternative splicing via comparative analysis of SMRT- and Illumina-based RNA-seq in strawberry. Plant J. 90, 164-176 (2017).

[18] Marquez, Y. et al. Transcriptome survey reveals increased complexity of the alternative splicing landscape in Arabidopsis. Genome Res. 22, 1184-1195 (2012).

[19] Campbell, M. A., Haas, B. J., Hamilton, J. P., Mount, S. M. & Buell, C. R. Comprehensive analysis of alternative splicing in rice and comparative analyses with Arabidopsis. BMC Genomics https://doi.org/10.1186/1471-2164-7-327 (2006).

[20] Wang, B. & Brendel, V. Genome wide comparative analysis of alternative splicing in plants. Proc. Natl Acad. Sci. USA 103, 7175-7180 (2006).

[21] Wang, T. T. et al. Comprehensive profiling of rhizome-associated alternative splicing and alternative polyadenylation in moso bamboo (Phyllostachys edulis). Plant J. 91, 684-699 (2017).

[22] Wang, M. J. et al. A global survey of alternative splicing in allopolyploid cotton: landscape, complexity and regulation. New Phytol. 217, 163-178 (2017).

[23] Cheng, B., Furtado, A. & Henry, R. J. Long-read sequencing of the coffee bean transcriptome reveals the diversity of full-length transcripts. Gigascience 6, 1-13 (2017).

[24] Shulaev, V. et al. The genome of woodland strawberry (Fragaria vesca). Nat. Genet. 43, 109-116 (2010).

[25] Edger, P. P. et al. Single-molecule sequencing and optical mapping yields an improved genome of woodland strawberry (Fragaria vesca) with chromosome-scale contiguity. Gigascience 7, 1-7 (2018).

[26] Li, Y. P. et al. Genome re-annotation of the wild strawberry Fragaria vesca using extensive Illumina- and SMRT-based RNA-seq datasets. DNA Res. 25, 61-70 (2018).

[27] Ehlers, K. et al. The MADS box genes ABS, SHP1, and SHP2 are essential for the coordination of cell divisions in ovule and seed coat development and for endosperm formation in Arabidopsis thaliana. PLoS One 11, 165075-165104 (2016).

[28] Maier, A. et al. Light and the E3 ubiquitin ligase COP1/SPA control the protein stability of the MYB transcriptions factors PAP1 and PAP2 involved in anthocyanin accumulation in Arabidopsis. Plant J. 74, 638-651 (2013).

[29] Kunieda, T. & Mitsuda, N. NAC family proteins NARS1/NAC2 and NARS2/NAM in the outer integument regulate embryogenesis in Arabidopsis. Plant Cell 20, 2631-2642 (2008).

[30] Dreher, T. W. & Poovaiah, B. W. Changes in auxin content during development in strawberry fruits. J. Plant Growth Regul. 1, 267-276 (1982).

[31] Archbold, D. D. & Dennis, F. G. J. Quantification of free ABA and free and conjugated IAA in strawberry achene and receptacle tissue during fruit development. J. Am. Soc. Hortic. Sci. 109, 330-335 (1984).

[32] Staswick, P. E. et al. Characterization of an Arabidopsis enzyme family that conjugates amino acids to Indole-3-Acetic Acid. Plant Cell 17, 616-627 (2005).

[33] Gray, W. M. et al. Auxin regulates SCFTIR1-dependent degradation of AUX/IAA proteins. Nature 414, 271-276 (2001).

[34] Tiwari, S. B., Hagen, G. & Guilfoyle, T. J. The roles of auxin response factor domains in auxin-responsive transcription. Plant Cell 15, 533-543 (2003).

[35] Zufall, R. A. & Rausher, M. D. Genetic changes associated with floral adaptation restrict future evolutionary potential. Nature 428, 847-850 (2004).

[36] Fragaria ananassa GDR RefTrans V1: https://www.rosaceae.org/analysis/230.

[37] Seymour, G. B. et al. Fruit development and ripening. Annu. Rev. Plant Biol. 64, 219-224 (2013).

[38] Medina-Puche, L. et al. MYB10 plays a major role in the regulation of flavonoid/phenylpropanoid metabolism during ripening of Fragaria x ananassa fruits. J. Exp. Bot. 65, 401-417 (2015).

[39] Hocher, V. et al. Changes in indole-3-acetic acid levels during tomato (Lycopersicon esculentum Mill.) seed development. Plant Cell Rep. 11, 253-256 (1992).

[40] Pattison, R. J. & Catalá, C. Evaluating auxin distribution in tomato (Solanum lycopersicum) through an analysis of the PIN and AUX/LAX gene families. Plant J. 70, 585-598 (2012).

[41] Huang, H. et al. Flower colour modification of chrysanthemum by suppression of f3’h and overexpression of the exogenous senecio cruentus f3'5’h gene. PLoS One 8, 74395 (2013).

[42] Han, Y. et al. Ectopic expression of apple F3′ H genes contributes to anthocyanin accumulation in the Arabidopsis tt7 mutant grown under nitrogen stress. Plant Physiol. 153, 806-820 (2010).

[43] Hackl, T., Hedrich, R., Schultz, J. & Förster, F. proovread: large-scale high-accuracy PacBio correction through iterative short read consensus. Bioinformatics 30, 3004-3011 (2014).

[44] Li, W. Z. & Godzik, A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22, 1658-1659 (2006).

[45] Kong, L. et al. CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 35, 345-349 (2007).

[46] Sun, L. et al. Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res. 41, 166 (2013).

[47] Wang, L. G. et al. CPAT: Coding-Potential Assessment Tool using an alignment-free logistic regression model. Nucleic Acids Res. 41, 74 (2013).

[48] Finn, R. D. et al. The Pfam protein families database: towards a more sustainable future. Nucleic Acids Res. 44, 279-285 (2016).

[49] Wu, T. D. & Watanabe, C. K. GMAP: a genomic mapping and alignment program form RNA and EST sequences. Bioinformatics 21, 1859-1875 (2005).

[50] Foissac, S. & Sammeth, M. ASTALAVISTA: dynamic and flexible analysis of alternative splicing events in custom gene datasets. Nucleic Acids Res. 35, 297 (2007).

[51] Altschul, S. F. et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389-3402 (1997).

[52] Waterhouse, R. M. et al. BUSCO applications from quality assessments to gene prediction and phylogenomics. Mol. Biol. Evol. 35, 543-548 (2017).

[53] Hirakawa, H. et al. Dissection of the octoploid strawberry genome by deep sequencing of the genomes of Fragaria species. DNA Res. 21, 169-181 (2013).

[54] Langmead, B. & Salzberg, S. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357-359 (2012).

[55] Anders, S. & Huber, W. Differential expression analysis for sequence count data. Genome Biol. 11, 106 (2010).

[56] Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate: a practical and powerful approach to multiple hypothesis testing. J. R. Stat. Soc. B 57, 289-300 (1995).

[57] Walvoort, D. J. J., Brus, D. J. & deGruijter, J. J. An r package for spatial coverage sampling and random sampling from compact geographical strata by k-means. Comput. Geosci. 36, 1261-1267 (2010).

[58] Jin, J. P. et al. PlantTFDB 4.0: toward a central hub for transcription factors and regulatory interactions in plants. Nucleic Acids Res. 45, 1040-1045 (2017).