Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Genome-wide analyses of four major histone modifications in Arabidopsis hybrids at the germinating seed stage

BMC Genomics volume 18, Article number: 137 (2017) Cite this article

Abstract

Background

Hybrid vigour (heterosis) has been used for decades in cropping agriculture, especially in the production of maize and rice, because hybrid varieties exceed their parents in plant biomass and seed yield. The molecular basis of hybrid vigour is not fully understood. Previous studies have suggested that epigenetic systems could play a role in heterosis.

Results

In this project, we investigated genome-wide patterns of four histone modifications in Arabidopsis hybrids in germinating seeds. We found that although hybrids have similar histone modification patterns to the parents in most regions of the genome, they have altered patterns at specific loci. A small subset of genes show changes in histone modifications in the hybrids that correlate with changes in gene expression. Our results also show that genome-wide patterns of histone modifications in geminating seeds parallel those at later developmental stages of seedlings.

Conclusion

Ler/C24 hybrids showed similar genome-wide patterns of histone modifications as the parents at an early germination stage. However, a small subset of genes, such as FLC, showed correlated changes in histone modification and in gene expression in the hybrids. The altered patterns of histone modifications for those genes in hybrids could be related to some heterotic traits in Arabidopsis, such as flowering time, and could play a role in hybrid vigour establishment.

Background

Hybrid vigour, or heterosis, is a phenomenon where the progeny derived from crosses between two accessions of a species have increased performance compared to their parents. Hybrids have been used in agriculture for over a century, achieving increased seed yields in crops, such as maize [1] and rice [2]. Besides seed yield, heterosis also applies to growth rate and biomass. In Arabidopsis, hybrids generated by crosses between different ecotypes show strong heterosis in many traits, especially in vegetative biomass and in seed yield [3–5].

The molecular bases of heterosis are still unclear. Genome-wide analyses of Arabidopsis hybrid transcriptomes indicate that thousands of genes have expression levels different from the average levels of the genes in the parents (non-additive expression) [6–9]. Epigenetic systems may be involved in heterosis (reviewed in [10]). In hybrids, the altered patterns of two epigenetic systems, small RNAs and DNA methylation, can lead to changed expression levels of some associated genes [8, 11–13].

Little is known about the role of another epigenetic system in hybrids, histone modification. Histone modification refers to post-translational modification of amino acids in histone proteins, including methylation, acetylation and phosphorylation of particular residues. The histone marks H3K4me3 and H3K9ac, are associated with active gene expression, and the histone marks H3K27me3 and H3K9me2 are associated with genes with low transcript levels [14]. The active marks and H3K27me3 are located along chromosome arms, while H3K9me2 is located in heterochromatin and pericentromeric regions, usually in transposable elements (TEs).

Genome-wide patterns of histone modification have been documented in hybrids of a number of plant species. Compared to parental lines, hybrids of rice and maize have changed patterns of histone modification in some regions of the genome, and the changes in histone modifications correlate with changes in levels of gene expression [15, 16]. In Arabidopsis hybrids, the genome-wide histone modification patterns of the parents are mostly inherited in the F1 genome [17, 18], but there are some genes with non-additive levels of gene expression which also have non-additive levels of histone modification. For example, in allotetraploids derived from crosses between A. thaliana and A. arenosa, the altered levels of histone marks H3K9ac and H3K4me2 are associated with non-additive expression levels of the key circadian clock genes which are involved in energy production and storage [19].

Non-additive gene expression has been reported in hybrids at early stages of seedling development, and could play a vital role in the establishment of biomass heterosis at later developmental stages [4, 6, 20, 21]. We investigated the patterns of four histone marks, H3K4me3, H3K9ac, H3K27me3 and H3K9me2 in the genomes of hybrids at the early stage of seed germination. The hybrids were produced in reciprocal crosses between two Arabidopsis ecotypes, Landsberg erecta (Ler) and C24. Our results show that although the hybrids have unchanged histone modification patterns compared to the parents in most regions of the genome, histone modification patterns are altered at specific loci. We compared the non-additive changes of histone modifications with the non-additive levels of gene expression, and found that for genes with changed histone modifications in hybrids, the majority do not have corresponding changes in gene expression. A correlation between changes in gene expression and histone modification was found in only a small subset of genes in the hybrids. Previous studies on histone modifications concentrated on mature developmental stages of seedlings; however, little is known about the histone modification patterns in germinating seeds. Our results suggest that geminating seeds have genome-wide patterns of histone modifications similar to those in seedlings at later developmental stages.

Results

Genome-wide patterns of histone modifications in germinating seeds

In germinating seeds, high levels of H3K4me3, H3K9ac and H3K27me3 were found along the arms of all five chromosomes but not in the pericentromeric regions (Additional file 1: Figure S1). In contrast, H3K9me2 preferentially targets the pericentromeric regions of chromosomes. The chromosomal patterns of the four histone marks in the hybrids are similar to those in the parents (Additional file 1: Figure S1).

Analyses of the distributions of histone marks in genes and their surrounding regions show that H3K4me3 and H3K9ac are enriched downstream of the transcription start sites of protein-coding genes, while H3K27me3 occurs along the gene body (Fig. 1a and c). These three marks do not occur in transposable elements (Fig. 1b and d). H3K9me2 has high levels in long TEs (≥1.5 kb) and lower levels in the short TEs (<1.5 kb; Fig. 1d) and in the gene body of protein-coding genes (Fig. 1a and c). H3K4me3 and H3K9ac occur at most protein-coding genes (Fig. 1c; Additional file 1: Figure S2a-c). High levels of H3K27me3 were present in approximately 20% (5627 genes) of protein-coding genes; these genes usually have low levels of the active marks (Fig. 1c). H3K9me2 is found only in a small number of protein-coding genes (Additional file 1: Figure S2c), in which the other three marks are not detected (Fig. 1c). The results show that histone modification patterns in germinating seeds parallel those in seedlings at later stages [22–24].

Fig. 1
figure 1

Global patterns of the four histone marks in C24. a Average levels of histone modification signals over gene body and the ±3 kb surrounding regions in protein-coding genes. TSS, transcription start site; TTS, transcription termination site. b Average levels of histone modification signals over long TE (≥1.5 kb) and the ±3 kb regions. c Heat maps show the enrichment levels of histone modifications in all the protein-coding genes in Arabidopsis. Genes are ordered from top to bottom with decreasing K4me3 levels over the gene body regions. The intensity of red represents the levels of modifications. d Heat maps show the enrichment levels of histone modifications in all the TE in Arabidopsis. TEs are ordered by length

Full size image

We compared the histone modification data with our transcriptome data collected at the same time point of seed germination [21]. The expression of protein-coding genes is positively correlated with the active marks, H3K4me3 and H3K9ac, and is negatively correlated with the repressive marks, H3K27me3 and H3K9me2 (Additional file 1: Figure S3a-d). There was a negative correlation between H3K9me2 and expression of TE genes (genes with TEs in the gene body; Additional file 1: Figure S3e).

Genes with altered histone modification levels in parents and hybrids

Statistical analyses of whole genome sequencing data showed that the majority of genes with histone marks (18261 for H3K4me3; 16368 for H3K9ac; 8801 for H3K27me3; 10177 for H3K9me2) have similar levels of histone modifications in the parents, Ler and C24 (Fig. 2). Approximately 6% of the genes are differentially modified in the parents for the two active marks, H3K4me3 and H3K9ac, and there is a higher proportion of differentially modified genes (DMGs) with the repressive marks (14% for H3K27me3 and 20% for H3K9me2). Of the DMGs, the majority have higher levels of the two active marks and H3K9me2 in Ler than in C24, and more genes have higher H3K27me3 levels in C24 than in Ler (Fig. 2).

Fig. 2
figure 2

The proportions of genes with differential histone modification levels between Ler and C24 for the four histone marks. Fold-change ≥ 1.5 and p-value ≤ 0.05

Full size image

We showed that hybrids have similar histone modification patterns to the parents at a chromosomal scale. However, at localised regions of the genome, there are hundreds of genes with altered histone modification levels in the hybrids relative to the average histone modification levels of the genes in the parents. There are 1,430 genes (16.5%) differentially modified with H3K27me3 in Ler × C24, whereas only 0.07% - 1.7% of genes have altered levels of the other marks (Fig. 3a). Almost all of H3K4me3-, H3K9ac and H3K27me3-associated DMGs have high modification levels relative to the average histone modification levels of the parents, whereas most H3K9me2-associated DMGs have decreased levels in C24 × Ler (Fig. 3a). In the reciprocal hybrids, Ler × C24 has significantly more genes differentially modified with the active marks and H3K27me3 than C24 × Ler, but has significantly fewer H3K9me2-associated DMGs (Fig. 3a). Of the DMGs in hybrids, only a small subset are also differentially modified in the parents (Fig. 3b; Additional file 1: Figure S4). For all four marks, there are small numbers of the DMGs found in both reciprocal hybrids, suggesting that in germinating seeds the changes in histone modifications occur at particular genes in both hybrids. In addition, there is a trend for Ler × C24 to have more genes with increased levels of each of the histone marks relative to the average levels of the parents than C24 × Ler (Fig. 4a), where the majority of genes have additive modification levels similar to the average modification levels of the parents.

Fig. 3
figure 3

The genes identified with differential histone modification levels in Ler/C24 hybrids relative to the average modification levels of the parents. a The numbers of differentially modified genes (fold-change ≥ 1.25 and p-value ≤ 0.05) for the four histone marks in the reciprocal hybrids. The percentages indicate proportions of DMGs in all the genes with each histone mark. Asterisks indicate significant differences of gene numbers between reciprocal hybrids (*p-value ≤ 0.05; **p-value ≤ 0.01; ***p-value ≤ 0.001). P-values were calculated using the Fisher’s Exact Test. b Venn-diagram showing the overlaps between the differentially modified genes by K27me3 in the parents and hybrids

Full size image
Fig. 4
figure 4

Box plots showing histone modification levels relative to the average modification levels of the parents (=1) in the reciprocal hybrids at the genes where the parents have differential modification levels. HP: high-parent; LP: Low-parent. Squares indicate the outliers

Full size image

Alterations in histone modifications correlate with alterations in gene expression in specific genes in hybrids

We used our previous transcriptome datasets to determine if gene expression in the hybrids is influenced by the changes in histone modifications. In the parents, we found that expression levels are positively correlated with the levels of active marks (Additional file 1: Figure S5). The expected negative correlation between repressive marks and gene expression was only observed for a subset of the analysed genes. A similar poor correlation between gene expression and H3K27me3 was previously observed in hybrids between rice subspecies [15].

To minimise the influence of genes that are not regulated by histone marks, only the DMGs in the parents having the expected direction of change in gene expression were selected for the analyses. There are 407 genes for H3K4me3, 457 genes for H3K9ac, 149 genes for H3K27me3 and 97 genes for H3K9me2 meeting these criteria. Unlike the DMGs between the parents, we found poor correlations between gene expression and histone modifications when comparing hybrids to the average levels of histone modifications and gene expression in the parents (Fig. 5a; Additional file 1: Figure S6), suggesting that histone modifications may regulate non-additive gene expression only in specific genes in hybrids. There are 47 genes in the hybrids showing non-additive changes in both gene expression and histone modification (Additional file 2: Data S2), including FLOWERING LOCUS C (FLC), a gene which negatively regulates flowering time and whose expression is repressed by high level of H3K27me3 [25]. At the germinating seed stage, the expression level of FLC is similar to the level in C24 and are up-regulated in the hybrids relative to the average levels of the parents, and is associated with decreased levels of H3K27me3 (Fig. 5b). Consistent with the high expression levels of FLC in the hybrids, the flowering times of Ler/C24 hybrids are later than Ler but similar to C24 [5].

Fig. 5
figure 5

Changes in gene expression are associated with changes in histone modification in hybrids. a Correlations between gene expression and histone modifications in Ler × C24 relative to the average levels of the parents. Red dots represent genes with significant changes in gene expression. Blue dots represent non-significant genes. To determine significant changes in gene expression: fold-change ≥ 1.3, fold-change ≤ 0.01. b K27me3 and transcript levels (read counts) of FLOWERING LOCUS C (FLC; AT5G10140) in the parents and hybrids at germinating seed stage. AVG. indicates the average levels of K27me3 and transcript of the parents. Asterisks indicate that hybrid levels significantly differentiate from the average levels of the parents (p-value ≤ 0.01)

Full size image

High histone modification levels in Ler × C24 are associated with Ler alleles

To reveal the consequences of interaction between the two parental epigenomes in the hybrids, we investigated the modification levels of the Ler and C24 alleles in the hybrids for genes with sufficient SNPs (single-nucleotide polymorphism). 4404, 3772, 366 and 153 genes were identified as being associated with H3K4me3, H3K9ac, H3K27me3 and H3K9me2, respectively. Of these genes, only a small proportion have ratios of Ler alleles and C24 alleles in the hybrids different from the ratios in the parents (Additional file 1: Figure S7), suggesting that hybrids retain allelic patterns of histone modification similar to those in the parents at this germination stage. However, different allelic patterns of histone modifications were found between the reciprocal hybrids for the genes where the parents show differential modification levels (288 for H3K4me3; 291 for H3K9ac; 95 for H3K27me3; 20 for H3K9me2). Except for the H3K9me2-targeted genes, the majority of genes tend to have increased modification levels of Ler alleles and unchanged levels of C24 alleles relative to the expected levels (half of the parental levels) in Ler × C24 (Fig. 6a; Additional file 1: Figure S8), indicating that the high modification levels in Ler × C24 are due to the increased modification levels of Ler alleles. Similar increased modification levels were not seen at either allele of the same genes in the C24 × Ler hybrid.

Fig. 6
figure 6

Allelic changes in histone modifications in the reciprocal hybrids. a box plot showing the allelic K4me3 levels relative to the expected levels in the reciprocal hybrids at the DMGs in the parents. Squares indicate the outliers. b Correlations between gene expression and histone modification at Ler and C24 alleles in the reciprocal hybrids. Dots represent individual genes. Expected levels of Ler and C24 alleles equal half of the corresponding parental levels. c The K4me3 and transcript levels (read counts) of AT4G11845 in hybrid (Ler_a + C24_b) and at the Ler and C24 alleles (Ler_a and C24_b) in the hybrids compared to the expected levels. Expected level of hybrid equals the average levels of the parents; expected levels of Ler and C24 alleles equal half of the corresponding parental levels

Full size image

The histone modification data were compared to the transcriptome data to determine if histone modifications regulate gene expression at the allele level. We found poor correlations between the alterations in gene expression and in histone modifications for both Ler and C24 alleles in the reciprocal hybrids (Fig. 6b; Additional file 1: Figure S9). Both hybrids have unchanged expression levels of Ler alleles and increased expression levels of C24 alleles. This is not consistent with the allelic patterns of histone modifications (Additional file 1: Figure S10), suggesting a limited role of histone modifications in the regulation of allelic gene expression. However, the changes in histone modifications are consistent with the changes in gene expression of alleles of some specific genes (Fig. 6c; Additional file 1: Figure S11). For example, one gene (AT4G11845) encoding an interleukin-1 receptor-associated kinase 4 protein shows up-regulated expression in the hybrids compared to the parents, due to the up-regulation of C24 alleles (Fig. 6c). Consistent with gene expression, we found that the H3K4me3 levels of C24 alleles at this gene are elevated in both reciprocal hybrids.

Discussion

Histone modification patterns at the germinating seed stage parallel those at later stages

Our findings suggest that geminating seeds have genome-wide patterns of histone modifications consistent with the patterns in later developmental stages of seedlings. H3K4me3, H3K9ac and H3K27me3 are located in the chromosome arms [22, 24, 26], whereas H3K9me2 is enriched in the pericentromeric regions [27]. In gene regions, the two active marks peak at the 3’ side of transcription start sites of protein-coding genes [23], and H3K27me3 targets the entire gene body [22]. H3K9me2 is not found in protein-coding genes, instead, it occurs in TEs [24, 28]. We found that there were higher H3K9me2 levels in long TEs than in short TEs. This could be due to the TEs of different sizes being targeted by different DNA methylation pathways [29]. The co-occurrence of the active marks with H3K27me3 has been found in a small proportion of genes [24, 30]. Our results showed a negative correlation between levels of the active marks and levels of H3K27me3 among protein-coding genes, indicating an incompletely mutually exclusive pattern between the two types of marks in the genome. We found that the gene expression is positively correlated with active marks and is negatively correlated with repressive marks, suggesting similar roles of histone modification in gene regulation at the germination stage and at later seedling stages.

Histone modification patterns in Arabidopsis hybrids

Genome-wide additive patterns of histone modification were previously reported for Arabidopsis hybrids at seedling stages [17, 31]. Our results in germinating seeds showed that, although the F1 hybrids of Ler and C24 have the histone modification patterns of their parents at chromosomal levels, they have altered levels of histone modifications in some localised regions of the genome. We identified hundreds of genes with altered histone modification levels (DMGs) in hybrids. Unlike the observations for small RNAs and DNA methylation in hybrids [12, 13], the changes in histone modifications are not dependent on differences between the parents. This may be because we examined histone modification patterns at an early developmental stage. The majority of DMGs in the hybrids have elevated modification levels relative to the average histone modification levels of the parents; similar increases for H3K4me3 and H3K27me3 were observed in rice hybrids at a later stage [15].

The reciprocal hybrids showed different patterns of histone modification in germinating seeds. Ler × C24 has significantly more DMGs than C24 × Ler, except for K9me2-targeted DMGs, and has higher levels of histone modifications than C24 × Ler. The increased levels of histone modifications in Ler × C24 are mainly contributed by increased modification of Ler alleles. The Ler alleles of the same genes in the other reciprocal hybrid, C24 × Ler, retain similar modification levels as those in the parents. The differences between the reciprocal hybrids are probably because the Ler alleles of specific genes in Ler × C24 (from maternal parent) undergo developmental pathways different to the Ler alleles in C24 × Ler (from paternal parent), leading to increased levels of histone modification in only one of the hybrids. The bias between reciprocal hybrids seems to only exist at this early germination stage, because similar histone modification patterns were found in reciprocal hybrids at later stages [15, 31]. Our previous work showed earlier changes in gene expression in Ler × C24 than in C24 × Ler at three days after sowing the seeds [21]. Similar earliness may apply to histone modification patterns in Ler × C24, suggesting that similar increases in histone modifications could occur in C24 × Ler at a later developmental stage.

Potential roles of histone modifications in non-additive gene expression

Our results show poor correlation between changes in histone modifications and changes in gene expression in hybrids relative to the average levels of the parents. At allelic levels, the hybrids display changes in histone modification not paralleled by changes in gene expression. These results suggest that non-additive gene expression in the hybrids during seed germination is not due to alterations of histone modifications for the majority of genes in the genome. Histone modifications may play vital roles in regulating some genes that are responsible for specific biological processes in hybrids. A previous study showed that H3K4me2 and H3K9ac are involved in changing circadian rhythm of allotetraploids derived from the crosses between A. thaliana and A. arenosa, by regulating the expression levels of key circadian clock genes [19]. In our results, there are 47 genes showing non-additive changes in both gene expression and histone modification in the hybrids at the germinating seed stage; some of these genes are involved in critical biological processes, such as root development (AT1G51190), membrane transportation (AT2G38940) and fatty acid degradation (AT4G14440). Our results also show that the loss of H3K27me3 on FLC leads to high FLC expression levels in Ler/C24 hybrids at the germination stage, consistent with the later flowering times of the hybrids [5].

Conclusions

Our genome-wide analyses of four major histone marks showed that Ler/C24 hybrids have similar histone modification patterns compared to the parents at an early germination stage. However, altered levels of histone modifications were found at hundreds of genes in the hybrids. A small number of these genes, such as FLC, show consistent changes in histone modification and in gene expression. The altered histone modifications in the hybrids at germination stage may be involved in heterosis establishment, or at least could be related to some heterotic traits in Arabidopsis hybrids, such as flowering time.

Methods

Plant materials and growth conditions

Seeds of Ler and C24 are available from The Arabidopsis Information Resource (TAIR). Seeds of Ler/C24 hybrid were obtained from hand-pollinated crosses between Ler and C24, and parental seeds were obtained from self-pollinated parents with restricted number of pollinated stigmas. Sterilized seeds of the parents and reciprocal hybrids were placed onto Murashige and Skoog (MS) medium for stratification at 4 °C in the dark for three days. Two replicates of 2000 seeds of each plant line were harvested for ChIP-Seq library preparation. Embryos were isolated from imbibed seeds following the protocol described in previous literature [32].

Chromatin immune-precipitation (ChIP)

Chromatin immune-precipitation was performed according to Helliwell et al., 2006. Micrococcal nuclease (MNase; ThermoFisher, 88216) was used for chromatin digestion. Samples were sonicated at 40 amplitude for two cycles of 15 s by using a UP400S sonicator (Hielscher). Antibodies for precipitating nucleosomes with the four histone marks were anti-H3K4me3 (Milipore,07-473), anti-H3K9ac (Milipore, 07–352), anti-H3K27me3 (Milipore, 07–449) and anti-H3K9me2 (diagenode, pAb-060-050). The precipitated ChIP DNA was purified using a MinElute Reaction Cleanup Kit (Qiagen). ChIP-Seq libraries were prepared using NEBNext® ChIP-Seq Library Prep Reagent Set for Illumina® (NEB, E6200) and NEBNext® Multiplex Oligos for Illumina® (NEB, E7335), following manufacturer’s manuals. Quality and quantity of ChIP-Seq library samples were measured by using a 2100 Bioanalyzer Instrument (Agilent Technologies). ChIP-Seq libraries were verified by using three control genes (Additional file 1: Figure S12), ACTIN7, AGAMOUS and TA3, which are known to be targeted by H3K4me3 (and H3K9ac), H3K27me3 and H3K9me2, respectively.

Bioinformatic analyses

Sequenced reads (Additional file 2: Data S1) of ChIP-Seq were mapped to the TAIR10 reference genome using Biokanga align (http://sourceforge.net/projects/biokanga/) with default settings and additionally applying parameters -M5 and –y10. For H3K9me2 libraries, instead of using –M5, parameters -r5 and -R500 were applied in mapping. Histone modification peaks were called using MACS2 callpeak [33] with p-value less or equal to 10−6 for K4me3 and K9ac and with p-values less or equal to 10−2 for K27me3 and K9me2. Peaks of K27me3 and K9me2 were called with additional parameter –broad, as these two marks usually form broad and continuous peaks. Duplicated reads (reads having same nucleotide sequences) generated were removed by using MACS2 filterdup. For corresponding peaks in the parents and hybrids, peak coordinates were merged using BEDTools merge [34]. Read counts were obtained for the merged peaks by using BEDTools intersect. Statistical analyses were applied on read counts of peaks between hybrids and parents by using DESEQ (http://bioconductor.org/packages/2.13/bioc/html/DESeq.html). Genes with non-additive levels of histone modifications in hybrids were identified by against the average modification levels of the parents (fold-change ≥ 1.25, p-value ≤ 0.05). Profiles and heat-maps of histone modification levels over genes were drawn by using ngsplot [35]. The histone modification targeted genes were identified by confirming the distances between genes and histone modification peaks. Genes were considered to be targeted by K4me3 and K9ac when the peaks overlap the 1.5 kb regions upstream of the TSS, while genes were considered to be targeted by K9me2 when the peaks overlap the 1 kb regions upstream or downstream of the TSS. For K27me3-targeted genes, at least 50% of peak regions and 50% of gene body regions overlap with each other. Heat maps of histone modification levels of genes in the parents and hybrids were drawn using GENE-E (http://www.broadinstitute.org/cancer/software/GENE-E/) and genes were clustered based on the Pearson correlation algorithm. Venn diagrams were made by using Venny (http://bioinfogp.cnb.csic.es/tools/venny/index.html).

SNP analyses and transcriptome analyses were performed following the methods described in [21].

Abbreviations

DMG:

Differentially modified gene

FLC :

FLOWERING LOCUS C

SNP:

Single-nucleotide polymorphism

TE:

Transposable element

References

  1. Crow JF. 90 years ago: the beginning of hybrid maize. Genetics. 1998;148(3):923–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  2. Yuan LP. Hybrid rice breeding in China. Manila; International Rice Research Institute; 1998.

    Google Scholar 

  3. Barth S, Busimi AK, Friedrich Utz H, Melchinger AE. Heterosis for biomass yield and related traits in five hybrids of Arabidopsis thaliana L. Heynh Heredity. 2003;91(1):36–42.

    Article  CAS  PubMed  Google Scholar 

  4. Meyer RC, Torjek O, Becher M, Altmann T. Heterosis of biomass production in Arabidopsis. Establishment during early development. Plant Physiol. 2004;134(4):1813–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Groszmann M, Gonzalez-Bayon R, Greaves IK, Wang L, Huen AK, Peacock WJ, Dennis ES. Intraspecific Arabidopsis hybrids show different patterns of heterosis despite the close relatedness of the parental genomes. Plant Physiol. 2014;166(1):265.

    Article  PubMed  PubMed Central  Google Scholar 

  6. 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 U S A. 2012;109(18):7109–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Meyer RC, Witucka-Wall H, Becher M, Blacha A, Boudichevskaia A, Dormann P, Fiehn O, Friedel S, von Korff M, Lisec J, et al. 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–83.

    Article  CAS  PubMed  Google Scholar 

  8. Shen H, He H, Li J, Chen W, Wang X, Guo L, Peng Z, He G, Zhong S, Qi Y, et al. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. 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 USA. 2015;112(46):E6397–406.

  10. Greaves IK, Gonzalez-Bayon R, Wang L, Zhu AY, Liu PC, Groszmann M, Peacock WJ, Dennis ES. Epigenetic changes in hybrids. Plant Physiol. 2015;168(4):1197–205.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Ha M, Lu J, Tian L, Ramachandran V, Kasschau KD, Chapman EJ, Carrington JC, Chen X, Wang X-J, Chen ZJ. Small RNAs serve as a genetic buffer against genomic shock in Arabidopsis interspecific hybrids and allopolyploids. Proc Natl Acad Sci. 2009;106(42):17835–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Greaves IK, Groszmann M, Ying H, Taylor JM, Peacock WJ, Dennis ES. Trans chromosomal methylation in Arabidopsis hybrids. Proc Natl Acad Sci U S A. 2012;109(9):3570–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Roudier F, Teixeira FK, Colot V. Chromatin indexing in Arabidopsis: an epigenomic tale of tails and more. Trends Genet. 2009;25(11):511–7.

    Article  CAS  PubMed  Google Scholar 

  15. He GM, Zhu XP, Elling AA, Chen LB, Wang XF, Guo L, Liang MZ, He H, Zhang HY, Chen FF, et al. Global epigenetic and transcriptional trends among Two rice subspecies and their reciprocal hybrids. Plant Cell. 2010;22(1):17–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. He G, Chen B, Wang X, Li X, Li J, He H, Yang M, Lu L, Qi Y, Wang X, et al. Conservation and divergence of transcriptomic and epigenomic variation in maize hybrids. Genome Biol. 2013;14(6):R57.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Moghaddam AM, Roudier F, Seifert M, Berard C, Magniette ML, Ashtiyani RK, Houben A, Colot V, Mette MF. Additive inheritance of histone modifications in Arabidopsis thaliana intra-specific hybrids. Plant J. 2011;67(4):691–700.

    Article  CAS  PubMed  Google Scholar 

  18. Dong X, Reimer J, Gobel U, Engelhorn J, He F, Schoof H, Turck F: Natural variation of H3K27me3 distribution between two Arabidopsis accessions and its association with flanking transposable elements. Genome Biology. 2012;13(12):R117.

  19. Ni ZF, Kim ED, Ha MS, Lackey E, Liu JX, Zhang YR, Sun QX, Chen ZJ. Altered circadian rhythms regulate growth vigour in hybrids and allopolyploids. Nature. 2009;457(7227):327–U327.

    Article  CAS  PubMed  Google Scholar 

  20. Meyer RC, Kusterer B, Lisec J, Steinfath M, Becher M, Scharr H, Melchinger AE, Selbig J, Schurr U, Willmitzer L, et al. QTL analysis of early stage heterosis for biomass in Arabidopsis. Theor Appl Genet. 2010;120(2):227–37.

    Article  PubMed  Google Scholar 

  21. Zhu A, Greaves IK, Liu PC, Wu L, Dennis ES, Peacock WJ: Early changes of gene activity in developing seedlings of Arabidopsis hybrids relative to parents may contribute to hybrid vigour. Plant J. 2016;88(4):597–607.

  22. Zhang XY, 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):1026–35.

    CAS  Google Scholar 

  23. Ha M, Ng DWK, Li WH, Chen ZJ. Coordinated histone modifications are associated with gene expression variation within and between species. Genome Res. 2011;21(4):590–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Roudier F, Ahmed I, Berard C, Sarazin A, Mary-Huard T, Cortijo S, Bouyer D, Caillieux E, Duvernois-Berthet E, Al-Shikhley L, et al. Integrative epigenomic mapping defines four main chromatin states in Arabidopsis. Embo J. 2011;30(10):1928–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Buzas DM, Robertson M, Finnegan EJ, Helliwell CA. Transcription-dependence of histone H3 lysine 27 trimethylation at the Arabidopsis polycomb target gene FLC. Plant J. 2011;65(6):872–81.

    Article  CAS  PubMed  Google Scholar 

  26. Bouyer D, Roudier F, Heese M, Andersen ED, Gey D, Nowack MK, Goodrich J, Renou JP, Grini PE, Colot V et al. Polycomb Repressive Complex 2 Controls the Embryo-to-Seedling Phase Transition. Plos Genet. 2011;7(3):e1002014.

  27. Bernatavichute YV, Zhang X, Cokus S, Pellegrini M, Jacobsen SE. Genome-wide association of histone H3 lysine nine methylation with CHG DNA methylation in Arabidopsis thaliana. Plos One. 2008;3(9):e3156.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Deleris A, Stroud H, Bernatavichute Y, Johnson E, Klein G, Schubert D, Jacobsen SE. Loss of the DNA Methyltransferase MET1 Induces H3K9 Hypermethylation at PcG Target Genes and Redistribution of H3K27 Trimethylation to Transposons in Arabidopsis thaliana. Plos Genet.2012;8(11):e1003062. 

  29. Stroud H, Do T, Du JM, Zhong XH, Feng SH, Johnson L, Patel DJ, Jacobsen SE. Non-CG methylation patterns shape the epigenetic landscape in Arabidopsis. Nat Struct Mol Biol. 2014;21(1):64–72.

    Article  CAS  PubMed  Google Scholar 

  30. Sequeira-Mendes J, Araguez I, Peiro R, Mendez-Giraldez R, Zhang XY, Jacobsen SE, Bastolla U, Gutierrez C. The functional topography of the Arabidopsis genome is organized in a reduced number of linear motifs of chromatin states. Plant Cell. 2014;26(6):2351–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Moghaddam AMB, Fuchs J, Czauderna T, Houben A, Mette MF. Intraspecific hybrids of Arabidopsis thaliana revealed no gross alterations in endopolyploidy, DNA methylation, histone modifications and transcript levels. Theor Appl Genet. 2010;120(2):215–26.

    Article  Google Scholar 

  32. Perry SE, Wang H. Rapid isolation of Arabidopsis thaliana developing embryos. BioTechniques. 2003;35(2):278–80. 282.

    CAS  PubMed  Google Scholar 

  33. Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nussbaum C, Myers RM, Brown M, Li W et al. Model-based Analysis ofChIP-Seq (MACS). Genome Biology. 2008;9(9):R137.

  34. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Shen L, Shao NY, Liu XC, Nestler E. ngs.plot: Quick mining and visualization of next-generation sequencing data by integrating genomic databases.BMC Genomics. 2014;15:284.

Download references

Acknowledgements

The authors wish to thank Ming-Bo Wang and Jean Finnegan for suggestions on the project and for reviewing the manuscript.

Funding

Not applicable.

Availability of data and materials

Sequence data from this article can be found in the GenBank data libraries under accession number SRX2024391.

Authors’ contributions

AZ, IKG, WJP, and ESD. designed research; AZ performed experiments; IKG supervised data analysis; AZ analysed data; and AZ, WJP, and ESD wrote the paper. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information

Authors and Affiliations

  1. Commonwealth Scientific and Industrial Research Organization, Agriculture and Food, Canberra, Australian Capital Territory, 2601, Australia

    Anyu Zhu, Ian K. Greaves, Elizabeth S. Dennis & W. James Peacock

  2. Faculty of Science, University of Technology, Sydney, New South Wales, 2007, Australia

    Anyu Zhu, Elizabeth S. Dennis & W. James Peacock

Authors
  1. Anyu Zhu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Ian K. Greaves
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Elizabeth S. Dennis
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. W. James Peacock
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to W. James Peacock.

Additional files

Additional file 1: Figure. S1.

The distributions of the four histone marks over the five chromosomes in Ler and C24 and their hybrids. Figure. S2. The correlation patterns between histone marks and between histone marks and genes in C24. Figure S3. Correlations between histone modification and gene expression in C24. Figure S4. Venn-diagram showing the overlaps between the differentially modified genes by K4me3, K9ac and K9me2, in the parents and hybrids. Figure S5. Correlations between changes in histone modifications and changes in gene expression between Ler and C24 in germinating seeds. Figure S6. Correlations between changes in histone modifications and changes in gene expression in the hybrids in germinating seeds. Figure S7. Proportions of genes with allelic ratios of histone modifications in the hybrids different from those between the parents. Figure S8. Box plots showing the allelic modification levels relative to the expected levels in the reciprocal hybrids at the DMGs in the parents. Figure S9. Correlations between gene expression and histone modification at Ler and C24 alleles in the reciprocal hybrids. Figure S10. Box plots showing allelic expression and modification levels in the reciprocal hybrids at the genes having corresponding changes in expression and histone modification between the parents. Figure S11. Examples of genes showing corresponding allelic changes in gene expression and modifications. Figure S12. Relative levels of four histone marks at ACTIN 7, AGAMOUS and TA3 using ChIP libraries for deep sequencing. (PDF 3608 kb)

Additional file 2: Data S1.

Number of sequenced reads mapped to reference genome. Data S2. List of the genes with non-additive expression and consistent non-additive histone modifications in LxC and CxL at germination stage. (XLSX 5193 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhu, A., Greaves, I.K., Dennis, E.S. et al. Genome-wide analyses of four major histone modifications in Arabidopsis hybrids at the germinating seed stage. BMC Genomics 18, 137 (2017). https://doi.org/10.1186/s12864-017-3542-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-017-3542-8

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us