The study of two barley Type I-like MADS-boxgenes as potential targets of epigenetic regulation during seed development
© Kapazoglou et al.; licensee BioMed Central Ltd. 2012
Received: 10 February 2012
Accepted: 30 August 2012
Published: 17 September 2012
MADS-box genes constitute a large family of transcription factors functioning as key regulators of many processes during plant vegetative and reproductive development. Type II MADS-box genes have been intensively investigated and are mostly involved in vegetative and flowering development. A growing number of studies of Type I MADS-box genes in Arabidopsis, have assigned crucial roles for these genes in gamete and seed development and have demonstrated that a number of Type I MADS-box genes are epigenetically regulated by DNA methylation and histone modifications. However, reports on agronomically important cereals such as barley and wheat are scarce.
Here we report the identification and characterization of two Type I-like MADS-box genes, from barley (Hordeum vulgare), a monocot cereal crop of high agronomic importance. Protein sequence and phylogenetic analysis showed that the putative proteins are related to Type I MADS-box proteins, and classified them in a distinct cereal clade. Significant differences in gene expression among seed developmental stages and between barley cultivars with varying seed size were revealed for both genes. One of these genes was shown to be induced by the seed development- and stress-related hormones ABA and JA whereas in situ hybridizations localized the other gene to specific endosperm sub-compartments. The genomic organization of the latter has high conservation with the cereal Type I-like MADS-box homologues and the chromosomal position of both genes is close to markers associated with seed quality traits. DNA methylation differences are present in the upstream and downstream regulatory regions of the barley Type I-like MADS-box genes in two different developmental stages and in response to ABA treatment which may be associated with gene expression differences.
Two barley MADS-box genes were studied that are related to Type I MADS-box genes. Differential expression in different seed developmental stages as well as in barley cultivars with different seed size was evidenced for both genes. The two barley Type I MADS-box genes were found to be induced by ABA and JA. DNA methylation differences in different seed developmental stages and after exogenous application of ABA is suggestive of epigenetic regulation of gene expression. The study of barley Type I-like MADS-box genes extends our investigations of gene regulation during endosperm and seed development in a monocot crop like barley.
KeywordsMADS-Box Epigenetic regulation Chromatin DNA methylation Histone methylation Seed development Endosperm Retrotransposon Barley
In angiosperms, the endosperm of the developing seeds is formed as a result of the double fertilization event. Fertilization of the egg cell by a sperm cell from the male gametophyte generates the diploid embryo from which the tissues, organs, and shoot meristems of the plant will be generated. Fertilization of the adjacent central cell by a second sperm cell forms a triploid endosperm which supports embryo growth and development by producing storage proteins, lipids and starch [1, 2]. During this process, a large number of genes are activated. Epigenetic regulatory considerations become important in relation to which parental allele will be expressed and in which reproductive tissue, a factor that ultimately governs among other things the size of the endosperm [3–6].
Developmental transition of plants from the vegetative to the reproductive stage and from the floral stage to the seed stage relies on the activity of MADS-box transcription factors. With the completion of the Arabidopsis genome, more than 100 genes encoding for MADS-box transcription factors were uncovered which can be phylogenetically classified in five clades, termed MIKC, Mα, Mβ, Mγ and Mδ (Arabidopsis Genome Initiative, 2000) [7, 8]. The Mα, Mβ, Mγ subfamilies comprise the Type I lineage, whereas the MIKC and Mδ constitute the Type II lineage. Similarly, phylogenetic analysis based on the completed rice genome revealed 44 Type II MADS-box genes, of the MIKC and Mδ clades, and 35 Type I MADS-box genes of the Mα, Mβ, Mγ clades . The two lineages are proposed to have arisen by a gene duplication that took place in the common eukaryotic ancestor more than a billion years ago [10–13]. MIKC genes harbor the characteristic MADS-box domain (M), the Intervening domain (I), the Keratin-like domain (K) and the C-terminal domain (C). In contrast, the Type I genes contain only the MADS-box domain (M). Although the Type I genes represent about 60% of the total MADS genes in Arabidopsis, they have only recently started to be studied, whereas the Type II genes have been extensively investigated both structurally and functionally [8, 14–17]. This is due mostly to the fact that Type I MADS genes function during early gametogenesis, embryogenesis and seed formation and consequently homozygous mutants are lethal. In addition their expression is very low during these developmental stages [8, 18]. Among the few Type I family members that have been studied is PHERES 1 (PHE1) [also known as AGAMOUS-LIKE37 (AGL37)], a Mγ-type gene with an important role in Arabidopsis gamete and seed development . Epigenetic regulatory mechanisms have been implicated in the transcriptional control of this gene. The PHE1 gene is expressed transiently at high levels immediately after fertilization in a parentally imprinted manner where the paternal allele is expressed whereas the maternal allele is silenced [19, 20]. Both DNA and histone methylation are responsible for these silencing events [19–21]. The PRC2 Polycomb group complex is partially responsible for keeping the maternal allele silent in the female gametophyte and in the seed after fertilization [19–21]. This is achieved through the histone methylation activity of the PRC2 complex, conferred by one of its subunits, the histone methyltransferase MEDEA (MEA) with histone 3 lysine 27 trimethylation activity (H3K27me3). The PRC2 complex also restricts the expression of PHE1 to the chalazal domain of the endosperm after fertilization. Arabidopsis mea mutants show upregulated PHE1 expression and form defective seed-like structures before fertilization and endosperm overproliferation after fertilization [19, 20, 22]. Furthermore, a distantly located region downstream of paternal PHE1 was found to have a DNA methylation requirement for PHE1 expression .
Additional Type I MADS-box genes have been implicated in gamete and seed development. The AGAMOUS-LIKE23 (AGL23) gene was found to regulate female gametophyte formation and normal embryo development . AGAMOUS-LIKE 80 (AGL80) was shown to be critical for central cell and endosperm development in Arabidopsis . In agl80 mutants the formation of the central cell is defective, and endosperm development fails to initiate after fertilization. Similarly, AGAMOUS-LIKE 61 (AGL61) also termed DIANA, plays an important role in central cell and endosperm formation, in Arabidopsis [16, 25]. Like agl80, agl61 mutants have aberrant central cell morphology which degenerates before fertilization occurs, and endosperm development does not take place post fertilization [16, 25]. Both genes are expressed exclusively in the central cell and early endosperm and this along with their similar mutant phenotype suggests that AGL61 and AGL80 proteins may function as heterodimers within the central cell. In agreement to this, AGL61 and AGL80 proteins were found to interact in yeast two-hybrid assays [16, 25]. Of equal interest is AGL62 with high expression in the early nuclear endosperm and sharp decline right before cellularization . In agl62 mutants, the endosperm cellularizes prematurely indicating that the AGL62 is required for repression of precocious cellularization during the syncytial phase. AGL62 is under the epigenetic control of PRC2 genes, as in prc2 mutants AGL62 fails to become silent and endosperm cellularization is arrested . Thus, AGL62 seems to regulate the timing of endosperm cellularization, which is triggered epigenetically by PRC2-mediated AGL62 silencing. Likewise, a recent study on AGL36 demonstrated that the expression of the AGL36 maternal allele is epigenetically controlled in a sequential manner, firstly upregulated by the function of DEMETER (DME), a DNA glycosylase enzyme responsible for demethylating MEDEA, and then downregulated during endosperm development by the PRC2 complex . Understanding the genetic and epigenetic processes controlling the timing of endosperm cellularization could be of particular interest in agriculture, as premature or delayed cellularization is associated with small and large seed size and weight, respectively.
An extensive genome-wide study of 60 Type I MADS-box genes in Arabidopsis uncovered a cell-type-specific expression pattern during female gametophyte and early endosperm development . Most genes are expressed in the central cell and antipodal cells in the female gametophyte and in the chalazal and peripheral endosperm of 1–2 days after fertilization (DAF) developing seeds. These data are in agreement with the results from earlier functional studies of individual genes like PHE1, AGL23, AGL61, AGL62 and AGL80 and propose a role for Type I MADS-box genes in female gametophyte and endosperm development.
Considerably less is known about Type I MADS-box genes in other plant species, especially monocots. In wheat, 42 MADS-box genes were identified in silico, of which 8 were classified as Type I. In rice, a global-scale microarray expression analysis identified an expression pattern for all Type II and Type I MADS-box genes during vegetative and reproductive development . Overall Type I MADS-box genes had lower expression levels than Type II MADS-box genes. Certain Type I genes were expressed throughout development, whereas others exhibited more specific expression in particular tissues. For example, four rice Type I genes, of the subclass Mα, were found predominately expressed in seeds 5–20 days after fertilization suggesting a role for these genes in seed development. Recently, a rice Type I MADS-box gene, OsMADS87, was reported to be maternally expressed in rice endosperm and associated to endosperm developmental transitions caused by interspecific hybridization .
Cereal crops account for about 50% of global human calorific intake with the endosperm of cereal seeds being one of the most important sources (faostat.fao.org). Chief among these monocots are barley, wheat, rice and maize. Besides human nutrition, cereal crops also represent major sources of feedstock, fiber and recently biofuel substrates. Contrary to Arabidopsis and other dicots where the endosperm is consumed during seed development, in monocots, such as barley and other cereals, the endosperm persists and constitutes the nutritional part of the seed containing storage proteins and starch.
Considering the implications of Type I MADS-box genes on seed development and the agronomic importance of cereal endosperm, we set out to identify and characterize Type I MADS-box genes in barley. At least three of the Arabidopsis Type I MADS-box genes, PHE1, AGL62, and AGL36 are under epigenetic regulation mediated, in part, by the chromatin repressive enzymatic complex PRC2. Our group has recently characterized barley genes encoding a putative PRC2 Polycomb group complex . To further focus our research efforts in this area and study potential barley targets of a PRC2 complex, we report here the identification and structural characterization of two barley Type I-like MADS-box genes. Their expression has been studied in different tissues and seed developmental stages, in two cultivars with varying seed size, and after exogenous application of the developmental- and stress-related phytohormone ABA. In addition, their genomic organization was examined and compared to their cereal homologues. Finally, the 5’upstream regions were analyzed for conserved cis regulatory elements and the DNA methylation patterns of upstream and downstream regions were investigated in two tissues with differential Type I-like MADS-box gene expression.
Commercial barley cultivars, Caresse, Byzantio, and Ippolytos differing in seed size and weight were planted in the field and were the source of total RNA for expression analysis. For Caresse, the weight of 1000 grains is 50–55 gr, and 98% of seeds have diameter longer than 2.5 mm, for Byzantio the weight of 1000 grains is 36–42 gr and 75% of seeds have diameter longer than 2.5 mm, whereas for Ippolytos, seeds weight 25–31 gr per 1000 grains and only 35–45% of seeds have diameter longer than 2.5 mm (http://www.cerealinstitute.gr).
Seven-day-old seedlings (Caresse) grown in a growth chamber (16 hours (h) light, 8 h darkness, at 22°C) were sprayed with 100 μM ABA, (abscisic acid +/− cis, trans-ABA, SIGMA) and 100 μM JA (methyl jasmonate, ALDRICH). Aerial parts of plants were collected at 6 h and 24 h after treatment and immediately stored in liquid nitrogen. Aerial parts from five plants were pooled together for RNA extraction for each time point. Control plants were sprayed with water plus 0.2% Tween.
RNA isolation and first strand cDNA synthesis
Total RNA was isolated from roots, shoots, apical meristems, first leaves of seedlings, flowers before fertilization (immature flowers), seeds 1–3, 3–5, 5–10, 10–15, 15–20 days after fertilization (DAF), and aerial parts after hormonal treatments, respectively, using TRI REAGENT (SIGMA) according to the instructions of the manufacturer. First strand cDNA synthesis was performed using 1.0 μg total RNA, 0.5 μg 3’ RACE Adapter primer, 5’-GGCCACGCGTCGACTAGTAC (T)17-3’ (Invitrogen), 1 mM dNTPs and 200U of Superscript II (Invitrogen) in 20 μL total volume, according to the specifications of the manufacturer.
Protein sequence analysis
Type I MADS-box sequences used for alignments and phylogenetic tree construction
Orysa sativa indica
Orysa sativa japonica
The genomic sequences of Brachypodium distachyon _Bradi2g59120, Brachypodium distachyon_ Bradi2g59190, Oryza sativa_OsMADS65_Q9XJ61 (Os01g0922800), and Zea mays_ZmMADS1_GRMZM2G171650_B4FML1, were downloaded from the Phytozome database (http://www.phytozome.net/). The sequences of the two full length barley Type-like MADS-box ESTs were used to interrogate the barley database http://webblast.ipk-gatersleben.de/barley/index.php for detection of genomic sequences. Genomic organization of exons and introns was obtained using the mRNA-to-genomic alignment Spidey tool, in NCBI (http://www.ncbi.nlm.nih.gov/spidey). Detection of retroelements was performed with the MASiVE and LTRharvester tools (http://tools.bat.ina.certh.gr/masive/) (http://tools.bat.ina.certh.gr/ltrharvester/) and homology was visualized with Circoletto (http://tools.bat.ina.certh.gr/circoletto/), all three tools developed in-house at INA, by the Bioinformatics Analysis Team (BAT).
Mapping in silico
An in silico visual comparative analysis was performed against the Barley, OPA 2009, Consensus-Hordeum-OPA-2009-3 H , the Barley, OPA123-2008, Consensus-Hordeum-OPA123-2008-3 H and the Barley, OWB, OPA2008-Hordeum-OWB-OPA2008-3 H  using the HarvEST tool and the comparative map viewer (Cmap) available at Gramene (http://www.gramene.org/cmap/).
Expression analysis of barley Type I-like MADS-box genes
Qualitative RT-PCR and quantitative real-time RT-PCR was performed with cDNA synthesized from 1 μg of total RNA from roots, stems, meristems, leaves, immature flowers, seeds 1–3 DAF, 3–5 DAF, 5–10 DAF, 10–15 DAF and 15–20 DAF and aerial parts of seedlings after ABA and JA treatment. For real-time PCR, each sample reaction was set up in a PCR reaction mix (20 μl) containing 5 μl of the 1:50 diluted cDNA, 0.25 μM of each primer and 1X Platinum SYBR Green qPCR Supermix-UDG (Invitrogen, Paisley, UK) and using the Corbett Rotor Gene 6000. Each reaction was performed in triplicates. General thermocycler conditions were 50°C for 2 min, 95°C for 2 min, then 42 cycles of 95°C for 15 sec, annealing [53°C for HvOS1 and 57°C for HvOS2 for 20 sec, extension 72°C for 20 sec, then 72°C for 5 min. To identify the PCR products a melting curve was performed from 65°C to 95°C with observations every 0.2°C and a 10-s hold between observations. Relative quantification was performed using actin as the reference gene and HvActinF / HvActinR as primers. The barley genes HVA22 and HvADC (arginine decarboxylase 2) which are known to be induced by ABA and JA, respectively [39, 40], were used as positive controls. All primers used in expression analysis correspond to non-conserved regions and are shown in (Additional file 1: Table S2).
In situ hybridization experiments were performed as it has been described previously . Briefly, seeds were fixed in 4% (w/v) paraformaldehyde supplemented with 0.25% (v/v) glutaraldehyde in 10 mM sodium phosphate buffer (pH 7.4) for 4 h in a vacuum aspirator. Fixed tissues were block-stained in 0.5% (w/v) safranin, dehydrated through ethanol series, embedded in paraffin and cut into 8 mm-thin sections. Antisense RNA probes labelled with digoxigenin-11-rUTP (Boehringer Mannheim, Mannheim, Germany) were originated from PCR-generated templates incorporating T3 polymerase sites. The probe was designed close to the 3’-UTR of the gene and its length was 202 bp. Sections were prepared for hybridization as described before  and hybridized overnight at 42°C in 50% (v/v) formamide, 300 mM NaCl, 10 mM Tris–HCl pH 7.5, 1 mM EDTA, 0.02% (w/v) Ficoll, 0.02% (w/v) polyvinylpyrrolidone, 0.025% (w/v) bovine serum albumin (BSA), 10% (v/v) dextran sulfate and 60 mM DTT. After hybridization, the sections were treated with a solution containing 500 mM NaCl, 1 mM EDTA, 10 mM Tris–HCl and 50 μg/ml RNase A. Finally, sections were washed several times in a 2xSSC solution. Hybridization signals were visualized with anti-digoxigenin antibodies conjugated with alkaline phosphatase. Images were processed using Photoshop 7 software (Adobe Systems Inc., San Jose, CA, USA).
Genome walking and 5’ upstream sequence analysis
Genome walking experiments were conducted on leaf genomic DNA obtained from the barley Caresse cultivar. The Clontech genome walking kit was used and procedures were conducted according to manufacturer’s specifications. Sequence information was obtained by doing three progressive rounds of genome walking using EST sequence data as starting points. Statistically significant prediction for CpG islands was performed using the online predictor, which is part of the sequence manipulation suite at http://www.bioinformatics.org/SMS/index.html. The prediction of the putative cis acting elements was accomplished using the TSSP /Prediction of PLANT Promoters algorithm (Using RegSite Plant DB, Softberry Inc.) in the SoftBerry database (http://linux1.softberry.com/cgi-bin/programs/promoter/tssp.pl) and PlantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/). The position of the putative 5’ upstream regulatory elements is indicated according to the start codon (ATG).
DNA methylation assays
Genomic DNA was prepared from immature flowers and from 1–3 DAF seeds (Caresse) with Qiagen columns following the protocol of the manufacturer (Qiagen Plant genomic DNA kit). Cytosine DNA methylation was analyzed by restricting 1 μg of genomic DNA from each sample with the methylation-dependent enzyme McrBC (NEB Biolabs), according to the manufacturer’s instructions, and PCR-amplifying equal quantities of McrBC-treated and untreated samples. Primers used are shown in (Additional file 1: Table S2).
Identification of Type I-like MADS-box genes in barley and protein sequence analyses
Initial efforts to identify orthologues of Type I MADS-box genes in barley EST databases via BLAST and other sequence similarity-based approaches failed to result in significant hits. Additional efforts using degenerate primers designed using Arabidopsis sequence data and barley cDNA from multiple cultivars did not result in target gene isolation either (data not shown). Hence a more rigorous bioinformatics approach was employed and a monocot-specific C-terminus probe was generated using existing sequence data for MADS-box putative proteins. The sequence of this probe was the following:
GAGXXVNGXOXXXNXDXXXXOOQXXLKEIAXWXXQNNAOXXDANOLEKLEOLLTOALRNTKXKKMLXOONXG, where X is any amino acid and O is a possible gap. This probe was used to isolate several cereal Type I-like MADS-box sequences from available EST databases and the GenBank non-redundant CDS translated database. Of particular interest were two hits from the barley EST database, chosen based on the presence of a MADS-box domain and the absence of a K-domain, characteristic of Type I MADS-box proteins. These candidates were: BG365393 and TC178280 (TIGR database, http://www.tigr.org) (see Materials and Methods). During the course of our investigations another study characterizing these sequences with respect to their role in the vernalization process was reported . These authors designated these sequences as ODDSOC1(HvOS1) and ODDSOC2(HvOS2), respectively. In our recent study on the identification and characterization of two SOC1-like gene homologues from barley  it was shown that two HvSOC1-like proteins containing the K-box domain are closely related to the Type II MADS-box proteins whereas they are more distantly related to ODDSOC1 and ODDSOC2. However, since Greenup et al. (2010)  were the first to report these genes we will utilize the nomenclature introduced by these authors and hereafter will be referring to the sequences BG365393 and TC178280 as HvOS1 and HvOS2, respectively.
Expression analysis of the Type I-like HvMADS-box genes in different tissues and during seed development
Real time PCR analysis was employed to examine and compare the expression of barley Type I-like HvMADS-box genes at different seed developmental stages and in different cultivars, Caresse (a large-seed cultivar) and Ippolytos (a small-seed cultivar) (Figure 4B).
For HvOS1 , an induction of approximately 8 fold in Caresse 1–3 DAF and 3–5 DAF seed, and a decrease thereafter, was observed. In Ippolytos, HvOS1 was induced in 1–3 and 3–5 DAF seeds by approximately 2 fold. Transcript levels dropped by about 2 fold in 5–10 DAF seeds whereas they increased slightly in 10–15 DAF as compared to immature flowers. In Caresse 15–20 DAF seeds transcript levels increased by about 2 fold, whereas in Ippolytos 15–20 DAF they decreased by about 5 fold (Figure 4B).
A marked increase of approximately 8 fold in HvOS2 transcript accumulation was observed in Caresse seeds 1–3 DAF and 3–5 DAF, as compared to immature unfertilized flowers. Expression levels dropped thereafter in 5–10, 10–15 and 15–20 DAF seeds to levels comparable to those of unfertilized flowers. In the small-seed cultivar, Ippolytos, HvOS2 exhibited substantial transcript accumulation in seeds 1–3 DAF and 3–5 DAF. In 5–10 DAF and 10–15 DAF seeds HvOS2 had a pronounced increase of approximately 10–15 fold, in contrast to Caresse. In Caresse seeds 15–20 DAF there was a slight increase of HvOS2 transcript. Conversely, a substantial decrease of about 10 fold was seen in Ippolytos 15–20 DAF seeds (Figure 4B). Expression of HvOS2 was examined in another large-seed cultivar, Byzantio. Similar to Caresse, HvOS2 transcript levels were significantly lower in 5–10, 10–15 DAF seeds as compared to Ippolytos (Additional File 2).
More detailed expression of the HvOS2 gene within the seed was examined by in situ localizations using cross sections of seeds and DIG-labelled antisense probe. Strong hybridization signal was observed in the aleurone layer and the first layer of endosperm cells adjacent to the aleurone layer that was absent from the sense-control (Figure 4C), suggesting cell-specific expression in the endosperm. In addition, strong hybridization signal was observed in embryo cells between sections hybridized with antisense as well as sense probes suggesting no specific expression in the embryo. HvOS1 analysis was not possible due to lack of appropriate sequences specific for the HvOS1 gene satisfying the particular needs of this experiment.
Expression analysis of the Type I-like HvMADS-box genes in response to JA and ABA
HvOS1 showed a significant induction of about 3 fold and 12 fold at 6 h and 24 h, respectively, after JA treatment, as compared to the control. The HvOS2 transcript was increased by about 3 fold at 24 h after JA treatment (Figure 5A). Two barley genes HVA22 and HvADC2 known to be induced by ABA and JA, respectively, were used as positive control genes for the ABA and JA treatments (Figure 5B).
Greenup et al. (2010) have assigned a position for HvOS1 and HvOS2 on barley chromosome 3 H, in syntenic regions to Brachypodium and rice, and proximal to the unigene 17348 . Analysis in HarvEST found unigene 17348 to reside close to the POPA3_0092 SNP marker on the comparative barley maps. This position is in proximity to important QTL markers for seed quality traits such as grain protein content (QGpc.StMo-3 H) and grain yield (QYld.StMo-3 H) (Additional File 4).
Presence of regulatory elements in the 5’ upstream regions of Type I-like HvMADS-box genes
5’ upstream sequences of the Type I-like HvMADS-box genes were isolated (542 bp 5’ from the ATG site for HvOS1, and 2185 bp 5’ from the ATG site for HvOS2) using genome walking in the cultivar Caresse. 5’ upstream regions of these genes were interrogated for known cis-acting regulatory elements and results were analyzed with special attention to elements identified and characterized in cereals. The 5’ upstream region of both HvOS1 and HvOS2 genes showed the presence of endosperm-specific-elements, and ABA and JA responsive elements (Figure 7 and Additional File 5). HvOS1 also contains a gibberellin-responsive-element. TATA and CAAT boxes were also detected. Further inspection of these sequences was carried out in order to identify regions which may be prone to cytosine-methylation. Predicted CpG islands which may play important roles in the regulation of gene expression via DNA methylation-directed silencing were identified in both sequences (Figure 7 and Additional File 5).
DNA methylation differences in the 5’ and 3’ flanking regions of Type I-like HvMADS-box genes
In order to examine the DNA methylation pattern of the Type I-like HvMADS-box regulatory regions and uncover potential links to gene expression differences among different developmental stages, we performed McrBC digestions of genomic DNA and subsequent PCR amplification analysis. McrBC is a methylation-dependent restriction enzyme which digests DNA that is methylated at two or more cytosines, thus reducing PCR amplification of a selected fragment proportionally to the methylated cytosines . Therefore the presence or absence (as well as reduction) of amplicons after PCR amplification of McrBC-digested genomic DNA with primers specific for a particular region suggests lower or higher degree of methylation within this region, respectively. We chose two developmental stages with the largest differences in expression, Caresse immature flowers and 1–3 DAF seeds.
For HvOS2, two non-coding regions, a 1170 bp region [located at −1525 bp upstream from the ATG site (region 1)], and a 264 bp region located 3’downstream from the coding region (region 2), and two coding regions (region 3 and region 4) were analyzed (Figure 8, middle panel). PCR amplification of immature flower genomic DNA digested with McrBC did not produce the expected 1170 bp upstream fragment contrary to the digested sample, suggesting the presence of methylated sites within region 1 in immature flower genomic DNA. The 1170 bp amplicon was produced from both McrBC-digested and undigested genomic DNA from 1–3 DAF seeds suggesting a lower degree of methylation within region 1 in 1–3 DAF genomic DNA. The presence of higher DNA methylation level in the promoter of HvOS2 in immature flowers may be associated with downregulation of this gene in this tissue as compared to the high levels of expression in 1–3 DAF seeds. Examination of the 3’ downstream region (region 2) revealed that the expected fragment was amplified after McrBC-PCR amplification of immature flower genomic DNA although to a lesser extent than the undigested sample (Figure 8, middle panel). Interestingly, amplification of McrBC-digested 1–3 DAF genomic DNA produced nearly no amplicon, implying a much higher degree of DNA methylation in the 3’downstream region of HvOS2 in 1–3 DAF than in immature flowers. Finally, examination of two sites within the coding region, in Exon 1 and Exon 5 (region 3 and region 4), did not show any decrease of PCR amplification in McrBC-digested as compared to undigested DNA, either in immature flowers or in 1–3 DAF (Figure 8, middle panel-inlet).
Two regions were also examined within the HvActin gene, one within the coding region (region 1) and one in the 3’region downstream from the coding region (region 2). There were no obvious differences in the amplification of these regions between McrBC-digested and non-digested genomic DNA, suggesting no methylation sites in these fragments (Figure 8, bottom panel). The larger than expected size of the PCR product after amplification of region 1 (300 bp instead of 235 bp) implies the presence of a small intron in that region.
In the current study we present the identification, protein sequence structural analysis, phylogenetic characterization, expression profiles, genomic organization and promoter analysis of two genes encoding Type I-like MADS-box proteins in barley. During the course of our studies Greenup et al. (2010) reported the identification of these genes in barley and named them HvOS1 and HvOS2. These authors focused their investigations on the role of these genes in the vernalization response and proposed a cereal-specific pathway controlling vernalization-induced flowering in temperate cereals different from the one reported for Arabidopsis.
We have focused our study on detailed in silico structural analysis of HvOS1 and HvOS2 protein sequences and on gene expression during seed development, in different cultivars, and in response to the seed developmental- and stress-related hormones ABA and JA. Moreover we have examined 5’ and 3’ regulatory regions of HvOS1 and HvOS2 for differential DNA methylation patterns that may be associated with differential gene expression.
Protein sequence analysis
Both HvOS1 and HvOS2 encode proteins possessing a MADS-box domain, lacking the typical K-box of Type II MADS-box proteins, and having highest resemblance to two Type I MADS-box proteins from wheat, TaAGL-33 and TaAGL-42, respectively. Phylogenetic analysis showed that the two barley Type I-like MADS-box proteins together with their putative orthologues from wheat (TaAGL-33 and TaAGL-42), brachypodium (Bradi2g59190, Bradi2g59120), rice (OsMADS65) and maize (ZmB4FML1), form a cluster which is closer related to the alpha clade of Type I MADS-box proteins (Mα). Interestingly this group does not contain sequences from Arabidopsis or other dicots.
Differential expression during vegetative, reproductive and seed development
Qualitative RT-PCR expression analysis in different tissues showed that the HvOS1 transcript was expressed in roots, apical meristem, young shoots, leaves, and immature flowers in Ippolytos but was present only in immature flowers in Caresse and was not detected in the other tissues examined. On the other hand, the HvOS2 transcript is present in roots, apical meristem, young shoots, leaves of seedlings and in immature flowers in two different barley cultivars, Caresse and Ippolytos. In that sense, HvOS1 expression in Ippolytos and HvOS2 expression in both cultivars resembles the Arabidopsis Type I MADS-box gene AGL28 gene which is expressed widely in vegetative and reproductive tissues .
Type I MADS-box genes have been shown to be expressed in the female gametophyte and the developing seeds in Arabidopsis and rice and have been assigned crucial roles in seed development through functional characterization in Arabidopsis. In order to start understanding the role of the two barley Type I-like MADS-box genes in seed development we conducted quantitative real-time PCR analysis to examine their expression in unfertilized flowers and during different seed developmental stages in two barley cultivars with varying seed size, Caresse (large seed) and Ippolytos (small seed). Quantitative real time PCR analysis demonstrated significant differences in gene expression among different seed developmental stages within a cultivar and between two cultivars with different seed size. In particular, both genes exhibited a remarkable induction post-fertilization, in 1–3 DAF and 3–5 DAF seeds in Caresse (approximately 8 fold) and a decline at the later stages of seed development. Likewise, an increase in transcript levels was evidenced for both genes in 1–3 DAF and 3–5 DAF seeds, in Ippolytos. However, in Ippolytos a remarkable increase was also observed for the HvOS2 gene at the later seed stages of 10–15 DAF and 15–10 DAF and a large reduction in 15–20 DAF, whereas such changes were not observed in Caresse.
The expression pattern of Type I-like HvMADS-box genes agrees with the expression of most Type I MADS-box genes in Arabidopsis, where they were found to be induced in early syncytial endosperm up to the 4 DAF seed developmental stage [8, 18, 48]. HvOS2 resembles its closest wheat homologue, TaAGL33, being expressed in 6–19 DAF seeds as well as in roots and leaves of seedlings . On the other hand, expression of the wheat TaAGL42, the closest homologue to HvOS1, is not detected either in immature flowers or in 6–19 DAF seeds of wheat, implying a different role for this homologue in wheat at least in the particular cultivar studied. The rice OsMADS65 homologue has been reported to be widely expressed . Microarray and quantitative PCR analysis has generated an expression profile for the Type I MADS-box genes in rice . Most genes are expressed, during the vegetative and reproductive stages and during seed development. Interestingly, quantitative PCR analysis showed that three genes of the Mα clade, OsMADS71, OsMADS78 and OsMADS79, have a substantial induction in 5–10 DAF and 10–15 DAF seeds resembling the expression pattern of HvOS2 in the Ippolytos cultivar. OsMADS71, OsMADS78 and OsMADS79 are the closest rice homologues to the barley Type I-like HvMADS-box genes, after OsMADS65. Protein sequence similarity together with similar expression patterns may reflect similar functional roles for these homologues in these two closely related cereal species.
Significant differences in gene expression were observed between Caresse, a large-seed cultivar and Ippolytos, a small-seed cultivar. A striking increase in HvOS2 transcript was found in Ippolytos 5–10 DAF seeds (~20fold) and 10–15 DAF seeds (~15fold) as compared to the same seed developmental stages in Caresse, suggesting cultivar-dependent differential gene expression. An increase in HvOS1 transcript was also evidenced in Ippolytos 10–15 DAF seeds (~4fold) as compared to the same seed stage in Caresse. Similar to Caresse, examination of HvOS2 expression in the large-seed cultivar, Byzantio, also found significant differences (10–15 fold) at the 5–10 and 10–15 DAF seed stages between Byzantio and the small-seed cultivar Ippolytos, supporting a possible association of HvOS2 expression and seed size.
Differential expression of HvOS1 and HvOS2 (pronounced at the 5–10 and 10–15 DAF seed stages) between Caresse (a large-seed cultivar) and Ippolytos (a small-seed cultivar) may suggest an association of these genes with the size of seed. In barley, endosperm cellularization begins at approximately 4 DAF and ends at 6–8 DAF, when the seed maturation process begins [1, 2, 49]. It may be possible that the differences in Type I-like HvMADS-box gene expression between the two cultivars during these critical stages of endosperm development are associated with the processes of cellularization and seed filling and ultimately with grain size. Certainly, a large number of genotypes with different seed size should be further tested to conclusively assess any association of these genes with the size of seed.
In situ hybridization experiments localized HvOS2 in the aleurone layer and in endosperm-specific cells situated adjacent to the aleurone layer, in 10 DAF seeds. Specific expression of this gene within the endosperm may be associated with particular gene expression programmes involved in seed maturation.
Analysis of the genomic organization of the HvOS2 revealed resemblance with its closest cereal homologues OsMADS65, Bradi2g59190, Bradi2g59120 and ZmB4FML1. Extensive conservation in the exon-intron organization is evident among all five genes with five exons and four introns at approximately the same relative positions and approximately the same size except from the fifth exon of the maize gene ZmB4FML1. Genomic organization for HvOS1 was not possible due to lack of barley sequence information. With the rapid advancement in the sequencing of the barley genome [50, 51] the complete genomic sequence of both HvOS1 and HvOS2 will become available soon for further analysis.
A strikingly large first intron is present in all five cereal homologues studied. In OsMADS65, Bradi2g59190, Bradi2g59120, ZmB4FML1 this intron accounts for about 85-90% of the entire length of the gene bodies (here exons plus introns). Notably, our analysis showed that the maize ZmB4FML1 intron has high similarity to a full length maize Sirevirus retrotransposon, Copia PZMAY_CS_U68408_OPIE. Likewise, a fragment of the first intron of the Bradi2g59120 gene exhibits high similarity to remnants of PBDIS_GX_ADDN010010, a brachypodium Gypsy retrotransposon. Homology with retrotransposons was not detected in the first intron of barley HvOS2, thus far.
Recently a sensitive comparative analysis of the genomes of the plant-specific Sireviruses retrotransposons was performed [45, 46], in an effort to understand their significance in the evolution of their plant hosts and their putative role in the epigenetic regulation of plant genes. In addition, our group reported a systematic analysis of >10,600 intact and ~28,000 degenerate maize Sireviruses  discovered by the MASiVE algorithm in maize, and showed that Sireviruses account for 90% of the Copia type population, comprising 21% of the maize genome, and reside in near-gene areas. The element PZMAY_CS_U68408_OPIE found within the intron of the maize ZmB4FML1 gene belongs to Opie, one of the most abundant families of the Copia genera and its location within an intron may suggest additional insertion sites of Sireviruses within genes.
It has become more and more evident that intact or remnants of transposable elements may be implicated in the transcriptional regulation of flowering plants especially in the control of gene imprinting [5, 53]. Most imprinted genes, mainly expressed only from the maternal allele, are found in the female gametophyte and in the endosperm of the developing seed . Massive demethylation of transposable elements and regions generating siRNAs was shown to take place in the Arabidopsis endosperm as opposed to the embryo [53, 55] and it has been hypothesized that fragments of transposable element inserted near genic regulatory elements may guide DNA methylation at specific alleles and govern imprinted gene expression (for example silencing of paternal alleles) [4, 56]. Presence of transposable elements within introns in genes of cereals may perform similar regulatory functions.
In silico mapping of HvOS1 and HvOS2 identified a corresponding SNP marker on barley chromosome 3 H close to regions associated with important agronomical traits such as grain protein content and grain yield. This could prove important for the development of functional markers to be used in marker assisted selection for breeding programmes.
Promoter sequence analysis revealed putative CpG regions within the 5’upstream regions of both HvOS1 and HvOS2. In addition several regulatory elements were identified, including elements related to endosperm and seed development such as endosperm-specific elements and ABA and gibberrellin-responsive elements. These observations together with the expression patterns at different seed developmental stages and in response to the seed development-related hormone ABA, mentioned above, reinforces the suggestion for functional roles for these genes in barley seed development.
JA responsive elements were also identified on the promoters of both HvOS1 and HvOS2. ABA and JA elements identified in the 5’upstream regions of HvOS1 and HvOS2 were also detected in the promoters of several of the rice homologues such as OsMADS65 and OsMADS84 and brachypodium homologues, Bradi2g59190 and Bradi2g59120 (data not shown). The hormones JA and ABA are known to mediate abiotic stress responses in plants. It is possible that Type I-like HvMADS-box genes are responsive to ABA- or JA- associated stress conditions. In support to this the TaAGL42 homologue is induced by cold stress in two winter wheat cultivars  and a number of Type I-MADS-box genes in rice were found to be induced by abiotic stress such as cold, drought, and increased salinity .
DNA methylation analysis was undertaken to reveal whether differentially methylated regions exist in the 5’upstream and 3’downstream regions of HvOS1 and HvOS2 between two different seed developmental stages with different gene expression levels. In particular, using the McrBC-PCR assay it was shown that the 5’upstream region of HvOS2 is highly methylated in Caresse immature flowers as compared to 1–3 DAF seeds where lower degree of methylation was observed. It is likely that the observed downregulation of HvOS2 in immature flowers may be due, in part, to a DNA methylation silencing effect. Similarly, a higher degree of DNA methylation in the 5’ upstream region of HvOS1 in immature flowers as opposed to low DNA methylation in 1–3 DAF seeds suggests that downregulation of HvOS1 in immature flower maybe be partially due to a DNA methylation repressive mechanism. Notably, a genome-wide analysis in Arabidopsis demonstrated extensive global DNA demethylation in 5’ upstream and 3’ downstream regions in endosperm but not in embryo tissue and this is associated with gene activation in endosperm and gene silencing in the embryo, respectively . A similar epigenetic mechanism could be operating in cereals seeds. The endosperm in monocot species (like cereals) is not consumed by the embryo during seed development as is the case in dicots, but rather develops to the mature seed which contains storage proteins and starch and constitutes the nutritional part of the seed. Although our study did not separate embryo from endosperm tissue, we may assume that the seed of 1–3 DAF in barley is mainly represented by endosperm syncytium with very small contribution from the embryo. Further studies will enable us to draw any parallels between dicots and monocots concerning DNA methylation and endosperm development.
DNA methylation in the 3’downstream region of HvOS2 was found to be increased in 1–3 DAF seeds where the gene is highly induced, than in immature flowers where the gene is downregulated. This resembles the methylation status of the Arabidopsis PHE1 where the 3’downstream region of the expressed paternal allele was found to be methylated, whereas the silenced maternal allele was unmethylated in the same region .
Since HvOS1 was found to be induced upon ABA treatment we investigated whether DNA differential methylation is evidenced in ABA treated and untreated seedlings. Notably, examination of the HvOS1 3’UTR region (region 2) in ABA-treated seedlings revealed higher degree of DNA methylation in leaves of ABA-treated seedlings than in untreated seedlings. This is the first time that DNA methylation has been investigated in MADS-box genes in association to hormonal treatment. It is likely that differential DNA methylation may be associated to differential expression of this barley gene in response to ABA.
Chromatin Immunoprecipitation (ChIP) experiments performed by Greenup et al. (2010) demonstrated the presence of the H3K4me3 gene activation mark, deposited by a Trithorax (TrX) methyltransferase, and found higher H3K4me3 in non vernalized than vernalized tissue, in agreement with HvOS2 decreased expression upon vernalization . They also demonstrated the presence of H3K27me3 mark (the repressive histone modification mark deposited by the PRC2 complex) on the 5’ UTR region of HvOS2, though levels were low both in vernalized and non vernalized tissue. This may be related to the fact that DNA methylation may exclude sites from H3K27me as mentioned in Weinhofer et al. 2010 . The observations by these authors together with our current results suggest that HvOS1 and HvOS2 are targets of epigenetic regulation by histone modifications or DNA methylation.
Epigenetic regulation of endosperm-related genes such as Type I MADS-box genes and PRC2 Polycomb group genes like MEA, through DNA methylation and/or histone modifications, has been well documented in Arabidopsis especially for PHE1, AGL62 and AGL36 and MEA [22, 26, 27, 59]. Our current results on barley Type I-like MADS-box genes point to a role for these genes in endosperm development. Moreover, differences in DNA methylation at different seed developmental stages with different gene expression suggest epigenetic regulation of endosperm-associated genes in barley. In support to this, our group has recently characterized barley epigenetic regulators such as histone deacetylases (HDACs), histone acetyltransferases (HATs), PRC2 Polycomb group genes, Trithorax (TrX) histone methyltransferase and histone demethylase (HDM) genes, and investigated their roles during seed development [30, 44, 60–64].
In the present work two Type I-like MADS-box gene homologues were studied during seed development in barley. HvOS1 and HvOS2 are related to Type I MADS-box genes and exhibit differential expression in particular seed developmental stages in a cultivar-specific manner. HvOS1 is induced by the phytohormone ABA and HvOS2 is detected in specific endosperm sub- compartments. Together these observations point to a role for HvOS1 and HvOS2 in seed development. DNA methylation differences associated with differences in gene expression are suggestive of epigenetic regulation of HvOS1 and HvOS2 in barley, in accordance to Arabidopsis Type I MADS genes. The current study provides further important knowledge for the understanding of gene regulation during seed development in barley and other agronomically important cereals.
We would like to thank Dr Konstantinos Bladenopoulos (NAGREF) for providing seed material and Dr Alexandros Bousios (INA/CERTH) for critical reading of the manuscript. This work was supported by a PENED grant (Ο3Ε_402/2003) and by the EU-COST ACTION 406. Continuous support for the Institute of Agrobiotechnology/CERTH from the General Secretariat of Research and Technology of Greece is also acknowledged.
- Olsen OA: Nuclear endosperm development in cereals and Arabidopsis thaliana. Plant Cell. 2004, 16: S214-S227. 10.1105/tpc.017111.PubMedPubMed CentralView ArticleGoogle Scholar
- Sabelli PA, Larkins BA: The Development of Endosperm in Grasses. Plant Physiol. 2009, 149 (1): 14-26. 10.1104/pp.108.129437.PubMedPubMed CentralView ArticleGoogle Scholar
- Bauer MJ, Fischer RL: Genome demethylation and imprinting in the endosperm. Curr Opin Plant Biol. 2011, 14 (2): 162-167. 10.1016/j.pbi.2011.02.006.PubMedPubMed CentralView ArticleGoogle Scholar
- Feng S, Jacobsen SE, Reik W: Epigenetic Reprogramming in Plant and Animal Development. Science. 2010, 330 (6004): 622-627. 10.1126/science.1190614.PubMedPubMed CentralView ArticleGoogle Scholar
- Gehring M, Choi Y, Fischer RL: Imprinting and seed development. Plant Cell. 2004, 16 (Suppl): S203-213. 10.1105/tpc.017988.PubMedPubMed CentralView ArticleGoogle Scholar
- Huh JH, Bauer MJ, Hsieh TF, Fischer RL: Cellular Programming of Plant Gene Imprinting. Cell. 2008, 132: 735-744. 10.1016/j.cell.2008.02.018.PubMedView ArticleGoogle Scholar
- Parenicova L, de Folter S, Kieffer M, Horner DS, Favalli C, Busscher J, Cook HE, Ingram RM, Kater MM, Davies B, et al: Molecular and Phylogenetic Analyses of the Complete MADS-Box Transcription Factor Family in Arabidopsis: New Openings to the MADS World. Plant Cell. 2003, 15 (7): 1538-1551. 10.1105/tpc.011544.PubMedPubMed CentralView ArticleGoogle Scholar
- Masiero S, Colombo L, Grini PE, Schnittger A, Kater MM: The Emerging Importance of Type I MADS Box Transcription Factors for Plant Reproduction. Plant Cell. 2011, 23 (3): 865-872. 10.1105/tpc.110.081737.PubMedPubMed CentralView ArticleGoogle Scholar
- Arora R, Agarwal P, Ray S, Singh A, Singh V, Tyagi A, Kapoor S: MADS-box gene family in rice: genome-wide identification, organization and expression profiling during reproductive development and stress. BMC Genomics. 2007, 8 (1): 242-10.1186/1471-2164-8-242.PubMedPubMed CentralView ArticleGoogle Scholar
- Gramzow L, Ritz MS, Theiben G: On the origin of MADS-domain transcription factors. Trends Genet. 2010, 26 (4): 149-153. 10.1016/j.tig.2010.01.004.PubMedView ArticleGoogle Scholar
- Alvarez-Buylla ER, Pelaz S, Liljegren SJ, Gold SE, Burgeff C, Ditta GS, Ribas de Pouplana L, Martinez-Castilla L, Yanofsky MF: An ancestral MADS-box gene duplication occurred before the divergence of plants and animals. PNAS. 2000, 97 (10): 5328-5333. 10.1073/pnas.97.10.5328.PubMedPubMed CentralView ArticleGoogle Scholar
- Liu C, Zhang J, Zhang N, Shan H, Su K, Zhang J, Meng Z, Kong H, Chen Z: Interactions among Proteins of Floral MADS-Box Genes in Basal Eudicots: Implications for Evolution of the Regulatory Network for Flower Development. Molecular Biology And Evolution. 2010, 27: 1598-1611. 10.1093/molbev/msq044.PubMedView ArticleGoogle Scholar
- Wu W, Huang X, Cheng J, Li Z, de Folter S, Huang Z, Jiang X, Pang H, Tao S: Conservation and Evolution in and among SRF- and MEF2-Type MADS Domains and Their Binding Sites. Molecular Biology And Evolution. 2011, 28 (1): 501-511. 10.1093/molbev/msq214.PubMedView ArticleGoogle Scholar
- Ferrario S, Immink RGH, Angenent GC: Conservation and diversity in flower land. Curr Opin Plant Biol. 2004, 7 (1): 84-91. 10.1016/j.pbi.2003.11.003.PubMedView ArticleGoogle Scholar
- Sd F, Angenent GC: trans meets cis in MADS science. Trends Plant Sci. 2006, 11 (5): 224-231. 10.1016/j.tplants.2006.03.008.View ArticleGoogle Scholar
- Bemer M, Wolters-Arts M, Grossniklaus U, Angenent GC: The MADS Domain Protein DIANA Acts Together with AGAMOUS-LIKE80 to Specify the Central Cell in Arabidopsis Ovules. Plant Cell. 2008, 20 (8): 2088-2101. 10.1105/tpc.108.058958.PubMedPubMed CentralView ArticleGoogle Scholar
- Dornelas MC, Patreze CM, Angenent GC, Immink RGH: MADS: the missing link between identity and growth?. Trends Plant Sci. 2010, 16 (2): 89-97.PubMedView ArticleGoogle Scholar
- Bemer M, Heijmans K, Airoldi C, Davies B, Angenent GC: An Atlas of Type I MADS Box Gene Expression during Female Gametophyte and Seed Development in Arabidopsis. Plant Physiol. 2010, 154 (1): 287-300. 10.1104/pp.110.160770.PubMedPubMed CentralView ArticleGoogle Scholar
- Kohler C, Hennig L, Spillane C, Pien S, Gruissem W, Grossniklaus U: The Polycomb-group protein MEDEA regulates seed development by controlling expression of the MADS-box gene PHERES1. Genes Dev. 2003, 17: 1540-1553. 10.1101/gad.257403.PubMedPubMed CentralView ArticleGoogle Scholar
- Kohler C, Hennig L, Page D, Gagliardini V, Gruissem W, Grossniklaus U: The Arabidopsis thaliana MEDEA Polycomb group protein controls expression of PHERES 1 by parental imprinting. Nat Genet. 2005, 37: 28-30.PubMedGoogle Scholar
- Makarevich G, Leroy O, Akinci U, Schubert D, Clarenz O, Goodrich J, Grossniklaus U, Kohler C: Different Polycomb group complexes regulate common target genes in Arabidopsis. EMBO Rep. 2006, 7: 947-952. 10.1038/sj.embor.7400760.PubMedPubMed CentralView ArticleGoogle Scholar
- Makarevich G, Villar CBR, Erilova A, Kohler C: Mechanism of PHERES1 imprinting in Arabidopsis. Journal of Cell Sciences. 2008, 121: 906-912. 10.1242/jcs.023077.View ArticleGoogle Scholar
- Colombo M, Masiero S, Vanzulli S, Lardelli P, Kater MM, Colombo L: AGL23, a type I MADS-box gene that controls female gametophyte and embryo development in Arabidopsis. Plant J. 2008, 54 (6): 1037-1048. 10.1111/j.1365-313X.2008.03485.x.PubMedView ArticleGoogle Scholar
- Portereiko MF, Lloyd A, Steffen JG, Punwani JA, Otsuga D, Drews GN: AGL80 Is Required for Central Cell and Endosperm Development in Arabidopsis. Plant Cell. 2006, 18 (8): 1862-1872. 10.1105/tpc.106.040824.PubMedPubMed CentralView ArticleGoogle Scholar
- Steffen JG, Kang I-H, Portereiko MF, Lloyd A, Drews GN: AGL61 Interacts with AGL80 and Is Required for Central Cell Development in Arabidopsis. Plant Physiol. 2008, 148 (1): 259-268. 10.1104/pp.108.119404.PubMedPubMed CentralView ArticleGoogle Scholar
- Kang I-H, Steffen JG, Portereiko MF, Lloyd A, Drews GN: The AGL62 MADS Domain Protein Regulates Cellularization during Endosperm Development in Arabidopsis. The Plant Cell Online. 2008, 20 (3): 635-647. 10.1105/tpc.107.055137.View ArticleGoogle Scholar
- Shirzadi R, Andersen ED, Bjerkan KN, Gloeckle BM, Heese M, Ungru A, Winge P, Koncz C, Aalen RB, Schnittger A, et al: Genome-Wide Transcript Profiling of Endosperm without Paternal Contribution Identifies Parent-of-Origin-Dependent Regulation of AGAMOUS-LIKE36. PLoS Genet. 2011, 7 (2): e1001303-10.1371/journal.pgen.1001303.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhao T, Ni Z, Dai Y, Yao Y, Nie X, Sun Q: Characterization and expression of 42 MADS-box genes in wheat (Triticum aestivum). Mol Genet Genomics. 2006, 276 (4): 334-350. 10.1007/s00438-006-0147-3.PubMedView ArticleGoogle Scholar
- Ishikawa R, Ohnishi T, Kinoshita Y, Eiguchi M, Kurata N, Kinoshita T: Rice interspecies hybrids show precocious or delayed developmental transitions in the endosperm without change to the rate of syncytial nuclear division. Plant J. 2011, 65: 798-806. 10.1111/j.1365-313X.2010.04466.x.PubMedView ArticleGoogle Scholar
- Kapazoglou A, Tondelli A, Papaefthimiou D, Ampatzidou H, Francia E, Stanca M, Bladenopoulos K, Tsaftaris A: Epigenetic chromatin modifiers in barley: IV. The study of barley Polycomb group (PcG) genes during seed development and in response to external ABA. BMC Plant Biol. 2010, 10 (1): 73-10.1186/1471-2229-10-73.PubMedPubMed CentralView ArticleGoogle Scholar
- Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, Gavin OL, Gunasekaran P, Ceric G, Forslund K, et al: The Pfam protein families database. Nucleic Acids Res. 2010, 38 (suppl 1): D211-D222.PubMedPubMed CentralView ArticleGoogle Scholar
- Overbeek R, Fonstein M, De Souza M, Pusch GD, Maltsev N: The use of gene clusters to infer functional coupling. PNAS. 1999, 96 (6): 2896-2901. 10.1073/pnas.96.6.2896.PubMedPubMed CentralView ArticleGoogle Scholar
- Katoh K, Kuma K-i, Toh H, Miyata T: MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucl Acids Res. 2005, 33 (2): 511-518. 10.1093/nar/gki198.PubMedPubMed CentralView ArticleGoogle Scholar
- Waterhouse AM, Procter JB, Martin DMA, Clamp M, Barton GJ: Jalview Version 2- a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009, 25 (9): 1189-1191. 10.1093/bioinformatics/btp033.PubMedPubMed CentralView ArticleGoogle Scholar
- Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 2004, 5 (2): 150-163. 10.1093/bib/5.2.150.PubMedView ArticleGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Bio Evol. 1987, 4 (4): 406-425.Google Scholar
- Close T, Bhat P, Lonardi S, Wu Y, Rostoks N, Ramsay L, Druka A, Stein N, Svensson J, Wanamaker S, et al: Development and implementation of high-throughput SNP genotyping in barley. BMC Genomics. 2009, 10 (1): 582-10.1186/1471-2164-10-582.PubMedPubMed CentralView ArticleGoogle Scholar
- Szucs P, Blake VC, Bhat PR, Chao S, Close TJ, Cuesta-Marcos A, Muehlbauer GJ, Ramsay L, Waugh R, Hayes PM: An Integrated Resource for Barley Linkage Map and Malting Quality QTL Alignment. The Plant Genome. 2009, 2 (2): 134-140. 10.3835/plantgenome2008.01.0005.View ArticleGoogle Scholar
- Shen Q, Chen C-N, Brands A, Pan S-M, Tuan-Hua DH: The stress- and abscisic acid-induced barley gene HVA22: developmental regulation and homologues in diverse organisms. Plant Mol Biol. 2001, 45 (3): 327-340. 10.1023/A:1006460231978.PubMedView ArticleGoogle Scholar
- Walia H, Wilson C, Condamine P, Liu X, Ismail A, Close TJ: Large-scale expression profiling and physiological characterization of jasmonic acid-mediated adaptation of barley to salinity stress. Plant Cell Environ. 2007, 30: 410-421. 10.1111/j.1365-3040.2006.01628.x.PubMedView ArticleGoogle Scholar
- Kouri ED, Labrou NE, Garbis SD, Kalliampakou KI, Stedel C, Dimou M, Udvardi MK, Katinakis P, Flemetakis E: Molecular and Biochemical Characterization of the Parvulin-Type PPIases in Lotus japonicus. Plant Physiol. 2009, 150 (3): 1160-1173. 10.1104/pp.108.132415.PubMedPubMed CentralView ArticleGoogle Scholar
- Scheres B, Van De Wiel C, Zalensky A, Horvath B, Spaink H, Van Eck H, Zwartkruis F, Wolters A-M, Gloudemans T, Van Kammen A, et al: The ENOD12 gene product is involved in the infection process during the pea-rhizobium interaction. Cell. 1990, 60 (2): 281-294. 10.1016/0092-8674(90)90743-X.PubMedView ArticleGoogle Scholar
- Greenup AG, Sasani S, Oliver SN, Talbot MJ, Dennis ES, Hemming MN, Trevaskis B: ODDSOC2 is a MADS box floral repressor that is down-regulated by vernalization in temperate cereals. Plant Physiol. 2010, 153: 1062-1073. 10.1104/pp.109.152488.PubMedPubMed CentralView ArticleGoogle Scholar
- Papaefthimiou D, Kapazoglou A, Tsaftaris A: Cloning and characterization of SOC1 homologues in barley (Hordeum vulgare) and their expression during seed development and in response to vernalization. Physiologia Plantarum. (in press)
- Darzentas N, Bousios A, Apostolidou V, Tsaftaris AS: MASiVE: Mapping and Analysis of SireVirus Elements in plant genome sequences. Bioinformatics. 2010, 26 (19): 2452-2454. 10.1093/bioinformatics/btq454.PubMedView ArticleGoogle Scholar
- Bousios A, Darzentas N, Tsaftaris A, Pearce SR: Highly conserved motifs in non-coding regions of Sirevirus retrotransposons: the key for their pattern of distribution within and across plants?. Bmc Genomics. 2010, 11:Google Scholar
- Rossi V, Locatelli S, Varotto S, Donn G, Pirona R, Henderson DA, Hartings H, Motto M: Maize histone deacetylase hda101 is involved in plant development, gene transcription, and sequence-specific modulation of histone modification of genes and repeats. Plant Cell. 2007, 19 (4): 1145-1162. 10.1105/tpc.106.042549.PubMedPubMed CentralView ArticleGoogle Scholar
- Day RC, Herridge RP, Ambrose BA, Macknight RC: Transcriptome Analysis of Proliferating Arabidopsis Endosperm Reveals Biological Implications for the Control of Syncytial Division, Cytokinin Signaling, and Gene Expression Regulation. Plant Physiol. 2008, 148 (4): 1964-1984. 10.1104/pp.108.128108.PubMedPubMed CentralView ArticleGoogle Scholar
- Agarwal P, Kapoor S, Tyagi AK: Transcription factors regulating the progression of monocot and dicot seed development. Bioessays. 2011, 33 (3): 189-202. 10.1002/bies.201000107.PubMedView ArticleGoogle Scholar
- Wicker T, Mayer KFX, Gundlach H, Martis M, Steuernagel B, Scholz U, Simkova H, Kubalakova M, Choulet F, Taudien S, et al: Frequent Gene Movement and Pseudogene Evolution Is Common to the Large and Complex Genomes of Wheat, Barley, and Their Relatives. Plant Cell. 2011, 23 (5): 1706-1718. 10.1105/tpc.111.086629.PubMedPubMed CentralView ArticleGoogle Scholar
- Mayer KFX, Martis M, Hedley PE, Simkova H, Liu H, Morris JA, Steuernagel B, Taudien S, Roessner S, Gundlach H, et al: Unlocking the Barley Genome by Chromosomal and Comparative Genomics. Plant Cell. 2011, 23 (4): 1249-1263. 10.1105/tpc.110.082537.PubMedPubMed CentralView ArticleGoogle Scholar
- Bousios A, Kourmpetis YA, Pavlidis P, Minga E, Tsaftaris A, Darzentas N: The turbulent life of Sirevirus retrotransposons and the evolution of the maize genome: more than ten thousand elements tell the story. Plant J. 2012, 69: 475-488. 10.1111/j.1365-313X.2011.04806.x.PubMedView ArticleGoogle Scholar
- Hsieh TF, Ibarra CA, Silva P, Zemach A, Eshed-Williams L, Fischer RL, Zilberman D: Genome-wide demethylation of Arabidopsis endosperm. Science (New York, NY. 2009, 324 (5933): 1451-1454. 10.1126/science.1172417.View ArticleGoogle Scholar
- Jullien PE, Berger F: Gamete-specific epigenetic mechanisms shape genomic imprinting. Curr Opin Plant Biol. 2009, 12 (5): 637-642. 10.1016/j.pbi.2009.07.004.PubMedView ArticleGoogle Scholar
- Gehring M, Bubb KL, Henikoff S: Extensive demethylation of repetitive elements during seed development underlies gene imprinting. Science New York, NY. 2009, 324 (5933): 1447-1451. 10.1126/science.1171609.View ArticleGoogle Scholar
- Mosher RA, Melnyk CW: siRNAs and DNA methylation: seedy epigenetics. Trends Plant Sci. 2010, 15 (4): 204-210. 10.1016/j.tplants.2010.01.002.PubMedView ArticleGoogle Scholar
- Winfield M, Lu C, Wilson I, Coghill J, Edwards K: Cold- and light-induced changes in the transcriptome of wheat leading to phase transition from vegetative to reproductive growth. BMC Plant Biol. 2009, 9 (1): 55-10.1186/1471-2229-9-55.PubMedPubMed CentralView ArticleGoogle Scholar
- Weinhofer I, Hehenberger E, Roszak P, Hennig L, Kohhler C: H3K27me3 Profiling of the Endosperm Implies Exclusion of Polycomb Group Protein Targeting by DNA Methylation. PLoS Genet. 2010, 6 (10): e1001152-10.1371/journal.pgen.1001152.PubMedPubMed CentralView ArticleGoogle Scholar
- Hehenberger E, Kradolfer D, Kohler C: Endosperm cellularization defines an important developmental transition for embryo development. Development. 2012, 139 (11): 2031-2039. 10.1242/dev.077057.PubMedView ArticleGoogle Scholar
- Demetriou K, Kapazoglou A, Bladenopoulos K, Tsaftaris A: Epigenetic Chromatin Modifiers in Barley: II. Characterization and Expression Analysis of the HDA1 Family of Barley Histone Deacetylases During Development and in Response to Jasmonic Acid. Plant Molecular Biology Reporter. 2010, 28 (1): 9-21. 10.1007/s11105-009-0121-4.View ArticleGoogle Scholar
- Demetriou K, Kapazolgou A, Tondelli A, Francia E, Stanca MS, Bladenopoulos K, Tsaftaris AS: Epigenetic chromatin modifiers in barley: I. Cloning, mapping and expression analysis of the plant specific HD2 family of histone deacetylases from barley, during seed development and after hormonal treatment. Physiol Plant. 2009, 136: 358-368. 10.1111/j.1399-3054.2009.01236.x.PubMedView ArticleGoogle Scholar
- Papaefthimiou D, Tsaftaris AS: Characterization of a drought inducible trithorax-like H3K4 methyltransferase from barley. Biologia Plantarum. 2012, 56 (4): 683-692. 10.1007/s10535-012-0125-z.View ArticleGoogle Scholar
- Papaefthimiou D, Likotrafiti E, Kapazoglou A, Bladenopoulos K, Tsaftaris A: Epigenetic chromatin modifiers in barley: III. Isolation and characterization of the barley GNAT-MYST family of histone acetyltransferases and responses to exogenous ABA. Plant Physiol Biochem. 2010, 48 (2–3): 98-107.PubMedView ArticleGoogle Scholar
- Papaefthimiou D, Tsaftaris A: Significant induction by drought of HvPKDM7–1, a gene encoding a jumonji-like histone demethylase homologue in barley (H. vulgare). Acta Physiologiae Plantarum. 2012, 34: 1187-1198. 10.1007/s11738-011-0915-5.View ArticleGoogle Scholar