Transposons, tandem repeats, and the silencing of imprinted genes

Heterochromatin repeat organization at an individual level: Rex1BD and the 14-3-3 protein coordinate to shape the epigenetic landscape within heterochromatin repeats

In eukaryotic cells, heterochromatin is typically composed of tandem DNA repeats and plays crucial roles in gene expression and genome stability. It has been reported that silencing at individual units within tandem heterochromatin repeats exhibits a position-dependent variation. However, how the heterochromatin is organized at an individual repeat level remains poorly understood. Using a novel genetic approach, our recent study identified a conserved protein Rex1BD required for position-dependent silencing within heterochromatin repeats. We further revealed that Rex1BD interacts with the 14-3-3 protein to regulate heterochromatin silencing by linking RNAi and HDAC pathways. In this review, we discuss how Rex1BD and the 14-3-3 protein coordinate to modulate heterochromatin organization at the individual repeat level, and comment on the biological significance of the position-dependent effect in heterochromatin repeats. We also identify the knowledge gaps that still need to be unveiled in the field.

H3K9 methylation regulates heterochromatin silencing through incoherent feedforward loops

Histone H3 lysine-9 methylation (H3K9me) is a hallmark of the condensed and transcriptionally silent heterochromatin. It remains unclear how H3K9me controls transcription silencing and how cells delimit H3K9me domains to avoid silencing essential genes. Here, using Arabidopsis genetic systems that induce H3K9me2 in genes and transposons de novo, we show that H3K9me2 accumulation paradoxically also causes the deposition of the euchromatic mark H3K36me3 by a SET domain methyltransferase, ASHH3. ASHH3-induced H3K36me3 confers anti-silencing by preventing the demethylation of H3K4me1 by LDL2, which mediates transcriptional silencing downstream of H3K9me2. These results demonstrate that H3K9me2 not only facilitates but orchestrates silencing by actuating antagonistic silencing and anti-silencing pathways, providing insights into the molecular basis underlying proper partitioning of chromatin domains and the creation of metastable epigenetic variation.

Rex1BD and the 14-3-3 protein control heterochromatin organization at tandem repeats by linking RNAi and HDAC

Tandem DNA repeats are often organized into heterochromatin that is crucial for genome organization and stability. Recent studies revealed that individual repeats within tandem DNA repeats can behave very differently. How DNA repeats are assembled into distinct heterochromatin structures remains poorly understood. Here, we developed a genome-wide genetic screen using a reporter gene at different units in a repeat array. This screen led to identification of a conserved protein Rex1BD required for heterochromatin silencing. Our structural analysis revealed that Rex1BD forms a four-helix bundle structure with a distinct charged electrostatic surface. Mechanistically, Rex1BD facilitates the recruitment of Clr6 histone deacetylase (HDAC) by interacting with histones. Interestingly, Rex1BD also interacts with the 14-3-3 protein Rad25, which is responsible for recruiting the RITS (RNA-induced transcriptional silencing) complex to DNA repeats. Our results suggest that coordinated action of Rex1BD and Rad25 mediates formation of distinct heterochromatin structure at DNA repeats via linking RNAi and HDAC pathways.

Transcriptional activity of PIF and Pong-like Class II transposable elements in Triticeae

Background

Transposable elements are major contributors to genome size and variability, accounting for approximately 70–80% of the maize, barley, and wheat genomes. PIF and Pong-like elements belong to two closely-related element families within the PIF/Harbinger superfamily of Class II (DNA) transposons. Both elements contain two open reading frames; one encodes a transposase (ORF2) that catalyzes transposition of the functional elements and their related non-autonomous elements, while the function of the second is still debated. In this work, we surveyed for PIF- and Pong-related transcriptional activity in 13 diploid Triticeae species, all of which have been previously shown to harbor extensive within-genome diversity of both groups of elements.

Results

The results revealed that PIF elements have considerable transcriptional activity in Triticeae, suggesting that they can escape the initial levels of plant cell control and are regulated at the post-transcriptional level. Phylogenetic analysis of 156 PIF cDNA transposase fragments along with 240 genomic partial transposase sequences showed that most, if not all, PIF clades are transcriptionally competent, and that multiple transposases coexisting within a single genome have the potential to act simultaneously. In contrast, we did not detect any transcriptional activity of Pong elements in any sample.

Conclusions

The lack of Pong element transcription shows that even closely related transposon families can exhibit wide variation in their transposase transcriptional activity within the same genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1028-6) contains supplementary material, which is available to authorized users.

Are all repeats created equal? Understanding DNA repeats at an individual level

Repetitive DNA sequences, comprising up to 50 % of the genome in all eukaryotes, play important roles in a wide range of cellular functions, such as transcriptional regulation, genome stability, and cellular differentiation. However, due to technical difficulties in differentiating their sequences, DNA repeats remain one of the most mysterious parts of eukaryotic genomes. Key questions, such as how repetitive entities behave at individual level and how the internal architecture of these repeats is organized, are still poorly understood. Recent advances from our group reveal unexpected position-dependent variation within tandem DNA repeats in fission yeast. Despite sharing identical DNA sequences, the peri-centromeric repeats are organized into diverse epigenetic states and chromatin structures. We demonstrate that this position-dependent variation requires key heterochromatin factors and condensin. Our works further suggest that the peri-centromeric repeats are organized into distinct higher order structures that ensure a proper positioning of CENP-A, the centromere-specific histone H3 variant, to centromeres. These most recent developments offer insights into the mechanisms underlying the position effect within tandem DNA arrays, and have broad implications in the field of epigenetics and chromatin biology.

Condensin Promotes Position Effects within Tandem DNA Repeats via the RITS Complex

Tandem repetitive DNA is highly abundant in eukaryotic genomes and contributes to transcription control and genome stability. However, how the individual sequences within tandem repeats behave remains largely unknown. Here we develop a collection of fission yeast strains with a reporter gene inserted at different units in a tandem repeat array. We show that, contrary to what is usually assumed, transcriptional silencing and replication timing among the individual repeats differ significantly. RNAi-mediated H3K9 methylation is essential for the silencing position effect. A short hairpin RNA of ura4(+) induces silencing in trans within the tandem array in a position-dependent manner. Importantly, the position effect depends on the condensin subunit, cut3(+). Cut3 promotes the position effect via interaction with the RNA-induced transcriptional silencing (RITS) complex. This study reveals variations in silencing within tandem DNA repeats and provides mechanistic insights into how DNA repeats at the individual level are regulated.

Control of Paternally Expressed Imprinted UPWARD CURLY LEAF1, a Gene Encoding an F-Box Protein That Regulates CURLY LEAF Polycomb Protein, in the Arabidopsis Endosperm

Genomic imprinting, an epigenetic process in mammals and flowering plants, refers to the differential expression of alleles of the same genes in a parent-of-origin-specific manner. In Arabidopsis, imprinting occurs primarily in the endosperm, which nourishes the developing embryo. Recent high-throughput sequencing analyses revealed that more than 200 loci are imprinted in Arabidopsis; however, only a few of these imprinted genes and their imprinting mechanisms have been examined in detail. Whereas most imprinted loci characterized to date are maternally expressed imprinted genes (MEGs), PHERES1 (PHE1) and ADMETOS (ADM) are paternally expressed imprinted genes (PEGs). Here, we report that UPWARD CURLY LEAF1 (UCL1), a gene encoding an E3 ligase that degrades the CURLY LEAF (CLF) polycomb protein, is a PEG. After fertilization, paternally inherited UCL1 is expressed in the endosperm, but not in the embryo. The expression pattern of a β-glucuronidase (GUS) reporter gene driven by the UCL1 promoter suggests that the imprinting control region (ICR) of UCL1 is adjacent to a transposable element in the UCL1 5′-upstream region. Polycomb Repressive Complex 2 (PRC2) silences the maternal UCL1 allele in the central cell prior to fertilization and in the endosperm after fertilization. The UCL1 imprinting pattern was not affected in paternal PRC2 mutants. We found unexpectedly that the maternal UCL1 allele is reactivated in the endosperm of Arabidopsis lines with mutations in cytosine DNA METHYLTRANSFERASE 1 (MET1) or the DNA glycosylase DEMETER (DME), which antagonistically regulate CpG methylation of DNA. By contrast, maternal UCL1 silencing was not altered in mutants with defects in non-CpG methylation. Thus, silencing of the maternal UCL1 allele is regulated by both MET1 and DME as well as by PRC2, suggesting that divergent mechanisms for the regulation of PEGs evolved in Arabidopsis.

Tandem repeat modification during double-strand break repair induced by an engineered TAL effector nuclease in zebrafish genome

Tandem repeats (TRs) are abundant and widely distributed in eukaryotic genomes. TRs are thought to have various functions in gene transcription, DNA methylation, nucleosome position and chromatin organization. Variation of repeat units in the genome is observed in association with a number of diseases, such as Fragile X Syndrome, Huntington's disease and Friedreich's ataxia. However, the underlying mechanisms involved are poorly understood, largely owing to the technical limitations in modification of TRs at definite sites in the genome in vivo. Transcription activator-like effector nucleases (TALENs) are widely used in recent years in gene targeting for their specific binding to target sequences when engineered in vitro. Here, we show that the repair of a double-strand break (DSB) induced by TALENs adjacent to a TR can produce serial types of mutations in the TR region. Sequencing analysis revealed that there are three types of mutations induced by the DSB repair, including indels only within the TR region or within the flanking TALEN target region or simutaneously within both regions. Therefore, desired TR mutant types can be conveniently obtained by using engineered TALENs. These results demonstrate that TALENs can serve as a convenient tool for modifying TRs in the genome in studying the functions of TRs.

Chromosomes carrying meiotic avoidance loci in three apomictic eudicot Hieracium subgenus Pilosella species share structural features with two monocot apomicts

The LOSS OF APOMEIOSIS (LOA) locus is one of two dominant loci known to control apomixis in the eudicot Hieracium praealtum. LOA stimulates the differentiation of somatic aposporous initial cells after the initiation of meiosis in ovules. Aposporous initial cells undergo nuclear proliferation close to sexual megaspores, forming unreduced aposporous embryo sacs, and the sexual program ceases. LOA-linked genetic markers were used to isolate 1.2 Mb of LOA-associated DNAs from H. praealtum. Physical mapping defined the genomic region essential for LOA function between two markers, flanking 400 kb of identified sequence and central unknown sequences. Cytogenetic and sequence analyses revealed that the LOA locus is located on a single chromosome near the tip of the long arm and surrounded by extensive, abundant complex repeat and transposon sequences. Chromosomal features and LOA-linked markers are conserved in aposporous Hieracium caespitosum and Hieracium piloselloides but absent in sexual Hieracium pilosella. Their absence in apomictic Hieracium aurantiacum suggests that meiotic avoidance may have evolved independently in aposporous subgenus Pilosella species. The structure of the hemizygous chromosomal region containing the LOA locus in the three Hieracium subgenus Pilosella species resembles that of the hemizygous apospory-specific genomic regions in monocot Pennisetum squamulatum and Cenchrus ciliaris. Analyses of partial DNA sequences at these loci show no obvious conservation, indicating that they are unlikely to share a common ancestral origin. This suggests convergent evolution of repeat-rich hemizygous chromosomal regions containing apospory loci in these monocot and eudicot species, which may be required for the function and maintenance of the trait.

Assembling a puzzle of dispersed retrotransposable sequences in the genome of chickpea (Cicer arietinum L.)

Several repetitive elements are known to be present in the genome of chickpea (Cicer arietinum L.) including satellite DNA and En/Spm transposons as well as two dispersed, highly repetitive elements, CaRep1 and CaRep2. PCR was used to prove that CaRep1, CaRep2, and previously isolated CaRep3 of C. arietinum represent different segments of a highly repetitive Ty3-gypsy-like retrotransposon (Metaviridae) designated CaRep that makes up large parts of the intercalary heterochromatin. The full sequence of this element including the LTRs and untranslated internal regions was isolated by selective amplification. The restriction pattern of CaRep was different within the annual species of the genus Cicer, suggesting its rearrangement during the evolution of the genus during the last 100 000 years. In addition to CaRep, another LTR and a non-LTR retrotransposon family were isolated, and their restriction patterns and physical localization in the chickpea genome were characterized. The LINE-like element CaLin is only of comparatively low abundance and reveals a considerable heterogeneity. The Ty1-copia-like element (Pseudoviridae) CaTy is located in the distal parts of the intercalary heterochromatin and adjacent euchromatic regions, but it is absent from the centromeric regions. These results together with earlier findings allow to depict the distribution of retroelements on chickpea chromosomes, which extensively resembles the retroelement landscape of the genome of the model legume Medicago truncatula Gaertn.

DNA methylation profiling in zebrafish

DNA methylation on cytosine in vertebrates such as zebrafish serves to silence gene expression by interfering with the binding of certain transcription factors and through the recruitment of repressive chromatin machinery. Cytosine DNA methylation is chemically stable and heritable through the germline - but also reversible through many modes, making it a useful and dynamic epigenetic modification. Virtually all of the enzymes and factors involved in the deposition, binding, and removal of cytosine methylation are conserved in zebrafish, and therefore the organism an excellent model for understanding the use of DNA methylation in the control of gene regulation and other processes. Here, we discuss the main approaches to quantifying DNA methylation levels genome-wide in zebrafish: one is an established method for revealing regional methylation (methylated DNA immunoprecipitation (MeDIP)), and the other is an emerging method that reveals DNA methylation at base-pair resolution (shotgun bisulphite sequencing). We also introduce some of the analytical methods that are useful for identifying regions of hypo- or hyper-methylation, and ways to identify differentially methylated regions.

Chromodomains direct integration of retrotransposons to heterochromatin

The enrichment of mobile genetic elements in heterochromatin may be due, in part, to targeted integration. The chromoviruses are Ty3/gypsy retrotransposons with chromodomains at their integrase C termini. Chromodomains are logical determinants for targeting to heterochromatin, because the chromodomain of heterochromatin protein 1 (HP1) typically recognizes histone H3 K9 methylation, an epigenetic mark characteristic of heterochromatin. We describe three groups of chromoviruses based on amino acid sequence relationships of their integrase C termini. Genome sequence analysis indicates that representative chromoviruses from each group are enriched in gene-poor regions of the genome relative to other retrotransposons, and when fused to fluorescent marker proteins, the chromodomains target proteins to specific subnuclear foci coincident with heterochromatin. The chromodomain of the fungal element, MAGGY, interacts with histone H3 dimethyl- and trimethyl-K9, and when the MAGGY chromodomain is fused to integrase of the Schizosaccharomyces pombe Tf1 retrotransposon, new Tf1 insertions are directed to sites of H3 K9 methylation. Repetitive sequences such as transposable elements trigger the RNAi pathway resulting in their epigenetic modification. Our results suggest a dynamic interplay between retrotransposons and heterochromatin, wherein mobile elements recognize heterochromatin at the time of integration and then perpetuate the heterochromatic mark by triggering epigenetic modification.

A Ds-insertion mutant of OSH6 (Oryza sativa Homeobox 6) exhibits outgrowth of vestigial leaf-like structures, bracts, in rice

OSH6 (Oryza sativa Homeobox6) is an ortholog of lg3 (Liguleless3) in maize. We generated a novel allele, termed OSH6-Ds, by inserting a defective Ds element into the third exon of OSH6, which resulted in a truncated OSH6 mRNA. The truncated mRNA was expressed ectopically in leaf tissues and encoded the N-terminal region of OSH6, which includes the KNOX1 and partial KNOX2 subdomains. This recessive mutant showed outgrowth of bracts or produced leaves at the basal node of the panicle. These phenotypes distinguished it from the OSH6 transgene whose ectopic expression led to a "blade to sheath transformation" phenotype at the midrib region of leaves, similar to that seen in dominant Lg3 mutants. Expression of a similar truncated OSH6 cDNA from the 35S promoter (35S::DeltaOSH6) confirmed that the ectopic expression of this product was responsible for the aberrant bract development. These data suggest that OSH6-Ds interferes with a developmental mechanism involved in bract differentiation, especially at the basal nodes of panicles.

A genome-wide transcriptional activity survey of rice transposable element-related genes

A genome-wide survey of the transcriptional activity of TE-related genes that were associated with fifteen developmental stages and stress conditions revealed clear, albeit low, general transcription of TE-related genes.

Background

Transposable element (TE)-related genes comprise a significant portion of the gene catalog of grasses, although their functions are insufficiently characterized. The recent availability of TE-related gene annotation from the complete genome sequence of rice (Oryza sativa) has created an opportunity to conduct a comprehensive evaluation of the transcriptional activities of these potentially mobile elements and their related genes.

Results

We conducted a genome-wide survey of the transcriptional activity of TE-related genes associated with 15 developmental stages and stress conditions. This dataset was obtained using a microarray encompassing 2,191 unique TE-related rice genes, which were represented by oligonucleotide probes that were free from cross-hybridization. We found that TE-related genes exhibit much lower transcriptional activities than do non-TE-related genes, although representative transcripts were detected from all superfamilies of both type I and II TE-related genes. The strongest transcriptional activities were detected in TE-related genes from among the MULE and CACTA superfamilies. Phylogenetic analyses suggest that domesticated TE-related genes tend to form clades with active transcription. In addition, chromatin-level regulations through histone and DNA modifications, as well as enrichment of certain cis elements in the promoters, appear to contribute to the transcriptional activation of representative TE-related genes.

Conclusion

Our findings reveal clear, albeit low, general transcription of TE-related genes. In combination with phylogenetic analysis, transcriptional analysis has the potential to lead to the identification of domesticated TEs with adapted host functions.

ROSINA (RSI) is part of a CACTA transposable element, TamRSI, and links flower development to transposon activity

ROSINA (RSI) was isolated as a DNA binding factor able to bind to the CArG-box present in the promoter of the MADS-box gene DEFICIENS of Antirrhinum majus. The mosaic nature of RSI and its multi-copy presence in the A. majus genome indicated that RSI could be a part of a mobile genetic element. Here we show that RSI is a part of a CACTA transposable element system of A. majus, named TamRSI, which has evolved and is still evolving within the terminal inverted repeats (TIRs) of this CACTA transposon. Interestingly, RSI is always found in opposite orientation with respect to the transcription of a second gene present within the CACTA transposon, which encodes a putative TRANSPOSASE (TNP). This structural configuration has not yet been described for any member of the CACTA transposons superfamily. Internal deletion derivatives of the TamRSI produce aberrant RSI transcripts (RSI-ATs) that carry parts of the RSI RNA fused to parts of the TNP RNA. In addition, an intriguing seed phenotype shown by RNAi transgenic lines generated to silence RSI, relate TamRSI to epigenetic mechanisms and associate the control of flower development to transposon activity.

Studies on the mechanism of RNAi-dependent heterochromatin assembly

Assembly of heterochromatin at centromeric DNA regions in the fission yeast Schizosaccharomyces pombe involves an intimate interplay between chromatin modifying complexes and components of the RNAi pathway. The RNA-induced transcriptional silencing (RITS) complex, containing Chp1, Ago1, Tas3, and centromeric siRNAs, localizes to centromeric DNA repeats and is required for the assembly and maintenance of heterochromatin. RITS brings together two types of molecular recognition modules: a chromodomain protein, which binds to lysine 9 methylated histone H3 (H3K9), and Argonaute, which binds to specific sequences by siRNA-directed base-pairing interactions. The RNA-directed RNA polymerase complex (RDRC), composed of Rdp1, the Hrr1 helicase, and the Cid12 Poly(A) polymerase family member, synthesizes double-stranded RNA and creates the substrate for Dicer to generate siRNAs. RDRC physically associates with RITS, and both complexes localize to noncoding centromeric RNAs and centromeric DNA repeats, suggesting that recognition of nascent RNA transcripts may be involved in localization of these complexes to specific chromosome regions. In support of this possibility, tethering of the RITS complex to the transcript of the normally euchromatic ura4 (+) gene results in siRNA generation and RNAi- and heterochromatin-dependent silencing of the ura4 (+) gene. Finally, silencing of a subset of endogenous and transgene promoters within heterochromatic DNA domains occurs by RNAi-dependent degradation of nascent transcripts by a mechanism that we have termed co-transcriptional gene silencing (CTGS).

MicroRNAs in the search for understanding human diseases

MiroRNAs (miRNAs) are double-stranded, noncoding RNA molecules (with an average size of 22bp) that serve as post-transcriptional regulators of gene expression in higher eukaryotes. miRNAs play an important role in development and other cellular processes by hybridizing with complementary target mRNA transcripts, preventing their translation and thereby destabilizing the target transcripts. Though hundreds of miRNAs have been discovered in a variety of organisms, little is known about their cellular function. They have been implicated in the regulation of developmental timing and pattern formation, restriction of differentiation potential, regulation of insulin secretion, resistance to viral infection, and in genomic rearrangements associated with carcinogenesis or other genetic disorders, such as fragile X syndrome. Recent evidence suggests that the number of unique miRNA genes in humans exceeds 1000, and may be as high as 20,000. It is estimated that 20-30% of all human mRNAs are miRNA targets. During the last few years, special attention has been given to miRNAs as candidate drug targets for cancer, diabetes mellitus, obesity, and viral diseases.

Noncoding RNAs and gene silencing

Noncoding RNA has long been proposed to control gene expression via sequence-specific interactions with regulatory regions. Here, we review the role of noncoding RNA in heterochromatic silencing and in the silencing of transposable elements (TEs), unpaired DNA in meiosis, and developmentally excised DNA. The role of cotranscriptional processing by RNA interference and by other mechanisms is discussed, as well as parallels with RNA silencing in imprinting, paramutation, polycomb silencing, and X inactivation. Interactions with regulatory sequences may well occur, but at the RNA rather than at the DNA level.

Zebra fish Dnmt1 and Suv39h1 regulate organ-specific terminal differentiation during development

DNA methylation and histone methylation are two key epigenetic modifications that help govern heterochromatin dynamics. The roles for these chromatin-modifying activities in directing tissue-specific development remain largely unknown. To address this issue, we examined the roles of DNA methyltransferase 1 (Dnmt1) and the H3K9 histone methyltransferase Suv39h1 in zebra fish development. Knockdown of Dnmt1 in zebra fish embryos caused defects in terminal differentiation of the intestine, exocrine pancreas, and retina. Interestingly, not all tissues required Dnmt1, as differentiation of the liver and endocrine pancreas appeared normal. Proper differentiation depended on Dnmt1 catalytic activity, as Dnmt1 morphants could be rescued by active zebra fish or human DNMT1 but not by catalytically inactive derivatives. Dnmt1 morphants exhibited dramatic reductions of both genomic cytosine methylation and genome-wide H3K9 trimethyl levels, leading us to investigate the overlap of in vivo functions of Dnmt1 and Suv39h1. Embryos lacking Suv39h1 had organ-specific terminal differentiation defects that produced largely phenocopies of Dnmt1 morphants but retained wild-type levels of DNA methylation. Remarkably, suv39h1 overexpression rescued markers of terminal differentiation in Dnmt1 morphants. Our results suggest that Dnmt1 activity helps direct histone methylation by Suv39h1 and that, together, Dnmt1 and Suv39h1 help guide the terminal differentiation of particular tissues.

Diverse chromatin remodeling genes antagonize the Rb-involved SynMuv pathways in C. elegans

In Caenorhabditis elegans, vulval cell-fate specification involves the activities of multiple signal transduction and regulatory pathways that include a receptor tyrosine kinase/Ras/mitogen-activated protein kinase pathway and synthetic multivulva (SynMuv) pathways. Many genes in the SynMuv pathways encode transcription factors including the homologs of mammalian Rb, E2F, and components of the nucleosome-remodeling deacetylase complex. To further elucidate the functions of the SynMuv genes, we performed a genome-wide RNA interference (RNAi) screen to search for genes that antagonize the SynMuv gene activities. Among those that displayed a varying degree of suppression of the SynMuv phenotype, 32 genes are potentially involved in chromatin remodeling (called SynMuv suppressor genes herein). Genetic mutations of two representative genes (zfp-1 and mes-4) were used to further characterize their positive roles in vulval induction and relationships with Ras function. Our analysis revealed antagonistic roles of the SynMuv suppressor genes and the SynMuv B genes in germline-soma distinction, RNAi, somatic transgene silencing, and tissue specific expression of pgl-1 and the lag-2/Delta genes. The opposite roles of these SynMuv B and SynMuv suppressor genes on transcriptional regulation were confirmed in somatic transgene silencing. We also report the identifications of ten new genes in the RNAi pathway and six new genes in germline silencing. Among the ten new RNAi genes, three encode homologs of proteins involved in both protein degradation and chromatin remodeling. Our findings suggest that multiple chromatin remodeling complexes are involved in regulating the expression of specific genes that play critical roles in developmental decisions.

Chromatin structure and gene regulation in T cell development and function

Transcription factors control gene expression programs in the context of the chromatin structure of their target genes. DNA methylation, post-translational histone modifications such as acetylation and methylation, and higher order chromatin organization allow the maintenance of gene expression patterns through mitosis, but how do they accommodate developmentally regulated changes in gene expression programs? Although histone acetylation and deacetylation are in dynamic equilibrium and mechanisms for the removal of methyl groups from histones are emerging, the extent to which there is active demethylation of DNA remains controversial. Looking at chromatin in the three-dimensional space of the nucleus, recent work demonstrates that gene regulation involves contacts between regulatory elements within genes or gene clusters on the same chromosome (in cis) and between different chromosomes (in trans). Finally, non-coding RNAs make a significant contribution to transcriptional and post-transcriptional gene silencing. Together, these advances contribute to an understanding of how gene expression programs are established, maintained and modified during development.