Variation burst during dedifferentiation and increased CHH-type DNA methylation after 30 years of in vitro culture of sweet orange

Abstract

Somaclonal variation arising from tissue culture may provide a valuable resource for the selection of new germplasm, but may not preserve true-to-type characteristics, which is a major concern for germplasm conservation or genome editing. The genomic changes associated with dedifferentiation and somaclonal variation during long-term in vitro culture are largely unknown. Sweet orange was one of the earliest plant species to be cultured in vitro and induced via somatic embryogenesis. We compared four sweet orange callus lines after 30 years of constant tissue culture with newly induced calli by comprehensively determining the single-nucleotide polymorphisms, copy number variations, transposable element insertions, methylomic and transcriptomic changes. We identified a burst of variation during early dedifferentiation, including a retrotransposon outbreak, followed by a variation purge during long-term in vitro culture. Notably, CHH methylation showed a dynamic pattern, initially disappearing during dedifferentiation and then more than recovering after 30 years of in vitro culture. We also analyzed the effects of somaclonal variation on transcriptional reprogramming, and indicated subgenome dominance was evident in the tetraploid callus. We identified a retrotransposon insertion and DNA modification alternations in the potential regeneration-related gene CLAVATA3/EMBRYO SURROUNDING REGION-RELATED 16. This study provides the foundation to harness in vitro variation and offers a deeper understanding of the variation introduced by tissue culture during germplasm conservation, somatic embryogenesis, gene editing, and breeding programs.

References

[1] Jiang F, Feng Z, Liu H et al. Involvement of plant stem cells or stem cell-like cells in dedifferentiation. Front Plant Sci. 2015; 6: 1-6.

[2] Sugiyama M . Historical review of research on plant cell dedifferentiation. J Plant Res. 2015; 128: 349-59.

[3] Larkin PJ, Scowcroft WR . Somaclonal variation-a novel source of variability from cell cultures for plant improvement. Theor Appl Genet. 1981; 60: 197-214.

[4] Krishna H, Alizadeh M, Singh D et al. Somaclonal variations and their applications in horticultural crops improvement. 3 Biotech. 2016; 6: 54-71.

[5] Li J, Manghwar H, Sun L et al. Whole genome sequencing reveals rare off-target mutations and considerable inherent genetic or/and somaclonal variations in CRISPR/Cas9-edited cotton plants. Plant Biotechnol J. 2019; 17: 858-68.

[6] Wei FJ, Kuang L-Y, Oung H-M et al. Somaclonal variation does not preclude the use of rice transformants for genetic screening. Plant J. 2016; 85: 648-59.

[7] Wang X, Xu Y, Zhang S et al. Genomic analyses of primitive, wild and cultivated citrus provide insights into asexual reproduction. Nat Genet. 2017; 49: 765-72.

[8] Spillane C, Curtis MD, Grossniklaus U . Apomixis technology development-virgin births in farmers’ fields? Nat Biotechnol. 2004; 22: 687-91.

[9] Wu GA, Prochnik S, Jenkins J et al. Sequencing of diverse mandarin, pummelo and orange genomes reveals complex history of admixture during citrus domestication. Nat Biotechnol. 2014; 32: 656-62.

[10] Xu Q, Chen L-L, Ruan X et al. The draft genome of sweet orange (Citrus sinensis). Nat Genet. 2013;45:59-66.

[11] Chaturvedi H, Mitra G . Clonal propagation of Citrus from somatic callus cultures. HortScience. 1974; 9: 118-20.

[12] Wang QM, Wang L . An evolutionary view of plant tissue culture: somaclonal variation and selection. Plant Cell Rep. 2012; 31: 1535-47.

[13] Rodriguez Lopez CM, Wetten AC, Wilkinson MJ . Progressive erosion of genetic and epigenetic variation in callus-derived cocoa (Theobroma cacao) plants . New Phytol. 2010; 186: 856-68.

[14] Alkan C, Coe BP, Eichler EE . Genome structural variation discovery and genotyping. Nat Rev Genet. 2011; 12: 363-76.

[15] Eid J, Fehr A, Gray J et al. Real-time DNA sequencing from single polymerase molecules. Science. 2009; 323: 133-8.

[16] Sedlazeck FJ, Rescheneder P, Smolka M et al. Accurate detection of complex structural variations using single-molecule sequencing. Nat Methods. 2018; 15: 461-8.

[17] Gaur U, Li K, Mei S et al. Research progress in allele-specific expression and its regulatory mechanisms. J Appl Genet. 2013; 54: 271-83.

[18] Skelly DA, Johansson M, Madeoy J et al. A powerful and flexible statistical framework for testing hypotheses of allele-specific gene expression from RNA-seq data. Genome Res. 2011; 21: 1728-37.

[19] Edger PP, Smith R, MR MK et al. Subgenome dominance in an interspecific hybrid, synthetic allopolyploid, and a 140-year-old naturally established neo-allopolyploid monkeyflower. Plant Cell. 2017; 29: 2150-67.

[20] Zhao B, Cao J-F, Hu G-J et al. Core cis-element variation confers subgenome-biased expression of a transcription factor that functions in cotton fiber elongation. New Phytol. 2018; 218: 1061-75.

[21] Shao L, Xing F, Xu C et al. Patterns of genome-wide allele-specific expression in hybrid rice and the implications on the genetic basis of heterosis. Proc Natl Acad Sci U S A. 2019; 116: 5653-8.

[22] Shen R, Wang L, Liu X et al. Genomic structural variation-mediated allelic suppression causes hybrid male sterility in rice. Nat Commun. 2017; 8: 1310-9.

[23] Paschold A, Larson NB, Marcon C et al. Nonsyntenic genes drive highly dynamic complementation of gene expression in maize hybrids. Plant Cell. 2014; 26: 3939-48.

[24] Waters AJ, Makarevitch I, Noshay J et al. Natural variation for gene expression responses to abiotic stress in maize. Plant J. 2017; 89: 706-17.

[25] Zhang X, Borevitz JO . Global analysis of allele-specific expression in Arabidopsis thaliana . Genetics. 2009; 182: 943-54.

[26] Todesco M, Balasubramanian S, Hu TT et al. Natural allelic variation underlying a major fitness trade-off in Arabidopsis thaliana . Nature. 2010; 465: 632-6.

[27] Peschke VM, Phillips RL, Gengenbach BG . Discovery of transposable element activity among progeny of tissue culture-derived maize plants. Science. 1987; 238: 804-7.

[28] Hirochika H, Sugimoto K, Otsuki Y et al. Retrotransposons of rice involved in mutations induced by tissue culture. Proc Natl Acad Sci U S A. 1996; 93: 7783-8.

[29] Lin C, Lin X, Hu L et al. Dramatic genotypic difference in, and effect of genetic crossing on, tissue culture-induced mobility of retrotransposon Tos17 in rice. Plant Cell Rep. 2012; 31: 2057-63.

[30] Kazuhiro K, Kazuki T, Masamitsu W et al. The plant MITE mPing is mobilized in anther culture. Nature. 2003; 421: 167-70.

[31] Komatsu M, Shimamoto K, Kyozuka J . Two-step regulation and continuous retrotransposition of the rice LINE-type retrotransposon karma . Plant Cell. 2003; 15: 1934-44.

[32] Kaeppler SM, Kaeppler HF, Rhee Y . Epigenetic aspects of somaclonal variation in plants. Plant Mol Biol. 2000; 43: 179-88.

[33] Miguel C, Marum L . An epigenetic view of plant cells cultured in vitro: somaclonal variation and beyond. J Exp Bot. 2011; 62: 3713-25.

[34] Stroud H, Ding B, Simon SA et al. Plants regenerated from tissue culture contain stable epigenome changes in rice. Elife. 2013; 2: e00354.

[35] Ong-Abdullah M, Ordway JM, Jian N et al. Loss of Karma transposon methylation underlies the mantled somaclonal variant of oil palm. Nature. 2015; 525: 533-7.

[36] Gent JI, Ellis NA, Guo L et al. CHH islands: de novo DNA methylation in near-gene chromatin regulation in maize. Genome Res. 2013; 23: 628-37.

[37] Cassells AC, Curry RF . Oxidative stress and physiological, epigenetic and genetic variability in plant tissue culture: implications for micropropagators and genetic engineers. Plant Cell Tissue Organ Cult. 2001; 64: 145-57.

[38] Etienne H, Bertrand B . Somaclonal variation in Coffea arabica: effects of genotype and embryogenic cell suspension age on frequency and phenotype of variants. Tree Physiol. 2003; 23: 419-26.

[39] Belles-Boix E, Hamant O, Witiak SM et al. KNAT6: an Arabidopsis homeobox gene involved in meristem activity and organ separation . Plant Cell. 2006; 18: 1900-7.

[40] Hung FY, Chen F-F, Li C et al. The Arabidopsis LDL1/2-HDA6 histone modification complex is functionally associated with CCA1/LHY in regulation of circadian clock genes. Nucleic Acids Res. 2018; 46: 10669-81.

[41] Luo M, Yu C-W, Chen F-F et al. Histone deacetylase HDA6 is functionally associated with AS1 in repression of KNOX genes in Arabidopsis . PLoS Genet. 2012; 8: e1003114.

[42] Hundertmark M, Hincha DK . LEA (late embryogenesis abundant) proteins and their encoding genes in Arabidopsis thaliana . BMC Genomics. 2008; 9: 118-39.

[43] Rodriguez-Leal D, Xu C, Kwon C-T et al. Evolution of buffering in a genetic circuit controlling plant stem cell proliferation. Nat Genet. 2019; 51: 786-92.

[44] Poetsch AR, Boulton SJ, Luscombe NM . Genomic landscape of oxidative DNA damage and repair reveals regioselective protection from mutagenesis. Genome Biol. 2018; 19: 215-37.

[45] Chuong EB, Elde NC, Feschotte C . Regulatory activities of transposable elements: from conflicts to benefits. Nat Rev Genet. 2017; 18: 71-86.

[46] Naito K, Zhang F, Tsukiyama T et al. Unexpected consequences of a sudden and massive transposon amplification on rice gene expression. Nature. 2009; 461: 1130-4.

[47] Ke L, Zhou Z, Xu X-W et al. Evolutionary dynamics of lincRNA transcription in nine citrus species. Plant J. 2019; 98: 912-27.

[48] Li H, Durbin R . Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009; 25: 1754-60.

[49] Li H, Handsaker B, Wysoker A et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009; 25: 2078-9.

[50] Alexander DH, Novembre J, Lange K . Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009; 19: 1655-64.

[51] Boeva V, Popova T, Bleakley K et al. Control-FREEC: a tool for assessing copy number and allelic content using next-generation sequencing data. Bioinformatics. 2012; 28: 423-5.

[52] Harris RS . Improved pairwise alignment of genomic DNA. Ph.D. Thesis. The Pennsylvania State University; 2007.

[53] Blanchette M, Kent WJ, Riemer C et al. Aligning multiple genomic sequences with the threaded blockset aligner. Genome Res. 2004; 14: 708-15.

[54] Hubisz MJ, Pollard KS, Siepel A . PHAST and RPHAST: phylogenetic analysis with space/time models. Brief Bioinformatics. 2010; 12: 41-51.

[55] Cingolani P, Platts A, Wang Le L et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin). 2012; 6: 80-92.

[56] Choi Y, Chan AP . PROVEAN web server: a tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics. 2015; 31: 2745-7.

[57] Hackl T, Hedrich R, Schultz J et al. proovread: large-scale high-accuracy PacBio correction through iterative short read consensus. Bioinformatics. 2014; 30: 3004-11.

[58] Robinson JT, Thorvaldsdóttir H, Winckler W et al. Integrative genomics viewer. Nat Biotechnol. 2011; 29: 24-6.

[59] Altschul SF, Gish W, Miller W et al. Basic local alignment search tool. J Mol Biol. 1990; 215: 403-10.

[60] Langmead B, Salzberg SL . Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012; 9: 357-9.

[61] Kim D, Pertea G, Trapnell C et al. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013; 14: 295-311.

[62] Trapnell C, Williams BA, Pertea G et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010; 28: 511-5.

[63] Du Z, Zhou X, Ling Y et al. agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010; 38: 64-70.

[64] Langfelder P, Horvath S . WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008; 9: 559-71.

[65] Su G, Morris JH, Demchak B et al. Biological network exploration with Cytoscape 3. Curr Protoc Bioinformatics. 2014; 47: 11-24.

[66] Romanel A, Lago S, Prandi D et al. ASEQ: fast allele-specific studies from next-generation sequencing data. BMC Medical Genom. 2015; 8: 9-20.

[67] Bolger AM, Lohse M, Usadel B . Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014; 30: 2114-20.

[68] Krueger F, Andrews SR . Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 2011; 27: 1571-2.

[69] Krzywinski M, Schein J, Birol I et al. Circos: an information aesthetic for comparative genomics. Genome Res. 2009; 19: 1639-45.

[70] Akalin A, Kormaksson M, Li S et al. methylKit: a comprehensive R package for the analysis of genome-wide DNA methylation profiles. Genome Biol. 2012; 13: 87-95.