- Research article
- Open Access
Transcriptome and DNA methylome reveal insights into yield heterosis in the curds of broccoli (Brassica oleracea L var. italic)
BMC Plant Biologyvolume 18, Article number: 168 (2018)
Curds are the main edible organs, which exhibit remarkable yield heterosis in F1 hybrid broccoli. However, the molecular basis underlying heterosis in broccoli remains elusive.
In the present study, transcriptome profiles revealed that the hybridization made most genes show additive expression patterns in hybrid broccoli. The differentially expressed genes including the non-additively expressed genes detected in the hybrid broccoli and its parents were mainly involved in light, hormone and hydrogen peroxide-mediated signaling pathways, responses to stresses, and regulation of floral development, which suggested that these biological processes should play crucial roles in the yield heterosis of broccoli. Among them, light and hydrogen peroxide-mediated signaling pathways represent two novel classes of regulatory processes that could function in yield or biomass heterosis of plants. Totally, 53 candidate genes closely involved in curd yield heterosis were identified. Methylome data indicated that the DNA methylation ratio of the hybrids was higher than that of their parents. However, the DNA methylation levels of most sites also displayed additive expression patterns. These sites with differential methylation levels were predominant in the intergenic regions. In most cases, the changes of DNA methylation levels in gene regions did not significantly affect their expression levels.
The differentially expressed genes, the regulatory processes and the possible roles of DNA methylation modification in the formation of curd yield heterotic trait were discovered. These findings provided comprehensive insights into the curd yield heterosis in broccoli, and were significant for breeding high-yield broccoli varieties.
Heterosis or hybrid vigor refers to the phenomenon in which the F1 hybrids have higher biomass, yields, growth rates and resistance to environmental stresses than their parents [1, 2]. In agriculture, heterosis has been used extensively in breeding high yield corps. The plantations of hybrid rice, maize, sorghum, and wheat all worldwide have significantly contributed to increasing the grain yield and ensuring global food security . Although heterosis has been successfully utilized in crop breeding practices, the molecular mechanism of heterosis remains poorly understood, and the explorations involved in how and why the F1 hybrids exhibit superior performance in plants never stopped since this phenomenon was discovered . Previous studies have proposed three main genetic hypotheses, namely, dominance , over-dominance [6, 7] and epistasis , to explain the formation of heterosis in diverse plant species. The dominance hypothesis posited that the inferior parental alleles in the hybrids are complemented by the superior or dominant alleles from the other parent. The over-dominance hypothesis suggested that heterosis arises from the allelic interactions within each of many genetic loci. The epistasis hypothesis emphasized that the interactions between different parental genes in hybrids lead to heterosis. However, these three models were mainly theoretical, and none of them could fully explain the heterosis in plants. With the development of quantitative genetics, quantitative trait locus (QTL) mapping studies were widely conducted to identify the gene loci that controlled the heterotic traits. These works confirmed the involvement of many QTLs in biomass or yield heterosis in several crops, such as rice and maize [9,10,11], thereby providing potential molecular markers for further cloning the genes regulating this phenomenon. In recent years, transcriptome analysis based on high-throughput sequencing technologies offered a powerful strategy for exploring genome-wide gene expression profiles and confirming differentially expressed genes (DEGs) in hybrids and their parents, which may be directly associated with the formation of heterosis. In model plants Arabidopsis and rice, a series of candidate genes regulating heterotic traits were identified by comparative transcriptome data [12,13,14,15]. The possible relationship and interactions among these genes were uncovered, thus providing comprehensive insights into the molecular basis of heterosis. In addition to genetic regulations, more and more evidence indicated that epigenetic modifications such as DNA methylation, small RNA-mediated gene regulation and histone modification also play crucial roles in the formation of heterosis [1, 16,17,18,19]. All of these investigations have broadened the knowledge on the heterosis in plants. However, most of these studies were focused on Arabidopsis and various crops, such as rice, maize and sorghum. As such, the molecular basis of heterosis in other plant species remains largely unknown.
Broccoli (Brassica oleracea L var. italic) is one of the major vegetables among the B. oleracea varieties, and it is planted worldwide. Unlike most B. oleracea varieties, broccoli forms a specialized organ named curd during floral development, which composed of many flower buds and shortened inflorescence branches. The curds are the main edible organs, and the yield of curds directly determines the economic value of broccoli. Heterosis is popular in broccoli. Compared with their parents, hybrid broccoli exhibits significantly increased biomass. Especially, the curd weight of hybrids is considerably heavier than those of the parents. However, the basis of heterosis in broccoli is scarce. In the present study, transcriptome analysis was conducted. A series of genes were confirmed to show significantly differential expression levels in the hybrid broccoli and their parents, and the functions and possible relationships of the DEGs were explored. In addition, the DNA methylation patterns of the hybrids and their parents were determined. Sites with differential DNA methylation levels were also identified. These findings provided comprehensive insights into the formation of yield heterosis, and gave novel clues that could potentially assist hybrid breeding in broccoli.
Remarkable biomass heterosis in the hybrid broccoli
Two hybrid combinations in which their F1 hybrids both showed obvious yield heterosis were analyzed in broccoli. Compared with their parental lines, the broccoli F1 hybrids displayed significantly larger curds, bigger leaves and stronger roots (Fig. 1). The total individual fresh weight in the hybrid NKR-06 was 3.6 times of the mid-parent value (MPV), and the curds of NKR-06 were 400% heavier than the MPV. Similarly, the total individual fresh weight in the hybrid Bro-12 was 120% heavier than the MPV. The weight of the curds in Bro-12 was 6.8 times of the MPV (Fig. 1). The results identified that the hybrids NKR-06 and Bro-12 demonstrated significant biomass heterosis. Specifically, the curds of both F1 hybrids exhibited obviously stronger heterosis in yield trait than those of their parental lines.
More genes were actively expressed in the hybrids than those in their parents
Over 50 million high-quality clean reads were detected in each sample, most of which were mapped to the Brassica oleracea var. oleracea reference genome and transcriptome (Additional file 1: Table S1). In the NKR-04/− 05/− 06 hybrid triad, 47, 962 genes showed active expression, among which 44, 040 and 44, 812 genes were detected in the parental lines NKR-04 and NKR-05, respectively. More genes (46, 442) were detected in the hybrid NKR-06 than those in the parental lines. Similarly, in the Bro-10/− 11/− 12 hybrid triad, the number of genes showing transcriptional expression in the hybrid Bro-12 (46, 810) was greater than those in the maternal line Bro-10 (45, 093) and paternal line Bro-11(45, 529). However, in each hybrid combinations, the number of actively expressed genes was similar between two parental lines (Fig. 2a, Additional file 1: Table S2). Pairwise comparison analysis further confirmed that the majority of genes (over 91%) showed common expression between the two parental lines in each hybrid triad, whereas genes showing specific expression in the paternal line were slightly greater than those in the material line (Fig. 2b). Compared with the parental lines, more genes showed specific expression in the hybrids. For example, 3.19% and 3.09% of the expressed genes specifically expressed in the hybrid NKR-06 and Bro-12, respectively. Only 0.98–1.59% of the expressed genes showed specific expression in the parental lines (Fig. 2b). In addition, although the two hybrids NKR-06 and Bro-12 generated from different hybrid combinations, more genes (93.4%) shared common expression in them than those in their parental lines (85.9%) (Fig. 2b). Moreover, correlation and clustering analysis of all of the expressed genes confirmed that the two hybrids achieved the highest correlation coefficient (R = 0.98) (Fig. 2c). These results indicated that although the parental lines are different, the hybridization in each broccoli hybrid combination activated more gene expression in the F1 hybrids compared with their parents.
DEGs with increasing trend of expression levels were predominant
Besides the count of actively expressed genes in two hybrid triads, the differential expression profiles of genes were further studied. In the NKR-04 and NKR-05 hybrid combination, 789 genes displayed significantly differential expression levels between the two parental lines (corrected p-value < 0.05). Among them, the proportion of genes showing significantly higher and lower expression levels in the maternal line NKR-04 than those in the paternal line NKR-05 was similar (Table 1). However, the DEGs between the hybrid NKR-06 and either of the parental lines were obviously less than those detected in two parents (Table 1). In total, 936 genes showed significantly differential expression levels between the hybrid NKR-06 and the parents (corrected p-value < 0.05), among which most of the DEGs (over 80%) showed higher expression levels in the hybrid NKR-06 than those in its parents. Moreover, the expression levels 115 of the 936 genes were significantly different from those in both parents. Compared with the parents, the majority of the 115 genes (79 out of 115) showed higher expression levels in the hybrid NKR-06 (Table 1, Fig. 3a and Additional file 2: Data set 1). Similarly, the DEGs between the two parental lines Bro-10 and Bro-11 were considerably greater than those detected between the hybrid Bro-12 and either of its parental lines. In total, 932 DEGs were detected between the hybrid Bro-12 and its parents, of which the majority of the DEGs showed higher expression levels in the hybrid than those in the parents. 153 of the 932 DEGs showed significantly higher expression levels in the hybrid Bro-12 than those in its parents, and most of them increased expression levels in the hybrid (Table 1, Fig. 3b, Additional file 2: Data set 2). Among the 115 DEGs detected in the NKR-04/− 05/− 06 hybrid triad and the 153 DEGs detected in the Bro-10/− 11/− 12 hybrid triad, 53 DEGs were shared in both hybrid triads. The majority of the 53 DEGs (50 out of 53) showed significantly higher expression levels in the hybrids than those in their parents (Fig. 3, Additional file 2: Data set 3). These results further confirmed that the hybridization not only activated more gene expression in the broccoli F1 hybrids, but also made more genes displayed significantly higher expression levels in the hybrids than those in their parents.
Few genes showed non-additive expression patterns
Genes in the hybrid having the expression levels significantly deviated from the MPV were designated as non-additively expressed genes. In the NKR-04/− 05/− 06 hybrid triad, among the 47, 962 genes only 60 genes (0.13%) showed non-additive expression patterns (Additional file 2: Data set 4), implying that the vast majority of the genes in this hybrid triad exhibited additive expression patterns. Further analysis indicated that 51 out of the 60 non-additively expressed genes (85%) exhibited higher expression levels than the MPV (Fig. 4a, Additional file 2: Data set 4). Similarly, a small proportion of genes (198/48, 545 = 0.4%) in the Bro-10/− 11/− 12 hybrid triad showed non-additive expression patterns (Additional file 2: Data set 4). Among the 198 non-additively expressed genes, the expression levels of 163 genes (83%) were higher than the MPV (Fig. 4b). In total, 258 genes displayed non-additive expression patterns, among which 37 genes were shared in the two hybrid triads. Moreover, the expression levels of 35 out of the 37 genes were significantly higher than those of their two parents (Additional file 2: Data set 4). These results demonstrated that the majority of genes exhibited additive expression patterns in the two hybrid triads. In the non-additively expressed genes, genes exhibiting higher expression levels than the MPV were predominant.
Gene regulatory networks involved in light, hormone and stress responses were altered in the hybrids
All DEGs were functional classification by gene ontology (GO) analysis (Additional file 2: Data set 5 and Data set 6). GO terms showed significantly enriched in biological processes were predominant, which were further clustered by ReviGO. The results indicated that GO terms involved in responses to environmental stimuli in vivo and vitro, such as “light response”, “response to hormone”, “response to biotic and abiotic stress”, “response to hydrogen peroxide” and “water and ion homeostasis”, were significantly overrepresented. In light response processes, the majority of GO terms were related to red light signaling pathways, followed by blue light signaling pathways, and response to light intensity and UV. Abscisic acid, jasmonic acid, salicylic acid, brassinosteroid and auxin-mediated signaling pathways were predominant in hormone response processes. Signaling pathways involved in response to various biotic and abiotic stresses were overrepresented in stress response processes. Responses to water deprivation and ion transport, especially iron, cadmium and manganese transport, were overrepresented in GO terms involved in water and ion homeostasis (Table 2, Additional file 2: Data set 5 and Data set 6). Gene expression profile assay indicated that the majority of the DEGs targeting these overrepresented response processes exhibited higher expression levels in the hybrids than those in both or either of the parents, except for the DGEs identified in the “response to biotic and abiotic stress” process (Fig. 5).
Gene regulation networks involved in organ growth and development were also altered in the hybrids
In addition to gene regulation networks involved in responses to environmental stimuli, which were significantly adjusted in the hybrids, a large proportion of the GO terms targeted by these DEGs displayed prominently enriched in pathways associated with growth and development, such as “leaf, root and shoot growth and development”, “floral development”, “embryo development”, “cell division and differentiation” and“photosynthesis and chloroplast development” (Table 2, Additional file 2: Data set 5 and Data set 6). Among them, GO terms involved in floral development were mainly targeted by the DEGs, which exhibited significantly differential expression levels in the hybrid compared with both parents, and/or non-additive expression patterns. The majority of these genes were annotated to cytochrome P450, type III polyketide synthases, glucan endo-1, 3-beta-glucosidase, fatty acyl-CoA reductase, non-specific lipid-transfer proteins, lipid-transfer proteins, AT-hook motif nuclear-localized proteins, Tetraketide alpha-pyrone reductase, 4-coumarate-CoA ligase, alpha-dioxygenase, monocopper oxidase, ubiquitin-related proteins, tetraketide alpha-pyrone reductase and tapetum-specific protein, which played important roles in floral development (Additional file 2: Data set 3). Another prominently overrepresented biological processes involved in growth and development were “biosynthesis and metabolism of sugar, fatty acid and carbohydrate” and “DNA/RNA and protein synthesis and modification” (Table 2, Additional file 2: Data set 5 and Data set 6), implying that the biosynthesis and metabolism processes in the hybrid broccoli were more active than those in their parents. Gene expression profile assay also indicated that the majority of the DEGs targeting these overrepresented growth and development processes exhibited higher expression levels in the hybrids than those in both or either of the parents (Additional file 3: Figure S1).
Expression patterns of DEGs were confirmed by qRT-PCR
Sixteen genes that showed significantly differential expression levels in the hybrids and their parents were selected to conduct quantitative real-time RT-PCR (qRT-PCR) identification. The expression patterns of 15 out of the 16 DEGs detected by qRT-PCR were highly consistent with the data produced by high-throughput sequencing (Fig. 6). In addition, the expression patterns of eight DEGs, including fatty acyl-CoA reductase, peroxidase 9, phytochrome A, cytochrome P450, glucan endo-1,3-beta-glucosidase A6, metal transporter Nramp1, alpha-dioxygenase 1, and beta-D-glucopyranosyl abscisate beta-glucosidase, were further confirmed in another hybrid triad of broccoli (i.e., BNR-M (male), BNR-F (female) and BNR-H (F1 hybrid)). The F1 hybrid BNR-H also exhibited remarkable curd yield heterosis. All eight genes displayed similar expression patterns in the BNR-M/-F/-H hybrid triad to those in the NKR-04/− 05/− 06 and the Bro-10/− 11/− 12 hybrid triads (Additional file 3: Figure S2).
F1 hybrids showed higher DNA methylation ratio than their parents
The DNA methylation at the two popular sites (i.e., CG and CHG. H = T or A) in the whole genomes of the two hybrids and their parents were analyzed by using MethylRAD method (Additional file 2: Data set 7). The total DNA methylation ratios (methylated CG and CHG sites/total CG and CHG sites) were 14.23%, 12.66% and 12.87% in NKR-06, NKR-04, and NKR-05, respectively. The total DNA methylation ratios were 16.60%, 15.70%, and 15.54% in Bro-12, Bro-10, and Bro-11, respectively (Fig. 7a, b, Additional file 1: Table S3). These data confirmed that the hybrids had slightly higher total DNA methylation ratio than their corresponding parents as well as the MPV (Fig. 7a, b). In addition, the data indicated that the DNA methylation ratio at the CG sites was higher than that at the CHG sites in each chromosome, except for the chromosome 9 of the hybrid Bro-12 (Additional file 1: Table S3). However, the distribution patterns of most methylated sites at the different elements of genomes were similar in the hybrids and their parents. Both the CG and CHG sites at the intergenic regions appeared to easily occur methylation, followed by the regions at the gene body, TSS1500 (upstream 1500 bp of transcription start sites (TSS)), intron, 1stexon and TSS200 (upstream 200 bp of TSS) (Additional file 3: Figure S3). Among them, the DNA methylation levels of the gene regions, including the TSS, the gene coding region, and the transcription termination site (TTS) were explored. The results confirmed that although the DNA methylation levels differed between the two parents of each hybrid triad, the DNA methylation levels of the hybrids at either the CG sites or the CHG sites in the gene regions were similar to the MPV (Fig. 8, Additional file 1: Table S4, Additional file 3: Figure S4).
Sites exhibiting differential methylation levels between the parents were greater than those between the hybrids and the parents
In the NKR-04/− 05/− 06 hybrid triad, compared with the parental lines NKR-04 and NKR-05, 1, 304 and 2, 732 CG sites showed significantly differential methylation levels in the hybrid NKR-06, respectively. These sites represented 464 and 783 genes with significantly methylation levels, respectively. By contrast, 4, 254 CG sites exhibited significantly differential methylation levels between two parents, which made 1, 322 genes show significantly differential methylation levels. At the CHG sequence regions, 2, 536 and 4, 492 sites exhibited significantly differential methylation levels between NKR-06 and NKR-04, and between NKR-06 and NKR-05, indicating that 296 and 514 genes showed significantly differential methylation levels, respectively. However, as many as 5, 355 CHG sites exhibited significantly differential methylation levels between the parents NKR-04 and NKR-05. These sites represented 664 genes showing significantly differential methylation levels (Additional file 2: Data set 8 and Data set 9). Similarly, in the Bro-10/− 11/− 12 hybrid triad, the sites and genes showing significantly differential methylation levels between two parents were greater than those between the hybrid and each of the parental line (Additional file 2: Data set 10 and Data set 11).
Sites with differential methylation levels were predominant in the intergenic regions
The distributions of the sites with differential methylation levels were explored. The data indicated that between the parents or between the hybrid and the parents, the sites with differential methylation levels were mostly enriched in the intergenic regions, followed by the gene coding regions (1stExon + other extrons), the TSS1500 (upstream 1500 bp of TSS), intron and TSS200 (upstream 200 bp of TSS) (Figs. 7c-f, Additional file 1: Table S5).
Few DEGs exhibited differential methylation levels
The relationships between the differential expression levels of genes and their DNA methylation levels were analyzed. In the NKR-04/− 05/− 06 hybrid triad, 753 genes showed significantly differential methylation levels between NKR-06 and NKR-04. Genes with significantly differential methylation levels between NKR-06 and NKR-05 were 1, 280. Between two parents, 1, 908 genes showed significantly differential methylation levels. Comparative transcriptome analysis indicated that 571 genes showed differential expression levels between NKR-06 and NKR-04. 480 and 798 genes showed differential expression levels between NKR-06 and NKR-05, and between two parents, respectively. However, compared with NRK-04 and NKR-05, only four and 12 genes simultaneously showed differential expression levels and DNA methylation levels in the hybrid NKR-06, respectively. Between NRK-04 and NKR-05, 35 DEGs showed differential methylation levels. Similarly, in the Bro-10/− 11/− 12 hybrid triad, only seven, 5 and 36 genes simultaneously exhibited differential methylation levels and expression levels between Bro-12 and Bro-10, between Bro-12 and Bro-11, and between two parents Bro-10 and Bro-11, respectively (Fig. 9, Additional file 1: Table S6). Among these DEGs which also exhibited differential methylation levels, a large proportion of them (60 DEGs) appeared to be negatively correlated with their DNA methylation levels. However, there are over 40 DEGs that their expression levels trended to show positively correlated with their DNA methylation levels (Fig. 9, Additional file 1: Table S6).
Transcriptome and other omics have provided powerful strategies for elucidating the molecular basis of heterosis in diverse plants [12, 15, 20]. In the present study, the transcriptome data confirmed that the hybridization activated more genes in the hybrid broccoli with obviously yield heterotic trait than those in their parents. However, the whole gene expression profiles of the hybrids and their parents were similar. Only a small proportion of the genes exhibited significantly differential expression levels in the hybrids and their parents. Investigations involved in the heterosis produced by intraspecific hybridization in Arabidopsis, rice and other plant species also confirmed that the majority of genes showed similar expression profiles in the hybrids and their parents, and that the expression levels of most genes were near the MPV [12, 21, 22]. These findings implied that extensive and dramatic changes of the whole gene expression profiles may not be required in the formation of certain heterotic traits through intraspecific hybridization, although the hybridization could make some genes activate or inhibit expression in the hybrids. The additive expression balance of most genes in the hybrids establishes an important foundation for the formation of heterosis. Nevertheless, the genes exhibiting significantly differential expression profiles in the hybrids and their parents, especially the non-additively expressed genes, provided valuable clues to elucidate the formation of heterotic traits. The functional annotation indicated that a large proportion of the DEGs were involved in response processes of environmental stimuli. Among them, “light response”, “response to hormone”, “response to biotic and abiotic stress” and “response to hydrogen peroxide” were overrepresented. The hormone-mediated signaling pathways targeted by the DEGs were mainly involved in abscisic acid, jasmonic acid, salicylic acid, brassinosteroid and auxin. All these hormones have been demonstrated to play important roles in the formation of heterosis [23,24,25,26,27]. Pathways involved in multiple biotic and abiotic stresses were also confirmed to involve in the heterosis [25, 28, 29]. These findings indicated that the DEGs involved in hormone-mediated regulation pathways and biotic and abiotic stresses should also play important roles in the formation of heterosis of broccoli. Interestingly, some DEGs were confirmed to function in light response and hydrogen peroxide response, which were few reported to be associated with the heterosis in plants. Light is not only the energy of photosynthesis but also an important signal regulator to function in photomorphogenesis and other growth and development processes [30,31,32]. A series of GO terms, such as “response to light intensity”, “response to red or far red light”, “detection of visible light” and “response to blue light”, displayed significantly enriched in light-mediated response processes in present study. Remarkably, a recent report indicated that a key factor of light signaling networks, PIF4, could contribute to heterosis by regulating the auxin pathway . In addition, hydrogen peroxide is an important signaling molecule. It can interplay with diverse phytohormone signals, such as abscisic acid, salicylic acid, jasmonic acid, auxin and brassinosteroid, to regulate plant developmental processes and stress responses [33,34,35]. Recent reports confirmed that hydrogen peroxide could induce the oxidation of the BZR1 transcription factor, a master regulator of brassinosteroid signal pathway. Oxidative modification enhances BZR1 transcriptional activity by promoting its interaction with key regulators in the auxin signal pathway and light signal pathway, including ARF6 and PIF4 . In the present study, multiple pathways involved in hydrogen peroxide response were significantly adjusted in the hybrid broccoli. These findings suggested that besides processes involved in hormone responses, and stress responses, light and hydrogen peroxide-mediated signaling pathways may represent two novel classes of regulatory processes that function in yield heterosis in broccoli. Among them, hydrogen peroxide as a key signaling molecule should play crucial roles in the formation of yield heterosis of broccoli by interplaying with diverse hormone-mediated signaling pathways as well as light signaling pathway.
Besides the response processes of environmental stimuli, several pathways involved in growth and development, such as “leaf, root and shoot growth and development”, “floral development”, “embryo development” and “cell division and differentiation”, also exhibited overrepresented in the GO terms targeted by the DEGs. Among them, the overrepresented pathways involved in floral development were mainly targeted by those DEGs, which exhibited significantly differential expression levels in the hybrid broccoli compared with both parents, and/or non-additive expression patterns. In total, 53 such DEGs, including 35 non-additively expressed genes, were identified in two broccoli hybrid triads. 13 out of the 53 DEGs were confirmed to play important roles in the anther and/or pollen development. For example, CYP703A2 and CYP704B1 are essential for the sporopollenin formation in diverse plant species [37,38,39]. Acyl-CoA Synthetase 5 (ACOS5) encodes an acyl-CoA synthetase protein that is required for sporopollenin monomer biosynthesis [28, 40]. Plant specific type III polyketide synthase (PSK) genes, PSKA and PSKB, encode hydroxyalkyl a-pyrone synthases to catalyze the condensation of malonyl-CoA units with various CoA ester starter molecules that are essential for the pollen development [41, 42]. Tetraketide alpha-pyrone reductase 1 (TKPR1) and TKPR2 are confirmed to function together with ACOS5 and PKSA/B to hydroxylate a-pyrone compounds that are important sporopollenin precursors of pollen development . ABC transporter G family member 26 (ABCG26), an ATP binding cassette transporter, plays crucial roles in Arabidopsis pollen exine formation and male fertility [44, 45]. The transposable element silencing via AT-hook (TEK) encodes an AT-hook nuclear localized (AHL) protein in Arabidopsis that is highly homologous with the AHL16 in broccoli. TEK is highly expressed in tapetum during the tetrad stage and determines nexine formation in the pollen wall . Moreover, the data indicated that all of these DEGs involved in floral development showed significantly over-parental expression patterns in the hybrids. The curd of broccoli mainly composed of flower buds. Accordingly, it should speculate that significantly increased expression levels of these genes could promote the growth and development of the flower organ, and therefore may function in the formation of yield heterosis of curds in the hybrid broccoli.
Genetic diversity even very small difference between two parents may be crucial but not the only factor affecting the formation of heterotic traits in F1 hybrids. Epigenetic modification, such as DNA methylation, has been confirmed to function in this phenomenon [15, 47, 48]. The data indicated that the hybridization made more sites facilitate cytosine methylation in the hybrid broccoli. However, the DNA methylation levels of most sites in the hybrids were near the MPV, and showed an additive character. The sites with differential methylation levels were predominant in the intergenic regions. Only few DEGs simultaneously exhibited differential methylation levels. In Arabidopsis, the heterotic hybrids showed increased DNA methylation ratios compared with the parents. The increased DNA methylation of the hybrid genomes also predominantly occurred in the non-protein-coding regions. The regions with different DNA methylation patterns were typically covered by small RNAs, implying that the RNA-directed DNA methylation (RdDM) pathway directed DNA methylation in the hybrids and play an important role in heterosis . The distributions of small RNAs were not explored in the present study. Nevertheless, it should speculate that the remodeling of DNA methylation directed by the non-coding RNAs should play roles for the formation of curd yield heterosis in broccoli. Further investigations are required to elucidate how these different DNA methylation sites and their relationship with gene expression function in this phenomenon.
The gene expression profiles and DNA methylation patterns of the hybrid broccoli and their parents were uncovered. The gene regulatory networks involved in signaling pathways of light, hormone and hydrogen peroxide, responses to stresses, and regulation of floral development were identified to play crucial roles in the yield heterosis of broccoli. Moreover, the present investigation indicated that the light and hydrogen peroxide-mediated signaling pathways represent two novel classes of regulatory processes that could function in yield or biomass heterosis of plants. Totally, 53 candidate genes closely involved in curd yield heterosis were identified. In addition, the data confirmed that the sites with differential methylation levels were predominant in the intergenic regions of hybrid broccoli. The changes of DNA methylation levels in gene regions did not significantly affect their expression levels. The remodeling of DNA methylation mainly directed by the non-coding RNAs then was speculated to also play roles in the formation of heterotic trait of the hybrid broccoli. These findings provided comprehensive insights into the curd yield heterosis in broccoli, and were significant for breeding high-yield broccoli varieties.
Homozygous broccoli parental lines (NKR-04, NKR-05, Bro-10 and Bro-11) and their respective F1 hybrids (NKR-06 and Bro-12) were used in this study. The F1 hybrid NKR-06 was generated by crossing NKR-04 (maternal line) with NKR-05 (paternal line). The F1 hybrid Bro-12 was generated by crossing Bro-10 (maternal line) with Bro-11 (paternal line). All of the parental materials were self-bred at least six generations to harvest homozygous seeds. All of the F1 hybrids exhibited remarkable curd yield heterosis. The seeds of F1 hybrid broccoli and their parental lines (stored in Tianjin Kernel Vegetable Research Institute, Tianjin, China) were planted in a seeding room under the controlled conditions with a 16/8 h light/dark cycle at 25 °C and 22 °C, respectively. The 25-day-old seedlings were transplanted in the soil of the greenhouse with natural light cycle at 15 °C–28 °C. At 70 d after transplanting, the leaves, curds, and roots of each material (at least 20 individual plants) were measured and then harvested, immediately frozen in liquid nitrogen, and stored at − 80 °C for subsequent use in RNA and DNA isolation.
RNA library construction and high-throughput sequencing
Total RNAs from the70-day-old curds of hybrid broccoli (NKR-06 and Bro-12) and their parents (NKR-04, NKR-05, Bro-10 and Bro-11) were isolated by TRIzol reagent (Invitrogen, CA, USA). RNAs with high purity (RNA integrity number, RIN > 8.0) were used to construct sequencing libraries by TruSeq Stranded mRNA LT Sample Prep Kit (Illumina, San Diego, CA, USA). These libraries were sequenced on the Illumina HiSeq™ 2500 sequencing platform and 150 bp/125 bp paired-end reads were generated. For each sample, at least 4 Gb raw data including over 50 million clean reads should be generated to get enough transcriptional information. In addition, to ensure the reliability of the sequencing data, three batches of independently isolated RNAs from a same sample were equally mixed to compose a RNA pool, which was used to construct sequencing library. In total, six RNA libraries were constructed and then sequenced with three technological replicates.
Transcriptome data analysis
Clean reads were generated by discarding the low-quality reads (more than 20% bases in a read which have the quality score < 15), adaptors, reads with unknown base (N) content > 5% and other contaminants from the raw reads by NGS QC Toolkit software. All clean reads were annotated and classified by mapping them to the reference genome and transcript sequences of Brassica oleracea var. oleracea (https://www.ncbi.nlm.nih.gov/genome/10901?genome_assembly_id=59537) by bowtie2 and Tophat (http://tophat.cbcb.umd.edu/) with the default value. The expressed genes were confirmed based on the annotation information of the clean reads. The specific/common expression profiles of the genes among these samples were visualized by Venn diagram (http://bioinformatics.psb.ugent.be/webtools/Venn/). The expression level of each gene was calculated and normalized by the fragments per kilobase of transcript sequence per million base pairs sequenced (FPKM). Differential expression analysis of the genes was performed by using the R package DESeq (R version 3.0.1). The significantly DEGs between the two arbitrary samples were identified based on the following thresholds: |log2 (fold-change (A/B))| > 1 and corrected p-value < 0.05. A and B represent the normalized expression level of genes in any two samples, respectively. To explore the additive and non-additive expression patterns of genes, the gene expression levels in the F1 hybrids were compared with MPV. Non-additively expressed genes were identified based on the following thresholds: |log2 (normalized expression levels of F1 hybird/MPV)| > 1 and corrected p-value < 0.05. A hierarchical cluster analysis of the DEGs was conducted by the R package “gplots” (R version 3.0.1). The Pearson correlation coefficient and hierarchical cluster of all samples were analyzed by the “cor” and “hclust” in the R project (R version 3.0.1) with the default value. In addition, GO analysis of the DEGs were performed by using the agriGO (http:// bioinfo.cau.edu.cn/agriGO/) platform, and hypergeometric test was conducted to identify the significantly enriched GO terms (corrected p-value < 0.05). To further visualize the statistically overrepresented GO terms, the GO terms were analyzed by ReviGO (http://revigo.irb.hr/).
The differential expression patterns of the genes detected by transcriptome data were validated by qRT-PCR analysis. The specific primer pairs were designed for the detection of corresponding genes (Additional file 1: Table S7). The Actin gene from broccoli was selected as internal control. Faststart Universal SYBR Green Master (Roche, Germany) was used in all experiments. The relative expression levels of the genes were calculated by the comparative 2−ΔΔCT method according to the manufacturer’s recommendations. Three batches of independently isolated RNAs from each sample were used, and three technological replicates were performed to ensure the reliability of quantitative analysis.
MethylRAD library construction and high-throughput sequencing
The 70-day-old curds of hybrid broccoli and their parents were used to isolate the genomic DNAs by cetyltrimethyl ammonium bromide (CTAB) method, respectively. 150–200 ng of DNAs were used to construct the MethylRAD library with slight modifications . Briefly, the total DNAs were digested by using the methylation-dependent restriction enzyme FspEI (NEB, USA). FspEI can recognize 5-methylcytosine (5-mC) and 5-hydroxymethylcytosine (5-hmC) in the CmCGG and mCHG sites, and generate a double-stranded DNA break on the 30 side of the modified cytosine at a fixed distance (N12/N16). Accordingly, symmetrical DNA methylated sites were bidirectionally cleaved by FspEI to generate 32-base-long fragments. Then, two adaptors were added to the digested DNAs by T4 DNA ligase (NEB, USA), and the ligation products were amplified by specific primers. The target DNA fragments were further purified and subjected to sequencing by Illumina HiSeq X-Ten sequencing platform . Three batches of independently isolated DNAs from a same sample were equally mixed to compose a DNA pool, which was used to construct DNA methylation sequencing library. In total, six libraries were constructed and then sequenced with three technological replicates.
DNA methylation data analysis
Adaptors and low-quality reads were discarded by NGS QC Toolkit software. The clean reads that did not contain the expected FspEI restriction site were further excluded, and the reads containing the methylated CG or CHG sites, named MethylRAD-tags, were identified. The MethylRAD-tags were mapped to the reference genome of Brassica oleracea var. oleracea by SOAP program (version 2.21) with two mismatches allowed. The distribution and density of methylated cytosine sites on chromosomes were calculated and visualized by the Integrative Genomics Viewer (IGV) . Furthermore, the distributions of the methylated cytosine sites on different elements of the genome, especially on the different regions of genes, were evaluated. The DNA methylation levels of the sites were determined by the FPKM value, and the DNA methylation levels of the genes were then evaluated by summing the methylation levels of sites that were localized in the gene regions. The differential DNA methylation levels of sites and genes were identified by using the R package DESeq (R version 3.0.1). The thresholds were |log2 (fold-change (C/D)) > 1| and corrected p-value < 0.05. C and D represented the FPKM of sites or genes in any two samples, respectively.
- ABCG26 :
ABC transporter G family member 26
Cetyltrimethyl ammonium bromide
Differentially expressed genes
Fragments per kilobase of transcript sequence per million base pairs
Quantitative real-time RT-PCR
Quantitative trait locus
Transposable element silencing via AT-hook
Transcription start sites
Transcription termination site
Groszmann M, Greaves IK, Fujimoto R, Peacock WJ, Dennis ES. The role of epigenetics in hybrid vigour. Trends Genet. 2013;29(12):684–90.
Schnable PS, Springer NM. Progress toward understanding heterosis in crop plants. Annu Rev Plant Biol. 2013;64:71–88.
Duvick DN. Heterosis: feeding people and protecting natural resources. In: Coors JG, Pandey S, editors. Genetics and exploitation of heterosis in crops. Madison: American Society of Agronomy; 1997. p. 19–29.
Darwin CR. The effects of cross and self-fertilization in the vegetable kingdom. London: John Murray; 1876.
Birchler JA, Yao H, Chudalayandi S. Unraveling the genetic basis of hybrid vigor. Proc Natl Acad Sci U S A. 2006;103(35):12957–8.
Li ZK, Luo LJ, Mei HW, Wang DL, Shu QY, Tabien R, Zhong DB, Ying CS, Stansel JW, Khush GS, Paterson AH. Overdominant epistatic loci are the primary genetic basis of inbreeding depression and heterosis in rice. I. Biomass and grain yield. Genetics. 2001;158(4):1737–53.
Birchler JA, Auger DL, Riddle NC. In search of the molecular basis of heterosis. Plant Cell. 2003;15(10):2236–9.
Hochholdinger F, Hoecker N. Towards the molecular basis of heterosis. Trends Plant Sci. 2007;12(9):427–32.
Gu L, Wu Y, Jiang M, Si W, Zhang X, Tian D, Yang S. Dissimilar manifestation of heterosis in superhybrid rice at early-tillering stage under nutrient-deficient and nutrient-sufficient condition. Plant Physiol. 2016;172(2):1142–53.
Samayoa LF, Malvar RA, Butrón A. QTL for maize midparent heterosis in the heterotic pattern American dent × European flint under corn borer pressure. Front Plant Sci. 2017;8:573.
Garcia AA, Wang S, Melchinger AE, Zeng ZB. Quantitative trait loci mapping and the genetic basis of heterosis in maize and rice. Genetics. 2008;180(3):1707–24.
Wei G, Tao Y, Liu G, Chen C, Luo R, Xia H, Gan Q, Zeng H, Lu Z, Han Y, Li X, Song G, Zhai H, Peng Y, Li D, Xu H, Wei X, Cao M, Deng H, Xin Y, Fu X, Yuan L, Yu J, Zhu Z, Zhu L. A transcriptomic analysis of superhybrid rice LYP9 and its parents. Proc Natl Acad Sci U S A. 2009;106(19):7695–701.
Shen H, He H, Li J, Chen W, Wang X, Guo L, Peng Z, He G, Zhong S, Qi Y, Terzaghi W, Deng XW. Genome-wide analysis of DNA methylation and gene expression changes in two Arabidopsis ecotypes and their reciprocal hybrids. Plant Cell. 2012;24(3):875–92.
Groszmann M, Greaves IK, Albertyn ZI, Scofield GN, Peacock WJ, Dennis ES. Changes in 24-nt siRNA levels in Arabidopsis hybrids suggest an epigenetic contribution to hybrid vigor. Proc Natl Acad Sci U S A. 2011;108(6):2617–22.
He G, Zhu X, Elling AA, Chen L, Wang X, Guo L, Liang M, He H, Zhang H, Chen F, Qi Y, Chen R, Deng XW. Global epigenetic and transcriptional trends among two rice subspecies and their reciprocal hybrids. Plant Cell. 2010;22(1):17–33.
Gao M, Huang Q, Chu Y, Ding C, Zhang B, Su X. Analysis of the leaf methylomes of parents and their hybrids provides new insight into hybrid vigor in Populus deltoides. BMC Genet. 2014;15:S8.
Zhang Q, Wang D, Lang Z, He L, Yang L, Zeng L, Li Y, Zhao C, Huang H, Zhang H, Zhang H, Zhu JK. Methylation interactions in Arabidopsis hybrids require RNA-directed DNA methylation and are influenced by genetic variation. Proc Natl Acad Sci U S A. 2016;113(29):E4248–56.
Greaves IK, Eichten SR, Groszmann M, Wang A, Ying H, Peacock WJ, Dennis ES. Twenty-four-nucleotide siRNAs produce heritable trans-chromosomal methylation in F1 Arabidopsis hybrids. Proc Natl Acad Sci U S A. 2016;113(44):E6895–902.
Zhang L, Peng Y, Wei X, Dai Y, Yuan D, Lu Y, Pan Y, Zhu Z. Small RNAs as important regulators for the hybrid vigour of super-hybrid rice. J Exp Bot. 2014;65(20):5989–6002.
Chen Y, Zhou Q, Tian R, Ma Z, Zhao X, Tang J, Fu Z. Proteomic analysis reveals that auxin homeostasis influences the eighth internode length heterosis in maize (Zea mays). Sci Rep. 2018;8(1):7159.
Hu X, Wang H, Diao X, Liu Z, Li K, Wu Y, Liang Q, Wang H, Huang C. Transcriptome profiling and comparison of maize ear heterosis during the spikelet and floret differentiation stages. BMC Genomics. 2016;17(1):959.
Song GS, Zhai HL, Peng YG, Zhang L, Wei G, Chen XY, Xiao YG, Wang L, Chen YJ, Wu B, Chen B, Zhang Y, Chen H, Feng XJ, Gong WK, Liu Y, Yin ZJ, Wang F, Liu GZ, Xu HL, Wei XL, Zhao XL, Ouwerkerk PB, Hankemeier T, Reijmers T, van der Heijden R, Lu CM, Wang M, van der Greef J, Zhu Z. Comparative transcriptional profiling and preliminary study on heterosis mechanism of super-hybrid rice. Mol Plant. 2010;3(6):1012–25.
Hu S, Wang C, Sanchez DL, Lipka AE, Liu P, Yin Y, Blanco M, Lübberstedt T. Gibberellins promote brassinosteroids action and both increase heterosis for plant height in maize (Zea mays L.). Front Plant Sci. 2017;8:1039.
Wang L, Wu LM, Greaves IK, Zhu A, Dennis ES, Peacock WJ. PIF4-controlled auxin pathway contributes to hybrid vigor in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2017;114(17):E3555–62.
Groszmann M, Gonzalez-Bayon R, Lyons RL, Greaves IK, Kazan K, Peacock WJ, Dennis ES. Hormone-regulated defense and stress response networks contribute to heterosis in Arabidopsis F1 hybrids. Proc Natl Acad Sci U S A. 2015;112(46):E6397–406.
Li Y, Wang C, Liu X, Song J, Li H, Sui Z, Zhang M, Fang S, Chu J, Xin M, Xie C, Zhang Y, Sun Q, Ni Z. Up-regulating the abscisic acid inactivation gene ZmABA8ox1b contributes to seed germination heterosis by promoting cell expansion. J Exp Bot. 2016;67(9):2889–900.
Zhang Q, Li Y, Xu T, Srivastava AK, Wang D, Zeng L, Yang L, He L, Zhang H, Zheng Z, Yang DL, Zhao C, Dong J, Gong Z, Liu R, Zhu JK. The chromatin remodeler DDM1 promotes hybrid vigor by regulating salicylic acid metabolism. Cell Discov. 2016;2:16027.
Wang L, Greaves IK, Groszmann M, Wu LM, Dennis ES, Peacock WJ. Hybrid mimics and hybrid vigor in Arabidopsis. Proc Natl Acad Sci U S A. 2015;112(35):E4959–67.
Miller M, Song Q, Shi X, Juenger TE, Chen ZJ. Natural variation in timing of stress-responsive gene expression predicts heterosis in intraspecific hybrids of Arabidopsis. Nat Commun. 2015;6:7453.
Castillon A, Shen H, Huq E. Phytochrome interacting factors: central players in phytochrome-mediated light signaling networks. Trends Plant Sci. 2007;12(11):514–21.
Foreman J, Johansson H, Hornitschek P, Josse EM, Fankhauser C, Halliday KJ. Light receptor action is critical for maintaining plant biomass at warm ambient temperatures. Plant J. 2011;65(3):441–52.
de Wit M, Keuskamp DH, Bongers FJ, Hornitschek P, Gommers CMM, Reinen E, Martínez-Cerón C, Fankhauser C, Pierik R. Integration of phytochrome and cryptochrome signals determines plant growth during competition for light. Curr Biol. 2016;26(24):3320–6.
Saxena I, Srikanth S, Chen Z. Cross talk between H2O2 and interacting signal molecules under plant stress response. Front Plant Sci. 2016;7:570.
Meinhard M, Rodriguez PL, Gril E. The sensitivity of ABI2 to hydrogen peroxide links the abscisic acid-response regulator to redox signalling. Planta. 2002;214(5):775–82.
Yuan HM, Liu WC, Lu YT. CATALASE2 coordinates SA-mediated repression of both auxin accumulation and JA biosynthesis in plant defenses. Cell Host Microbe. 2017;21(2):143–55.
Tian Y, Fan M, Qin Z, Lv H, Wang M, Zhang Z, Zho W, Zhao N, Li X, Han C, Ding Z, Wang W, Wang ZY, Bai MY. Hydrogen peroxide positively regulates brassinosteroid signaling through oxidation of the BRASSINAZOLE-RESISTANT1 transcription factor. Nat Commun. 2018;9(1):1063.
Xiong SX, Lu JY, Lou Y, Teng XD, Gu JN, Zhang C, Shi QS, Yang ZN, Zhu J. The transcription factors MS188 and AMS form a complex to activate the expression of CYP703A2 for sporopollenin biosynthesis in Arabidopsis thaliana. Plant J. 2016;88(6):936–46.
Yang X, Wu D, Shi J, He Y, Pinot F, Grausem B, Yin C, Zhu L, Chen M, Luo Z, Liang W, Zhang D. Rice CYP703A3, a cytochrome P450 hydroxylase, is essential for development of anther cuticle and pollen exine. J Integr Plant Biol. 2014;56(10):979–94.
Singh M, Kumar M, Thilges K, Cho MJ, Cigan AM. MS26/CYP704B is required for anther and pollen wall development in bread wheat (Triticum aestivum L.) and combining mutations in all three homeologs causes male sterility. PLoS One. 2017;12(5):e0177632.
de Azevedo SC, Kim SS, Koch S, Kienow L, Schneider K, McKim SM, Haughn GW, Kombrink E, Douglas CJ. A novel fatty acyl-CoA Synthetase is required for pollen development and sporopollenin biosynthesis in Arabidopsis. Plant Cell. 2009;21(2):507–25.
Kim SS, Grienenberger E, Lallemand B, Colpitts CC, Kim SY, Souza Cde A, Geoffroy P, Heintz D, Krahn D, Kaiser M, Kombrink E, Heitz T, Suh DY, Legrand M, Douglas CJ. LAP6/POLYKETIDE SYNTHASE A and LAP5/POLYKETIDE SYNTHASE B encode hydroxyalkyl α-pyrone synthases required for pollen development and sporopollenin biosynthesis in Arabidopsis thaliana. Plant Cell. 2010;22(12):4045–66.
Zhu X, Yu J, Shi J, Tohge T, Fernie AR, Meir S, Aharoni A, Xu D, Zhang D, Liang W. The polyketide synthase OsPKS2 is essential for pollen exine and Ubisch body patterning in rice. J Integr Plant Biol. 2017;59(9):612–28.
Grienenberger E, Kim SS, Lallemand B, Geoffroy P, Heintz D, Souza Cde A, Heitz T, Douglas CJ, Legrand M. Analysis of TETRAKETIDE α-PYRONE REDUCTASE function in Arabidopsis thaliana reveals a previously unknown, but conserved, biochemical pathway in sporopollenin monomer biosynthesis. Plant Cell. 2010;22(12):4067–83.
Choi H, Jin JY, Choi S, Hwang JU, Kim YY, Suh MC, Lee Y. An ABCG/WBC-type ABC transporter is essential for transport of sporopollenin precursors for exine formation in developing pollen. Plant J. 2011;65(2):181–93.
Quilichini TD, Samuels AL, Douglas CJ. ABCG26-mediated polyketide trafficking and hydroxycinnamoyl spermidines contribute to pollen wall exine formation in Arabidopsis. Plant Cell. 2014;26(11):4483–98.
Lou Y, Xu XF, Zhu J, Gu JN, Blackmore S, Yang ZN. The tapetal AHL family protein TEK determines nexine formation in the pollen wall. Nat Commun. 2014;5:3855.
Ng DW, Miller M, Yu HH, Huang TY, Kim ED, Lu J, Xie Q, McClung CR, Chen ZJ. A role for CHH methylation in the parent-of-origin effect on altered circadian rhythms and biomass heterosis in Arabidopsis intraspecific hybrids. Plant Cell. 2014;26(6):2430–40.
Kawanabe T, Ishikura S, Miyaji N, Sasaki T, Wu LM, Itabashi E, Takada S, Shimizu M, Takasaki-Yasuda T, Osabe K, Peacock WJ, Dennis ES, Fujimoto R. Role of DNA methylation in hybrid vigor in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2016;113(43):E6704–11.
Shen G, Hu W, Zhang B, Xing Y. The regulatory network mediated by circadian clock genes is related to heterosis in rice. J Integr Plant Biol. 2015;57(3):300–12.
Wang S, Lv J, Zhang L, Dou J, Sun Y, Li X, Fu X, Dou H, Mao J, Hu X, Bao Z. MethylRAD: a simple and scalable method for genome-wide DNA methylation profiling using methylation-dependent restriction enzymes. Open Biol. 2015;5(11).
Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP. Integrative genomics viewer. Nat Biotechnol. 2011;29:24–6.
We thank Dr. Chengguo Duan of Shanghai Center for Plant Stress Biology and Center of Excellence in Molecular Plant Sciences, Chinese Academy of Sciences, Shanghai, China, for his valuable suggestions in data analysis.
This work was funded by grants from the Natural Science Foundation of China (No. 31470669 and No. 31401889) and Tianjin (No.14JCZDJC34000 and No. 15JCQNJC15100), and Graduate education quality improvement projects of Tianjin Agricultural University (2017YPY004). The funders had no role in study design, data collection, data analysis, data interpretation, the writing of the manuscript or decision to publish.
Availability of data and materials
The data sets supporting the results of the present study are included within this article (and its additional files).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. View of the transcriptome data in two broccoli hybrid combinations. Table S2. Normalized expression levels and functional annotation of genes detected in two broccoli hybrid combinations. Table S3. DNA methylation ratio in two hybrid broccoli and their parents. Table S4. Relative DNA methylation levels at different regions of genes. Table S5. Distributions of CG and CHG sites with differential methylation levels. Table S6. Genes simultaneously showing differential methylation levels and differential expression levels. Table S7. primers used in the study. (ZIP 1392 kb)
Data set 1. DEGs detected by pairwise comparison in NKR-06 and two parents. Data set 2. DEGs detected by pairwise comparison in Bro-12 and two parents. Data set 3. DEGs in the hybrids compared with both parents. Data set 4. Non-additively expressed genes in two hybrid triads. Data set 5. GO functional annotation of the DEGs in the NKR-04/− 05/− 06 hybrid triad. Data set 6. GO functional annotation of the DEGs in the Bro-10/− 11/− 12 hybrid triad. Data set 7. DNA methylation of CG and CHG sites detected in the hybrid broccoli and their parents. Data set 8. Differential DNA methylation levels of CG and CHG sites in NKR-06 hybrid triad. Data set 9. Genes with different methylation levels in NKR-06 hybrid triad. Data set 10 Differential DNA methylation levels of CG and CHG sites in Bro-12 hybrid triad. Data set 11. Genes with differential methylation levels in Bro-12 hybrid triad. (ZIP 5944 kb)
Figure S1. Expression levels of genes involved in several overrepresented biological processes. Figure S2. Expression profiles of several genes in BNR-H broccoli hybrid triad detected by qRT-PCR. Figure S3. Distributions of CG and CHG methylation sites at the different regions of genomes in the hybrids and their parents. Figure S4. DNA methylation levels at CHG sites in different regions of genes. (ZIP 34711 kb)