Search

Chronic lymphocytic leukemia (CLL) is characterized by disordered DNA methylation, suggesting these epigenetic changes might play a critical role in disease onset and progression. The methyltransferase DNMT3A is a key regulator of DNA methylation. Although DNMT3A somatic mutations in CLL are rare, we found that low DNMT3A expression is associated with more aggressive disease. A conditional knockout mouse model showed that homozygous depletion of Dnmt3a from B cells results in the development of CLL with 100% penetrance at a median age of onset of 5.3 months, and heterozygous Dnmt3a depletion yields a disease penetrance of 89% with a median onset at 18.5 months, confirming its role as a haploinsufficient tumor suppressor. B1a cells were confirmed as the cell of origin of disease in this model, and Dnmt3a depletion resulted in focal hypomethylation and activation of Notch and Myc signaling. Amplification of chromosome 15 containing the Myc gene was detected in all CLL mice tested, and infiltration of high-Myc–expressing CLL cells in the spleen was observed. Notably, hyperactivation of Notch and Myc signaling was exclusively observed in the Dnmt3a CLL mice, but not in three other CLL mouse models tested (Sf3b1-Atm, Ikzf3, and MDR), and Dnmt3a-depleted CLL were sensitive to pharmacologic inhibition of Notch signaling in vitro and in vivo. Consistent with these findings, human CLL samples with lower DNMT3A expression were more sensitive to Notch inhibition than those with higher DNMT3A expression. Altogether, these results suggest that Dnmt3a depletion induces CLL that is highly dependent on activation of Notch and Myc signaling. Significance: Loss of DNMT3A expression is a driving event in CLL and is associated with aggressive disease, activation of Notch and Myc signaling, and enhanced sensitivity to Notch inhibition.

Genetic and epigenetic intra-tumoral heterogeneity cooperate to shape the evolutionary course of cancer1. Chronic lymphocytic leukaemia (CLL) is a highly informative model for cancer evolution as it undergoes substantial genetic diversification and evolution after therapy2,3. The CLL epigenome is also an important disease-defining feature4,5, and growing populations of cells in CLL diversify by stochastic changes in DNA methylation known as epimutations6. However, previous studies using bulk sequencing methods to analyse the patterns of DNA methylation were unable to determine whether epimutations affect CLL populations homogeneously. Here, to measure the epimutation rate at single-cell resolution, we applied multiplexed single-cell reduced-representation bisulfite sequencing to B cells from healthy donors and patients with CLL. We observed that the common clonal origin of CLL results in a consistently increased epimutation rate, with low variability in the cell-to-cell epimutation rate. By contrast, variable epimutation rates across healthy B cells reflect diverse evolutionary ages across the trajectory of B cell differentiation, consistent with epimutations serving as a molecular clock. Heritable epimutation information allowed us to reconstruct lineages at high-resolution with single-cell data, and to apply this directly to patient samples. The CLL lineage tree shape revealed earlier branching and longer branch lengths than in normal B cells, reflecting rapid drift after the initial malignant transformation and a greater proliferative history. Integration of single-cell bisulfite sequencing analysis with single-cell transcriptomes and genotyping confirmed that genetic subclones mapped to distinct clades, as inferred solely on the basis of epimutation information. Finally, to examine potential lineage biases during therapy, we profiled serial samples during ibrutinib-associated lymphocytosis, and identified clades of cells that were preferentially expelled from the lymph node after treatment, marked by distinct transcriptional profiles. The single-cell integration of genetic, epigenetic and transcriptional information thus charts the lineage history of CLL and its evolution with therapy. A single-cell approach is used to follow the heritable stochastic changes to DNA methylation that occur in primary chronic lymphocytic leukaemia and healthy B cells, allowing the tracing of cell lineage histories and evolution during treatment with ibrutinib.

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 preneoplastic 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 significant differentially methylated regions in the high-count MBL-to-CLL transition. Patient methylomes showed remarkable stability with natural disease and posttherapy 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 with 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 disease onset. Significance: DNA methylation data from a large cohort of patients with MBL and CLL show that epigenetic transformation emerges early and persists throughout disease stages with limited subsequent changes. Our results indicate an early role for this aberrant landscape in the normal-to-preneoplastic transition that may reflect a pan-cancer mechanism. See related commentary by Rossi, p. 6. This article is highlighted in the In This Issue feature, p. 1

Leukemic relapse remains a major barrier to successful allogeneic hematopoietic stem cell transplantation (allo-HSCT) for aggressive hematologic malignancies. The basis for relapse of advanced lymphoid malignancies remains incompletely understood and may involve escape from the graft-versus-leukemia (GvL) effect. We hypothesized that for patients with chronic lymphocytic leukemia (CLL) treated with allo-HSCT, leukemic cell-intrinsic features influence transplant outcomes by directing the evolutionary trajectories of CLL cells. Integrated genetic, transcriptomic, and epigenetic analyses of CLL cells from 10 patients revealed that the clinical kinetics of post-HSCT relapse are shaped by distinct molecular dynamics. Early relapses after allo-HSCT exhibited notable genetic stability; single CLL cell transcriptional analysis demonstrated a cellular heterogeneity that was static over time. In contrast, CLL cells relapsing late after allo-HSCT displayed notable genetic evolution and evidence of neoantigen depletion, consistent with marked single-cell transcriptional shifts that were unique to each patient. We observed a greater rate of epigenetic change for late relapses not seen in early relapses or relapses after chemotherapy alone, suggesting that the selection pressures of the GvL bottleneck are unlike those imposed by chemotherapy. No selective advantage for human leukocyte antigen (HLA) loss was observed, even when present in pretransplant subpopulations. Gain of stem cell modules was a common signature associated with leukemia relapse regardless of posttransplant relapse kinetics. These data elucidate the biological pathways that underlie GvL resistance and posttransplant relapse.

In diploid cells, maternal and paternal copies of genes usually have similar transcriptional activity. Mammalian allele-specific epigenetic mechanisms such as X-chromosome inactivation (XCI) and imprinting were historically viewed as rare exceptions to this rule. The discovery of mitotically stable monoallelic autosomal expression (MAE) a decade ago revealed an additional allele-specific mode regulating thousands of mammalian genes. However, despite its prevalence, the mechanistic basis of MAE remains unknown. To uncover the mechanism of MAE maintenance, we devised a small-molecule screen for reactivation of silenced alleles across multiple loci using targeted RNA sequencing. Contrary to previous reports, we identified DNA methylation as a key mechanism of MAE mitotic maintenance. In contrast with the binary choice of the active allele in XCI, stringent transcriptome-wide analysis revealed MAE as a regulatory mode with tunable control of allele-specific expression, dependent on the extent of DNA methylation. In a subset of MAE genes, allelic imbalance was insensitive to changes in DNA methylation, implicating additional mechanisms in MAE maintenance in these loci. Our findings identify a key mechanism of MAE maintenance, reveal tunability of this mode of gene regulation, and provide the essential platform for probing the biological role of MAE in development and disease.

Transcription factors (TFs) direct developmental transitions by binding to target DNA sequences, influencing gene expression and establishing complex gene-regultory networks. To systematically determine the molecular components that enable or constrain TF activity, we investigated the genomic occupancy of FOXA2, GATA4 and OCT4 in several cell types. Despite their classification as pioneer factors, all three TFs exhibit cell-type-specific binding, even when supraphysiologically and ectopically expressed. However, FOXA2 and GATA4 can be distinguished by low enrichment at loci that are highly occupied by these factors in alternative cell types. We find that expression of additional cofactors increases enrichment at a subset of these sites. Finally, FOXA2 occupancy and changes to DNA accessibility can occur in G1-arrested cells, but subsequent loss of DNA methylation requires DNA replication.

Concomitant activation of the Wnt pathway and suppression of Mapk signalling by two small molecule inhibitors (2i) in the presence of leukaemia inhibitory factor (LIF) (hereafter termed 2i/L) induces a naive state in mouse embryonic stem (ES) cells that resembles the inner cell mass (ICM) of the pre-implantation embryo. Since the ICM exists only transiently in vivo, it remains unclear how sustained propagation of naive ES cells in vitro affects their stability and functionality. Here we show that prolonged culture of male mouse ES cells in 2i/L results in irreversible epigenetic and genomic changes that impair their developmental potential. Furthermore, we find that female ES cells cultured in conventional serum plus LIF medium phenocopy male ES cells cultured in 2i/L. Mechanistically, we demonstrate that the inhibition of Mek1/2 is predominantly responsible for these effects, in part through the downregulation of DNA methyltransferases and their cofactors. Finally, we show that replacement of the Mek1/2 inhibitor with a Src inhibitor preserves the epigenetic and genomic integrity as well as the developmental potential of ES cells. Taken together, our data suggest that, although short-term suppression of Mek1/2 in ES cells helps to maintain an ICM-like epigenetic state, prolonged suppression results in irreversible changes that compromise their developmental potential.

Mammalian cells stably maintain high levels of DNA methylation despite expressing both positive (DNMT3A/B) and negative (TET1–3) regulators. Here, we used wildtype and TET triple knockout human embryonic stem cells (ESCs), generated DNMT3-null as well as TET and DNMT3-null pentuple knockouts, and compared methylation patterns using whole genome bisulfite sequencing (WGBS). The greatest impact on global methylation levels was observed in DNMT3-deficient cells, including reproducible focal demethylation at thousands of normally methylated loci. This demethylation depends upon TET expression and only occurs when both DNMT3s are absent. Dynamic loci are enriched for hydroxymethylcytosine and overlap with subsets of putative somatic enhancers that are methylated in ESCs and can be activated upon differentiation. We observe similar dynamics in mouse ESCs that were less frequenct in epiblast stem cells (EpiSCs) and scarce in somatic tissues, suggesting a conserved pluripotency-linked mechanism. Taken together, our data reveal tightly regulated competition between DNMT3s and TETs at thousands of somatic regulatory sequences within pluripotent cells.

Mucin-1 kidney disease, previously described as medullary cystic kidney disease type 1 (MCKD1, OMIM 174000), is an autosomal dominant tubulointerstitial kidney disease recently shown to be caused by a single-base insertion within the variable number tandem repeat region of the MUC1 gene. Because of variable age of disease onset and often subtle signs and symptoms, clinical diagnosis of mucin-1 kidney disease and differentiation from other forms of hereditary kidney disease have been difficult. The causal insertion resides in a variable number tandem repeat region with high GC content, which has made detection by standard next-generation sequencing impossible to date. The inherently difficult nature of this mutation required an alternative method for routine detection and clinical diagnosis of the disease. We therefore developed and validated a mass spectrometry–based probe extension assay with a series of internal controls to detect the insertion event using 24 previously characterized positive samples from patients with mucin-1 kidney disease and 24 control samples known to be wild type for the variant. Validation results indicate an accurate and reliable test for clinically establishing the molecular diagnosis of mucin-1 kidney disease with 100% sensitivity and specificity across 275 tests called., Carlos Slim Foundation. Slim Initiative for Genomic Medicine

Hi-Culfite, a protocol combining Hi-C and whole-genome bisulfite sequencing (WGBS), determines chromatin contacts and DNA methylation simultaneously. Hi-Culfite also reveals relationships that cannot be seen when the two assays are performed separately. For instance, we show that loci associated with open chromatin exhibit context-sensitive methylation: when their spatial neighbors lie in closed chromatin, they are much more likely to be methylated.

DNA methylation plays an essential role in mammalian genomes and expression of the responsible enzymes is tightly controlled. Deregulation of the de novo DNA methyltransferase DNMT3B is frequently observed across cancer types, yet little is known about its ectopic genomic targets. Here, we used an inducible transgenic mouse model to delineate rules for abnormal DNMT3B targeting, as well as the constraints of its activity across different cell types. Our results explain the preferential susceptibility of certain CpG islands to aberrant methylation and point to transcriptional state and the associated chromatin landscape as the strongest predictors. Although DNA methylation and H3K27me3 are usually non-overlapping at CpG islands, H3K27me3 can transiently co-occur with DNMT3B-induced DNA methylation. Our genome-wide data combined with ultra-deep locus-specific bisulfite sequencing suggest a distributive activity of ectopically expressed Dnmt3b that leads to discordant CpG island hypermethylation and provides new insights for interpreting the cancer methylome.

DNA methylation is a key epigenetic modification involved in regulating gene expression and maintaining genomic integrity. Here we inactivated all three catalytically active DNA methyltransferases (DNMTs) in human embryonic stem cells (ESCs) using CRISPR/Cas9 genome editing to further investigate the roles and genomic targets of these enzymes. Disruption of DNMT3A or DNMT3B individually as well as of both enzymes in tandem results in viable, pluripotent cell lines with distinct effects on the DNA methylation landscape, as assessed by whole-genome bisulfite sequencing. Surprisingly, in contrast to findings in mouse, deletion of DNMT1 resulted in rapid cell death in human ESCs. To overcome this immediate lethality, we generated a doxycycline-responsive tTA-DNMT1* rescue line and readily obtained homozygous DNMT1-mutant lines. However, doxycycline-mediated repression of exogenous DNMT1* initiates rapid, global loss of DNA methylation, followed by extensive cell death. Our data provide a comprehensive characterization of DNMT-mutant ESCs, including single-base genome-wide maps of the targets of these enzymes.

Notch-mediated expansion of cord blood progenitors: maintenance of transcriptional and epigenetic fidelity

Summary Pluripotent stem cells provide a powerful system to dissect the underlying molecular dynamics that regulate cell fate changes during mammalian development. Here we report the integrative analysis of genome wide binding data for 38 transcription factors with extensive epigenome and transcriptional data across the differentiation of human embryonic stem cells to the three germ layers. We describe core regulatory dynamics and show the lineage specific behavior of selected factors. In addition to the orchestrated remodeling of the chromatin landscape, we find that the binding of several transcription factors is strongly associated with specific loss of DNA methylation in one germ layer and in many cases a reciprocal gain in the other layers. Taken together, our work shows context-dependent rewiring of transcription factor binding, downstream signaling effectors, and the epigenome during human embryonic stem cell differentiation.

In mammals, the canonical somatic DNA methylation landscape is established upon specification of the embryo proper and subsequently disrupted within many cancer types1–4. However, the underlying mechanisms that direct this genome-scale transformation remain elusive, with no clear model for its systematic acquisition or potential developmental utility5,6. Here, we analysed global remethylation from the mouse preimplantation embryo into the early epiblast and extraembryonic ectoderm. We show that these two states acquire highly divergent genomic distributions with substantial disruption of bimodal, CpG density-dependent methylation in the placental progenitor7,8. The extraembryonic epigenome includes specific de novo methylation at hundreds of embryonically protected CpG island promoters, particularly those that are associated with key developmental regulators and are orthologously methylated across most human cancer types9. Our data suggest that the evolutionary innovation of extraembryonic tissues may have required co-option of DNA methylation-based suppression as an alternative to regulation by Polycomb-group proteins, which coordinate embryonic germ-layer formation in response to extraembryonic cues10. Moreover, we establish that this decision is made deterministically, downstream of promiscuously used—and frequently oncogenic—signalling pathways, via a novel combination of epigenetic cofactors. Methylation of developmental gene promoters during tumorigenesis may therefore reflect the misappropriation of an innate trajectory and the spontaneous reacquisition of a latent, developmentally encoded epigenetic landscape.

DNA methylation is a defining feature of mammalian cellular identity and is essential for normal development. Most cell types, except germ cells and pre-implantation embryos, display relatively stable DNA methylation patterns, with 70-80% of all CpGs being methylated. Despite recent advances, we still have a limited understanding of when, where and how many CpGs participate in genomic regulation. Here we report the in-depth analysis of 42 whole-genome bisulphite sequencing data sets across 30 diverse human cell and tissue types. We observe dynamic regulation for only 21.8% of autosomal CpGs within a normal developmental context, most of which are distal to transcription start sites. These dynamic CpGs co-localize with gene regulatory elements, particularly enhancers and transcription-factor-binding sites, which allow identification of key lineage-specific regulators. In addition, differentially methylated regions (DMRs) often contain single nucleotide polymorphisms associated with cell-type-related diseases as determined by genome-wide association studies. The results also highlight the general inefficiency of whole-genome bisulphite sequencing, as 70-80% of the sequencing reads across these data sets provided little or no relevant information about CpG methylation. To demonstrate further the utility of our DMR set, we use it to classify unknown samples and identify representative signature regions that recapitulate major DNA methylation dynamics. In summary, although in theory every CpG can change its methylation state, our results suggest that only a fraction does so as part of coordinated regulatory programs. Therefore, our selected DMRs can serve as a starting point to guide new, more effective reduced representation approaches to capture the most informative fraction of CpGs, as well as further pinpoint putative regulatory elements.

Recent molecular studies have revealed that, even when derived from a seemingly homogenous population, individual cells can exhibit substantial differences in gene expression, protein levels, and phenotypic output 1–5 , with important functional consequences 4,5 . Existing studies of cellular heterogeneity, however, have typically measured only a few pre-selected RNAs 1,2 or proteins 5,6 simultaneously because genomic profiling methods 3 could not be applied to single cells until very recently 7–10 . Here, we use single-cell RNA-Seq to investigate heterogeneity in the response of bone marrow derived dendritic cells (BMDCs) to lipopolysaccharide (LPS). We find extensive, and previously unobserved, bimodal variation in mRNA abundance and splicing patterns, which we validate by RNA-fluorescence in situ hybridization (RNA-FISH) for select transcripts. In particular, hundreds of key immune genes are bimodally expressed across cells, surprisingly even for genes that are very highly expressed at the population average. Moreover, splicing patterns demonstrate previously unobserved levels of heterogeneity between cells. Some of the observed bimodality can be attributed to closely related, yet distinct, known maturity states of BMDCs; other portions reflect differences in the usage of key regulatory circuits. For example, we identify a module of 137 highly variable, yet co-regulated, antiviral response genes. Using cells from knockout mice, we show that variability in this module may be propagated through an interferon feedback circuit involving the transcriptional regulators Stat2 and Irf7. Our study demonstrates the power and promise of single-cell genomics in uncovering functional diversity between cells and in deciphering cell states and circuits.

Although genetic lesions responsible for some mendelian disorders can be rapidly discovered through massively parallel sequencing of whole genomes or exomes, not all diseases readily yield to such efforts. We describe the illustrative case of the simple mendelian disorder medullary cystic kidney disease type 1 (MCKD1), mapped more than a decade ago to a 2-Mb region on chromosome 1. Ultimately, only by cloning, capillary sequencing and de novo assembly did we find that each of six families with MCKD1 harbors an equivalent but apparently independently arising mutation in sequence markedly under-represented in massively parallel sequencing data: the insertion of a single cytosine in one copy (but a different copy in each family) of the repeat unit comprising the extremely long (~1.5–5 kb), GC-rich (>80%) coding variable-number tandem repeat (VNTR) sequence in the MUC1 gene encoding mucin 1. These results provide a cautionary tale about the challenges in identifying the genes responsible for mendelian, let alone more complex, disorders through massively parallel sequencing., National Institutes of Health (U.S.) (Intramural Research Program), National Human Genome Research Institute (U.S.), Charles University (program UNCE 204011), Charles University (program PRVOUK-P24/LF1/3), Czech Republic. Ministry of Education, Youth, and Sports (grant NT13116-4/2012), Czech Republic. Ministry of Health (grant NT13116-4/2012), Czech Republic. Ministry of Health (grant LH12015), National Institutes of Health (U.S.) (Harvard Digestive Diseases Center, grant DK34854)

Background The ability to measure DNA methylation precisely and efficiently continues to drive our understanding of this modification in development and disease. Whole genome bisulfite sequencing has the advantage of theoretically capturing all cytosines in the genome at single-nucleotide resolution, but it has a number of significant practical drawbacks that become amplified with increasing sample numbers. All other technologies capture only a fraction of the cytosines that show dynamic regulation across cell and tissue types. Results Here, we present a novel hybrid selection design focusing on loci with dynamic methylation that captures a large number of differentially methylated gene-regulatory elements. We benchmarked this assay against matched whole genome data and profiled 25 human tissue samples to explore its ability to detect differentially methylated regions. Conclusions Our target capture design fills a major gap left by all other assays that exist to map DNA methylation. It maintains the ability to link cytosine methylation to genetic differences, the single-base resolution and the analysis of neighboring cytosines while notably reducing the cost per sample by focusing the sequencing effort on the most informative and relevant regions of the genome. Electronic supplementary material The online version of this article (doi:10.1186/s13072-016-0105-1) contains supplementary material, which is available to authorized users.

We recently used in situ Hi-C to create kilobase-resolution 3D maps of mammalian genomes. Here, we combine these maps with new Hi-C, microscopy, and genome-editing experiments to study the physical structure of chromatin fibers, domains, and loops. We find that the observed contact domains are inconsistent with the equilibrium state for an ordinary condensed polymer. Combining Hi-C data and novel mathematical theorems, we show that contact domains are also not consistent with a fractal globule. Instead, we use physical simulations to study two models of genome folding. In one, intermonomer attraction during polymer condensation leads to formation of an anisotropic "tension globule." In the other, CCCTC-binding factor (CTCF) and cohesin act together to extrude unknotted loops during interphase. Both models are consistent with the observed contact domains and with the observation that contact domains tend to form inside loops. However, the extrusion model explains a far wider array of observations, such as why loops tend not to overlap and why the CTCF-binding motifs at pairs of loop anchors lie in the convergent orientation. Finally, we perform 13 genome-editing experiments examining the effect of altering CTCF-binding sites on chromatin folding. The convergent rule correctly predicts the affected loops in every case. Moreover, the extrusion model accurately predicts in silico the 3D maps resulting from each experiment using only the location of CTCF-binding sites in the WT. Thus, we show that it is possible to disrupt, restore, and move loops and domains using targeted mutations as small as a single base pair.

Eliminating the bacterial cloning step has been a major factor in the vastly improved efficiency of massively parallel sequencing approaches. However, this also has made it a technical challenge to produce the modern equivalent of the Fosmid- or BAC-end sequences that were crucial for assembling and analyzing complex genomes during the Sanger-based sequencing era. To close this technology gap, we developed Fosill, a method for converting Fosmids to Illumina-compatible jumping libraries. We constructed Fosmid libraries in vectors with Illumina primer sequences and specific nicking sites flanking the cloning site. Our family of pFosill vectors allows multiplex Fosmid cloning of end-tagged genomic fragments without physical size selection and is compatible with standard and multiplex paired-end Illumina sequencing. To excise the bulk of each cloned insert, we introduced two nicks in the vector, translated them into the inserts, and cleaved them. Recircularization of the vector via coligation of insert termini followed by inverse PCR generates a jumping library for paired-end sequencing with 101-base reads. The yield of unique Fosmid-sized jumps is sufficiently high, and the background of short, incorrectly spaced and chimeric artifacts sufficiently low, to enable applications such as mapping of structural variation and scaffolding of de novo assemblies. We demonstrate the power of Fosill to map genome rearrangements in a cancer cell line and identified three fusion genes that were corroborated by RNA-seq data. Our Fosill-powered assembly of the mouse genome has an N50 scaffold length of 17.0 Mb, rivaling the connectivity (16.9 Mb) of the Sanger-sequencing based draft assembly.

DNA methylation is highly dynamic during mammalian embryogenesis. It is broadly accepted that the paternal genome is actively depleted of 5-methylcytosine at fertilization, followed by passive loss that reaches a minimum at the blastocyst stage. However, this model is based on limited data, and so far no base-resolution maps exist to support and refine it. Here we generate genome-scale DNA methylation maps in mouse gametes and from the zygote through post-implantation. We find that the oocyte already exhibits global hypomethylation, particularly at specific families of long interspersed element 1 and long terminal repeat retroelements, which are disparately methylated between gametes and have lower methylation values in the zygote than in sperm. Surprisingly, the oocyte contributes a unique set of differentially methylated regions (DMRs)--including many CpG island promoters--that are maintained in the early embryo but are lost upon specification and absent from somatic cells. In contrast, sperm-contributed DMRs are largely intergenic and become hypermethylated after the blastocyst stage. Our data provide a genome-scale, base-resolution timeline of DNA methylation in the pre-specified embryo, when this epigenetic modification is most dynamic, before returning to the canonical somatic pattern.

An improved understanding of enhancers in mammalian genomes could facilitate the design of new regulatory elements. Melnikov et al. synthesize thousands of ~90 nt enhancer variants, assay their activity in human cells and use the data to rationally optimize synthetic enhancers.

Massively parallel sequencing of cDNA has enabled deep and efficient probing of transcriptomes. Current approaches for transcript reconstruction from such data often rely on aligning reads to a reference genome, and are thus unsuitable for samples with a partial or missing reference genome. Here we present the Trinity method for de novo assembly of full-length transcripts and evaluate it on samples from fission yeast, mouse and whitefly, whose reference genome is not yet available. By efficiently constructing and analyzing sets of de Bruijn graphs, Trinity fully reconstructs a large fraction of transcripts, including alternatively spliced isoforms and transcripts from recently duplicated genes. Compared with other de novo transcriptome assemblers, Trinity recovers more full-length transcripts across a broad range of expression levels, with a sensitivity similar to methods that rely on genome alignments. Our approach provides a unified solution for transcriptome reconstruction in any sample, especially in the absence of a reference genome.

SummaryDespite rapid progress in characterizing transcription factor-driven reprogramming of somatic cells to an induced pluripotent stem cell (iPSC) state, many mechanistic questions still remain. To gain insight into the earliest events in the reprogramming process, we systematically analyzed the transcriptional and epigenetic changes that occur during early factor induction after discrete numbers of divisions. We observed rapid, genome-wide changes in the euchromatic histone modification, H3K4me2, at more than a thousand loci including large subsets of pluripotency-related or developmentally regulated gene promoters and enhancers. In contrast, patterns of the repressive H3K27me3 modification remained largely unchanged except for focused depletion specifically at positions where H3K4 methylation is gained. These chromatin regulatory events precede transcriptional changes within the corresponding loci. Our data provide evidence for an early, organized, and population-wide epigenetic response to ectopic reprogramming factors that clarify the temporal order through which somatic identity is reset during reprogramming.

Genetic, epigenetic and transcriptional heterogeneity cooperate to fuel cancer's ability to evolve and adapt to therapy. In chronic lymphocytic leukemia (CLL), we have shown through bulk DNA methylation (DNAme) sequencing that the growing CLL populations diversify through stochastic DNAme changes (epimutations), impacting transcriptional heterogeneity, clonal evolution and clinical outcome. To directly integrate across epigenetic, genetic and transcriptional intra-leukemic cell-to-cell variation, we developed a high-throughput multi-modality platform that jointly interrogates the methylome, transcriptome and genetic driver mutations from the same single cell (Fig. A). We applied it to >2,000 B cells from 6 healthy donors and 12 CLL samples. We found that the common clonal CLL origin resulted in an elevated but uniform epimutation rate (i.e., low cell-to-cell epimutation variability). In contrast, in index sorted B cell subsets, ranging in maturity from naïve B cells (CD27-IgM+IgD+++IgG-) to memory B cells (CD27+IgG+), variable epimutation rates reflect cells with diverse evolutionary ages across the B cell differentiation trajectory (Fig. B). Thus, we posited that epimutation can serve as a "molecular clock", enabling high-resolution lineage reconstruction, applicable directly to patient samples. CLL lineage tree topology revealed earlier branching, and longer branch lengths than normal B cells, consistent with rapid drift after transformation, and a greater proliferative history (Fig. C). In contrast to CLL topologies, non-clonal normal B cell trees provided a smaller increase in clustering accuracy compared with parsimony-based trees (P < 0.001; Fig. D). To validate the inferred tree topology, we leveraged our multi-modal capture of DNAme and genetic drivers in single cells. Indeed, genetic subclones mapped accurately to distinct clades inferred solely based on epimutation information [e.g., a clade composed of SF3B1 mutants and another clade composed of SF3B1 wild-type cells (P = 7 x 10-9; Fig. E, F)]. Using the joint single-cell transcriptional profiling, we found that cells in SF3B1 mutated clade also displayed higher alternative 3' splicing than their wild-type counterparts (P = 0.01526; Fig. G). Cells belonging to SF3B1 mutated clade were marked by expression changes in genes related to DNA damage (e.g., KLF8) and Notch signaling (e.g., DTX4) (P < 0.05; Fig. H, I). The direct linking of single-cell transcriptional data with lineage identity also showed that transcriptional similarity between cells decreases as a function of their lineage distance (P < 0.05; Fig. K). Notably, the molecular clock feature of epimutations enabled precise timing of subclonal divergence event in the CLL's evolutionary history, estimated to have occurred 2180±219 days after the emergence of the parental clone (Fig. J). To examine potential lineage biases during therapy, we performed serial multimodal single-cell profiling of a CLL patient without subclonal genetic drivers, prior to and during ibrutinib-associated lymphocytosis. The lineage trees revealed a distinct clade of cells preferentially expelled from the lymph node, marked by a distinct transcriptional profile, likely representing ibrutinib-sensitive cells. Lastly, we hypothesized that frequent epigenetic modifier mutations seen in hematological malignancies may increase the epimutation rate promoting intra-leukemic cellular diversity. To test this hypothesis, we applied our multi-modality single-cell platform to cells from TET2 knock-out (KO) mouse models, and observed higher epimutation rates, closely associated with higher cell-to-cell transcriptional heterogeneity compared with wild-type cells (P < 0.0001; Fig. L, M). In summary, we revealed that CLL show uniformly elevated epimutation rates, reflecting the common evolutionary age of these cells. By leveraging the heritable information captured through epimutation data, we provided a native lineage tracing applicable directly to patient samples. With this approach, we showed that CLL lineage topology exhibit early branching and long-branch length, consistent with exponential growth. Finally, we demonstrated that multimodal single-cell profiling enables projection of genetic and transcriptional identity onto the lineage tree, providing a direct measurement of the heritability of transcriptional profiles as a function of lineage distance. Figure. Figure. Disclosures Wu: Neon Therapeutics: Equity Ownership.

available in PMC 2010 November 2., RNA-Seq provides an unbiased way to study a transcriptome, including both coding and noncoding genes. To date, most RNA-Seq studies have critically depended on existing annotations, and thus focused on expression levels and variation in known transcripts. Here, we present Scripture, a method to reconstruct the transcriptome of a mammalian cell using only RNA-Seq reads and the genome sequence. We apply it to mouse embryonic stem cells, neuronal precursor cells, and lung fibroblasts to accurately reconstruct the full-length gene structures for the vast majority of known expressed genes. We identify substantial variation in protein-coding genes, including thousands of novel 5′-start sites, 3′-ends, and internal coding exons. We then determine the gene structures of over a thousand lincRNA and antisense loci. Our results open the way to direct experimental manipulation of thousands of non-coding RNAs, and demonstrate the power of ab initio reconstruction to render a comprehensive picture of mammalian transcriptomes., Merkin Family Foundation for Stem Cell Research, Howard Hughes Medical Institute, National Human Genome Research Institute (U.S.), Burroughs Wellcome Fund, Broad Institute

Chromosomal Mapping The conformation of the genome in the nucleus and contacts between both proximal and distal loci influence gene expression. In order to map genomic contacts, Lieberman-Aiden et al. (p. 289 , see the cover) developed a technique to allow the detection of all interactions between genomic loci in the eukaryotic nucleus followed by deep sequencing. This technology was used to map the organization of the human genome and to examine the spatial proximity of chromosomal loci at one megabase resolution. The map suggests that the genome is partitioned into two spatial compartments that are related to local chromatin state and whose remodeling correlates with changes in the chromatin state.

Maternal-Effect Dominant Embryonic Arrest (“ Medea ”) factors are selfish nuclear elements that combine maternal-lethal and zygotic-rescue activities to gain a postzygotic survival advantage. We show that Medea 1 activity in Tribolium castaneum is associated with a composite Tc1 transposon inserted just downstream of the neurotransmitter reuptake symporter bloated tubules ( blot ), whose Drosophila ortholog has both maternal and zygotic functions. The 21.5-kb insertion contains defective copies of elongation initiation factor-3 , ATP synthase subunit C , and an RNaseD -related gene, as well as a potentially intact copy of a prokaryotic DUF1703 gene. Sequence comparisons suggest that the current distribution of Medea 1 reflects global emanation after a single transpositional event in recent evolutionary time. The Medea system in Tribolium represents an unusual type of intragenomic conflict and could provide a useful vehicle for driving desirable genes into populations.

Conserved noncoding elements (CNEs) constitute the majority of sequences under purifying selection in the human genome, yet their function remains largely unknown. Experimental evidence suggests that many of these elements play regulatory roles, but little is known about regulatory motifs contained within them. Here we describe a systematic approach to discover and characterize regulatory motifs within mammalian CNEs by searching for long motifs (12–22 nt) with significant enrichment in CNEs and studying their biochemical and genomic properties. Our analysis identifies 233 long motifs (LMs), matching a total of ≈60,000 conserved instances across the human genome. These motifs include 16 previously known regulatory elements, such as the histone 3′-UTR motif and the neuron-restrictive silencer element, as well as striking examples of novel functional elements. The most highly enriched motif (LM1) corresponds to the X-box motif known from yeast and nematode. We show that it is bound by the RFX1 protein and identify thousands of conserved motif instances, suggesting a broad role for the RFX family in gene regulation. A second group of motifs (LM2*) does not match any previously known motif. We demonstrate by biochemical and computational methods that it defines a binding site for the CTCF protein, which is involved in insulator function to limit the spread of gene activation. We identify nearly 15,000 conserved sites that likely serve as insulators, and we show that nearby genes separated by predicted CTCF sites show markedly reduced correlation in gene expression. These sites may thus partition the human genome into domains of expression.

Although RNA-seq is a powerful tool, the considerable time and cost associated with library construction has limited its utilization for various applications. RNAtag-Seq, an approach to generate multiple RNA-seq libraries in a single reaction, lowers time and cost per sample, and it produces data on prokaryotic and eukaryotic samples that are comparable to those generated by traditional strand-specific RNA-seq approaches.

Human pluripotent stem cell derived models that accurately recapitulate neural development in vitro and allow for the generation of specific neuronal subtypes are of major interest to the stem cell and biomedical community. Notch signaling, particularly through the Notch effector HES5, is a major pathway critical for the onset and maintenance of neural progenitor cells (NPCs) in the embryonic and adult nervous system1-3. This can be exploited to isolate distinct populations of human embryonic stem (ES) cell derived NPCs4. Here, we report the transcriptional and epigenomic analysis of six consecutive stages derived from a HES5-GFP reporter ES cell line5 differentiated along the neural trajectory aimed at modeling key cell fate decisions including specification, expansion and patterning during the ontogeny of cortical neural stem and progenitor cells. In order to dissect the regulatory mechanisms that orchestrate the stage-specific differentiation process, we developed a computational framework to infer key regulators of each cell state transition based on the progressive remodeling of the epigenetic landscape and then validated these through a pooled shRNA screen. We were also able to refine our previous observations on epigenetic priming at transcription factor binding sites and show here that they are mediated by combinations of core and stage- specific factors. Taken together, we demonstrate the utility of our system and outline a general framework, not limited to the context of the neural lineage, to dissect regulatory circuits of differentiation.

Deep mutational scanning has emerged as a promising tool for mapping sequence-activity relationships in proteins, RNA and DNA. In this approach, diverse variants of a sequence of interest are first ranked according to their activities in a relevant pooled assay, and this ranking is then used to infer the shape of the fitness landscape around the wild-type sequence. Little is currently know, however, about the degree to which such fitness landscapes are dependent on the specific assay conditions from which they are inferred. To explore this issue, we performed deep mutational scanning of APH(3')II, a Tn5 transposon-derived kinase that confers resistance to aminoglycoside antibiotics, in E. coli under selection with each of six structurally diverse antibiotics at a range of inhibitory concentrations. We found that the resulting fitness landscapes showed significant dependence on both antibiotic structure and concentration. This shows that the notion of essential amino acid residues is context-dependent, but also that this dependence can be exploited to guide protein engineering. Specifically, we found that differential analysis of fitness landscapes allowed us to generate synthetic APH(3')II variants with orthogonal substrate specificities.

YAC-based and bacterial-clone based STS-content maps were constructed that served as the framework physical maps for the positional cloning of a candidate gene for hereditary hemochromatosis. The YAC-based map comprises 43 YACs and 86 STS and spans ∼8 Mb of DNA between the class I region of the major histocompatibility complex on human chromosome 6p21.3 and D6S276 in 6p22. Comparison with published maps revealed a hole in the MIT/Whitehead and CEPH YAC maps that includes the immediate region around the hemochromatosis gene itself. Approximately 3 Mb of DNA was covered by a bacterial clone contig that consists of 38 BACs, 45 PACs, 26 P1 clones and one λ phage. The bacterial clone-based STS map comprises 153 STSs. A contiguous block of 8 STSs could be amplified from both human chromosome 6 and 5. Further characterization of selected STSs and bacterial clones by radiation hybrid mapping and fluorescence in situ hybridization, respectively, revealed the presence of a multicopy DNA segment, more than one bacterial clone length in size, which is duplicated near the chromosome-6 centromere and part of which is present in multiple copies on chromosome 5. Possible implications of the incomplete public YAC–contig map and of the multicopy segment for physical mapping and linkage disequilibrium studies of the hemochromatosis candidate region are discussed.[The sequence data described in this paper have been submitted to GenBank under accession nos. G31205–G31325.]

International audience; Understanding the extent of genomic transcription and its functional relevance is a central goal in genomics research. However, detailed genome-wide investigations of transcriptome complexity in major mammalian organs have been scarce. Here, using extensive RNA-seq data, we show that transcription of the genome is substantially more widespread in the testis than in other organs across representative mammals. Furthermore, we reveal that meiotic spermatocytes and especially postmeiotic round spermatids have remarkably diverse transcriptomes, which explains the high transcriptome complexity of the testis as a whole. The widespread transcriptional activity in spermatocytes and spermatids encompasses protein-coding and long noncoding RNA genes but also poorly conserves intergenic sequences, suggesting that it may not be of immediate functional relevance. Rather, our analyses of genome-wide epigenetic data suggest that this prevalent transcription, which most likely promoted the birth of new genes during evolution, is facilitated by an overall permissive chromatin in these germ cells that results from extensive chromatin remodeling.

SummaryDifferentiation of human embryonic stem cells (hESCs) provides a unique opportunity to study the regulatory mechanisms that facilitate cellular transitions in a human context. To that end, we performed comprehensive transcriptional and epigenetic profiling of populations derived through directed differentiation of hESCs representing each of the three embryonic germ layers. Integration of whole-genome bisulfite sequencing, chromatin immunoprecipitation sequencing, and RNA sequencing reveals unique events associated with specification toward each lineage. Lineage-specific dynamic alterations in DNA methylation and H3K4me1 are evident at putative distal regulatory elements that are frequently bound by pluripotency factors in the undifferentiated hESCs. In addition, we identified germ-layer-specific H3K27me3 enrichment at sites exhibiting high DNA methylation in the undifferentiated state. A better understanding of these initial specification events will facilitate identification of deficiencies in current approaches, leading to more faithful differentiation strategies as well as providing insights into the rewiring of human regulatory programs during cellular transitions.

RNA-seq is an effective method for studying the transcriptome, but it can be difficult to apply to scarce or degraded RNA from fixed clinical samples, rare cell populations or cadavers. Recent studies have proposed several methods for RNA-seq of low-quality and/or low-quantity samples, but the relative merits of these methods have not been systematically analyzed. Here we compare five such methods using metrics relevant to transcriptome annotation, transcript discovery and gene expression. Using a single human RNA sample, we constructed and sequenced ten libraries with these methods and compared them against two control libraries. We found that the RNase H method performed best for chemically fragmented, low-quality RNA, and we confirmed this through analysis of actual degraded samples. RNase H can even effectively replace oligo(dT)-based methods for standard RNA-seq. SMART and NuGEN had distinct strengths for measuring low-quantity RNA. Our analysis allows biologists to select the most suitable methods and provides a benchmark for future method development.

Hereditary haemochromatosis (HH), which affects some 1 in 400 and has an estimated carrier frequency of 1 in 10 individuals of Northern European descent, results in multi-organ dysfunction caused by increased iron deposition, and is treatable if detected early. Using linkage-disequilibrium and full haplotype analysis, we have identified a 250-kilobase region more than 3 megabases telomeric of the major histocompatibility complex (MHC) that is identical-by-descent in 85% of patient chromosomes. Within this region, we have identified a gene related to the MHC class I family, termed HLA-H, containing two missense alterations. One of these is predicted to inactivate this class of proteins and was found homozygous in 83% of 178 patients. A role of this gene in haemochromatosis is supported by the frequency and nature of the major mutation and prior studies implicating MHC class I-like proteins in iron metabolism.

Genetic and epigenetic heterogeneity of cancer cells fundamentally shapes cancer progression and relapse. In chronic lymphocytic leukemia (CLL), we previously reported that intra-leukemic epigenetic diversity in DNA methylation (DNAme) follows a stochastic pattern reminiscent of genetic 'trial and error' in cancer evolution. We measured stochastic DNAme from bulk sequencing by observing the proportion of discordantly methylated sequencing reads (PDR), and found that higher PDR was associated with greater cell-to-cell transcriptional heterogeneity and adverse clinical outcome. However, bulk DNAme sequencing does not allow to phase stochastic DNAme across distant genomic loci for a single cell, as it is limited to length of a short sequencing read. Thus, it can only provide the average PDR for a population of cells, rather than the cell-to-cell variation in this important epigenetic feature. Therefore, to define stochastic DNAme changes at the single cell level, we optimized a multiplexed single cell reduced representation bisulfite sequencing (MscRRBS) protocol to allow high-throughput single cell DNAme sequencing. This protocol significantly improves scalability by multiplexing cells with shorter inline-barcodes at an early stage and utilizing SPRI beads purification to eliminate adapter-dimers. We applied MscRRBS to 393 single CD19pos B cells from two healthy volunteers and 111 single cells from a CLL sample. 88% of cells were evaluable with greater than 100,000 covered CpGs (average of 436,230 CpGs per cell). We achieved bisulfite conversion rates of 99.7%+/-0.0001 (mean+/- SD), without a significant reduction in coverage. A downsampling analysis showed that 2.1 million reads per cell provided 85% of CpG coverage with only marginal increase in coverage with further sequencing. Biallelic coverage was observed in 4.6+/-2% of germline SNPs. With MscRRBS, we measure the PDR of each individual cell. As expected from our prior bulk RRBS analysis, we found that the average PDR across cells was higher in CLL compared with B cells from healthy adult volunteers (0.39+/-0.01 vs. 0.26+/-0.08, P The higher uniformity of PDR in CLL may reflect the relationship between epi-mutation rate (measured through PDR) and the evolutionary age of the cells. As additional stochastic DNAme changes are generated with each generation, we hypothesize that PDR estimates the number of generations in a cells' evolutionary history. Thus, CLL cells have uniformly high PDR reflecting a high but uniform number of generations in their history, consistent with a single common cell of origin. In contrast, normal B cells have diverse histories with newly formed naive cells intermixed with long-lived memory cells. Consistent with this hypothesis, healthy donor CD27neg cells showed lower PDR and less PDR cell-to-cell variation compared with CD27pos cells (mean absolute PDR difference 0.008+/-0.01 vs. 0.03+/-0.03, respectively, P As MscRRBS allows complete phasing of DNAme across distant genomic loci, we further calculated the odds ratio of concordance in methylation state between any two neighboring CpG as a function of their genomic distance, compared with any two randomly paired CpGs. Through this procedure we quantified the properties of DNAme concordance decay and uncovered different genomic scales of stochastically disordered methylation. Finally, the MscRRBS also allowed to reconstruct phylo(epi)genetic relationship between the cells, and provided accurate estimate of the rates of stochastic epi-mutation across the genome. Thus, single cell DNAme is a powerful novel tool to define epigenetic diversification and its impact on CLL evolution. Disclosures No relevant conflicts of interest to declare.

Understanding the principles governing mammalian gene regulation has been hampered by the difficulty in measuring in vivo binding dynamics of large numbers of transcription factors (TF) to DNA. Here, we develop a high-throughput Chromatin ImmunoPrecipitation (HT-ChIP) method to systematically map protein-DNA interactions. HT-ChIP was applied to define the dynamics of DNA binding by 25 TFs and 4 chromatin marks at 4 time-points following pathogen stimulus of dendritic cells. Analyzing over 180,000 TF-DNA interactions we find that TFs vary substantially in their temporal binding landscapes. This data suggests a model for transcription regulation whereby TF networks are hierarchically organized into cell differentiation factors, factors that bind targets prior to stimulus to prime them for induction, and factors that regulate specific gene programs. Overlaying HT-ChIP data on gene-expression dynamics shows that many TF-DNA interactions are established prior to the stimuli, predominantly at immediate-early genes, and identified specific TF ensembles that coordinately regulate gene-induction.

Cross-talk between DNA methylation and histone modifications drives the establishment of composite epigenetic signatures and is traditionally studied using correlative rather than direct approaches. Here, we present sequential ChIP-bisulfite-sequencing (ChIP-BS-seq) as an approach to quantitatively assess DNA methylation patterns associated with chromatin modifications or chromatin-associated factors directly. A chromatin-immunoprecipitation (ChIP)-capturing step is used to obtain a restricted representation of the genome occupied by the epigenetic feature of interest, for which a single-base resolution DNA methylation map is then generated. When applied to H3 lysine 27 trimethylation (H3K27me3), we found that H3K27me3 and DNA methylation are compatible throughout most of the genome, except for CpG islands, where these two marks are mutually exclusive. Further ChIP-BS-seq-based analysis in Dnmt triple-knockout (TKO) embryonic stem cells revealed that total loss of CpG methylation is associated with alteration of H3K27me3 levels throughout the genome: H3K27me3 in localized peaks is decreased while broad local enrichments (BLOCs) of H3K27me3 are formed. At an even broader scale, these BLOCs correspond to regions of high DNA methylation in wild-type ES cells, suggesting that DNA methylation prevents H3K27me3 deposition locally and at a megabase scale. Our strategy provides a unique way of investigating global interdependencies between DNA methylation and other chromatin features.

Clone-based genome maps can be constructed by determining the presence or absence of sequence-tagged sites (STSs) in a redundant collection of yeast artificial chromosome clones (YACs). While STS-content mapping has proven to be an effective means of ordering clone ends and STSs along chromosomes, the exact physical map positions of these landmarks are not determined. This fundamental weakness can be overcome by RecA-assisted restriction endonuclease (RARE) cleavage, a method that exploits the binding specificity on duplex DNA of a RecA-protein-oligodeoxynucleotide complex to enhance the cleavage specificity of a restriction endonuclease. This technique allows selective cleavage at individual members of a large set of restriction sites. RARE-cleavage mapping was applied to a contig comprising 5 overlapping YACs spanning 580 kb on human chromosome 14. An STS-content map comprising 10 YAC-end specific STSs and one internal STS was constructed. RARE cleavage was performed on 2 YACs that span the entire contig at the EcoRI sites defining the vector-insert junctions of all 5 YACs, as well as at a HhaI site within the STS that was initially used to screen the YAC library for the clones in the contig. The sizes of the RARE-cleavage fragments were measured by pulsed-field gel electrophoresis and used to convert the STS-content map into a true physical map that indicates precise positions of clone ends and STSs.

available in PMC 2011 November 01., Cellular RNA levels are determined by the interplay of RNA production, processing and degradation. However, because most studies of RNA regulation do not distinguish the separate contributions of these processes, little is known about how they are temporally integrated. Here we combine metabolic labeling of RNA at high temporal resolution with advanced RNA quantification and computational modeling to estimate RNA transcription and degradation rates during the response of mouse dendritic cells to lipopolysaccharide. We find that changes in transcription rates determine the majority of temporal changes in RNA levels, but that changes in degradation rates are important for shaping sharp 'peaked' responses. We used sequencing of the newly transcribed RNA population to estimate temporally constant RNA processing and degradation rates genome wide. Degradation rates vary significantly between genes and contribute to the observed differences in the dynamic response. Certain transcripts, including those encoding cytokines and transcription factors, mature faster. Our study provides a quantitative approach to study the integrative process of RNA regulation., Human Frontier Science Program (Strasbourg, France), Howard Hughes Medical Institute, Burroughs Wellcome Fund (Career Award at the Scientific Interface)

Analysis of DNA methylation patterns relies increasingly on sequencing-based profiling methods. The four most frequently used sequencing-based technologies are the bisulfite-based methods MethylC-seq and reduced representation bisulfite sequencing (RRBS), and the enrichment-based techniques methylated DNA immunoprecipitation sequencing (MeDIP-seq) and methylated DNA binding domain sequencing (MBD-seq). We applied all four methods to biological replicates of human embryonic stem cells to assess their genome-wide CpG coverage, resolution, cost, concordance and the influence of CpG density and genomic context. The methylation levels assessed by the two bisulfite methods were concordant (their difference did not exceed a given threshold) for 82% for CpGs and 99% of the non-CpG cytosines. Using binary methylation calls, the two enrichment methods were 99% concordant and regions assessed by all four methods were 97% concordant. We combined MeDIP-seq with methylation-sensitive restriction enzyme (MRE-seq) sequencing for comprehensive methylome coverage at lower cost. This, along with RNA-seq and ChIP-seq of the ES cells enabled us to detect regions with allele-specific epigenetic states, identifying most known imprinted regions and new loci with monoallelic epigenetic marks and monoallelic expression.

This unit describes a protocol for the targeted enrichment of exons from randomly sheared genomic DNA libraries using an in-solution hybrid selection approach for sequencing on an Illumina Genome Analyzer II. The steps for designing and ordering a hybrid selection oligo pool are reviewed, as are critical steps for performing the preparation and hybrid selection of an Illumina paired-end library. Critical parameters, performance metrics, and analysis workflow are discussed.

The three-dimensional folding of chromosomes compartmentalizes the genome and and can bring distant functional elements, such as promoters and enhancers, into close spatial proximity (2-6). Deciphering the relationship between chromosome organization and genome activity will aid in understanding genomic processes, like transcription and replication. However, little is known about how chromosomes fold. Microscopy is unable to distinguish large numbers of loci simultaneously or at high resolution. To date, the detection of chromosomal interactions using chromosome conformation capture (3C) and its subsequent adaptations required the choice of a set of target loci, making genome-wide studies impossible (7-10). We developed Hi-C, an extension of 3C that is capable of identifying long range interactions in an unbiased, genome-wide fashion. In Hi-C, cells are fixed with formaldehyde, causing interacting loci to be bound to one another by means of covalent DNA-protein cross-links. When the DNA is subsequently fragmented with a restriction enzyme, these loci remain linked. A biotinylated residue is incorporated as the 5' overhangs are filled in. Next, blunt-end ligation is performed under dilute conditions that favor ligation events between cross-linked DNA fragments. This results in a genome-wide library of ligation products, corresponding to pairs of fragments that were originally in close proximity to each other in the nucleus. Each ligation product is marked with biotin at the site of the junction. The library is sheared, and the junctions are pulled-down with streptavidin beads. The purified junctions can subsequently be analyzed using a high-throughput sequencer, resulting in a catalog of interacting fragments. Direct analysis of the resulting contact matrix reveals numerous features of genomic organization, such as the presence of chromosome territories and the preferential association of small gene-rich chromosomes. Correlation analysis can be applied to the contact matrix, demonstrating that the human genome is segregated into two compartments: a less densely packed compartment containing open, accessible, and active chromatin and a more dense compartment containing closed, inaccessible, and inactive chromatin regions. Finally, ensemble analysis of the contact matrix, coupled with theoretical derivations and computational simulations, revealed that at the megabase scale Hi-C reveals features consistent with a fractal globule conformation.

Strand-specific, massively-parallel cDNA sequencing (RNA-Seq) is a powerful tool for novel transcript discovery, genome annotation, and expression profiling. Despite multiple published methods for strand-specific RNA-Seq, no consensus exists as to how to choose between them. Here, we developed a comprehensive computational pipeline to compare library quality metrics from any RNA-Seq method. Using the well-annotated Saccharomyces cerevisiae transcriptome as a benchmark, we compared seven library construction protocols, including both published and our own novel methods. We found marked differences in strand-specificity, library complexity, evenness and continuity of coverage, agreement with known annotations, and accuracy for expression profiling. Weighing each method’s performance and ease, we identify the dUTP second strand marking and the Illumina RNA ligation methods as the leading protocols, with the former benefitting from the current availability of paired-end sequencing. Our analysis provides a comprehensive benchmark, and our computational pipeline is applicable for assessment of future protocols in other organisms.

Background Recent studies in budding yeast have shown that antisense transcription occurs at many loci. However, the functional role of antisense transcripts has been demonstrated only in a few cases and it has been suggested that most antisense transcripts may result from promiscuous bi-directional transcription in a dense genome. Results Here, we use strand-specific RNA sequencing to study anti-sense transcription in Saccharomyces cerevisiae. We detect 1,103 putative antisense transcripts expressed in mid-log phase growth, ranging from 39 short transcripts covering only the 3' UTR of sense genes to 145 long transcripts covering the entire sense open reading frame. Many of these antisense transcripts overlap sense genes that are repressed in mid-log phase and are important in stationary phase, stress response, or meiosis. We validate the differential regulation of 67 antisense transcripts and their sense targets in relevant conditions, including nutrient limitation and environmental stresses. Moreover, we show that several antisense transcripts and, in some cases, their differential expression have been conserved across five species of yeast spanning 150 million years of evolution. Divergence in the regulation of antisense transcripts to two respiratory genes coincides with the evolution of respiro-fermentation. Conclusions Our work provides support for a global and conserved role for antisense transcription in yeast gene regulation., Canadian Friends of the Hebrew University, Howard Hughes Medical Institute, Human Frontier Science Program (Strasbourg, France), Burroughs Wellcome Fund (Career Award at the Scientific Interface), National Institutes of Health (U.S.). Pioneer Award, Broad Institute of MIT and Harvard, U.S.-Israel Binational Science Foundation (BSF), National Human Genome Research Institute (U.S.), Alfred P. Sloan Foundation

Global studies of transcript structure and abundance in cancer cells enable the systematic discovery of aberrations that contribute to carcinogenesis, including gene fusions, alternative splice isoforms, and somatic mutations. We developed a systematic approach to characterize the spectrum of cancer-associated mRNA alterations through integration of transcriptomic and structural genomic data, and we applied this approach to generate new insights into melanoma biology. Using paired-end massively parallel sequencing of cDNA (RNA-seq) together with analyses of high-resolution chromosomal copy number data, we identified 11 novel melanoma gene fusions produced by underlying genomic rearrangements, as well as 12 novel readthrough transcripts. We mapped these chimeric transcripts to base-pair resolution and traced them to their genomic origins using matched chromosomal copy number information. We also used these data to discover and validate base-pair mutations that accumulated in these melanomas, revealing a surprisingly high rate of somatic mutation and lending support to the notion that point mutations constitute the major driver of melanoma progression. Taken together, these results may indicate new avenues for target discovery in melanoma, while also providing a template for large-scale transcriptome studies across many tumor types.