The reference genome and full-length transcriptome of pakchoi provide insights into cuticle formation and heat adaption

Abstract

Brassica rapa includes various vegetables with high economic value. Among them, green petiole type pakchoi (B. rapa ssp. chinensis) is one of the major vegetables grown in southern China. Compared with other B. rapa varieties, green petiole type pakchoi shows a higher level of heat resistance, which is partially derived from the rich epicuticular wax. Here we sequence a high-quality genome of green petiole type pakchoi, which has been widely used as the parent in breeding. Our results reveal that long terminal repeat retrotransposon insertion plays critical roles in promoting the genome expansion and transcriptional diversity of pakchoi genes through preferential insertions, particularly in cuticle biosynthetic genes. After whole-genome triplication, over-retained pakchoi genes escape stringent selection pressure, and among them a set of cuticle-related genes are retained. Using bulked-segregant analysis of a heat-resistant pakchoi cultivar, we identify a frame-shift deletion across the third exon and the subsequent intron of BrcCER1 in candidate regions . Using Nanopore long-read sequencing, we analyze the full-length transcriptome of two pakchoi cultivars with opposite sensitivity to high temperature. We find that the heat-resistant pakchoi cultivar can mitigate heat-caused leaf damage by activating an unfolded protein response, as well as by inhibiting chloroplast development and energy metabolism, which are presumably mediated by both transcriptional regulation and splicing factors. Our study provides valuable resources for Brassica functional genomics and breeding research, and deepens our understanding of plant stress resistance.

References

[1] Warwick SI, Francis A, Al-Shehbaz IA . Brassicaceae: species checklist and database on CD-rom. Plant Syst and Evol. 2006; 259: 249-25.

[2] Cheng F, Wu J, Wang X . Genome triplication drove the diversification of Brassica plants. Hortic Res. 2014; 1: 14024.

[3] Song X, Wei Y, Xiao D et al. Brassica carinata genome characterization clarifies U’s triangle model of evolution and polyploidy in Brassica . Plant Physiol. 2021; 186: 388-406.

[4] Nagaharu U . Genome analysis in brassica with special reference to the experimental formation of B. napus and peculiar mode of fertilization. Jpn J Bot. 1935; 7: 389-452.

[5] Bowers JE, Chapman BA, Rong J et al. Unravelling angiosperm genome evolution by phylogenetic analysis of chromosomal duplication events. Nature. 2003; 422: 433-8.

[6] Jiao Y, Wickett NJ, Ayyampalayam S et al. Ancestral polyploidy in seed plants and angiosperms. Nature. 2011; 473: 97-100.

[7] Yang YW, Lai KN, Tai PY et al. Rates of nucleotide substitution in angiosperm mitochondrial DNA sequences and dates of divergence between Brassica and other angiosperm lineages. J Mol Evol. 1999; 48: 597-604.

[8] Town CD, Cheung F, Maiti R et al. Comparative genomics of Brassica oleracea and Arabidopsis thaliana reveal gene loss, fragmentation, and dispersal after polyploidy . Plant Cell. 2006; 18: 1348-59.

[9] Wang X, Wang H, Wang J et al. The genome of the mesopolyploid crop species Brassica rapa . Nat Genet. 2011; 43: 1035-9.

[10] Cheng F, Mandáková T, Wu J et al. Deciphering the diploid ancestral genome of the mesohexaploid Brassica rapa . Plant Cell. 2013; 25: 1541-54.

[11] Zhao J, Wang X, Deng B et al. Genetic relationships within Brassica rapa as inferred from AFLP fingerprints. Theor Appl Genet. 2005; 110: 1301-14.

[12] Xiao D, Wang H, Basnet RK et al. Genetic dissection of leaf development in Brassica rapa using a genetical genomics approach. Plant Physiol. 2014; 164: 1309-25.

[13] Zhang X, Liu Z, Wang P et al. Fine mapping of BrWax1, a gene controlling cuticular wax biosynthesis in Chinese cabbage (Brassica rapa L. ssp. pekinensis) . Mol Breed. 2013; 32: 867-74.

[14] Sturaro M, Motto M, Hemadaranjan A . Plant cuticular waxes: biosynthesis and functions. Adv Plant Physiology. 2006; 9: 229-51.

[15] Pu Y, Gao J, Guo Y et al. A novel dominant glossy mutation causes suppression of wax biosynthesis pathway and deficiency of cuticular wax in Brassica napus . Plant Biol. 2013; 13: 215-5.

[16] Li P, Su T, Zhao X et al. Assembly of the non-heading pak choi genome and comparison with the genomes of heading Chinese cabbage and the oilseed yellow sarson. Plant Biotechnol J. 2021; 19: 966-76.

[17] Li Y, Liu GF, Ma LM et al. A chromosome-level reference genome of non-heading Chinese cabbage [ Brassica campestris (syn. Brassica rapa) ssp. chinensis] . Hortic Res 2020; 7: 212

[18] Zhang L et al. Improved Brassica rapa reference genome by single-molecule sequencing and chromosome conformation capture technologies. Hortic Res. 2018; 5: 50.

[19] Simão FA, Waterhouse RM, Ioannidis P et al. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015; 31: 3210-2.

[20] Campbell MS, Holt C, Moore B et al. Genome annotation and curation using MAKER and MAKER-P. Curr Protoc Bioinformatics. 2014; 48: 4.11.1-39.

[21] Liu S, Liu Y, Yang X et al. The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes. Nat Commun. 2014; 5: 3930.

[22] Zheng N, Schulman BA, Song L et al. Structure of the Cul1-Rbx1- Skp1-F box^Skp2 SCF ubiquitin ligase complex . Nature. 2002; 416: 703-9.

[23] Xu G, Ma H, Nei M et al. Evolution of F-box genes in plants: different modes of sequence divergence and their relationships with functional diversification. Proc Natl Acad Sci USA 2009; 106: 835-40.

[24] Calderón-Villalobos LIA, Nill C, Marrocco K et al. The evolutionarily conserved Arabidopsis thaliana F-box protein AtFBP7 is required for efficient translation during temperature stress. Gene. 2007; 392: 106-16.

[25] Calderón Villalobos LIA, Lee S, de Oliveira C et al. A combinatorial TIR1/AFB-aux/IAA co-receptor system for differential sensing of auxin. Nat Chem Biol. 2012; 8: 477-85.

[26] Cai X, Chang L, Zhang T et al. Impacts of allopolyploidization and structural variation on intraspecific diversification in Brassica rapa . Genome Biol. 2021; 22: 166.

[27] Kuroda H, Yanagawa Y, Takahashi N et al. A comprehensive analysis of interaction and localization of Arabidopsis SKP1-LIKE (ASK) and F-box (FBX) proteins. PLoS One. 2012; 7: e50009.

[28] Zhang Y, Xu W, Li Z et al. F-box protein DOR functions as a novel inhibitory factor for abscisic acid-induced stomatal closure under drought stress in Arabidopsis . Plant Physiol. 2008; 148: 2121-33.

[29] Pélissier T, Mathieu O . Glue for jumping elements: epigenetic means for controlling transposable elements in plants. In: Grandbastien M-A, Casacuberta JM, eds. Plant Transposable Elements: Impact on Genome Structure and Function.Springer: Berlin, 2012, 125-45.

[30] Lisch D . How important are transposons for plant evolution? Nature Rev Genet. 2013; 14: 49-61.

[31] Feschotte C . Transposable elements and the evolution of regulatory networks. Nat Rev Genet. 2008; 9: 397-405.

[32] Bernard A, Joubès JM . Arabidopsis cuticular waxes: advances in synthesis, export and regulation . Prog Lipid Res. 2013; 52: 110-29.

[33] Bourdenx B, Bernard A, Domergue F et al. Overexpression of Arabidopsis ECERIFERUM1 promotes wax very-long-chain alkane biosynthesis and influences plant response to biotic and abiotic stresses. Plant Physiol. 2011; 156: 29-45.

[34] Owen Rowland HZ, Hepworth SR, Lam P et al. CER4 encodes an alcohol-forming fatty acyl-coenzyme a reductase involved in cuticular wax production in Arabidopsis . Plant Physiol. 2006; 142: 866-77.

[35] Hickey DA, Benkel B . Introns as relict retrotransposons: implications for the evolutionary origin of eukaryotic mRNA splicing mechanisms. J Theor Biol. 1986; 121: 283-91.

[36] Jaligot E, Hooi WY, Debladis E et al. DNA methylation and expression of the EgDEF1 gene and neighboring retrotransposons in mantled somaclonal variants of oil palm. PLoS One. 2014; 9: e91896.

[37] Galindo-González L, Mhiri C, Deyholos MK et al. LTR-retrotransposons in plants: engines of evolution. Gene. 2017; 626: 14-25.

[38] Bhanot V, Fadanavis SV, Panwar J . Revisiting the architecture, biosynthesis and functional aspects of the plant cuticle: there is more scope. Environ Exp Bot. 2021; 183: 104364.

[39] Bueno A, Alfarhan A, Arand K et al. Effects of temperature on the cuticular transpiration barrier of two desert plants with water-spender and water-saver strategies. J Exp Bot. 2019; 70: 1613-25.

[40] Fedoroff NV, Battisti DS, Beachy RN et al. Radically rethinking agriculture for the 21st century. Science. 2010; 327: 833-4.

[41] Cheng F, Wu J, Fang L et al. Biased gene fractionation and dominant gene expression among the subgenomes of Brassica rapa . PLoS One. 2012; 7: e36442.

[42] Krishnamurthy P, Kim JA, Jeong MJ et al. Gene loss/retention and evolutionary pattern of ascorbic acid biosynthesis and recycling genes in Brassica rapa following whole genome triplication. Genes Genomics. 2016; 38: 1129-43.

[43] Ramírez-González RH, Borrill P, Lang D et al. The transcriptional landscape of polyploid wheat. Science. 2018; 361: eaar6089.

[44] Schnable JC, Springer NM, Freeling M . Differentiation of the maize subgenomes by genome dominance and both ancient and ongoing gene loss. Proc Natl Acad Sci USA 2011; 108: 4069-74.

[45] Bergh EVD, Hofberger JA, Schranz ME . Flower power and the mustard bomb: comparative analysis of gene and genome duplications in glucosinolate biosynthetic pathway evolution in Cleomaceae and Brassicaceae. Am J Bot. 2016; 103: 1212-22.

[46] Weike D, Song X, Liu T et al. Patterns of evolutionary conservation of ascorbic acid-related genes following whole-genome triplication in Brassica rapa . Genome Biol Evol. 2014; 7: 299-313.

[47] Pratt WB, Toft DO . Regulation of signaling protein function and trafficking by the hsp90/hsp70-based chaperone machinery. Exp Biol Med (Maywood). 2003; 228: 111-33.

[48] Arsova B, Hoja U, Wimmelbacher M et al. Plastidial thioredoxin z interacts with two fructokinase-like proteins in a thiol-dependent manner: evidence for an essential role in chloroplast development in Arabidopsis and Nicotiana benthamiana . Plant Cell. 2010; 22: 1498-515.

[49] Wimmelbacher M, Börnke F . Redox activity of thioredoxin z and fructokinase-like protein 1 is dispensable for autotrophic growth of Arabidopsis thaliana . J Exp Bot. 2014; 65: 2405-13.

[50] Lv Y, Shao G, Qiu J et al. White leaf and panicle 2, encoding a PEP-associated protein, is required for chloroplast biogenesis under heat stress in rice. J Exp Bot. 2017; 68: 5147-60.

[51] Sirpiö S, Allahverdiyeva Y, Holmström M et al. Novel nuclear-encoded subunits of the chloroplast NAD(P)H dehydrogenase complex. J Biol Chem. 2009; 284: 905-12.

[52] Sun Y, Geng Q, Du Y et al. Induction of cyclic electron flow around photosystem I during heat stress in grape leaves. Plant Sci. 2017; 256: 65-71.

[53] Wang S, Bai G, Wang S et al. Chloroplast RNA-binding protein RBD1 promotes chilling tolerance through 23S rRNA processing in Arabidopsis . PLoS Genet. 2016; 12: e1006027.

[54] Asakura Y, Barkan A . Arabidopsis orthologs of maize chloroplast splicing factors promote splicing of orthologous and species-specific group II introns . Plant Physiol. 2006; 142: 1656-63.

[55] Tanabe N, Yoshimura K, Kimura A et al. Differential expression of alternatively spliced mRNAs of Arabidopsis SR protein homologs, atSR30 and atSR45a, in response to environmental stress. Plant Cell Physiol. 2007; 48: 1036-49.

[56] Soltis DE, Smith SA, Cellinese N et al. Angiosperm phylogeny: 17 genes, 640 taxa. Am J Bot. 2011; 98: 704-30.

[57] Song X-M, Wang JP, Sun PC et al. Preferential gene retention increases the robustness of cold regulation in Brassicaceae and other plants after polyploidization. Hortic Res. 2020; 7: 20.

[58] Peterson-Burch BD, Nettleton D, Voytas DF . Genomic neighborhoods for Arabidopsis retrotransposons: a role for targeted integration in the distribution of the Metaviridae. Genome Biol. 2004; 5: R78.

[59] Jin W, Melo JR, Nagaki K et al. Maize centromeres: organization and functional adaptation in the genetic background of oat. Plant Cell. 2004; 16: 571-81.

[60] Miller WJ, Capy P . Applying mobile genetic elements for genome analysis and evolution. Mol Biotechnol. 2006; 33: 161-74.

[61] Todorovska E. Retrotransposons and their role in plant-genome evolution. Biotechnol Biotechnol Equip. 2014; 21: 294-305.

[62] Beulé T, Agbessi MDT, Dussert S et al. Genome-wide analysis of LTR-retrotransposons in oil palm. BMC Genomics. 2015; 16: 795.

[63] Koren S, Walenz BP, Berlin K et al. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 2017; 27: 722-36.

[64] Vaser R, Sović I, Nagarajan N et al. Fast and accurate de novo genome assembly from long uncorrected reads. Genome Res. 2017; 27: 737-46.

[65] Walker BJ, Abeel T, Shea T et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One. 2014; 9: e112963.

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

[67] Burton JN, Adey A, Patwardhan RP et al. Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Nat Biotechnol. 2013; 31: 1119-25.

[68] Marchler-Bauer A, Lu S, Anderson JB et al. CDD: a conserved domain database for the functional annotation of proteins. Nucleic Acids Res. 2011; 39: D225-9.

[69] Koonin EV, Fedorova ND, Jackson JD et al. A comprehensive evolutionary classification of proteins encoded in complete eukaryotic genomes. Genome Biol. 2004; 5: R7.

[70] Ogata H, Goto S, Sato K et al. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 1999; 27: 29-34.

[71] Boeckmann B, Bairoch A, Apweiler R et al. The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 2003; 31: 365-70.

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

[73] Dimmer EC, Huntley RP, Alam-Faruque Y et al. The UniProt-GO annotation database in 2011. Nucleic Acids Res. 2012; 40: D565-70.

[74] Harris MA, Clark J, Ireland A et al. The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res. 2004; 32: D258-61.

[75] Conesa A, Götz S, García-Gómez JM et al. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005; 21: 3674-6.

[76] Zheng Y, Jiao C, Sun H et al. iTAK: a program for genome-wide prediction and classification of plant transcription factors, transcriptional regulators, and protein kinases. Mol Plant. 2016; 9: 1667-70.

[77] Letunic I, Bork P . 20 years of the SMART protein domain annotation resource. Nucleic Acids Res. 2018; 46: D493-d496.

[78] Prakash A, Jeffryes M, Bateman A et al. The HMMER web server for protein sequence similarity search. Curr Protoc Bioinformatics. 2017; 60: 3.15.11-13.15.23.

[79] Finn RD, Bateman A, Clements J et al. Pfam: the protein families database. Nucleic Acids Res. 2014; 42: D222-30.

[80] Tarailo-Graovac M, Chen N . Using RepeatMasker to identify repetitive elements in genomic sequences. Curr Protoc Bioinformatics. 2009; Chapter 4: Unit 4.10.

[81] Claire H, Arnoux S, Moisset M et al. PASTEC: an automatic transposable element classification tool. PLoS One. 2014; 9: e91929.

[82] Xu Z, Wang H . LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Res. 2007; 35: W265-8.

[83] Ou S, Jiang N . LTR_retriever: a highly accurate and sensitive program for identification of long terminal repeat retrotransposons. Plant Physiol. 2018; 176: 1410-22.

[84] Katoh K, Standley DM . MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013; 30: 772-80.

[85] Price MN, Dehal PS, Arkin AP . FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009; 26: 1641-50.

[86] Emms DM, Kelly S . OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 2019; 20: 238.

[87] Han MV, Thomas GW, Lugo-Martinez J et al. Estimating gene gain and loss rates in the presence of error in genome assembly and annotation using CAFE 3. Mol Biol Evol. 2013; 30: 1987-97.

[88] Sanderson MJ . r8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. Bioinformatics. 2003; 19: 301-2.

[89] Tang H, Bowers JE, Wang X et al. Synteny and collinearity in plant genomes. Science. 2008; 320: 486-8.

[90] Chen S, Zhou Y, Chen Y et al. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018; 34: i884-90.

[91] Kim D, Langmead B, Salzberg SL . HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015; 12: 357-60.

[92] Trapnell C, Roberts A, Goff L et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks. Nat Protoc. 2012; 7: 562-78.

[93] McKenna A, Hanna M, Banks E et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010; 20: 1297-303.

[94] Wang K, Li M, Hakonarson H . ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010; 38: e164.

[95] Li H . Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018; 34: 3094-100.

[96] Anders S, Huber W . Differential expression analysis for sequence count data. Genome Biol. 2010; 11: R106.

[97] Ernst J, Bar-Joseph Z . STEM: a tool for the analysis of short time series gene expression data. BMC Bioinformatics. 2006; 7: 191.

[98] Trincado JL, Entizne JC, Hysenaj G et al. SUPPA2: fast, accurate, and uncertainty-aware differential splicing analysis across multiple conditions. Genome Biol. 2018; 19: 40.