- Research article
- Open Access
The phenotypic predisposition of the parent in F1 hybrid is correlated with transcriptome preference of the positive general combining ability parent
BMC Genomicsvolume 15, Article number: 297 (2014)
Sprague and Tatum (1942) introduced the concepts of general combining ability (GCA) and specific combining ability (SCA) to evaluate the breeding parents and F1 hybrid performance, respectively. Since then, the GCA was widely used in cross breeding for elite parent selection. However, the molecular basis of GCA remains to unknown.
We studied the transcriptomes of three varieties and three F1 hybrids using RNA-Sequencing. Transcriptome sequence analysis revealed that the transcriptome profiles of the F1s were similar to the positive GCA-effect parent. Moreover, the expression levels of most differentially expressed genes (DEGs) were equal to the parent with a positive GCA effect. Analysis of the gene expression patterns of gibberellic acid (GA) and flowering time pathways that determine plant height and flowering time in rice validated the preferential transcriptome expression of the parents with positive GCA effect. Furthermore, H3K36me3 modification bias in the Pseudo-Response Regulators (PRR) gene family was observed in the positive GCA effect parents and demonstrated that the phenotype and transcriptome bias in the positive GCA effect parents have been epigenetically regulated by either global modification or specific signaling pathways in rice.
The results revealed that the transcriptome profiles and DEGs in the F1s were highly related to phenotype bias to the positive GCA-effect parent. The transcriptome bias toward high GCA parents in F1 hybrids attributed to H3K36me3 modification both on global modification level and specific signaling pathways. Our results indicated the transcriptome profile and epigenetic modification level bias to high GCA parents could be the molecular basis of GCA.
Selecting elite parents is of paramount importance in cross- and hybrid-breeding programs. The selection of parents from a phalanx of inbred lines, however, is extremely laborious and time-consuming and can be random. Adding to this complexity, parents with excellent agronomic traits do not always pass those traits on to their progeny. To evaluate breeding parents, Sprague and Tatum (1942) introduced the concepts of general combining ability (GCA) and specific combining ability (SCA), which allow the study and comparison of the performances of inbred lines in hybrid combination. GCA is used to designate the average performance of an inbred line in hybrid combination, and SCA is used to designate those cases in which certain combinations do relatively better or worse than would be expected on the basis of the average performance of the lines involved . GCA and SCA effects have been successfully used as criteria to evaluate elite parents in conventional crossbreeding and the performance of hybrid combinations.
Since the introduction of the GCA concept in 1942, very limited genetic studies on GCA have been conducted even though GCA is widely used in breeding programs for evaluation of the parents in early generations [2–5]. A few quantitative genetic loci with GCA have been identified recently. Qu et al. analyzed the QTLs of 10 agronomic traits for GCA using recombinant inbred line (RIL) populations with three testers in three testcross populations and a backcross recombinant inbred line (BCRIL) population of rice . They detected a large number of additive effects of QTLGCA loci. Qi et al. found that several genetic loci responding for GCA and SCA for five yield-related traits using a set of testcrosses with introgression lines (ILs) of maize under different environmental conditions. Total of 56 significant QTLGCA loci have been mapped . These studies have revealed that GCA effects, like the traits, are genetically controlled. At the molecular level, however, how the phenotypes associated with GCA effects are passed on to the F1 hybrids remains unknown.
High-throughput genome-wide analysis approaches such as microarray analysis and next generation sequencing have been used in the study of phenotypes such as hybrid vigor [8–10]. Genome-wide gene expression profiles related to heterosis have been studied [11–14]. Stupar et al. studied the gene expression profiles between maize F1 hybrids and their parents, and approximately 75% of the differentially expressed genes showed additive expression patterns between F1 hybrids and parents . They studied the genetic diversity and transcriptional variation with different maize hybrids and found that the genetic diversity was correlated with transcriptional variation, and little (less than 1%) of the gene expression in F1 hybrids was outside the parental range . Wei et al. studied the gene expression profile between the super hybrid rice LYP9 and its parents and suggested that the differentially expressed genes might correlate with heterosis. Moreover, Riedelsheimer et al. have analyzed metabolic pathways in maize hybrids and were able to predict GCA scores using genome-wide association analysis . They have used 285 crosses derived from the diverse inbred lines of maize with two testers and predicted their combining ability for seven biomass- and bioenergy related traits using 56,110 single nucleotide polymorphisms (SNPs) and 130 metabolites. The prediction accuracies were from 0.72 to 0.81 for SNPs and 0.60 to 0.80 for metabolites. A genome-wide analysis of the association between GCA and gene expression profiles has not been investigated. The contribution of the elite parents with GCA effects to their offspring also remains unexplored.
In the present study, we constructed a diallele crossing population with five rice parents and 10 F1 hybrids. We evaluated the GCA effects of three agronomic traits related to grain yield, heading date, plant height and grain number from five elite rice varieties. The results showed that 93–11 and Teqing (TQ) have positive GCA effects, and Guangluai 4# (GL), Aijiaonante (AJ) and Zhenshan 97 (ZS) showed negative GCA effects. The analysis of the transcriptome profiles of the leaves from three F1 hybrids, GL × 93-11, GL × TQ and 93-11 × TQ, and their parents revealed that transcriptome profiles were correlated to the positive-GCA-effects parent, showing obvious parental bias. Further analysis found that the expression levels of most of the DEGs were obviously biased towards to the positive GCA effect parent. Analysis of the gene expression patterns of gibberellic acid (GA) and flowering-time signaling pathways for plant height and flowering time validated the transcriptome bias to the positive GCA effect parents. Our results also indicated that the H3K36me3 modification showed bias to the positive GCA effect parent and demonstrated that phenotypes of the positive GCA effect parents were attributed to transcriptome bias. These results will be helpful in understanding the parental gene contribution to F1 hybrids and the molecular basis of GCA.
Evaluation of general combining ability of elite rice varieties
Five inbred rice varieties representing different breeding objectives from the 1970s to the present were chosen. Ten F1 hybrids were obtained from five varieties using a diallele crossing design (Additional file 1: Figure S1 and Additional file 2: Figure S2). Three important agronomic traits related to grain yield (that is, plant height, heading date and grains per panicle) were evaluated based on the phenotypes of F1 hybrids from the diallele crossing populations (Figure 1A, B and C). The results revealed that the GCA scores for the two rice varieties 93–11 and TQ were positive, whereas negative GCA scores were observed for the varieties GL, ZS and AJ (Figure 1D). Further analysis revealed that the phenotypes of the three agronomic traits in the F1 hybrids were always similar to or higher than 93–11 and TQ, when either was used as one parent in the cross (Figure 1A-C and Additional file 1: Figure S1). In contrast, the phenotypes of F1 hybrids derived from GL (Figure 1), AJ and ZS (Additional file 1: Figure S1 and Additional file 2: Figure S2) were significantly different compared to 93–11 and TQ. The parental phenotypes were significantly correlated with the F1 hybrids crossed to 93–11 and TQ (r > 0.97, P < 3.80E-07) but less significantly with the F1 hybrids resulting from crosses to GL, ZS and AJ (r = 0.81, P < 1.33E-03) (Table 1). These results indicated that the GCAs for the three agronomic traits relevant to grain yield in parents 93–11 and TQ were significantly higher than those of the varieties GL, AJ and ZS.
The phenotype bias toward the positive GCA effect parent is attributable to transcriptome bias
As shown in Figure 1 and Additional file 3: Figure S3, GL had an earlier heading date (68.4 days on average) and had a shorter plant height (88.9 cm on average), whereas TQ and 93–11 had later heading dates (90.1 and 93.7 days on average) and had a taller plant height (113.3 and 117.2 cm on average, respectively). The heading dates (90.5 and 109.5 days on average, respectively) and plant heights (113.9 and 117.0 cm on average, respectively) in the F1 hybrids derived from TQ and 93–11 were obviously biased toward TQ or 93–11. To understand the relationship between the F1 phenotype and the positive GCA effect parent, we conducted a transcriptome profile analysis of the three parents, GL, 93–11 and TQ, and their three F1 hybrids, GL × 93-11, GL × TQ and 93-11 × TQ, using RNA sequencing technology (Additional file 3: Figure S3) . A cluster analysis indicated that the transcriptomes of GL × TQ and GL × 93-11 are similar to those of the positive GCA effect parents, TQ and 93–11, respectively, and are significantly different from those of the negative GCA effect parent GL (Figure 2A and 2B). The transcriptome profile of 93-11 × TQ was more similar to 93–11 that had the higher positive GCA effect than to TQ (Figure 2C). The transcriptome similarity of the F1 hybrids to either the 93–11 or TQ parent was consistent with the phenotypes of the three traits in 93–11 and TQ. To explore the transcriptome profiles of the parents in the F1 hybrids, we analyzed the gene expression level in F1 hybrids. The results showed that 76.0% of the gene expression profiles in GL × TQ were similar to GL and 83.7% were similar to TQ (Figure 2D). Analogous results were found for the GL × 93-11 and 93-11 × TQ (i.e., 84% of the gene expression profiles in GL × 93-11 were similar to those in 93–11 and 67.1% of genes bias to GL) (Figure 2E). Of 90.2% expressed genes in 93-11 × TQ were similar to those in 93–11 and 82.6% of the genes bias to TQ (Figure 2F). Therefore, the phenotypes of the F1s derived from the positive GCA effect parents are correlated with a transcriptome bias in the F1s.
Differentially expressed genes in the F1 hybrids are those preferentially expressed in the positive GCA effect parent
Differentially expressed genes (DEGs) have been recognized to play important biological functions in heterosis [16–18]. To explore whether DEGs contribute to the positive GCA effect parent in the F1 hybrids, we analyzed the DEG profiles from three F1 hybrids. A total of 34,486-35,718 genes were expressed in the F1 hybrids and their parents (Figure 3A). We found that 22.1% to 37.0% genes were differentially expressed in the F1 hybrids comparing GL × TQ, GL × 93-11 and 93-11 × TQ with their parents, respectively (fold change >2.0, P < 0.05) Moreover, 75.5% to 84.1% of the DEGs between the F1 hybrids and the parents had similar expression levels to one parent (Figure 3B). Of these, 63.4%, 79.5% and 69.6% of the gene expression levels were similar to the positive GCA effect parents TQ and 93–11 in GL × TQ, GL × 93-11 and 93-11 × TQ, respectively (Figure 3C). The expression levels of the remaining genes were similar to the negative or lower GCA effect parents GL and TQ (Figure 3C). 7.9% to 19.7% of the DEGs expression level in the F1 hybrids were out of the parental ranges, whereas 4.8% to 10.2% of the DEGs showed mid-parent expression levels, but showing significant difference to both parents (Figure 3B). These results reveal that the expression levels of the majority of the DEGs were similar to those of the positive GCA parent.
Expression patterns of the genes in the regulatory pathways controlling flowering time and plant height in rice
To further support the observed correlation between the positive GCA effect phenotype in the parent and transcriptome bias in the F1, we chose the genes in two well-known signaling pathways controlling flowering time (Additional file 4: Table S1) and plant height in rice. The expression of EARLY HEADING DATE1 (Ehd1) and EARLY HEADING DATE2 (Ehd2) and repression of HEADING DATE1 (Hd1), HEADING DATE3a (Hd3a) and RICE FLOWERING LOCUS T1 (RFT1) are activated in a cascade to initiate rice flowering [19–25]. Expression of the Rice Pseudo-Response Regulators 1 (OsPRR1) and Grain number, plant height and heading date7 (Ghd7) involves a delay in the flowering time and an increase in the plant height and grain number [26–28]. We found that the expression levels of OsPRR1, Hd1 and Ghd7 were very high in TQ, 93–11, GL × TQ, GL × 93-11 and 93-11 × TQ versus GL (Table 2). By contrast, the expression levels of RFT1, Hd3a, Ehd1 and Ehd2 were higher in GL (Table 2). These results showed that the expression patterns of the flowering regulation genes were consistent with transcriptome bias towards the positive GCA effect parents in the F1 hybrids.
The expression patterns of the other set of genes involved in gibberellic acid (GA) metabolism and the signaling pathway that controls plant height in rice were also analyzed between F1s and parents [29–32] (Additional file 5: Table S2). In plants, bioactive GAs are synthesized from the precursor geranylgeranyl diphosphate by ent-copalyl diphosphate synthase (CPS) and ent-kaurene synthase (KS), followed by ent-kaurene oxidase (KO) and ent-kaurenoic acid oxidase (KAO). At the final stage, GA20-oxidase (GA20ox) and GA3-oxidase (GA3ox) catalyze the conversion of GA53/GA12 and GA1/GA4. The bioactive GAs and the precursors elongated the uppermost internode and were deactivated by GA2-oxidase (GA2ox) and (EUI) . As shown in Table 2, three genes, ent-copalyl diphosphate synthase (OsCPS1), ent-kaurenoic acid oxidase (OsKAO) and GA2-oxidase 6 (OsGA2ox6) had the higher expression levels in TQ, 93–11, GL × TQ, GL × 93-11 than in GL. In contrast, expression levels of the other three genes, GA3ox2 and GA20ox2 were less in TQ, 93–11, GL × TQ, GL × 93-11 and 93-11 × TQ than in GL. The expression profiles of those genes exactly matched the feedback and feed-forward regulation mechanism of GAs synthesis [31, 32]. Our results are consistence with the results reported in wheat hybrids [33, 34]. The results again demonstrated that the expression patterns of the genes corresponding to plant height exhibited a bias toward the positive GCA effect parents. Taken together, the data validated the transcriptome bias toward the positive GCA effect parents through individual metabolism pathways for plant height and heading date in rice and further demonstrated that the phenotype bias to the positive GCA effect parent in F1 hybrids is due to the transcriptome bias toward the positive GCA effect parents.
Transcriptome bias toward the positive GCA effect parent regulated by H3K36me3 modifications
Previous studies have described the histone modifications involved in the regulation of the transcriptome [35–38]. Previous studies have indicated that the trimethylation of histone H3 on lysine 27 (H3K27me3) in gene body represses gene expression [10, 39, 40] and trimethylated histone H3 on lysine 36 (H3K36me3) in the gene body activated gene expression [41–43]. To further explore the mechanism of transcriptome bias to the positive GCA effect parent in F1 hybrids, we analyzed the patterns of methylation at H3K36me3 and H3K27me3 in F1 hybrids. Genome-wide histone modifications of H3K36me3 and H3K27me3 in GL × 93-11, GL × TQ and their parents, GL, 93–11 and TQ, were analyzed with the aid of a DNA library prepared after chromatin immunoprecipitation (ChIP). The results showed 18.2%-18.3% genes were overlapped between H3K36me3 and H3K27me3 modifications. However, an obvious parental bias of H3K36me3 modification was detected. In GL × 93-11, 82.0% of the genes with H3K36me3 modifications were similar to 93–11 (FDR > 0.001), whereas 73.4% were similar to GL (FDR > 0.001) (Figure 4A). 93.0% of H3K36me3 modifications in GL × TQ were found to prefer TQ (FDR > 0.001), and 75.0% of H3K36me3 preferred GL (FDR > 0.001). For H3K27me3 modifications, 69.8% of the genes in GL × 93-11 with modifications were similar to GL (FDR > 0.001), and 65.2% were biased to 93–11 (FDR > 0.001) (Figure 4B). In GL × TQ, 90.7% genes with the H3K27me3 modifications preferred TQ (FDR > 0.001), whereas only 62.5% of H3K27me3 modifications were biased to GL (FDR > 0.001) (Figure 4B). There was no significant correlation of H3K27me3 between positive GCA effect and negative GCA effect parents in two F1 hybrids. Taken together, the results indicated that H3K36me3 modifications exhibited a bias toward the positive GCA effect parent in the F1 hybrids, which suggested that epigenetic regulation could involve regulation of the phenotypes and transcriptome bias toward positive GCA effect parents.
Furthermore, we analyzed the modification patterns of Cryptochrome 2 (CRY2), Pseudo-Response Regulator 37 and 95 (PRR37 and PRR95) in F1 hybrids (Figure 5). These genes are involved in the regulation of rice flowering time and adaptability. CRY2 is a photolyase-like blue-light receptor that mediates light responses in plants via interaction with the CIB1 (cryptochrome-interacting basic-helix-loop-helix) protein to promote CRY2-dependent floral initiation . The PRR37 gene is involved in the down-regulation of Hd3a gene expression to suppress flowering under long-day conditions. The varieties harboring nonfunctional alleles of PRR37 flower extremely early under natural, long-day conditions . We found that the H3K36me3 modification levels of three genes in GL × 93-11 or GL × TQ biased toward the positive GCA effect parents, TQ or 93–11 (Figure 5). The results are consistent with the heading date phenotype bias toward high GCA effect parents in the F1. Again, our data demonstrated that the phenotype and transcriptome biased toward the positive GCA effect parents were epigenetically regulated by both global modification and specific signaling pathways in rice.
The GCA concept has been studied for more than 70 years . It has being successfully applied in crop and livestock breeding for the evaluation of parent performance [3–5]. However, in spite of this, the genetic analysis of the GCA is very limited [6, 7], and the molecular mechanism of the GCA has not been documented. In the present study, we found that certain phenotypes in rice F1 hybrids were always biased toward the positive GCA effect parents in the traits tested (r > 0.97, P < 3.80E-07). Further studies showed that the transcriptome profiles in F1 hybrids were the same as phenotypes that are biased to the parents with positive GCA effect. The expression patterns of individual genes corresponding to plant height and flowering time demonstrated that the phenotype bias toward the positive GCA effect parent is attributed to the transcriptome profile and the specific epigenetic modification. Our data revealed that the phenotype bias in F1 hybrids could be attributed to transcriptome and H3K36me3 modification bias. Our findings provide molecular clues as to why and how the positive GCA effect parent has been widely used for crop and livestock improvement in breeding programs for half a century.
Although GCA is widely recognized and applied to crop and livestock improvement by breeders, it is based largely on breeding experiences, and therefore, it highly is unpredictable. Our findings provide molecular evidence of GCA through transcriptome and epi-genome analysis. This study first revealed that the transcriptome and epigenome in rice F1 hybrids are similar to the parent with the positive GCA performance, but significantly different from the negative GCA effect parents. We also found that the positive GCA effect parents showed better performance with regard to agronomic traits in F1 hybrids compared to the negative GCA effect parents. In the previous studies, a large number of genetic loci with dominance effect were detected in F1 hybrids [46–48]. eQTL analysis revealed that the gene expression regulation was an complicated regulation networks [49–51]. So, the phenomena of the phenotypes and the transcriptomes biased to high GCA parents in F1 hybrids might be consequence of the accumulation of the loci with dominance effect in the elite parents through artificial selection. Our results suggested that the transcriptomes profile in F1 hybrids could be optimized during elite parent selection, which give rise to elite agronomic traits, such as increase of biomass and grain yield. During the parent improvement process, these alleles corresponding to the agronomic traits satisfying different breeding objectives had been selected or kept. For instance, a representative elite variety of GL that met the breeding objectives of more panicle and less grain number per panicle had been developed in the 1970s in China. Later, an ideal plant type that has less panicle numbers and more spikelet per panicle was raised, and the breeding objectives of higher grain yields required high biomass that was due to later heading date and higher plant height. To meet the requirements of the breeding objective, the alleles corresponding to late heading date and higher plant height were accumulated and optimized through artificial selection pressure. Therefore, the modern elite parent 93–11 has a longer flowering time and taller plant height that produces higher biomass resulting in higher grain yield. Consequently, the alleles corresponding to later flowering time and taller plant height were selected and maintained in the genome of the elite parent. Our results again demonstrated the significance of selecting elite parents with excellent GCA performance in crossbreeding program.
Several studies have indicated that the mechanism of the allelic-specific expression could be regulated by epigenetic modification, including DNA methylation and various histone modifications [52, 53]. Previous studies have detected maternal alleles for some imprinted genes were hypomethylated in the endosperm, which contributed to the parent-of-origin expression pattern in reciprocal hybrids in maize[52, 53]. However, very few parent-of-origin effects were detected in plant vegetative tissues . Our global gene expression analysis has not found any parent-of-origin genes, but almost of genes expression in rice leaves showed as genotype-dependent fashion . In present work, the specific genes involved in plant height and circadian rhythms are matched to phenotype bias as well, suggesting that epigenetic regulation is highly selective and complicated. Limited data obtained in the present work may have led to finding the H3K36me3 modification only on the partial genes that were associated with the GCA effect, and these modifications did not match the expression level of the genes analyzed. The results of H3K27me3 modification profile biased to TQ in GL × TQ, but not bias in another F1 hybrids GL × 93-11 suggested that epigenetic modifications could be dynamic changes depending on genetic backgrounds, developmental stage, environment etc. Therefore, more research into epigenetic modifications, such as other histone modifications, acetylation and DNA methylation, are needed to fully understand the genome-wide regulation network.
Extensive studies have revealed that a large number of genes exhibited differential expression level between F1 hybrids and their parents [8, 10, 12, 13, 55, 56]. However, the regulation of gene expression profiles in the heterozygous state remains unknown. Our data revealed that the positive GCA effect parent in the F1 hybrids could determine the phenotypes of the important agronomic traits, such as plant height and heading date. Plant height and heading date are very important agronomic traits that affect rice grain yield [27, 57]. Regulation of the genes in the signal transduction pathway and syntheses of metabolic products have been reported to occur through SNPs in these gene sequences that could be optimized to the gene expression profiles in elite varieties . Our results indicated that the phenotypes of parental bias in F1 hybrids could be the consequence of transcriptome profile or some specific gene expression level.
One of the most important characteristics of the elite parent with high GCA effect is its wide adaptability. The wide adaptability of the elite parent could be due to both photoperiod and temperature neutral in order to adapt different ecological environments. Therefore, rice plant could maintain a constant growth period so that it can develop higher biomass and grain yield in different ecological environments. Previous study has reported the OsPRR37 had involved in the regulation network of the adaptability to grow in different latitudes [45, 57]. In present work, we found the H3K36me3 modification patterns of the OsPRR gene family, OsPRR37 and OsPRR95, showed similar patterns in F1 hybrids, high GCA parents 93–11 and TQ. The results indicated that the specific modification patterns of H3K36me3 in high GCA parents might contribute to the adaptability for elite parents and their offsprings. It is reasonable to assume that the PRR gene family could play an important role in controlling the adaptability of an elite parent that is responsible for GCA effects in rice.
The results presented here revealed that the transcriptome profiles and DEGs in the F1s were highly related to phenotype bias to the positive GCA-effect parent. The transcriptome bias was also demonstrated by analyzing the genes that controlling the specific phenotypes of plant height and flowering and H3K36me3 modification on global modification or specific signaling pathways. Our results indicated the molecular basis of GCA is both on transcriptome and epigenetic level.
Plant material and phenotype analysis
Five rice elite varieties, GL, TQ, 93–11, AJ and ZS, were chosen for this study. Their pedigree as showed in Additional file 6. GCA evaluation was conducted by a diallele crossing design. Briefly, 10 parallel crosses derived from five parents and the phenotypes were collected from a field in Wuhan in 2008 and 2009 and from the Hainan Island province in 2009. The plots were established in triplicate with 30 plants per plot in every season. Three agronomic traits, plant height (PH), heading date (HD) and grain number per panicle (GNP), were studied. The general combining ability (GCA) for each parent was calculated according to Griffing .
Nuclear RNA extraction
The second fully expanded leaves were harvested at the secondary branch differentiation stage, immediately frozen in liquid nitrogen and stored at -80°C. The leaves from the triplicate plots were pooled for RNA extraction. Nuclei were isolated from ~10 g of frozen leaves using the Plant Nuclei Isolation/Extraction Kit (Sigma, St. Louis, MO, USA). Total hnRNA was extracted from nuclei using Trizol (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions, and then treated with RNase-free DNase I (New England Biolabs, Ipswich, MA, USA) to remove any contaminating genomic DNA.
RNA-Seq Library construction
The Illumina mRNA-Seq Sample Prep Kit (Illumina, San Diego, CA USA) was used to prepare the sequencing library with 3 μg of nuclear RNA. Fragmentation buffer in the kit was added directly to hnRNA to produce short fragments of 200–700 bp, which served as the templates for first-strand cDNA synthesis using random hexamers. Second-strand cDNA was synthesized followed the protocol described in the kit and was purified using a QIAquick PCR Extraction Kit (Qiagen, Valencia,CA USA) and eluted in elution buffer (EB). The short fragments were then ligated to sequencing adapters. Suitable fragments of approximately 200 bp were selected as templates for amplification in a MyCycler PCR instrument (Bio-Rad, Hercules, CA USA) with the following program: denaturation at 98°C for 30 s followed by 15 cycles of 98°C for 10 s, 65°C for 30 s, and 72°C for 30 s plus a terminal hold at 72°C for 5 min. The samples were then purified using the QIAquick PCR Purification Kit according to the manufacturer’s protocol and eluted in 30 μL of EB. One μL aliquot of the construct was loaded onto an Agilent Technologies 2100 Bioanalyzer using the Agilent DNA 1000 Chip Kit (Agilent, Santa Clara, CA USA). After verifying the size and purity of the DNA fragments, the library was sequenced using an Illumina GA II x platform by BGI (Shenzhen, China).
ChIP-Seq library generation
Chromatin immunoprecipitation (ChIP) was performed with antibodies against trimethylated histone H3 on lysine 27 (H3K27me3, Abcam Cat. #ab6002, Cambridge, MA, USA) and trimethylated histone H3 on lysine 36 (H3K36me3) (Abcam, Cat. #ab9050, Cambridge, MA, USA) as described by Saleh et al. . The DNA was extracted by adding equal volumes of phenol/chloroform/isoamyl alcohol to each tube and vortexing briefly. The DNA was precipitated with 2.5 volumes of 100% EtOH, 1/10 volume of 3 M sodium acetate (pH 5.2) and 2 μl of glycogen (20 mg ml-1) at -80°C. Next, the ChIP DNA was used to generate Illumina sequencing libraries following the manufacturer’s protocol of Paired-End DNA Sample Prep Kit (Catalog #: PE-102-1001, Illumina, San Diego, CA, USA): appropriate fragments approximately 200 base pairs (bp) in length were selected as templates for amplification involving denaturation at 98°C for 30 s; 15 cycles of 98°C for 10 s, 65°C for 30 s and 72°C for 30 s; and a final incubation at 72°C for 5 min. The samples were then purified using the QIAquick PCR Purification Kit according to the manufacturer’s protocol and were eluted in 30 μl of elution buffer. One microliter of the library was loaded on an Agilent Technologies 2100 Bioanalyzer using the Agilent DNA 1000 chip kit (Agilent, part #5067–1504). After its size and purity were verified, the library was sequenced using the Illumina GAIIx platform.
Analysis of ChIP-Seq reads with the Reads Per Kb per Million reads (RPKM) method
A threshold of more than a 2-fold difference in gene expression levels and a False Discovery Rate (FDR) of less than 0.05 were used to identify the DEGs. An FDR less than 0.001 were used to identify the genes differentially modifying histones. The evaluation of the P value and FDR of DEGs was performed as described by Audic et al.  and Benjamini et al. , respectively. The t-test and correlation analysis were conducted using Microsoft Office Excel 2010.
Quantitation real time PCR
General combining ability
Specific combining ability
Differentially expressed gene
False discovery rate
Reads Per Kb per Million reads.
Sprague GF, Tatum LA: General Vs specific combining ability in single crosses of Corn. J Ame Soc Agron. 1942, 34 (10): 923-932. 10.2134/agronj1942.00021962003400100008x.
Viana JM, Delima RO, Mundim GB, Conde AB, Vilarinho AA: Relative efficiency of the genotypic value and combining ability effects on reciprocal recurrent selection. Theor Appl Genet. 2012, 126 (4): 889-899.
Moterle LM, Braccini AL, Scapim CA, Pinto RJ, Goncalves LS, do Amaral Junior AT, Silva TR: Combining ability of tropical maize lines for seed quality and agronomic traits. Genet Mol Res. 2011, 10 (3): 2268-2278.
Ahangar L, Ranjbar GA, Nouroozi M: Estimation of combining ability for yield and yield components in rice (Oryza sativa L.) cultivars using diallel cross. Pak J Biol Sci. 2008, 11 (9): 1278-1281. 10.3923/pjbs.2008.1278.1281.
Joshi SK, Sharma SN, Singhania DL, Sain RS: Combining ability in the F1 and F2 generations of diallel cross in hexaploid wheat (Triticum aestivum L. em. Thell). Hereditas. 2004, 141 (2): 115-121. 10.1111/j.1601-5223.2004.01730.x.
Qu Z, Li L, Luo J, Wang P, Yu S, Mou T, Zheng X, Hu Z: QTL mapping of combining ability and heterosis of agronomic traits in rice backcross recombinant inbred lines and hybrid crosses. PLoS One. 2012, 7 (1): e28463-10.1371/journal.pone.0028463.
Qi H, Huang J, Zheng Q, Huang Y, Shao R, Zhu L, Zhang Z, Qiu F, Zhou G, Zheng Y, Yue B: Identification of combining ability loci for five yield-related traits in maize using a set of testcrosses with introgression lines. Theor Appl Genet. 2012, 126 (2): 369-377.
Meyer RC, Witucka-Wall H, Becher M, Blacha A, Boudichevskaia A, Dormann P, Fiehn O, Friedel S, von Korff M, Lisec J, Melzer M, Repsilber D, Schmidt R, Scholz M, Selbig J, Willmitzer L, Altmann T: Heterosis manifestation during early Arabidopsis seedling development is characterized by intermediate gene expression and enhanced metabolic activity in the hybrids. Plant J. 2012, 71 (4): 669-683. 10.1111/j.1365-313X.2012.05021.x.
Fujimoto R, Taylor JM, Shirasawa S, Peacock WJ, Dennis ES: Heterosis of Arabidopsis hybrids between C24 and Col is associated with increased photosynthesis capacity. Proc Natl Acad Sci USA. 2012, 109 (18): 7109-7114. 10.1073/pnas.1204464109.
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 X: Global epigenetic and transcriptional trends among two rice subspecies and their reciprocal hybrids. Plant Cell. 2010, 22 (1): 17-33. 10.1105/tpc.109.072041.
Song G, Guo Z, Liu Z, Cheng Q, Qu X, Chen R, Jiang D, Liu C, Wang W, Sun Y, Zhang L, Zhu Y, Yang D: Global RNA sequencing reveals that genotype-dependent allele-specific expression contributes to differential expression in rice F1 hybrids. BMC Plant Biol. 2013, 13: 221-10.1186/1471-2229-13-221.
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 USA. 2009, 106 (19): 7695-7701. 10.1073/pnas.0902340106.
Stupar RM, Gardiner JM, Oldre AG, Haun WJ, Chandler VL, Springer NM: Gene expression analyses in maize inbreds and hybrids with varying levels of heterosis. BMC Plant Biol. 2008, 8: 33-10.1186/1471-2229-8-33.
Stupar RM, Springer NM: Cis-transcriptional variation in maize inbred lines B73 and Mo17 leads to additive expression patterns in the F1 hybrid. Genetics. 2006, 173 (4): 2199-2210. 10.1534/genetics.106.060699.
Riedelsheimer C, Czedik-Eysenberg A, Grieder C, Lisec J, Technow F, Sulpice R, Altmann T, Stitt M, Willmitzer L, Melchinger AE: Genomic and metabolic prediction of complex heterotic traits in hybrid maize. Nat Genet. 2012, 44 (2): 217-220. 10.1038/ng.1033.
Birchler JA, Yao H, Chudalayandi S, Vaiman D, Veitia RA: Heterosis. Plant Cell. 2010, 22 (7): 2105-2112. 10.1105/tpc.110.076133.
Birchler JA, Veitia RA: The gene balance hypothesis: implications for gene regulation, quantitative traits and evolution. New Phytol. 2010, 186 (1): 54-62. 10.1111/j.1469-8137.2009.03087.x.
Hochholdinger F, Hoecker N: Towards the molecular basis of heterosis. Trends Plant Sci. 2007, 12 (9): 427-432. 10.1016/j.tplants.2007.08.005.
Inigo S, Alvarez MJ, Strasser B, Califano A, Cerdan PD: PFT1, the MED25 subunit of the plant Mediator complex, promotes flowering through CONSTANS dependent and independent mechanisms in Arabidopsis. Plant J. 2012, 69 (4): 601-612. 10.1111/j.1365-313X.2011.04815.x.
Endo-Higashi N, Izawa T: Flowering time genes Heading date 1 and Early heading date 1 together control panicle development in rice. Plant Cell Physiol. 2011, 52 (6): 1083-1094. 10.1093/pcp/pcr059.
Takahashi Y, Teshima KM, Yokoi S, Innan H, Shimamoto K: Variations in Hd1 proteins, Hd3a promoters, and Ehd1 expression levels contribute to diversity of flowering time in cultivated rice. Proc Natl Acad Sci USA. 2009, 106 (11): 4555-4560. 10.1073/pnas.0812092106.
Matsubara K, Yamanouchi U, Wang ZX, Minobe Y, Izawa T, Yano M: Ehd2, a rice ortholog of the maize INDETERMINATE1 gene, promotes flowering by up-regulating Ehd1. Plant Physiol. 2008, 148 (3): 1425-1435. 10.1104/pp.108.125542.
Komiya R, Ikegami A, Tamaki S, Yokoi S, Shimamoto K: Hd3a and RFT1 are essential for flowering in rice. Development. 2008, 135 (4): 767-774. 10.1242/dev.008631.
Doi K, Izawa T, Fuse T, Yamanouchi U, Kubo T, Shimatani Z, Yano M, Yoshimura A: Ehd1, a B-type response regulator in rice, confers short-day promotion of flowering and controls FT-like gene expression independently of Hd1. Genes Dev. 2004, 18 (8): 926-936. 10.1101/gad.1189604.
Kojima S, Takahashi Y, Kobayashi Y, Monna L, Sasaki T, Araki T, Yano M: Hd3a, a rice ortholog of the Arabidopsis FT gene, promotes transition to flowering downstream of Hd1 under short-day conditions. Plant Cell Physiol. 2002, 43 (10): 1096-1105. 10.1093/pcp/pcf156.
Farre EM, Liu T: The PRR family of transcriptional regulators reflects the complexity and evolution of plant circadian clocks. Curr Opin Plant Biol. 2013, 16 (5): 621-629. 10.1016/j.pbi.2013.06.015.
Xue W, Xing Y, Weng X, Zhao Y, Tang W, Wang L, Zhou H, Yu S, Xu C, Li X, Zhang Q: Natural variation in Ghd7 is an important regulator of heading date and yield potential in rice. Nat Genet. 2008, 40 (6): 761-767. 10.1038/ng.143.
Murakami M, Ashikari M, Miura K, Yamashino T, Mizuno T: The evolutionarily conserved OsPRR quintet: rice pseudo-response regulators implicated in circadian rhythm. Plant Cell Physiol. 2003, 44 (11): 1229-1236. 10.1093/pcp/pcg135.
Ma Q, Hedden P, Zhang Q: Heterosis in rice seedlings: its relationship to gibberellin content and expression of gibberellin metabolism and signaling genes. Plant Physiol. 2011, 156 (4): 1905-1920. 10.1104/pp.111.178046.
Yamaguchi S: Gibberellin metabolism and its regulation. Annu Rev Plant Biol. 2008, 59: 225-251. 10.1146/annurev.arplant.59.032607.092804.
Olszewski N, Sun TP, Gubler F: Gibberellin signaling: biosynthesis, catabolism, and response pathways. Plant Cell. 2002, 14 (Suppl): S61-S80.
Hedden P, Phillips AL: Gibberellin metabolism: new insights revealed by the genes. Trends Plant Sci. 2000, 5 (12): 523-530. 10.1016/S1360-1385(00)01790-8.
Wang X, Yao Y, Peng H, Zhang Y, Lu L, Ni Z, Sun Q: The relationship of differential expression of genes in GA biosynthesis and response pathways with heterosis of plant height in a wheat diallel cross. Chin Sci Bull. 2009, 54: 3029-3034. 10.1007/s11434-009-0518-3.
Zhang Y, Ni Z, Yao Y, Nie X, Sun Q: Gibberellins and heterosis of plant height in wheat (Triticum aestivum L.). BMC Genet. 2007, 8: 40-
Jaenisch R, Bird A: Epigenetic regulation of gene expression: how the genome integrates intrinsic and environmental signals. Nat Genet. 2003, 33 (Suppl): 245-254.
Stancheva I, El-Maarri O, Walter J, Niveleau A, Meehan RR: DNA methylation at promoter regions regulates the timing of gene activation in Xenopus laevis embryos. Dev Biol. 2002, 243 (1): 155-165. 10.1006/dbio.2001.0560.
Futscher BW, Oshiro MM, Wozniak RJ, Holtan N, Hanigan CL, Duan H, Domann FE: Role for DNA methylation in the control of cell type specific maspin expression. Nat Genet. 2002, 31 (2): 175-179. 10.1038/ng886.
Monk M: Epigenetic programming of differential gene expression in development and evolution. Dev Genet. 1995, 17 (3): 188-197. 10.1002/dvg.1020170303.
Roudier F, Ahmed I, Berard C, Sarazin A, Mary-Huard T, Cortijo S, Bouyer D, Caillieux E, Duvernois-Berthet E, Al-Shikhley L, Giraut L, Despres B, Drevensek S, Barneche F, Derozier S, Brunaud V, Aubourg S, Schnittger A, Bowler C, Martin-Magniette ML, Robin S, Caboche M, Colot V: Integrative epigenomic mapping defines four main chromatin states in Arabidopsis. EMBO J. 2011, 30: 1928-1938. 10.1038/emboj.2011.103.
Zhang X, Clarenz O, Cokus S, Bernatavichute YV, Pellegrini M, Goodrich J, Jacobsen SE: Whole-genome analysis of histone H3 lysine 27 trimethylation in Arabidopsis. PLoS Biol. 2007, 5 (5): e129-10.1371/journal.pbio.0050129.
Pokholok DK, Harbison CT, Levine S, Cole M, Hannett NM, Lee TI, Bell GW, Walker K, Rolfe PA, Herbolsheimer E, Zeitlinger J, Lewitter F, Gifford DK, Young RA: Genome-wide map of nucleosome acetylation and methylation in yeast. Cell. 2005, 122 (4): 517-527. 10.1016/j.cell.2005.06.026.
Bannister AJ, Schneider R, Myers FA, Thorne AW, Crane-Robinson C, Kouzarides T: Spatial distribution of di- and tri-methyl lysine 36 of histone H3 at active genes. J Biol Chem. 2005, 280 (18): 17732-17736. 10.1074/jbc.M500796200.
Strahl BD, Grant PA, Briggs SD, Sun ZW, Bone JR, Caldwell JA, Mollah S, Cook RG, Shabanowitz J, Hunt DF, Allis CD: Set2 is a nucleosomal histone H3-selective methyltransferase that mediates transcriptional repression. Mol Cell Biol. 2002, 22 (5): 1298-1306. 10.1128/MCB.22.5.1298-1306.2002.
Liu H, Yu X, Li K, Klejnot J, Yang H, Lisiero D, Lin C: Photoexcited CRY2 interacts with CIB1 to regulate transcription and floral initiation in Arabidopsis. Science. 2008, 322 (5907): 1535-1539. 10.1126/science.1163927.
Koo BH, Yoo SC, Park JW, Kwon CT, Lee BD, An G, Zhang Z, Li J, Li Z, Paek NC: Natural variation in OsPRR37 regulates heading date and contributes to rice cultivation at a wide range of latitudes. Mol Plant. 2013, 6 (6): 1877-1888. 10.1093/mp/sst088.
Zhou G, Chen Y, Yao W, Zhang C, Xie W, Hua J, Xing Y, Xiao J, Zhang Q: Genetic composition of yield heterosis in an elite rice hybrid. Proc Natl Acad Sci USA. 2012, 109 (39): 15847-15852. 10.1073/pnas.1214141109.
Lu H, Romero-Severson J, Bernardo R: Genetic basis of heterosis explored by simple sequence repeat markers in a random-mated maize population. Theor Appl Genet. 2003, 107 (3): 494-502. 10.1007/s00122-003-1271-7.
Hua J, Xing Y, Wu W, Xu C, Sun X, Yu S, Zhang Q: Single-locus heterotic effects and dominance by dominance interactions can adequately explain the genetic basis of heterosis in an elite rice hybrid. Proc Natl Acad Sci USA. 2003, 100 (5): 2574-2579. 10.1073/pnas.0437907100.
Fu J, Cheng Y, Linghu J, Yang X, Kang L, Zhang Z, Zhang J, He C, Du X, Peng Z, Wang B, Zhai L, Dai C, Xu J, Wang W, Li X, Zheng J, Chen L, Luo L, Liu J, Qian X, Yan J, Wang J, Wang G: RNA sequencing reveals the complex regulatory network in the maize kernel. Nat Commun. 2013, 4: 2832-
Holloway B, Luck S, Beatty M, Rafalski JA, Li B: Genome-wide expression quantitative trait loci (eQTL) analysis in maize. BMC Genomics. 2011, 12: 336-10.1186/1471-2164-12-336.
Swanson-Wagner RA, DeCook R, Jia Y, Bancroft T, Ji T, Zhao X, Nettleton D, Schnable PS: Paternal dominance of trans-eQTL influences gene expression patterns in maize hybrids. Science. 2009, 326 (5956): 1118-1120. 10.1126/science.1178294.
Zhang M, Zhao H, Xie S, Chen J, Xu Y, Wang K, Guan H, Hu X, Jiao Y, Song W, Lai J: Extensive, clustered parental imprinting of protein-coding and noncoding RNAs in developing maize endosperm. Proc Natl Acad Sci USA. 2011, 108 (50): 20042-20047. 10.1073/pnas.1112186108.
Waters AJ, Makarevitch I, Eichten SR, Swanson-Wagner RA, Yeh CT, Xu W, Schnable PS, Vaughn MW, Gehring M, Springer NM: Parent-of-origin effects on gene expression and DNA methylation in the maize endosperm. Plant Cell. 2011, 23 (12): 4221-4233. 10.1105/tpc.111.092668.
Nodine MD, Bartel DP: Maternal and paternal genomes contribute equally to the transcriptome of early plant embryos. Nature. 2012, 482 (7383): 94-97. 10.1038/nature10756.
Meyer S, Pospisil H, Scholten S: Heterosis associated gene expression in maize embryos 6 days after fertilization exhibits additive, dominant and overdominant pattern. Plant Mol Biol. 2007, 63 (3): 381-391. 10.1007/s11103-006-9095-x.
Bao J, Lee S, Chen C, Zhang X, Zhang Y, Liu S, Clark T, Wang J, Cao M, Yang H, Wang SM, Yu J: Serial analysis of gene expression study of a hybrid rice strain (LYP9) and its parental cultivars. Plant Physiol. 2005, 138 (3): 1216-1231. 10.1104/pp.105.060988.
Yan W, Liu H, Zhou X, Li Q, Zhang J, Lu L, Liu T, Zhang C, Zhang Z, Shen G, Yao W, Chen H, Yu S, Xie W, Xing Y: Natural variation in Ghd7.1 plays an important role in grain yield and adaptation in rice. Cell Res. 2013, 23 (7): 969-971. 10.1038/cr.2013.43.
Griffing B: Concept of general and specific combining ability in relation to diallel crossing systems. Aust J Biol Sci. 1956, 9: 463-493.
Saleh A, Alvarez-Venegas R, Avramova Z: An efficient chromatin immunoprecipitation (ChIP) protocol for studying histone modifications in Arabidopsis plants. Nat Protoc. 2008, 3 (6): 1018-1025. 10.1038/nprot.2008.66.
Barski A, Cuddapah S, Cui K, Roh TY, Schones DE, Wang Z, Wei G, Chepelev I, Zhao K: High-resolution profiling of histone methylations in the human genome. Cell. 2007, 129 (4): 823-837. 10.1016/j.cell.2007.05.009.
Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 7 (10): 986-995.
Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I: Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001, 125 (1–2): 279-284.
Wang W, Liu Z, Guo Z, Song G, Cheng Q, Jiang D, Zhu Y, Yang D: Comparative transcriptomes profiling of photoperiod-sensitive male sterile rice Nongken 58S during the male sterility transition between short-day and long-day. BMC Genomics. 2011, 12: 462-10.1186/1471-2164-12-462.
We are grateful for BGI at Shenzhen carried out RNA sequencing and data analysis. This work was supported by National Basic Research Program of China (grant No. 2011CB100102). We are thankful to Prof. Penny von Wettstein-Knowles (Copenhagen University) for critical reading of the manuscript.
The authors declare that they have no competing interests.
GY, S performed the most of the experiments; ZB, G performed the epigenetics experiments; DM, J performed all field experiments; ZW, L, XF, Q, W, W, helped GY, S performed crossing and analysis of phenotypes; DC, Y and YG, Z supervised this study; DC, Y and GY, S designed the experiments and wrote the manuscript. All the authors discussed the results and contributed to the manuscript. All authors read and approved the final manuscript.