Epigenetic regulator function through mouse gastrulation.

During ontogeny, proliferating cells become restricted in their fate through the combined action of cell-type-specific transcription factors and ubiquitous epigenetic machinery, which recognizes universally available histone residues or nucleotides in a context-dependent manner1,2. The molecular functions of these regulators are generally well understood, but assigning direct developmental roles to them is hampered by complex mutant phenotypes that often emerge after gastrulation3,4. Single-cell RNA sequencing and analytical approaches have explored this highly conserved, dynamic period across numerous model organisms5-8, including mouse9-18. Here we advance these strategies using a combined zygotic perturbation and single-cell RNA-sequencing platform in which many mutant mouse embryos can be assayed simultaneously, recovering robust  morphological and transcriptional information across a panel of ten essential regulators. Deeper analysis of central Polycomb repressive complex (PRC) 1 and 2 components indicates substantial cooperativity, but distinguishes a dominant role for PRC2 in restricting the germline. Moreover, PRC mutant phenotypes emerge after gross epigenetic and transcriptional changes within the initial conceptus prior to gastrulation. Our experimental framework may eventually lead to a fully quantitative view of how cellular diversity emerges using an identical genetic template and from a single totipotent cell.

Pharmacological Improvement of Cystic Fibrosis Transmembrane Conductance Regulator Function Rescues Airway Epithelial Homeostasis and Host Defense in Children with Cystic Fibrosis.

Rationale: Pharmacological improvement of cystic fibrosis transmembrane conductance regulator (CFTR) function with elexacaftor/tezacaftor/ivacaftor (ETI) provides unprecedented improvements in lung function and other clinical outcomes in patients with cystic fibrosis (CF). However, ETI effects on impaired mucosal homeostasis and host defense at the molecular and cellular levels in the airways of patients with CF remain unknown. Objectives: To investigate effects of ETI on the transcriptome of nasal epithelial and immune cells from children with CF at the single-cell level. Methods: Nasal swabs from 13 children with CF and at least one F508del allele aged 6 to 11 years were collected at baseline and 3 months after initiation of ETI, subjected to single-cell RNA sequencing, and compared with swabs from 12 age-matched healthy children. Measurements and Main Results: Proportions of CFTR-positive cells were decreased in epithelial basal, club, and goblet cells, but not in ionocytes, from children with CF at baseline and were restored by ETI therapy to nearly healthy levels. Single-cell transcriptomics revealed an impaired IFN signaling and reduced expression of major histocompatibility complex classes I and II encoding genes in epithelial cells of children with CF at baseline, which was partially restored by ETI. In addition, ETI therapy markedly reduced the inflammatory phenotype of immune cells, particularly of neutrophils and macrophages. Conclusions: Pharmacological improvement of CFTR function improves innate mucosal immunity and reduces immune cell inflammatory responses in the upper airways of children with CF at the single-cell level, highlighting the potential to restore epithelial homeostasis and host defense in CF airways by early initiation of ETI therapy.

Non-classical circulating monocytes expressing high levels of microsomal prostaglandin E2 synthase-1 tag an aberrant IFN-response in systemic sclerosis.

Systemic sclerosis (SSc) is a complex disease that affects the connective tissue, causing fibrosis. SSc patients show altered immune cell composition and activation in the peripheral blood (PB). PB monocytes (Mos) are recruited into tissues where they differentiate into macrophages, which are directly involved in fibrosis. To understand the role of CD14+ PB Mos in SSc, a single-cell transcriptome analysis (scRNA-seq) was conducted on 8 SSc patients and 8 controls. Using unsupervised clustering methods, CD14+ cells were assigned to 11 clusters, which added granularity to the known monocyte subsets: classical (cMos), intermediate (iMos) and non-classical Mos (ncMos) or type 2 dendritic cells. NcMos were significantly overrepresented in SSc patients and showed an active IFN-signature and increased expression levels of PTGES, in addition to monocyte motility and adhesion markers. We identified a SSc-related cluster of IRF7+ STAT1+ iMos with an aberrant IFN-response. Finally, a depletion of M2 polarised cMos in SSc was observed. Our results highlighted the potential of PB Mos as biomarkers for SSc and provided new possibilities for putative drug targets for modulating the innate immune response in SSc.

An Unbiased Molecular Approach Using 3'-UTRs Resolves the Avian Family-Level Tree of Life.

Presumably, due to a rapid early diversification, major parts of the higher-level phylogeny of birds are still resolved controversially in different analyses or are considered unresolvable. To address this problem, we produced an avian tree of life, which includes molecular sequences of one or several species of ∼90% of the currently recognized family-level taxa (429 species, 379 genera) including all 106 family-level taxa of the nonpasserines and 115 of the passerines (Passeriformes). The unconstrained analyses of noncoding 3-prime untranslated region (3'-UTR) sequences and those of coding sequences yielded different trees. In contrast to the coding sequences, the 3'-UTR sequences resulted in a well-resolved and stable tree topology. The 3'-UTR contained, unexpectedly, transcription factor binding motifs that were specific for different higher-level taxa. In this tree, grebes and flamingos are the sister clade of all other Neoaves, which are subdivided into five major clades. All nonpasserine taxa were placed with robust statistical support including the long-time enigmatic hoatzin (Opisthocomiformes), which was found being the sister taxon of the Caprimulgiformes. The comparatively late radiation of family-level clades of the songbirds (oscine Passeriformes) contrasts with the attenuated diversification of nonpasseriform taxa since the early Miocene. This correlates with the evolution of vocal production learning, an important speciation factor, which is ancestral for songbirds and evolved convergent only in hummingbirds and parrots. As 3'-UTR-based phylotranscriptomics resolved the avian family-level tree of life, we suggest that this procedure will also resolve the all-species avian tree of life.

The salmon louse genome: Copepod features and parasitic adaptations.

Copepods encompass numerous ecological roles including parasites, detrivores and phytoplankton grazers. Nonetheless, copepod genome assemblies remain scarce. Lepeophtheirus salmonis is an economically and ecologically important ectoparasitic copepod found on salmonid fish. We present the 695.4 Mbp L. salmonis genome assembly containing ≈60% repetitive regions and 13,081 annotated protein-coding genes. The genome comprises 14 autosomes and a ZZ-ZW sex chromosome system. Assembly assessment identified 92.4% of the expected arthropod genes. Transcriptomics supported annotation and indicated a marked shift in gene expression after host attachment, including apparent downregulation of genes related to circadian rhythm coinciding with abandoning diurnal migration. The genome shows evolutionary signatures including loss of genes needed for peroxisome biogenesis, presence of numerous FNII domains, and an incomplete heme homeostasis pathway suggesting heme proteins to be obtained from the host. Despite repeated development of resistance against chemical treatments L. salmonis exhibits low numbers of many genes involved in detoxification.

Genomic insights into the physiology and ecology of the marine filamentous cyanobacterium Lyngbya majuscula.

Filamentous cyanobacteria of the genus Lyngbya are important contributors to coral reef ecosystems, occasionally forming dominant cover and impacting the health of many other co-occurring organisms. Moreover, they are extraordinarily rich sources of bioactive secondary metabolites, with 35% of all reported cyanobacterial natural products deriving from this single pantropical genus. However, the true natural product potential and life strategies of Lyngbya strains are poorly understood because of phylogenetic ambiguity, lack of genomic information, and their close associations with heterotrophic bacteria and other cyanobacteria. To gauge the natural product potential of Lyngbya and gain insights into potential microbial interactions, we sequenced the genome of Lyngbya majuscula 3L, a Caribbean strain that produces the tubulin polymerization inhibitor curacin A and the molluscicide barbamide, using a combination of Sanger and 454 sequencing approaches. Whereas ∼ 293,000 nucleotides of the draft genome are putatively dedicated to secondary metabolism, this is far too few to encode a large suite of Lyngbya metabolites, suggesting Lyngbya metabolites are strain specific and may be useful in species delineation. Our analysis revealed a complex gene regulatory network, including a large number of sigma factors and other regulatory proteins, indicating an enhanced ability for environmental adaptation or microbial associations. Although Lyngbya species are reported to fix nitrogen, nitrogenase genes were not found in the genome or by PCR of genomic DNA. Subsequent growth experiments confirmed that L. majuscula 3L is unable to fix atmospheric nitrogen. These unanticipated life history characteristics challenge current views of the genus Lyngbya.

Nephrocystin-5, a ciliary IQ domain protein, is mutated in Senior-Loken syndrome and interacts with RPGR and calmodulin.

Nephronophthisis (NPHP) is the most frequent genetic cause of chronic renal failure in children. Identification of four genes mutated in NPHP subtypes 1-4 (refs. 4-9) has linked the pathogenesis of NPHP to ciliary functions. Ten percent of affected individuals have retinitis pigmentosa, constituting the renal-retinal Senior-Loken syndrome (SLSN). Here we identify, by positional cloning, mutations in an evolutionarily conserved gene, IQCB1 (also called NPHP5), as the most frequent cause of SLSN. IQCB1 encodes an IQ-domain protein, nephrocystin-5. All individuals with IQCB1 mutations have retinitis pigmentosa. Hence, we examined the interaction of nephrocystin-5 with RPGR (retinitis pigmentosa GTPase regulator), which is expressed in photoreceptor cilia and associated with 10-20% of retinitis pigmentosa. We show that nephrocystin-5, RPGR and calmodulin can be coimmunoprecipitated from retinal extracts, and that these proteins localize to connecting cilia of photoreceptors and to primary cilia of renal epithelial cells. Our studies emphasize the central role of ciliary dysfunction in the pathogenesis of SLSN.

Transcriptome sequencing and microarray design for functional genomics in the extremophile Arabidopsis relative Thellungiella salsuginea (Eutrema salsugineum).

BACKGROUND: Most molecular studies of plant stress tolerance have been performed with Arabidopsis thaliana, although it is not particularly stress tolerant and may lack protective mechanisms required to survive extreme environmental conditions. Thellungiella salsuginea has attracted interest as an alternative plant model species with high tolerance of various abiotic stresses. While the T. salsuginea genome has recently been sequenced, its annotation is still incomplete and transcriptomic information is scarce. In addition, functional genomics investigations in this species are severely hampered by a lack of affordable tools for genome-wide gene expression studies. RESULTS: Here, we report the results of Thellungiella de novo transcriptome assembly and annotation based on 454 pyrosequencing and development and validation of a T. salsuginea microarray. ESTs were generated from a non-normalized and a normalized library synthesized from RNA pooled from samples covering different tissues and abiotic stress conditions. Both libraries yielded partially unique sequences, indicating their necessity to obtain comprehensive transcriptome coverage. More than 1 million sequence reads were assembled into 42,810 unigenes, approximately 50% of which could be functionally annotated. These unigenes were compared to all available Thellungiella genome sequence information. In addition, the groups of Late Embryogenesis Abundant (LEA) proteins, Mitogen Activated Protein (MAP) kinases and protein phosphatases were annotated in detail. We also predicted the target genes for 384 putative miRNAs. From the sequence information, we constructed a 44 k Agilent oligonucleotide microarray. Comparison of same-species and cross-species hybridization results showed superior performance of the newly designed array for T. salsuginea samples. The developed microarrays were used to investigate transcriptional responses of T. salsuginea and Arabidopsis during cold acclimation using the MapMan software. CONCLUSIONS: This study provides the first comprehensive transcriptome information for the extremophile Arabidopsis relative T. salsuginea. The data constitute a more than three-fold increase in the number of publicly available unigene sequences and will greatly facilitate genome annotation. In addition, we have designed and validated the first genome-wide microarray for T. salsuginea, which will be commercially available. Together with the publicly available MapMan software this will become an important tool for functional genomics of plant stress tolerance.

Comparative genomic analysis of fruiting body formation in Myxococcales.

Genetic programs underlying multicellular morphogenesis and cellular differentiation are most often associated with eukaryotic organisms, but examples also exist in bacteria such as the formation of multicellular, spore-filled fruiting bodies in the order Myxococcales. Most members of the Myxococcales undergo a multicellular developmental program culminating in the formation of spore-filled fruiting bodies in response to starvation. To gain insight into the evolutionary history of fruiting body formation in Myxococcales, we performed a comparative analysis of the genomes and transcriptomes of five Myxococcales species, four of these undergo fruiting body formation (Myxococcus xanthus, Stigmatella aurantiaca, Sorangium cellulosum, and Haliangium ochraceum) and one does not (Anaeromyxobacter dehalogenans). Our analyses show that a set of 95 known M. xanthus development-specific genes--although suffering from a sampling bias--are overrepresented and occur more frequently than an average M. xanthus gene in S. aurantiaca, whereas they occur at the same frequency as an average M. xanthus gene in S. cellulosum and in H. ochraceum and are underrepresented in A. dehalogenans. Moreover, genes for entire signal transduction pathways important for fruiting body formation in M. xanthus are conserved in S. aurantiaca, whereas only a minority of these genes are conserved in A. dehalogenans, S. cellulosum, and H. ochraceum. Likewise, global gene expression profiling of developmentally regulated genes showed that genes that upregulated during development in M. xanthus are overrepresented in S. aurantiaca and slightly underrepresented in A. dehalogenans, S. cellulosum, and H. ochraceum. These comparative analyses strongly indicate that the genetic programs for fruiting body formation in M. xanthus and S. aurantiaca are highly similar and significantly different from the genetic program directing fruiting body formation in S. cellulosum and H. ochraceum. Thus, our analyses reveal an unexpected level of plasticity in the genetic programs for fruiting body formation in the Myxococcales and strongly suggest that the genetic program underlying fruiting body formation in different Myxococcales is not conserved. The evolutionary implications of this finding are discussed.

Identification of genes essential for prey-independent growth of Bdellovibrio bacteriovorus HD100.

Bdellovibrio bacteriovorus HD100 is an obligate predatory bacterium that attacks and invades Gram-negative bacteria. The predator requires living bacteria to survive as growth and replication take place inside the bacterial prey. It is possible to isolate mutants that grow and replicate outside prey bacteria. Such mutants are designated host or prey independent, and their nutritional requirements vary. Some mutants are saprophytic and require prey extracts for extracellular growth, whereas other mutants grow axenically, which denotes the formation of colonies on complete medium in the absence of any prey components. The initial events leading to prey-independent growth are still under debate, and several genes may be involved. We selected new mutants by three different methods: spontaneous mutation, transposon mutagenesis, and targeted gene knockout. By all approaches we isolated mutants of the hit (host interaction) locus. As the relevance of this locus for the development of prey independence has been questioned, we performed whole-genome sequencing of five prey-independent mutants. Three mutants were saprophytic, and two mutants could grow axenically. Whole-genome analysis revealed that the mutation of a small open reading frame of the hit locus is sufficient for the conversion from predatory to saprophytic growth. Complementation experiments were performed by introduction of a plasmid carrying the wild-type hit gene into saprophytic mutants, and predatory growth could be restored. Whole-genome sequencing of two axenic mutants demonstrated that in addition to the hit mutation the colony formation on complete medium was shown to be influenced by the mutations of two genes involved in RNA processing. Complementation experiments with a wild-type gene encoding an RNA helicase, RhlB, abolished the ability to form colonies on complete medium, indicating that stability of RNA influences axenic growth.

Genome of Streptococcus oralis strain Uo5.

Streptococcus oralis, a commensal species of the human oral cavity, belongs to the Mitis group of streptococci, which includes one of the major human pathogens as well, S. pneumoniae. We report here the first complete genome sequence of this species. S. oralis Uo5, a high-level penicillin- and multiple-antibiotic-resistant isolate from Hungary, is competent for genetic transformation under laboratory conditions. Comparative and functional genomics of Uo5 will be important in understanding the evolution of pathogenesis among Mitis streptococci and their potential to engage in interspecies gene transfer.

Deep sequencing of naupliar-, cyprid- and adult-specific normalised Expressed Sequence Tag (EST) libraries of the acorn barnacle Balanus amphitrite.

In order to improve the genetic characterisation of the barnacle Balanus amphitrite, normalised EST libraries for the developmental stages, viz. nauplius (a mix of instars I and II), cyprid and adult, were generated. The libraries were sequenced independently using 454 technologies and 575,666 reads were generated. For adults, 4843 unique isotigs were estimated and 6754 and 7506 in the cyprid and naupliar stage, respectively. It was found that some of the previously proposed cyprid-specific bcs genes were also expressed during the naupliar and adult stage. Furthermore, as lectins have been hypothesised to influence settlement cue recognition in barnacles, the database was searched for lectin-like isotigs. Two proteins, uniquely expressed in either the cyprid or the adult stage, matched a mannose receptor, and their nucleotide sequences were 33% and 31% identical to a lectin (BRA-3) isolated from Megabalanus rosa. Further characterisation of these genes may suggest their involvement in settlement.

Preneoplastic Alterations Define CLL DNA Methylome and Persist through Disease Progression and Therapy.

Most human cancers converge to a deregulated methylome with reduced global levels and elevated methylation at select CpG islands. To investigate the emergence and dynamics of the cancer methylome, we characterized genome-wide DNA methylation in pre-neoplastic monoclonal B cell lymphocytosis (MBL) and chronic lymphocytic leukemia (CLL), including serial samples collected across disease course. We detected the aberrant tumor-associated methylation landscape at CLL diagnosis and found no significantly differentially methylated regions in the high-count MBL-to-CLL transition. Patient methylomes showed remarkable stability with natural disease and post-therapy progression. Single CLL cells were consistently aberrantly methylated, indicating a homogeneous transition to the altered epigenetic state, and a distinct expression profile together with MBL cells compared to normal B cells. Our longitudinal analysis reveals the cancer methylome to emerge early, which may provide a platform for subsequent genetically-driven growth dynamics and together with its persistent presence suggests a central role in the normal-to-cancer transition.

Analysis of a normalised expressed sequence tag (EST) library from a key pollinator, the bumblebee Bombus terrestris.

BACKGROUND: The bumblebee, Bombus terrestris (Order Hymenoptera), is of widespread importance. This species is extensively used for commercial pollination in Europe, and along with other Bombus spp. is a key member of natural pollinator assemblages. Furthermore, the species is studied in a wide variety of biological fields. The objective of this project was to create a B. terrestris EST resource that will prove to be valuable in obtaining a deeper understanding of this significant social insect. RESULTS: A normalised cDNA library was constructed from the thorax and abdomen of B. terrestris workers in order to enhance the discovery of rare genes. A total of 29'428 ESTs were sequenced. Subsequent clustering resulted in 13'333 unique sequences. Of these, 58.8 percent had significant similarities to known proteins, with 54.5 percent having a "best-hit" to existing Hymenoptera sequences. Comparisons with the honeybee and other insects allowed the identification of potential candidates for gene loss, pseudogene evolution, and possible incomplete annotation in the honeybee genome. Further, given the focus of much basic research and the perceived threat of disease to natural and commercial populations, the immune system of bumblebees is a particularly relevant component. Although the library is derived from unchallenged bees, we still uncover transcription of a number of immune genes spanning the principally described insect immune pathways. Additionally, the EST library provides a resource for the discovery of genetic markers that can be used in population level studies. Indeed, initial screens identified 589 simple sequence repeats and 854 potential single nucleotide polymorphisms. CONCLUSION: The resource that these B. terrestris ESTs represent is valuable for ongoing work. The ESTs provide direct evidence of transcriptionally active regions, but they will also facilitate further functional genomics, gene discovery and future genome annotation. These are important aspects in obtaining a greater understanding of this key pollinator species.

A 454 sequencing approach for large scale phylogenomic analysis of the common emperor scorpion (Pandinus imperator).

In recent years, phylogenetic tree reconstructions that rely on multiple gene alignments that had been deduced from expressed sequence tags (ESTs) have become a popular method in molecular systematics. Here, we present a 454 pyrosequencing approach to infer the transcriptome of the Emperor scorpion Pandinus imperator. We obtained 428,844 high-quality reads (mean length=223+/-50 b) from total cDNA, which were assembled into 8334 contigs (mean length 422+/-313 bp) and 26,147 singletons. About 1200 contigs were successfully annotated by BLAST and orthology search. Specific analyses of eight distinct hemocyanin sequences provided further proof for the quality of the 454 reads and the assembly process. The P. imperator sequences were included in a concatenated alignment of 149 orthologous genes of 67 metazoan taxa that covers 39,842 amino acids. After removal of low-quality regions, 11,168 positions were employed for phylogenetic reconstructions. Using Bayesian and maximum likelihood methods, we obtained strongly supported monophyletic Ecdysozoa, Arthropoda (excluding Tardigrada), Euarthropoda, Pancrustacea and Hexapoda. We also recovered the Myriochelata (Chelicerata+Myriapoda). Within the chelicerates, Pycnogonida form the sister group of Euchelicerata. However, Arachnida were found paraphyletic because the Acari (mites and ticks) were recovered as sister group of a clade comprising Xiphosura, Scorpiones and Araneae. In summary, we have shown that 454 pyrosequencing is a cost-effective method that provides sufficient data and coverage depth for gene detection and multigene-based phylogenetic analyses.

Ice-Age Climate Adaptations Trap the Alpine Marmot in a State of Low Genetic Diversity.

Some species responded successfully to prehistoric changes in climate [1, 2], while others failed to adapt and became extinct [3]. The factors that determine successful climate adaptation remain poorly understood. We constructed a reference genome and studied physiological adaptations in the Alpine marmot (Marmota marmota), a large ground-dwelling squirrel exquisitely adapted to the "ice-age" climate of the Pleistocene steppe [4, 5]. Since the disappearance of this habitat, the rodent persists in large numbers in the high-altitude Alpine meadow [6, 7]. Genome and metabolome showed evidence of adaptation consistent with cold climate, affecting white adipose tissue. Conversely, however, we found that the Alpine marmot has levels of genetic variation that are among the lowest for mammals, such that deleterious mutations are less effectively purged. Our data rule out typical explanations for low diversity, such as high levels of consanguineous mating, or a very recent bottleneck. Instead, ancient demographic reconstruction revealed that genetic diversity was lost during the climate shifts of the Pleistocene and has not recovered, despite the current high population size. We attribute this slow recovery to the marmot's adaptive life history. The case of the Alpine marmot reveals a complicated relationship between climatic changes, genetic diversity, and conservation status. It shows that species of extremely low genetic diversity can be very successful and persist over thousands of years, but also that climate-adapted life history can trap a species in a persistent state of low genetic diversity.

Increasing genomic information in bivalves through new EST collections in four species: development of new genetic markers for environmental studies and genome evolution.

The generation of EST information is an essential step in the genomic characterisation of species. In the context of the European Network Marine Genomics, a common goal was to significantly increase the amount of ESTs in commercial marine mollusk species and more specifically in the less studied but ecologically and commercially important groups, such as mussel and clam genera. Normalized cDNA libraries were constructed for four different relevant bivalves species (Crassostrea gigas, Mytilus edulis, Ruditapes decussatus and Bathymodiolus azoricus), using numerous tissues and physiological conditions. In this paper, we present the analysis of the 13,013 expressed sequence tags (ESTs) generated. Each EST library was independently assembled and 1300-3000 unique sequences were identified in each species. For the different species, functional categories could be assigned to only about 16 to 27% of ESTs using the GO annotation tool. All sequences have been incorporated into a publicly available database and form the basis for subsequent microarray design, SNP detection and polymorphism analysis, and the placement of novel markers on genetic linkage maps.

Comparative analysis of expressed sequence tag (EST) libraries in the seagrass Zostera marina subjected to temperature stress.

Global warming is associated with increasing stress and mortality on temperate seagrass beds, in particular during periods of high sea surface temperatures during summer months, adding to existing anthropogenic impacts, such as eutrophication and habitat destruction. We compare several expressed sequence tag (EST) in the ecologically important seagrass Zostera marina (eelgrass) to elucidate the molecular genetic basis of adaptation to environmental extremes. We compared the tentative unigene (TUG) frequencies of libraries derived from leaf and meristematic tissue from a control situation with two experimentally imposed temperature stress conditions and found that TUG composition is markedly different among these conditions (all P < 0.0001). Under heat stress, we find that 63 TUGs are differentially expressed (d.e.) at 25 degrees C compared with lower, no-stress condition temperatures (4 degrees C and 17 degrees C). Approximately one-third of d.e. eelgrass genes were characteristic for the stress response of the terrestrial plant model Arabidopsis thaliana. The changes in gene expression suggest complex photosynthetic adjustments among light-harvesting complexes, reaction center subunits of photosystem I and II, and components of the dark reaction. Heat shock encoding proteins and reactive oxygen scavengers also were identified, but their overall frequency was too low to perform statistical tests. In all conditions, the most abundant transcript (3-15%) was a putative metallothionein gene with unknown function. We also find evidence that heat stress may translate to enhanced infection by protists. A total of 210 TUGs contain one or more microsatellites as potential candidates for gene-linked genetic markers. Data are publicly available in a user-friendly database at http://www.uni-muenster.de/Evolution/ebb/Services/zostera .

Global delay in nascent strand DNA methylation.

Cytosine methylation is widespread among organisms and essential for mammalian development. In line with early postulations of an epigenetic role in gene regulation, symmetric CpG methylation can be mitotically propagated over many generations with extraordinarily high fidelity. Here, we combine BrdU labeling and immunoprecipitation with genome-wide bisulfite sequencing to explore the inheritance of cytosine methylation onto newly replicated DNA in human cells. Globally, we observe a pronounced lag between the copying of genetic and epigenetic information in embryonic stem cells that is reconsolidated within hours to accomplish faithful mitotic transmission. Populations of arrested cells show a global reduction of lag-induced intermediate CpG methylation when compared to proliferating cells, whereas sites of transcription factor engagement appear cell-cycle invariant. Alternatively, the cancer cell line HCT116 preserves global epigenetic heterogeneity independently of cell-cycle arrest. Taken together, our data suggest that heterogeneous methylation largely reflects asynchronous proliferation, but is intrinsic to actively engaged cis-regulatory elements and cancer.

Genome-wide tracking of dCas9-methyltransferase footprints.

In normal mammalian development cytosine methylation is essential and is directed to specific regions of the genome. Despite notable advances through mapping its genome-wide distribution, studying the direct contribution of DNA methylation to gene and genome regulation has been limited by the lack of tools for its precise manipulation. Thus, combining the targeting capability of the CRISPR-Cas9 system with an epigenetic modifier has attracted interest in the scientific community. In contrast to profiling the genome-wide cleavage of a nuclease competent Cas9, tracing the global activity of a dead Cas9 (dCas9) methyltransferase fusion protein is challenging within a highly methylated genome. Here, we report the generation and use of an engineered, methylation depleted but maintenance competent mouse ES cell line and find surprisingly ubiquitous nuclear activity of dCas9-methyltransferases. Subsequent experiments in human somatic cells refine these observations and point to an important difference between genetic and epigenetic editing tools that require unique experimental considerations.