- Research article
- Open Access
Genome-wide analysis of DNA methylation in photoperiod- and thermo-sensitive male sterile rice Peiai 64S
© Hu et al.; licensee BioMed Central. 2015
- Received: 19 October 2014
- Accepted: 3 February 2015
- Published: 19 February 2015
Epigenetic modifications play important roles in the regulation of plant development. DNA methylation is an important epigenetic modification that dynamically regulates gene expression during developmental processes. However, little studies have been reported about the methylation profiles of photoperiod- and thermo-sensitive genic male sterile (PTGMS) rice during the fertility transition.
In this study, using methylated DNA immunoprecipitation sequencing (MeDIP-seq), the global DNA methylation patterns were compared in the rice PTGMS line PA64S under two different environments (different temperatures and day lengths). The profiling of the DNA methylation under two different phenotypes (sterility and fertility) revealed that hypermethylation was observed in PA64S (sterility), and 1258 differentially methylated regions (DMRs) were found between PA64S (sterility) and PA64S (fertility). Twenty differentially methylated genes of them were further validated through bisulfite sequencing, and four of these genes were analyzed by qRT-PCR. Especially, a differentially methylated gene (LOC_Os08g38210), which encoded transcription factor BIM2, is a component of brassinosteroid signaling in rice. The hypermethylated BIM2 gene may suppress some downstream genes in brassinosteroid signaling pathway, and thus affect the male fertility in PA64S.
The results presented here indicated that hypermethylation was observed in PA64S (sterility). Gene Ontology (GO) analysis and KEGG analysis revealed that flavone and flavonol biosynthrsis, circadian rhythm, photosynthesis and oxidative phosphorylation pathways were involved in sterility-fertility transition of PA64S.
- Oryza sativa
- DNA methylation
- Fertility transition
- Gene expression
DNA methylation is an indispensable epigenetic modification of plant genomes. It occurs predominantly in the CG context in both plants and animals [1-3]. The state of CpG methylation regulates and stabilizes chromatin structure, and possibly regulates the accessibility of these DNA regions to the transcriptional machinery . Studies have shown that DNA methylation plays an important role in many plant processes, including transposon silencing, virus defence, and gene imprinting [5-7].
In plant, DNA methylation is mainly distributed in heterochromatic regions that are composed of tandem or inverted repeats, i.e., transposons. Transposable elements (TEs) are usually highly methylated along their entire length, but DNA methylation within genes is distributed away from the 5′ and 3′ ends. It has been proposed that DNA methylation is primarily used to inhibit the transcription and movement of TEs to protect genomes against uncontrolled insertions. Some studies have reported that transposons or retrotransposons are activated in response to environmental stress and may generate true genetic or epigenetic changes for adapting to stresses . In rice, epigenetic regulation has been reported to affect the transposition activity of Tos17 and modulates the activity of the neighboring genes . Highly methylated transposons may affect transcription of neighboring genes. For example, the patterns of DNA methylation in long terminal repeat (LTR) transposable elements differed between rice leaves and roots and affected the transcription of the neighboring genes .
Recently, three DNA treatment techniques have been typically employed for DNA methylation detection: (i) methylation sensitive amplified polymorphism (MSAP), (ii) methylated DNA enrichment (MeDIP), (iii) bisulfite-sequencing. With the advent of new technologies, the elucidation of the genome-wide methylation profiles is possible. In this regard, the immunoprecipitation of methylated DNA by monoclonal antibodies specific to 5-methylcytidine (5mC) combined with high-throughput sequencing (MeDIP-seq) has been used as a valuable tool to map methylated DNA at the genomic scale [11,12]. The MeDIP-seq approach employs an antibody against 5-methylcytosine or methyl-binding domain proteins to capture methylated DNA, which are subsequently subjected to next generation sequencing . After methylated DNA enrichment, the unmethylated DNA fragments are removed from the genome samples for reducing the redundancy. Thus, a relatively lower sequencing throughput is required, compared with bisulfite-sequencing. In addition, MeDIP is more unbiased and less limited than MSAP. Due to these advantages, MeDIP-seq is considered to have great potential for the development of cost-effective and unbiased strategies for whole genome DNA methylation profiling.
Rice (Oryza sativa L.) is one of the most important grain crops worldwide and provides a staple food for almost half of the world’s population . Recently, hybrid rice has made a tremendous contribution to food security both in China and many other countries which are losing arable land. Hybrid rice technologies are mainly based on two well-known male sterility systems, namely CMS (Cytoplasmic Male Sterile) and EGMS (Environmentally sensitive Genic Male Sterile). Compared with the three-line system (CMS system), the two-line system, is based on the discovery and application of EGMS lines, which serve as both the male sterile lines and maintainer lines under different environmental conditions . PTGMS (Photoperiod- Thermo-Sensitive Genic Male Sterile) is the major type of EGMS germplasm resources and has been widely used for breeding two-line hybrid rice. The PTGMS system is advantageous for broad restoration ability, easy maintenance, and multiplication [15-17]. Thus, the utilization of PTGMS lines is an important approach for the better exploitation of heterosis in rice.
Peiai 64S (PA64S) is one of the most important indica rice genic male sterile lines in two-line system. Its paternal line is Peiai64, and its maternal line is Nongken58, which is the first photoperiod- sensitive genic male sterile line discovered by Mingsong Shi in Hubei, China in 1973 . Using PA64S as the maternal line, more than 10 two-line hybrid rice varieties have been developed since 1996 in China. Among these new varieties, LiangYouPei9 (LYP9) is an elite super hybrid rice and increases grain yield by approximately 20% per hectare compared with other hybrid rice lines in the past years [7,19]. The fertility transition of PA64S was controlled by the photoperiod and temperature: the sterility trait exhibited at temperatures higher than 23.5°C during the anther development and long-day (LD) (≥14 h) conditions can suppress the fertility transition. Therefore, under lower temperatures (~21-23°C) and short-day (SD) (<14 h) conditions during anther development, PA64S can transform from male sterility to fertility [20,21].
To date, some studies have reported that a substitution of C-to-G in the male-fertility-associated noncoding RNA p/tms12-1 in PA64S produced a mutant small RNA, namely osa-smR5864m. This mutated noncoding small RNA gene may lead to PTGMS in PA64S . Furthermore, the transcriptome of the PGMS rice Nongken 58S was significantly suppressed under LD condition and the circadian rhythm pathway had been shown to be involved in the male sterility transition . Another noncoding RNA called LDMAR was also found in Nongken 58S and was required for the fertility transition. A siRNA (Psi–LDMAR) was detected to associate with the DNA methylation pattern of the LDMAR promoter in Nongken 58S. And increasing methylation in the promoter of LDMAR of Nongken 58S also reduced the expression of LDMAR, leading to male sterility . These results suggested that DNA methylation or RNA-dependent DNA methylation (RdDM) might be involved in the regulation of PTGMS.
Brassinosteroids (BRs), a class of steroid hormones, have been reported to control male fertility by regulating the expression of several key genes involved in anther and pollen development, such as SPL/NZZ, TDF1, AMS, AtMYB103 and MS1 [24,25]. The basic helix-loop-helix (bHLH) protein BIM1 (BES1-interacting Myc-like1) is a brassinosteroid signaling component involved in regulating BR-induced genes in Arabidopsis . BRI1- EMSSUPRESSOR1 (BES1) is shown to interact with BIM1 and together bind to the E-box of a BR-induced gene .
In this study, we investigated the global DNA methylation alterations in young panicles of the PTGMS line PA64S under two different environmental conditions using the MedIP-seq method. The aim of this work was to explore the DNA methylation patterns and their influences on gene expression during the transition from sterile to fertile in PA64S. The results will provide some information for a better understanding of the role of DNA methylation in PTGMS rice.
Cytological observation of pollen morphology
Genome-wide DNA methylation profiles of PA64S (S) and PA64S (F)
Data generated by MeDIP-seq in PA64S (S) and PA64S (F)
Total MeDIP-Seq data
Percentage of mapped reads in total reads
Unique mapped reads
Unique mapped bases
Percentage of unique mapped reads
To get insight into the DNA methylation status of PA64S (S) and PA64S (F), RPKMs (Reads per 1 kb window per million mapped reads) were calculated for each sample and different elements (Figure 4). The results showed that PA64S(S) had higher methylation level than that of PA64S (F). And the DNA methylation levels between PA64S(S) and PA64S (F) were observed to be significantly different on the upstream 2k, 5’UTR, intron, 3’UTR, and downstream 2k regions of the gene body (p < 0.01, p < 0.001). However, there were not significant differences on CDS (coding DNA sequence) of the differential methylated regions (DMRs).
DNA methylation patterns in genes and TE regions
Numbers of differentially methylated genes between PA64S(S) and PA64S (F) in different gene regions
PA64S(S) VS. PA64S(F)- Hyper
PA64S(S) VS. PA64S(F)- Hypo
Significant GO enrichment of hypermethylated genes in PA64S (S)
Upstream 2k/ Intron
NB-ARC domain containing protein
Early nodulin 75-like protein
zinc finger DHHC domain-containing protein
phosphate carrier protein, mitochondrial precursor
Amino acid permease family protein
pumilio-family RNA binding repeat containing protein
transposon protein, putative, CACTA, En/Spm sub-class
Probable calcium-binding protein CML16
transposon protein, putative, unclassified
Long chain base biosynthesis protein 2b
transferase family protein
pumilio-family RNA binding repeat containing protein
5'-UTR/ CDS/ Intron
MATE efflux family protein
prenylated rab acceptor
transposon protein, CACTA, En/Spm sub-class
transposon protein, CACTA, En/Spm sub-class
leucine-rich repeat family protein
Leucine Rich Repeat family protein
impaired sucrose induction 1
retrotransposon protein, Ty1-copia subclass
retrotransposon protein, Ty3-gypsy subclass
PE-PGRS family protein precursor
inorganic phosphate transporter
telomerase reverse transcriptase
FAM10 family protein
impaired sucrose induction 1
amino acid permease
DHHC zinc finger domain containing protein
exostosin family domain containing protein
transposon protein, Pong sub-class
homeobox domain containing protein
auxin response factor
cysteine proteinase EP-B 1 precursor
DEK C terminal domain containing protein
protein phosphatase inhibitor 2 containing protein
Rapid ALkalinization Factor RALF family protein
ent-kaurene synthase, chloroplast precursor
transmembrane amino acid transporter protein
helix-loop-helix DNA-binding domain containing protein
calcium-binding mitochondrial carrier
serine palmitoyltransferase 2
bile acid sodium symporter family protein
aspartic proteinase oryzasin-1 precursor
retrotransposon protein, Ty3-gypsy subclass
Transposable element methylation of PA64S(S) and PA64S (F) based on 10 kb tiled genome windows
Validation of MedIP data by bisulfite sequencing
Association of DNA methylation and gene expression
Analysis of the relationship between DNA methylation and gene expression revealed that hypermethylation in some regions were related to the expression of neighboring genes. DNA methylation level of PA64S(S) on the upstream 2k, 5’UTR, intron, 3’UTR, and downstream 2k regions were hypermethylated compared with PA64S (F) (p < 0.001) (Figure 4b and Table 2). Correspondingly, many DMR-associated genes may be down-regulated in PA64S (S).
The DNA methylation patterns and gene expression of four DMR-associated genes, namely LOC_Os03g51200 (H2A), LOC_Os08g38210 (transcription factor BIM2), LOC_Os06g16270 (heat shock factor binding protein 2) and LOC_Os03g27770 (heme oxygenase 2), were further examined in this study (Figure 6). Three of these had lower DNA methylation levels, but higher expression in PA64S (S). Furthermore, The DMR-associated gene LOC_Os08g38210 showed a higher DNA methylation level and a lower expression in PA64S (S) (Figure 6d).
Expression of the genes involved in brassinosteroid signaling pathway
Gene ontology (GO) and KEGG pathway enrichment analysis
In this study, genes that overlapped with the methylation peaks in the upstream 2k, downstream 2k or gene body regions were termed as methylated genes. A total of 36,321 methylated genes were found in the two samples: 32,678 in PA64S (S) and 32,990 in PA64S (F). Of them, 1258 methylated genes were identified to be differentially expressed genes between PA64S (S) and PA64S (F), and these included 725 genes that were hypermethylated and 533 genes that were hypomethylated in PA64S (S) (Table 2).
Overrepresented Gene Ontology (GO) categories in the hypermethylated genes in PA64S (S)
Number in input list
Number in BG/Ref
mitochondrial respiratory chain complex III
mitochondrial respiratory chain
mitochondrial inner membrane
hydroxymethyl-, formyl- and related transferase activity
transferase activity, transferring one-carbon groups
transcription regulator activity
electron transporter, transferring electrons within CoQH2-cytochrome c reductase complex activity
tetrahydrofolate metabolic process
folic acid and derivative metabolic process
pteridine and derivative metabolic process
cellular metabolic compound salvage
coenzyme metabolic process
ribonucleotide metabolic process
cofactor metabolic process
nucleotide metabolic process
cellular aromatic compound metabolic process
monocarboxylic acid metabolic process
heterocycle metabolic process
mitochondrial electron transport, ubiquinol to cytochrome c
mitochondrial ATP synthesis coupled electron transport
ATP synthesis coupled electron transport
To generate further insight view of pathway perturbation in PA64S (S) under long day length and high temperature conditions, we performed KEGG pathway analysis of the differentially methylated genes. As a result, 58 predicted pathways were enriched, including metabolism, circadian rhythm, photosynthesis, porphyrin and chlorophyll metabolism, peroxisome and oxidative phosphorylation (Additional files 6). Through the GO term enrichment and KEGG pathway analysis of the differentially regulated genes, we observed that certain genes were hypermethylated in PA64S (S) and that they were involved in photorespiration and mitochondrial ATP synthesis pathways. These pathways were closely connected to photosynthesis and energy synthesis. For instance, cytochrome c1-1 (LOC_Os05g23620) was involved in mitochondrial electron transport.
Although global DNA methylation surveys have been performed in some plants [3,31,32], little studies have been reported about the methylation profiles of PTGMS rice (for example PA64S) during the fertility transition. In this study, the objective was to perform a genome-wide identification of the methylated genes that affect the sterility-fertility transition in PTGMS rice PA64S. To confirm the results from MeDIP-seq, the methylation patterns of twenty regions in each sample were analyzed by bisulfite sequencing. And the methylation levels obtained by the two methods were generally in accord with each other. The read distribution analysis performed in this study found that uniquely mapped reads were enriched in the upstream 2k, 5’UTR, 3’UTR and downstream 2k regions. In addition, the rice sterile line PA64S(S) had higher methylation level than that of PA64S (F), and many DMR-associated genes were down-regulated in PA64S (S).
Transposable elements (TEs) are abundant features in plant genomes that are frequently marked by cytosine methylation. In this study, we examined many categories of transposable elements from the RepBase transposable element annotation database. LTR retroelements are the predominant type of transposable element sequences in plant genomes . Our results showed that the LTR categories had higher 5mC content in PA64S (F) than in PA64S (S). And the LINE and SINE retroelement categories had an increased in 5mC content in PA64S (S) compared with PAS64 (F), although this difference was not significant. Previous studies using methylation-deficient mutants have demonstrated that hypomethylation can result in transposable element mobilization . The increased 5mC content of TEs found in the present study suggested that the operation of adaptive genomic mechanisms might help protect genomes against transposable element spread. It would be of interest to measure the relationship of DNA methylation and TE elements during the fertility transition of PA64S.
In this study, the DNA methylation patterns of twenty DMR-associated genes were validated through bisulfite sequencing and we found that hypermethylation in some regions affected the expression of neighboring genes (Figures 5 and 6, Additional file 3). Although hypermethylated DMR-associated genes did not show significant enrichment of down-regulated genes, there were still some down-regulated genes with higher methylation levels (Figures 5 and 6). Some flower/anther development related genes were hypermethylated in PA64S (S), for example, LOC_Os06g23980 (putative MADS-box protein AGL16) and LOC_Os10g10560 (IPMI) (Additional file 3). Previous studies of a male sterile mutant (Ms-cd1) in Brassica oleracea have demonstrated that IPMI was specifically expressed in anther and filament and was not expressed in CMS NILs . The increased 5mC content in IPMI of PA64S(S) found in the present study might lead to a lower expression of IPMT and most likely disturbed the anther cell wall metabolism. Furthermore, photosystem and mitochondrial ATP synthesis associated genes, including LOC_Os12g08770 (photosystem I reaction center subunit N), LOC_Os07g17390 (PsbP) and LOC_Os03g01820 (vacuolar ATP synthase subunit E) were also hypermethyalted in PA64S (S) (Additional file 2). High temperature (>23.5°C) and LD (14 h) conditions increased the methylation levels of these genes and down-regulated their expression during anther development in PA64S, which might in turn lead to PA64S sterility.
In contrast, some heat shock proteins, such as LOC_Os08g43490 (heat shock protein DnaJ), LOC_Os06g16270 (heat shock factor binding protein 2) were found to be hypomethylated in PA64S (S) in our study, and the expression level of LOC_Os06g16270 was also higher in PA64S (S) (Figure 6c). Because heat shock proteins are always expressed in response to stress, their higher expression in PA64S (S) suggested that the adverse environment might affect anther development. These findings suggested that different environments might regulate the gene expression through DNA methylation in the gene body or flanking sequence of fertility associated genes.
In Arabidopsis and rice, the BR signaling pathway is similar. Generally, BRs bind the extracellular domain of BRI1, and then BRI1 heteromerizes with the BRI1-associated kinase 1 (BAK1). The receptor complex then activates the BR-signaling kinases (BSKs) and inactivates the glycogen synthase kinase-3-like kinase BIN2, resulting in nuclear accumulation of unphosphorylated forms of two transcription factors BZR1 and BES1, which activate the expression of BR-responsive genes [36-38]. In this study, some of the components in BR signaling pathway, such as OsBRI, OsBAK1, OsBZR and OsBIM2 were dramatically reduced in male sterile line PA64S (S) (Figure 8). Our results indicated that BR signaling involved in regulating the sterility-fertility transition of PA64S under different environments.
On the other hand, OsBIM2 (LOC_Os08g38210) showed a higher DNA methylation level in PA64S(S) and its expression dramatically reduced (Figure 6d). In Arabidopsis, BIM1 and SPL8 act cooperatively in a common complex or via synergistic pathways to promote male fertility . In our study, both OsBIM2 and OsSPL8 were downregulated in PA64S(S) (Figure 8). Since OsBIM2 is homolog with BIM1, the lower expression level will affect BR signaling pathway and then regulated the downstream genes which involved in anther and pollen development. Therefore, DNA methylation of BR signaling related genes may affect the expression of key genes (for example, OsMS2) for anther development, leading to male sterility under long-day and high temperature conditions. However, whether the hypermethylated OsBIM2 directly suppress key genes need to be further investigated.
The GO term significant enrichment analysis showed that many genes, including the transposons and retrotransposons, were hypermethylated in PA64S (S) (Tables 3 and 5 and Additional files 3 and 4). For example, LOC_Os09g38100 (phosphate carrier protein mitochondrial precursor) and LOC_Os06g40200 (calcium-binding mitochondrial carrier) were observed higher methylation levels in PA64S (S) than PA64S (F) (Table 3 and Additional file 3). These two genes participate in mitochondrial respiratory chain and facilitate mitochondrial ATP synthesis . Thus, their hypermethylated would affect the energy synthesis during pollen development. Another hypermethylated gene LOC_Os02g17780 (ent-kaurene synthase, chloroplast precursor) was an early intermediate in GA biosynthesis . Other significantly differentially expressed DMR-associated genes included LOC_Os12g40130 (phosphoribosyl transferase), LOC_Os04g30490 (MATE efflux family protein) and LOC_Os05g12210 (chalcone synthase). Some studies have reported that the adenine phosphoribosyl transferase activity deficient mutants of Arabidopsis thaliana were male sterile due to abortion of pollen development after the meiotic divisions of pollen mother cells . The MATE (Multidrug And Toxic Compound Extrusion) efflux protein family has also been reported to be differentially expressed in sterile and fertile NILs of Brassica oleracea . Chalcone synthase catalyzes the initial step of the branch of phenylpropanoid pathway which leads to flavonoids. The absence of chalcone synthase activity was found to have a pleiotropic effect in maize and petunia mutants, and pollen fertility/pollen germination and flavonoid synthesis was disrupted . Therefore, these differentially expressed genes might be involved in fertility transition.
To uncover the regulatory mechanisms of the fertility transition, the regulatory network of differentially methylated genes was examined in this study. Several important pathways were found, including flavone and flavonol biosynthesis, circadian rhythm, photosynthesis, porphyrin and chlorophyll metabolism, and oxidative phosphorylation (Additional files 6). Previous studies have reported that flavone and flavonol biosyntheses are associated with pollen development. Flavonoids are required for pollen viability in many species [42-44]. In addition, flavonols can promote pollen germination frequency and pollen tube growth in tobacco . Many higher plants use the day length as an environmental cue to switch from vegetative to reproductive growth and plants with defects in their circadian regulation cannot properly regulate the timing of the floral transition . In PGMS rice, the circadian rhythm genes OsPRR37 in inflorescences have been reported to directly affect the male sterility transition in Nongken 58S . In this study, some DMR-associated genes, such as LOC_Os05g12210 (chalcone synthase) and LOC_Os08g38210 (transcription factor BIM2), were found to be involved in the circadian rhythm. The other three pathways, namely photosynthesis, porphyrin and chlorophyll metabolism, and oxidative phosphorylation were related to carbon or energy metabolism [47,48]. The down-regulated genes of these pathways in PA64S (S) may lead to insufficient energy and materials for pollen development. Therefore, these five pathways were regarded as potentially related to the fertility transition of PA64S in the present study.
We systematically investigated the complete DNA methylome of photoperiod- and thermo-sensitive male sterile rice Peiai 64S under two different environments (different temperatures and day lengths). From the whole-genome DNA metylation map, hypermethylation was observed in PA64S (S) and many pathways such as circadian rhythm, photosynthesis and oxidative phosphorylation were participated in the fertility transition of PA64S. Furthermore, the hypermethylated OsBIM2 may be involved in brassinosteroid signaling pathway and affected the male fertility in PA64S.
The rice photo-thermo-sensitive genic male sterile line PA64S, which was maintained at Wuhan University, was used in this study. The seeds of PA64S were sown under two different natural ecological conditions: in Hainan (18°48′ N, 110°02′ E) from December, 2010 to March, 2011 (fertility, F), and in Wuhan (30°30′ N, 114°18′ E) from May 2011 to August 2011 (sterility, S). Young panicles of PA64S of the two phenotypes at the meiosis stage were collected, and the pollen fertility was assayed by staining with potassium iodide (1% I2-KI).
DNA preparation and methylated DNA immunoprecipitation
Total genomic DNA was extracted using a Genomic DNA Miniprep Kit (Axygen) following the manufacturer’s instructions. The DNA quality was then evaluated by agarose gel electrophoresis and a BioPhotometer Plus spectrophotometer (Eppendorf, Germany).
Subsequently, the two samples were sonicated to produce DNA fragments ranging from 100 to 500 bp. After end repairing, phosphorylating and A-tailing with Paired-End DNA Sample Prep kit (Illumina, USA), the DNA was ligated to an Illumina sequencing primer adaptor. Then the fragments were used for methylated DNA immunoprecipitation (MeDIP) enrichment using a Magnetic Methylated DNA Immunoprecipitation kit (Diagenod, Belgium) following the manufacturer’s recommendations and the qualifying DNA was used for PCR amplification. Then bands between 220 and 320 bp were excised from the gel and purified with a QIAquick Gel Extraction Kit (Qiagen, Germany). The products were quantified with a Quant-iT™ dsDNA HS Assay Kit (Invitrogen, USA) using an Agilent 2100 Analyzer (Agilent Technologies, USA). Following qPCR qualification, DNA libraries were sequenced on the Illumina HiSeq 2000 (Illumina, CA, USA) to generate paired-end 50-bp reads by the Beijing Genomics Institute (BGI, China).
The raw data obtained from Illumina sequencing were first processed to filter out reads containing adapters, unknown or low quality bases and were then were mapped to the rice reference genome (ftp://ftp.plantbiology.msu.edu/pub/data/Eukaryotic_Projects/o_sativa) by SOAPaligner v 2.21 (http://soap.genomics.org.cn/) with no more than 2 bp mismatches . The uniquely mapped data were retained for read distribution analysis, including the distribution in rice chromosomes and the distribution in different components of the genome.
Subsequently, the regions of differential methylation (DMRs) with DNA methylation peaks were employed for differential DNA methylation analysis. Methylated regions were deemed significant differentially methylated between PA64S (S) and PA64S (F) with a p-value < 0.05, a false discovery rate (FDR) < 0.05 and at least a 2.0-fold change in sequence counts. The base data manipulations and statistical analysis were performed using the R package (http://www.r-project.org). Significant differences of DNA methylation level between PA64S (S) and PA64S (F) were evaluated using Student’s t test.
Gene Ontology and KEGG pathway analysis
where N is the number of all genes with a GO annotation in the rice genome, n is the number of those N genes that are differentially expressed, M is the number of all of the genes that are annotated to certain GO terms, and m is the number of differentially expressed genes in M. The hypergeometric test was also performed using a user written program in R (R Core Team, Vienna, Austria). The KEGG pathway enrichment analysis identifies significantly enriched metabolic pathways or signal transduction pathways in differentially expressed genes comparing with the whole genome background. The calculating formula and the applied program were the same as that used in the GO analysis. Then, FDR (Q-value) ≤0.05 was used as the threshold to determine the most significantly enriched pathways in DEGs.
Genomic DNA was modified using an EpiTect Bisulfite kit (Qiagen, Germany) according to the manufacturer’s instructions. Briefly, 1 μg of DNA was first treated with bisulfite containing a C-T conversion reagent and then incubated at 95°C for 5 min, and 64°C for 5 h. The modified DNA was purified using a MinElute DNA spin columns (Qiagen, Germany) and stored at -20°C until use. Primers for selected genic sequences for bisulfite sequencing were designed by the Meth-Primer Program (http://www.urogene.org/methprimer/) and were given in Additional file 7. For each PCR reaction, 1.0 μl of bisulfite treated DNA was used in a 50 μl reaction system. These PCR products were gel-purified using a Gel Extraction Kit (Omega, USA), and then cloned into the pMD18-T vector (Takara, Dalian, China) and sequenced. At least fifteen clones were sequenced for each sample. The methylation levels which are expressed as percentage (%) per site for each of the three types of cytosines, CG, CHG and CHH, were calculated by dividing the number of non-converted (methylated) cytosines by the total number of cytosines within the assay. The analyses of the bisulfite sequencing results were conducted using the Kismeth website (http://katahdin.mssm.edu/kismeth).
RNA isolation and real time RT-PCR
Total RNA was extracted from young panicles using an RNA isolation kit (BioDev, Beijing, China) and treated with DNase I (Fermentas, Canada) to remove the genomic DNA contamination. Approximately 2 μg of RNA from each sample was used for reverse transcription using the RevertAid First Strand cDNA Synthesis Kit (Fermentas, Canada). The reaction was carried out in a total volume of 20 μl, including 2 μg of RNA, 1 μl of oligo (dN)18 primers (100 μmol/L), 4 μl 5× reaction buffer, 2 μl 10 mmol/L dNTP mix, 1 μl of RiboLock RNase Inhibitor (20 U/μl), and 1 μl of M-MuLV Reverse Transcriptase (Fermentas, Canada). The procedure was as follows: the mixture was first denatured at 70°C for 10 min and then incubated at 37°C for 5 min, followed by 42 °C for 60 min and at last incubated at 70°C for 15 min.
Quantitative real-time PCR (qPCR) was performed using a StepOne Plus real-time PCR system (Applied Biosystems) in a total volume of 10 μl reaction mixture, which included 1 μl cDNA, 0.4 μl each primer (5 μmol/L each), 0.2 μL ROX Reference Dye (50×) and SYBR premix Ex Taq (Tli RNaseH plus) (2×) (Takara, Dalian, China). The primers used for the qPCR experiments are given in Additional file 7. The qPCR reactions were repeated three times for each sample from three independent experiments. After the amplification steps, the melting curve was determined for each primer pair to verify that only one specific product had been amplified. Relative quantification was performed using the 2△△Ct method . Rice actin1 was used as a reference for the mRNA expression analyses.
Availability of supporting data
All MeDIP-seq raw data supporting the results of this article were deposited in the sequence read archive repository (http://www.ncbi.nlm.nih.gov/sra/) under Sequence Read Archive (SRA) database (accession # SRP050401).
This work was supported by the “973” Program of China (2013CB126900) and the National Nature Science Foundation of China (31471464).
- Feng S, Jacobsen SE, Reik W. Epigenetic reprogramming in plant and animal development. Science. 2010;330:622–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Xiang H, Zhu J, Chen Q, Dai FY, Li X, Li M, et al. Single base-resolution methylome of the silkworm reveals a sparse epigenomic map. Nat Biotechnol. 2010;28(5):516–20.View ArticlePubMedGoogle Scholar
- Zhang X, Yazaki J, Sundaresan A, Cokus S, Chan SWL, Chen HM, et al. Genome wide High-Resolution Mapping and Functional Analysis of DNA Methylation in Arabidopsis. Cell. 2006;126(6):1189–201.View ArticlePubMedGoogle Scholar
- Ruike Y, Imanaka Y, Sato F, Shimizu K, Tsujimoto G. Genome-wide analysis of aberrant methylation in human breast cancer cells using methyl-DNA immune- precipitation combined with high-throughput sequencing. BMC genomics. 2010;11:137.View ArticlePubMed CentralPubMedGoogle Scholar
- Cokus SJ, Feng S, Zhang X, Chen Z, Merriman B, Haudenschild CD, et al. Shotgun bisulphitesequencing of the Arabidopsis genome reveals DNA methylation patterning. Nature. 2008;452:215–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Lister R, O’Malley RC, Tonti-Filippini J, Gregory BD, Berry CC, Millar AH, et al. Highly integrated single-base resolution maps of the epigenome in Arabidopsis. Cell. 2008;133:523–36.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhang CJ, Chu HJ, Chen GX, Shi DW, Zuo M, Wang J, et al. Photosynthetic and biochemical activities in flag leaves of a newly developed super high-yield hybrid rice (Oryza sativa) and its parents during the reproductive stage. J Plant Res. 2007;120:209–17.View ArticlePubMedGoogle Scholar
- Reinders J, Wulff BB, Mirouze M, Mari-Ordonez A, Dapp M, Rozhon W, et al. Compromised stability of DNA methylation and transposon immobilization inmosaic Arabidopsis epigenomes. Genes Dev. 2009;23:939–50.View ArticlePubMed CentralPubMedGoogle Scholar
- Cheng C, Daigen M, Hirochika H. Epigenetic regulation of the rice retrotransposon Tos17. Mol Genet Genomics. 2006;276:378–90.View ArticlePubMedGoogle Scholar
- Kashkush K, Khasdan V. Large-scale survey of cytosine methylation of retrotransposons and the impact of readout transcription from long terminal repeats on expression of adjacent rice genes. Genetics. 2007;177:1975–85.View ArticlePubMed CentralPubMedGoogle Scholar
- Down TA, Rakyan VK, Turner DJ, Flicek P, Li H, Kulesha E, et al. A Bayesian deconvolution strategy for immunoprecipitation-based DNA methylome analysis. Nat Biotechnol. 2008;26:779–85.View ArticlePubMed CentralPubMedGoogle Scholar
- Tomazou EM, Rakyan VK, Lefebvre G, Andrews R, Ellis P, Jackson DK, et al. Generation of a genomic tiling array of the human major histocompatibility complex (MHC) and its application for DNA methylation analysis. BMC Med Genomics. 2008;1:19.View ArticlePubMed CentralPubMedGoogle Scholar
- Delseny M, Salses J, Cooke R, Sallaud C, Regad F, Lagoda P, et al. Rice genomics: Present and future. Plant Physiol Bioch. 2001;39:323–34.View ArticleGoogle Scholar
- Yuan LP. Purification and production of foundation seed of rice PGMS and TGMS lines. Hybrid Rice. 1994;6:1–3.Google Scholar
- Xu JJ, Wang BH, Wu YH, Du PN, Wang J, Wang M, et al. Fine mapping and candidate gene analysis of ptgms2-1, the photoperiod- thermo- sensitive genic male sterile gene in rice (Oryza sativa L.). Theor Appl Genet. 2011;122:365–72.View ArticlePubMedGoogle Scholar
- Yuan LP. Progress of two-line system hybrid rice breeding. Sci Agric Sin. 1990;23(3):1–6.Google Scholar
- Yang Q, Liang C, Zhuang W, Li J, Deng HB, Deng QY, et al. Characterization and identification of the candidate gene of rice thermo-sensitive genic male sterile gene tms5 by mapping. Planta. 2007;225:321–30.View ArticlePubMedGoogle Scholar
- Shi MS. The discovery and study of the photosensitive recessive male sterile rice. Sci Agric Sin. 1985;2:44–8.Google Scholar
- Wang WW, Meng B, Ge XM, Song SH, Yang Y, Yu XM, et al. Proteomic profiling of rice embryos from a hybrid rice cultivar and its parental lines. Proteomics. 2008;8:4808–21.View ArticlePubMedGoogle Scholar
- Xu ML, Zhou GQ, Chen LB. Response of fertility of Pei’ai 64S to temperature and photoperiod in rice. Sci Agric Sin. 1999;25:772–6.Google Scholar
- Zhou H, Liu QX, Li J, Jiang DG, Zhou LY, Wu P, et al. Photoperiod and thermo-sensitive genic male sterility in rice are caused by a point mutation in a novel noncoding RNA that produces a small RNA. Cell Res. 2012;22:649–60.View ArticlePubMed CentralPubMedGoogle Scholar
- Wang W, Liu ZW, Guo ZB, Song GY, Cheng Q, Jiang DM, et al. 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(1):462.View ArticlePubMed CentralPubMedGoogle Scholar
- Ding JH, Shen JQ, Mao HL, Xie WB, Li XH, Zhang QF. RNA-directed DNA methylation is involved in regulating photoperiod sensitive male sterility in rice. Mol Plant. 2012;5:1210–6.View ArticlePubMedGoogle Scholar
- Kutschera U, Wang ZY. Brassinosteroid action in flowering plants: A darwinian perspective. J Exp Bot. 2012;63:3511–22.View ArticlePubMed CentralPubMedGoogle Scholar
- Ye QQ, Zhu WJ, Li L, Zhang SS, Yin YH, Ma H, et al. Brassinosteroids control male fertility by regulating the expression of key genes involved in Arabidopsis anther and pollen development. Proc Natl Acad Sci USA. 2010;107:6100–5.View ArticlePubMed CentralPubMedGoogle Scholar
- Yin YH, Vafeados D, Tao Y, Yoshida SG, Asami T, Chory J. A new class of transcription factors mediates brassinosteroid-regulated gene expression in Arabidopsis. Cell. 2005;120:249–59.View ArticlePubMedGoogle Scholar
- Nemhauser JL, Mockler TC, Chory J. Interdependency of brassinosteroid and auxin signaling in Arabidopsis. PLoS Biol. 2004;2(9):e258.View ArticlePubMed CentralPubMedGoogle Scholar
- Elango N, Hunt BG, Goodisman MA, Yi SV. DNA methylation is widespread and associated with differential gene expression in castes of the honeybee, Apis mellifera. Proc Natl Acad Sci USA. 2009;106:11206–11.View ArticlePubMed CentralPubMedGoogle Scholar
- Suzuki MM, Kerr AR, De Sousa D, Bird A. CpG methylation is targeted to transcription units in an invertebrate genome. Genome Res. 2007;17:625–31.View ArticlePubMed CentralPubMedGoogle Scholar
- Xing SP, Quodt V, Chandler J, Hohmann S, Berndtgen R, Huijser P. SPL8 acts together with the brassinosteroid-signaling component BIM1 in controlling Arabidopsis thaliana male fertility. Plants. 2013;2:416–28.View ArticleGoogle Scholar
- Li X, Zhu JD, Hu FY, Ge S, Ye MZ, Xiang H, et al. Single-base resolution maps of cultivated and wild rice methylomes and regulatory roles of DNA methylation in plant gene expression. BMC Genomics. 2012;13:300.View ArticlePubMed CentralPubMedGoogle Scholar
- Vining K, Pomraning KR, Wilhelm LJ, Priest HD, Pellegrini M, Mockler TC, et al. Dynamic DNA cytosine methylation in the Populus trichocarpa genome: tissue-level variation and relationship to gene expression. BMC Genomics. 2012;13:27.View ArticlePubMed CentralPubMedGoogle Scholar
- Kumar A, Bennetzen JL. Plant retrotransposons. Annu Rev Genet. 1999;33:479–532.View ArticlePubMedGoogle Scholar
- Vining K, Pomraning KR, Wilhelm LJ, Ma C, Pellegrini M, Di YM, et al. Methylome reorganization during in vitro dedifferentiation an regeneration of Populus trichocarpa. BMC Plant Biol. 2013;13:92.View ArticlePubMed CentralPubMedGoogle Scholar
- Lou P, Kang J, Zhang G, Bonnema G, Fang Z, Wang X. Transcript profiling of a dominant male sterile mutant (Ms-cd1) in cabbage during flower bud development. Plant Sci. 2007;172:111–9.View ArticleGoogle Scholar
- Clouse SD. Brassinosteroid signal transduction: from receptor kinase activation to transcriptional networks regulating plant development. Plant Cell. 2011;23:1219–30.View ArticlePubMed CentralPubMedGoogle Scholar
- He JX, Gendron JM, Yang YL, Li JM, Wang ZY. The GSK3-like kinase BIN2 phosphorylates and destabilizes BZR1, a positive regulator of the brassinosteroid signaling pathway in Arabidopsis. Proc Natl Acad Sci USA. 2002;99:10185–90.View ArticlePubMed CentralPubMedGoogle Scholar
- Wang XL, Chory J. Brassinosteroids regulate dissociation of BKI1, a negative regulator of BRI1 signaling, from the plasma membrane. Science. 2006;313:1118–22.View ArticlePubMedGoogle Scholar
- Stael S, Rocha AG, Robinson AJ, Kmiecik P, Vothknecht UC, Teige M. Arabidopsis calcium-binding mitochondrial carrier proteins as potential facilitators of mitochondrial ATP-import and plastid SAM-import. FEBS Lett. 2011;585:3935–40.View ArticlePubMedGoogle Scholar
- Yamaguchi S, Sun TP, Kawaide H, Kamiya Y. The GA2 Locus of Arabidopsis thaliana Encodes ent-Kaurene Synthase of Gibberellin Biosynthesis. Plant physiol. 1998;116:1271–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Moffatt B, Somerville C. Positive selection for male-sterile mutants of Arabidopsis lacking adenine phosphoribosyl transferase activity. Plant physiol. 1988;86:1150–4.View ArticlePubMed CentralPubMedGoogle Scholar
- Mo YY, Nagel C, Taylor LP. Biochemical complementation of chalcone synthase mutants defines a role for flavonols in functional pollen. Proc Natl Acad Sci USA. 1992;89:7213–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Taylor LP, Jorgensen R. Conditional male-fertility in chalconesynthase -deficient petunia. J Hered. 1992;83:11–7.Google Scholar
- Buer CS, Imin N, Djordjevic MA. Flavonoids: new roles for old molecules. J Integr Plant Biol. 2010;52(1):98–111.View ArticlePubMedGoogle Scholar
- Ylstra B, Touraev A, Moreno RMB, Stöger E, van Tunen AJ, Vicente O, et al. Flavonols stimulate development, germination, and tube growth of tobacco pollen. Plant physiol. 1992;100(2):902–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Doyle MR, Davis SJ, Bastow RM, McWatters HG, Kozma-Bognár L, Nagy F, et al. The ELF4 gene controls circadian rhythms and flowering time in Arabidopsis thaliana. Nature. 2002;419(6902):74–7.View ArticlePubMedGoogle Scholar
- Bazzaz FA, Carlson RW, Harper JL. Contribution to reproductive effort by photosynthesis of flowers and fruits. Nature. 1979;279:554–5.View ArticleGoogle Scholar
- Hatefi Y. The mitochondrial electron transport and oxidative phosphorylation system. Annu Rev Biochem. 1985;54(1):1015–69.View ArticlePubMedGoogle Scholar
- Li R, Li Y, Kristiansen K, Wang J. SOAP: short oligonucleotide alignment program. Bioinformatics. 2008;24:713–4.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(T) (-Delta Delta C) method. Methods. 2001;25:402–8.View ArticlePubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.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.