Epigenetic profiling of heterochromatic satellite DNA

Sugar beet (Beta vulgaris) chromosomes consist of large heterochromatic blocks in pericentromeric, centromeric, and intercalary regions comprised of two different highly abundant DNA satellite families. To investigate DNA methylation at single base resolution at heterochromatic regions, we applied a method for strand-specific bisulfite sequencing of more than 1,000 satellite monomers followed by statistical analyses. As a result, we uncovered diversity in the distribution of different methylation patterns in both satellite families. Heavily methylated CG and CHG (H=A, T, or C) sites occur more frequently in intercalary heterochromatin, while CHH sites, with the exception of CAA, are only sparsely methylated, in both intercalary and pericentromeric/centromeric heterochromatin. We show that the difference in DNA methylation intensity is correlated to unequal distribution of heterochromatic histone H3 methylation marks. While clusters of H3K9me2 were absent from pericentromeric heterochromatin and restricted only to intercalary heterochromatic regions, H3K9me1 and H3K27me1 were observed in all types of heterochromatin. By sequencing of a small RNA library consisting of 6.76 million small RNAs, we identified small interfering RNAs (siRNAs) of 24 nucleotides in size which originated from both strands of the satellite DNAs. We hypothesize an involvement of these siRNAs in the regulation of DNA and histone methylation for maintaining heterochromatin.     

Endosperm-based incompatibilities in hybrid monkeyflowers

Endosperm is an angiosperm innovation central to their reproduction whose development, and thus seed viability, is controlled by genomic imprinting, where expression from certain genes is parent-specific. Unsuccessful imprinting has been linked to failed inter-specific and inter-ploidy hybridization. Despite their importance in plant speciation, the underlying mechanisms behind these endosperm-based barriers remain poorly understood. Here, we describe one such barrier between diploid Mimulus guttatus and tetraploid Mimulus luteus. The two parents differ in endosperm DNA methylation, expression dynamics, and imprinted genes. Hybrid seeds suffer from underdeveloped endosperm, reducing viability, or arrested endosperm and seed abortion when M. guttatus or M. luteus is seed parent, respectively, and transgressive methylation and expression patterns emerge. The two inherited M. luteus subgenomes, genetically distinct but epigenetically similar, are expressionally dominant over the M. guttatus genome in hybrid embryos and especially their endosperm, where paternal imprints are perturbed. In aborted seeds, de novo methylation is inhibited, potentially owing to incompatible paternal instructions of imbalanced dosage from M. guttatus imprints. We suggest that diverged epigenetic/regulatory landscapes between parental genomes induce epigenetic repatterning and global shifts in expression, which, in endosperm, may uniquely facilitate incompatible interactions between divergent imprinting schemes, potentially driving rapid barriers.

Diverged epigenetic/regulatory landscapes between parental genomes result in epigenetic repatterning in hybrids that drive global shifts in endosperm gene expression patterns.     

Lysine 27 dimethylation of Drosophila linker histone dH1 contributes to heterochromatin organization independently of H3K9 methylation

Post-translational modifications (PTMs) of core histones are important epigenetic determinants that correlate with functional chromatin states. However, despite multiple linker histone H1s PTMs have been identified, little is known about their genomic distribution and contribution to the epigenetic regulation of chromatin. Here, we address this question in Drosophila that encodes a single somatic linker histone, dH1. We previously reported that dH1 is dimethylated at K27 (dH1K27me2). Here, we show that dH1K27me2 is a major PTM of Drosophila heterochromatin. At mitosis, dH1K27me2 accumulates at pericentromeric heterochromatin, while, in interphase, it is also detected at intercalary heterochromatin. ChIPseq experiments show that >98% of dH1K27me2 enriched regions map to heterochromatic repetitive DNA elements, including transposable elements, simple DNA repeats and satellite DNAs. Moreover, expression of a mutated dH1K27A form, which impairs dH1K27me2, alters heterochromatin organization, upregulates expression of heterochromatic transposable elements and results in the accumulation of RNA:DNA hybrids (R-loops) in heterochromatin, without affecting H3K9 methylation and HP1a binding. The pattern of dH1K27me2 is H3K9 methylation independent, as it is equally detected in flies carrying a H3K9R mutation, and is not affected by depletion of Su(var)3–9, HP1a or Su(var)4–20. Altogether these results suggest that dH1K27me2 contributes to heterochromatin organization independently of H3K9 methylation.     

Organization and evolution of heterochromatin in malaria mosquitoes

The African malaria mosquito Anopheles gambiae was the first disease vector chosen for genome sequencing. Although its genome assembly has been facilitated by physical mapping, large gaps still pose a serious problem for accurate annotation and genome analysis. The majority of the gaps are located in regions of pericentromeric and intercalary heterochromatin. Genomic analysis has identified protein-coding genes and various classes of repetitive elements in the Anopheles heterochromatin. Molecular and cytogenetic studies have demonstrated that heterochromatin is a structurally heterogeneous and rapidly evolving part of the malaria mosquito genome.     

Bypassing reproductive barriers in hybrid seeds using chemically induced epimutagenesis

The triploid block, which prevents interploidy hybridizations in flowering plants, is characterized by a failure in endosperm development, arrest in embryogenesis, and seed collapse. Many genetic components of triploid seed lethality have been successfully identified in the model plant Arabidopsis thaliana, most notably the paternally expressed genes (PEGs), which are upregulated in tetraploid endosperm with paternal excess. Previous studies have shown that the paternal epigenome is a key determinant of the triploid block response, as the loss of DNA methylation in diploid pollen suppresses the triploid block almost completely. Here, we demonstrate that triploid seed collapse is bypassed in Arabidopsis plants treated with the DNA methyltransferase inhibitor 5-Azacytidine during seed germination and early growth. We identified strong suppressor lines showing stable transgenerational inheritance of hypomethylation in the CG context, as well as normalized expression of PEGs in triploid seeds. Importantly, differentially methylated loci segregate in the progeny of “epimutagenized” plants, which may allow epialleles involved in the triploid block response to be identified in future studies. Finally, we demonstrate that chemically induced epimutagenesis facilitates hybridization between different Capsella species, thus potentially emerging as a strategy for producing triploids and interspecific hybrids with high agronomic interest.

Genome-wide loss of DNA methylation induced by 5-Azacytidine allows interploidy and interspecific hybridization barriers to be bypassed in Arabidopsis and Capsella.     

Mapping of days to flower and seed yield in spring oilseed <Emphasis Type="Italic">Brassica napus</Emphasis> carrying genome content introgressed from <Emphasis Type="Italic">Brassica oleracea</Emphasis>

Earliness of flowering and maturity and high seed yield are important objectives of breeding spring Brassica napus canola. Previously, we have introgressed earliness of flowering from Brassica oleracea into spring B. napus canola through interspecific crossing between these two species. In this paper, we report quantitative trait locus (QTL) mapping of days to flower and seed yield by use of publicly available markers and markers designed based on flowering time genes and a doubled haploid population, derived from crossing of the spring canola parent and an early flowering line developed from a B. napus × B. oleracea cross, tested in nine field trials for over 5 years. Five genomic regions associated with days to flower were identified on C1, C2, C3, and C6 of which the single QTL of C1 was detected in all trials; in all cases, the allele introgressed from B. oleracea reduced the number of days to flower. BLASTn search in the Brassica genomes located the physical position of the QTL markers and identified putative flowering time genes in these regions. In the case of seed yield, ten QTL from eight linkage groups were detected; however, none could be consistently detected in all trials. The QTL region of C1 associated with days to flower did not show significant association with seed yield in more than 80% of the field trials; however, in a single trial, the allele introgressed from B. oleracea exerted a negative effect on seed yield. Thus, the genomic regions and molecular markers identified in this research could potentially be used in breeding for the development of early flowering B. napus canola cultivars without affecting seed yield in a majority of the environments.     

Young but not relatively old retrotransposons are preferentially located in gene‐rich euchromatic regions in tomato (Solanum lycopersicum) plants

Long terminal repeat (LTR) retrotransposons are the major DNA components of flowering plants. They are generally enriched in pericentromeric heterochromatin regions of their host genomes, which could result from the preferential insertion of LTR retrotransposons and the low effectiveness of purifying selection in these regions. To estimate the relative importance of the actions of these two factors on their distribution pattern, the LTR retrotransposons in Solanum lycopersicum (tomato) plants were characterized at the genome level, and then the distribution of young elements was compared with that of relatively old elements. The current data show that old elements are mainly located in recombination‐suppressed heterochromatin regions, and that young elements are preferentially located in the gene‐rich euchromatic regions. Further analysis showed a negative correlation between the insertion time of LTR retrotransposons and the recombination rate. The data also showed there to be more solo LTRs in genic regions than in intergenic regions or in regions close to genes. These observations indicate that, unlike in many other plant genomes, the current LTR retrotransposons in tomatoes have a tendency to be preferentially located into euchromatic regions, probably caused by their severe suppression of activities in heterochromatic regions. These elements are apt to be maintained in heterochromatin regions, probably as a consequence of the pericentromeric effect in tomatoes. These results also indicate that local recombination rates and intensities of purifying selection in different genomic regions are largely responsible for structural variation and non‐random distribution of LTR retrotransposons in tomato plants.     

A genome-wide association study of seed protein and oil content in soybean

Background

Association analysis is an alternative to conventional family-based methods to detect the location of gene(s) or quantitative trait loci (QTL) and provides relatively high resolution in terms of defining the genome position of a gene or QTL. Seed protein and oil concentration are quantitative traits which are determined by the interaction among many genes with small to moderate genetic effects and their interaction with the environment. In this study, a genome-wide association study (GWAS) was performed to identify quantitative trait loci (QTL) controlling seed protein and oil concentration in 298 soybean germplasm accessions exhibiting a wide range of seed protein and oil content.

Results

A total of 55,159 single nucleotide polymorphisms (SNPs) were genotyped using various methods including Illumina Infinium and GoldenGate assays and 31,954 markers with minor allele frequency >0.10 were used to estimate linkage disequilibrium (LD) in heterochromatic and euchromatic regions. In euchromatic regions, the mean LD (r ² ) rapidly declined to 0.2 within 360 Kbp, whereas the mean LD declined to 0.2 at 9,600 Kbp in heterochromatic regions. The GWAS results identified 40 SNPs in 17 different genomic regions significantly associated with seed protein. Of these, the five SNPs with the highest associations and seven adjacent SNPs were located in the 27.6-30.0 Mbp region of Gm20. A major seed protein QTL has been previously mapped to the same location and potential candidate genes have recently been identified in this region. The GWAS results also detected 25 SNPs in 13 different genomic regions associated with seed oil. Of these markers, seven SNPs had a significant association with both protein and oil.

Conclusions

This research indicated that GWAS not only identified most of the previously reported QTL controlling seed protein and oil, but also resulted in narrower genomic regions than the regions reported as containing these QTL. The narrower GWAS-defined genome regions will allow more precise marker-assisted allele selection and will expedite positional cloning of the causal gene(s).     

FISH to meiotic pachytene chromosomes of tomato locates the root-knot nematode resistance gene Mi-1 and the acid phosphatase gene Aps-1 near the junction of euchromatin and pericentromeric heterochromatin of chromosome arms 6S and 6L, respectively

 The root-knot nematode resistance gene Mi-1 in tomato has long been thought to be located in the pericentromeric heterochromatin region of the long arm of chromosome 6 because of its very tight genetic linkage (approx. 1 cM) to the markers Aps-1 (Acid phosphatase 1) and yv (yellow virescent). Using Mi-BAC clones and an Aps-1 YAC clone in fluorescence in situ hybridisation (FISH) to pachytene chromosomes we now provide direct physical evidence showing that Mi-1 is located at the border of the euchromatin and heterochromatin regions in the short arm (6S) and Aps-1 in the pericentromeric heterochromatin of the long arm (6L) close to the euchromatin. Taking into account both the estimated DNA content of hetero- and euchromatin regions and the compactness of the tomato chromosomes at pachytene (2 Mb/μm), our data suggest that Mi-1 and Aps-1 are at least 40 Mb apart, a base pair-to-centiMorgan relationship that is more than 50-fold higher than the average value of 750 kb/cM of the tomato genome. An integrated cytogenetic-molecular map of chromosome 6 is presented that provides a framework for physical mapping. Received: 24 July 1998 / Accepted: 14 August 1998     

Locus- and Site-Specific DNA Methylation of 19 kDa Zein Genes in Maize

An interesting question in maize development is why only a single zein gene is highly expressed in each of the 19-kDa zein gene clusters (A and B types), z1A2-1 and z1B4, in the immature endosperm. For instance, epigenetic marks could provide a structural difference. Therefore, we investigated the DNA methylation of the arrays of gene copies in both promoter and gene body regions of leaf (non-expressing tissue as a control), normal endosperm, and cultured endosperm. Although we could show that expressed genes have much lower methylation levels in promoter regions than silent ones in both leaf and normal endosperm, there was surprisingly also a difference in the pattern of the z1A and z1B gene clusters. The expression of z1B gene is suppressed by increased DNA methylation and activated with reduced DNA methylation, whereas z1A gene expression is not. DNA methylation in gene coding regions is higher in leaf than in endosperm, whereas no significant difference is observed in gene bodies between expressed and non-expressed gene copies. A median CHG methylation (25–30%) appears to be optimal for gene expression. Moreover, tissue-cultured endosperm can reset the DNA methylation pattern and tissue-specific gene expression. These results reveal that DNA methylation changes of the 19-kDa zein genes is subject to plant development and tissue culture treatment, but varies in different chromosomal locations, indicating that DNA methylation changes do not apply to gene expression in a uniform fashion. Because tissue culture is used to produce transgenic plants, these studies provide new insights into variation of gene expression of integrated sequences.     

Structure and DNA methylation pattern of partially heterochromatinised endosperm nuclei in Gagea lutea (Liliaceae)

Pentaploid endosperm nuclei in certain Gagea species exhibit large masses of sticky and dense chromatin, not observed in somatic nuclei. These heterochromatin masses most probably stem from the triploid chalasal polar nucleus of the embryo sac, thus representing an example of facultative heterochromatinisation in plants. In the present investigation, we studied the nuclei in Gagea lutea (L.) Ker-Gawl. endosperm tissue. The position of the heterochromatin in interphase nuclei was observed by confocal laser scanning microscopy (CLSM) and the DNA methylation status of the euchromatin and heterochromatin was analysed by immunolabelling with an antibody raised against 5-methylcytosine (anti-5-mC). In young endosperms, heterochromatin was relatively dispersed, occupying some peripheral and inner parts of the nuclei. In a later endosperm development, the nuclei became smaller and more pycnotic, and the heterochromatin masses were placed predominantly near the nuclear periphery. The distribution of anti-5-mC labelling on the heterochromatic regions was unequal: some parts appeared hypermethylated while other parts were, like the euchromatin, not labelled. During mitosis, the labelling intensity of all the chromosomes was approximately the same, thus indicating that there are no cytologically detectable methylation differences among the individual sets of chromosomes. However, differences in the anti-5-mC signal intensity along individual chromosomes were observed, resulting in banding patterns with highly positive bands apparently representing constitutive heterochromatic regions. From these results it is obvious that facultative heterochromatinisation, in contrast to constitutive heterochromatinisation, need not be strictly accompanied by a prominent DNA hypermethylation. Received: 24 April 1997 / Accepted: 28 July 1997     

Drosophila Dalmatian combines sororin and shugoshin roles in establishment and protection of cohesion

Sister chromatid cohesion is crucial to ensure chromosome bi‐orientation and equal chromosome segregation. Cohesin removal via mitotic kinases and Wapl has to be prevented in pericentromeric regions in order to protect cohesion until metaphase, but the mechanisms of mitotic cohesion protection remain elusive in Drosophila. Here, we show that dalmatian (Dmt), an ortholog of the vertebrate cohesin‐associated protein sororin, is required for protection of mitotic cohesion in flies. Dmt is essential for cohesion establishment during interphase and is enriched on pericentromeric heterochromatin. Dmt is recruited through direct association with heterochromatin protein‐1 (HP1), and this interaction is required for cohesion. During mitosis, Dmt interdependently recruits protein phosphatase 2A (PP2A) to pericentromeric regions, and PP2A binding is required for Dmt to protect cohesion. Intriguingly, Dmt is sufficient to protect cohesion upon heterologous expression in human cells. Our findings of a hybrid system, in which Dmt exerts both sororin‐like establishment functions and shugoshin‐like heterochromatin‐based protection roles, provide clues to the evolutionary modulation of eukaryotic cohesion regulation systems.     

Role for rodent Smc6 in pericentromeric heterochromatin domains during spermatogonial differentiation and meiosis

Chromatin structure and function are for a large part determined by the six members of the structural maintenance of chromosomes (SMC) protein family, which form three heterodimeric complexes: Smc1/3 (cohesin), Smc2/4 (condensin) and Smc5/6. Each complex has distinct and important roles in chromatin dynamics, gene expression and differentiation. In yeast and Drosophila, Smc6 is involved in recombinational repair, restarting collapsed replication forks and prevention of recombination in repetitive sequences such as rDNA and pericentromeric heterochromatin. Although such DNA damage control mechanisms, as well as highly dynamic changes in chromatin composition and function, are essential for gametogenesis, knowledge on Smc6 function in mammalian systems is limited. We therefore have investigated the role of Smc6 during mammalian spermatogonial differentiation, meiosis and subsequent spermiogenesis. We found that, during mouse spermatogenesis, Smc6 functions as part of meiotic pericentromeric heterochromatin domains that are initiated when differentiating spermatogonia become irreversibly committed toward meiosis. To our knowledge, we are the first to provide insight into how commitment toward meiosis alters chromatin structure and dynamics, thereby setting apart differentiating spermatogonia from the undifferentiated spermatogonia, including the spermatogonial stem cells. Interestingly, Smc6 is not essential for spermatogonial mitosis, whereas Smc6-negative meiotic cells appear unable to finish their first meiotic division. Importantly, during meiosis, we find that DNA repair or recombination sites, marked by γH2AX or Rad51 respectively, do not co-localize with the pericentromeric heterochromatin domains where Smc6 is located. Considering the repetitive nature of these domains and that Smc6 has been previously shown to prevent recombination in repetitive sequences, we hypothesize that Smc6 has a role in the prevention of aberrant recombination events between pericentromeric regions during the first meiotic prophase that would otherwise cause chromosomal aberrations leading to apoptosis, meiotic arrest or aneuploidies.