The phenotypic predisposition of the parent in F1 hybrid is correlated with transcriptome preference of the positive general combining ability parent
© Song et al.; licensee BioMed Central Ltd. 2014
Received: 20 January 2014
Accepted: 10 April 2014
Published: 22 April 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.
KeywordsGeneral combining ability (GCA) Transcriptome profiling Hybrid rice Transcriptome and phenotype bias Diallelic crosses
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
The correlation of phenotypes between each parent and their F1s
The phenotype bias toward the positive GCA effect parent is attributable to transcriptome bias
Differentially expressed genes in the F1 hybrids are those preferentially expressed in the positive GCA effect parent
Expression patterns of the genes in the regulatory pathways controlling flowering time and plant height in rice
Expression level of flowering- and GA metabolism-related genes
Expression level (RPKM)
GA metabolism pathway
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
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.
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.
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Birchler JA, Yao H, Chudalayandi S, Vaiman D, Veitia RA: Heterosis. Plant Cell. 2010, 22 (7): 2105-2112. 10.1105/tpc.110.076133.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Yamaguchi S: Gibberellin metabolism and its regulation. Annu Rev Plant Biol. 2008, 59: 225-251. 10.1146/annurev.arplant.59.032607.092804.PubMedView ArticleGoogle Scholar
- Olszewski N, Sun TP, Gubler F: Gibberellin signaling: biosynthesis, catabolism, and response pathways. Plant Cell. 2002, 14 (Suppl): S61-S80.PubMed CentralPubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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-PubMed CentralPubMedView ArticleGoogle Scholar
- Jaenisch R, Bird A: Epigenetic regulation of gene expression: how the genome integrates intrinsic and environmental signals. Nat Genet. 2003, 33 (Suppl): 245-254.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Monk M: Epigenetic programming of differential gene expression in development and evolution. Dev Genet. 1995, 17 (3): 188-197. 10.1002/dvg.1020170303.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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-PubMedGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Griffing B: Concept of general and specific combining ability in relation to diallel crossing systems. Aust J Biol Sci. 1956, 9: 463-493.Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 7 (10): 986-995.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.