Genome-wide analysis of changes in miRNA and target gene expression reveals key roles in heterosis for Chinese cabbage biomass

Abstract

Heterosis is a complex phenomenon in which hybrids show better phenotypic characteristics than their parents do. Chinese cabbage (Brassica rapa L. spp. pekinensis) is a popular leafy crop species, hybrids of which are widely used in commercial production; however, the molecular basis of heterosis for biomass of Chinese cabbage is poorly understood. We characterized heterosis in a Chinese cabbage F₁ hybrid cultivar and its parental lines from the seedling stage to the heading stage; marked heterosis of leaf weight and biomass yield were observed. Small RNA sequencing revealed 63 and 50 differentially expressed microRNAs (DEMs) at the seedling and early-heading stages, respectively. The expression levels of the majority of miRNA clusters in the F₁ hybrid were lower than the mid-parent values (MPVs). Using degradome sequencing, we identified 1,819 miRNA target genes. Gene ontology (GO) analyses demonstrated that the target genes of the MPV-DEMs and low parental expression level dominance (ELD) miRNAs were significantly enriched in leaf morphogenesis, leaf development, and leaf shaping. Transcriptome analysis revealed that the expression levels of photosynthesis and chlorophyll synthesis-related MPV-DEGs (differentially expressed genes) were significantly different in the F₁ hybrid compared to the parental lines, resulting in increased photosynthesis capacity and chlorophyll content in the former. Furthermore, expression of genes known to regulate leaf development was also observed at the seedling stage. Arabidopsis plants overexpressing BrGRF4.2 and bra-miR396 presented increased and decreased leaf sizes, respectively. These results provide new insight into the regulation of target genes and miRNA expression patterns in leaf size and heterosis for biomass of B. rapa.

References

[1] Chen, Z.J. Genomic and epigenetic insights into the molecular bases of heterosis. Nat. Rev. Genet. 14, 471-482 (2013).

[2] Schnable, P. & Springer, N. Progress toward understanding heterosis in crop plants. Annu. Rev. Plant Biol. 64, 71-88 (2013).

[3] Birchler, J., Yao, H., Chudalayandi, S., Vaiman, D. & Veitia, R . Heterosis. Plant Cell 22, 2105-2112 (2010).

[4] Bruce, A.B. The mendelian theory of heredity and the augmentation of vigor. Science 32, 627-628 (1910).

[5] East, E.M. Heterosis. Genetics 21, 375 (1936).

[6] Schnell, F. & Cockerham, C. Multiplicative vs. arbitrary gene action in heterosis. Genetics 131, 461-469 (1992).

[7] Hochholdinger, F. & Höcker, N. Towards the molecular basis of heterosis. Trends Plant Sci. 12, 427-432 (2007).

[8] He, G., Elling, A. & Deng, X. The epigenome and plant development. Annu. Rev. Plant Biol. 62, 411-435 (2010).

[9] Lippman, Z. & Zamir, D. Heterosis: revisiting the magic. Trends Genet. 23, 60-66 (2007).

[10] Jones-Rhoades, M.W., Bartel, D.P. & Bartel, B. MicroRNAs and their regulatory roles in plants. Annu. Rev. Plant Biol. 57, 19-53 (2006).

[11] Voinnet, O. Origin, biogenesis, and activity of plant microRNAs. Cell 136, 669-687 (2009).

[12] Shen, Y. et al. Analysis of transcriptional and epigenetic changes in hybrid vigor of allopolyploid Brassica napus uncovers key roles for small RNAs. Plant J. 91, 874-893 (2017).

[13] Ha, M. et al. Small RNAs serve as a genetic buffer against genomic shock in Arabidopsis interspecific hybrids and allopolyploids. Proc. Natl Acad. Sci. USA 106, 17835-17840 (2009).

[14] Li, A. et al. mRNA and small RNA transcriptomes reveal insights into dynamic homoeolog regulation of allopolyploid heterosis in Nascent Hexaploid Wheat. Plant Cell 26, 1878-1900 (2014).

[15] Xu, M. et al. Developmental functions of miR156-regulated SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes in Arabidopsis thaliana. PLoS Genet. 12, e1006263 (2016).

[16] Ben-Gera, H. et al. Auxin-mediated lamina growth in tomato leaves is restricted by two parallel mechanisms. Plant J. 86, 443-457 (2016).

[17] Bresso, E.G., Chorostecki, U., Rodriguez, R.E., Palatnik, J.F. & Schommer, C. Spatial control of gene expression by miR319-regulated TCP transcription factors in leaf development. Plant Physiol. 176, 1694-1708 (2018).

[18] Omidbakhshfard, M.A., Proost, S., Fujikura, U. & Mueller-Roeber, B. Growth-regulating factors (GRFs): a small transcription factor family with important functions in plant biology. Mol. Plant 8, 998-1010 (2015).

[19] Groszmann, M. et al. Intraspecific Arabidopsis hybrids show different patterns of heterosis despite the close relatedness of the parental genomes. Plant Physiol. 166, 265-280 (2014).

[20] Guo, M. et al. Maize ARGOS1 (ZAR1) transgenic alleles increase hybrid maize yield. J. Exp. Bot. 65, 249-260 (2013).

[21] Fujimoto, R., Taylor, J.M., Shirasawa, S., Peacock, W.J. & Dennis, E.S. Heterosis of Arabidopsis hybrids between C24 and Col is associated with increased photosynthesis capacity. Proc. Natl Acad. Sci. USA 109, 7109-7114 (2012).

[22] Zhang, S.F., Song, Z.H. & Zhao, X.Y. Breeding of interactive genic sterile line in chinese cabbage (Brassica pekinensis Rupr) and utilization model. Acta. Hort. Sin. 17, 117-125 (1990).

[23] Cheng, F. et al. Subgenome parallel selection is associated with morphotype diversification and convergent crop domestication in Brassica rapa and Brassica oleracea. Nat. Genet. 48, 1218-1224 (2016).

[24] Li, H. et al. Transcriptome and DNA methylome reveal insights into yield heterosis in the curds of broccoli (Brassica oleracea L var. italic). BMC Plant Biol. 18, 168 (2018).

[25] Saeki, N. et al. Molecular and cellular characteristics of hybrid vigour in a commercial hybrid of Chinese cabbage. BMC Plant Biol. 16, 45 (2016).

[26] Kong, X. et al. Transcriptome analysis of biological pathways associated with heterosis in Chinese cabbage. Genomics 6, 4732-4741 (2020).

[27] Yi, H. et al. Genome-wide analysis of heterosis-related genes in non-heading Chinese cabbage. J. Plant Biotechnol. 44, 208-219 (2017).

[28] Liu, T. et al. Enhanced photosynthetic activity in pak choi hybrids is associated with increased grana thylakoids in chloroplasts. Plant J. 103, 2211-2224 (2020).

[29] Su, T. et al. Development of cost-effective single nucleotide polymorphism marker assays for genetic diversity analysis in Brassica rapa. Mol. Breed. 38, 42 (2018).

[30] Addo-Quaye, C., Miller, W. & Axtell, M.J. CleaveLand: a pipeline for using degradome data to find cleaved small RNA targets. Bioinformatics 25, 130-131 (2008).

[31] Gonzalez, N. & Beemster, G. David and Goliath: what can the tiny weed Arabidopsis teach us to improve biomass production in crops? Curr. Opin. Plant Biol. 12, 157-164 (2009).

[32] Breuninger, H. & Lenhard, M. Control of tissue and organ growth in plants. Curr. Top. Dev. Biol. 91, 185-220 (2010).

[33] Ni, Z. et al. Altered circadian rhythms regulate growth vigour in hybrids and allopolyploids. Nature 457, 327-331 (2009).

[34] Darwin, C. The Effects of Cross and Self-Fertilisation in the Vegetable Kingdom (John Murray, 1876).

[35] Wang, C. et al. Dissecting a heterotic gene through GradedPool-Seq mapping informs a rice-improvement strategy. Nat. Commun. 10, 2982 (2019).

[36] Huang, X. et al. Genomic architecture of heterosis for yield traits in rice. Nature 537, 629-633 (2016).

[37] Liu, H. et al. Genome-wide identification and analysis of heterotic loci in three maize hybrids. Plant Biotechnol. J. 18, 185-194 (2019).

[38] Riedelsheimer, C. et al. Genomic and metabolic prediction of complex heterotic traits in hybrid maize. Nat. Genet. 44, 217-220 (2012).

[39] Kenan-Eichler, M. et al. Wheat hybridization and polyploidization results in deregulation of small RNAs. Genetics 188, 263-272 (2011).

[40] Krieger, U. et al. The flowering gene single flower truss drives heterosis for yield in tomato. Nat. Genet. 42, 459-463 (2010).

[41] Ng, D. et al. A role for CHH methylation in the parent-of-origin effect on altered circadian rhythms and biomass heterosis in Arabidopsis intraspecific hybrids. Plant Cell 26, 2430-2440 (2014).

[42] Han, J., Hou, X., Shi, G., Geng, J. & Deng, X. Genetic model analysis of leaf-weight traits in non-heading Chinese cabbage (Brassica campestris ssp. chinensis Makino). Hereditas 29, 1149-1153 (2007).

[43] Li, D. et al. Integrated analysis of high-throughput sequencing data shows abscisic acid-responsive genes and miRNAs in strawberry receptacle fruit ripening. Hortic. Res. 6, 26 (2019).

[44] Rodriguez, R.E., Debernardi, J.M. & Palatnik, J.F. Morphogenesis of simple leaves: regulation of leaf size and shape. Wires. Dev. Biol. 3, 41-57 (2014).

[45] Debernardi, J. et al. Post-transcriptional control of GRF transcription factors by microRNA miR396 and GIF co-activator affects leaf size and longevity. Plant J. 79, 413-426 (2014).

[46] Fhdiana, I., Meehan, L., Harkins, K., Chory, J. & Rodermel, S. Lhcb transcription is coordinated with cell size and chlorophyll accumulation. Plant Physiol. 112, 953-963 (1996).

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

[48] Love, M.I., Wolfgang, H. & Simon, A. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).

[49] German, M.A., Luo, S., Schroth, G., Meyers, B.C. & Green, P.J. Construction of parallel analysis of RNA ends (PARE) libraries for the study of cleaved miRNA targets and the RNA degradome. Nat. Protoc. 4, 356-362 (2009).

[50] Patel, R.K. & Jain, M. NGS QC Toolkit: a toolkit for quality control of next generation sequencing data. PLoS ONE 7, e30619 (2012).

[51] Kim, D., Langmead, B. & Salzberg, S.L. HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357-360 (2015).

[52] Mihaela, P. et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290-295 (2015).

[53] Trapnell, C. et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28, 511-515 (2010).

[54] Young, M.D., Wakefield, M.J., Smyth, G.K. & Oshlack, A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 11, 1-12 (2010).

[55] Li, P. et al. Genome-wide analysis of mRNA and lncRNA expression and mitochondrial genome sequencing provide insights into the mechanisms underlying a novel cytoplasmic male sterility system, BVRC-CMS96, in Brassica rapa. Theor. Appl. Genet. 133, 2157-2170 (2020).

[56] Livak, K.J. & Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 25, 402-408 (2001).

[57] Wang, X. et al. The genome of the mesopolyploid crop species Brassica rapa. Nat. Genet. 43, 1035-1039 (2011).

[58] Clough, S. & Bent, A. Floral dip: a simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana. Plant J. 16, 735-743 (1998).

[59] Li, P. et al. BrLAS, a GRAS transcription factor from Brassica rapa, is involved in drought stress tolerance in transgenic Arabidopsis. Front. Plant Sci. 9, 1792 (2018).