STAT3 gain-of-function mutations connect leukemia with autoimmune disease by pathological NKG2Dhi CD8+ T cell dysregulation and accumulation.

The association between cancer and autoimmune disease is unexplained, exemplified by T cell large granular lymphocytic leukemia (T-LGL) where gain-of-function (GOF) somatic STAT3 mutations correlate with co-existing autoimmunity. To investigate whether these mutations are the cause or consequence of CD8+ T cell clonal expansions and autoimmunity, we analyzed patients and mice with germline STAT3 GOF mutations. STAT3 GOF mutations drove the accumulation of effector CD8+ T cell clones highly expressing NKG2D, the receptor for stress-induced MHC-class-I-related molecules. This subset also expressed genes for granzymes, perforin, interferon-γ, and Ccl5/Rantes and required NKG2D and the IL-15/IL-2 receptor IL2RB for maximal accumulation. Leukocyte-restricted STAT3 GOF was sufficient and CD8+ T cells were essential for lethal pathology in mice. These results demonstrate that STAT3 GOF mutations cause effector CD8+ T cell oligoclonal accumulation and that these rogue cells contribute to autoimmune pathology, supporting the hypothesis that somatic mutations in leukemia/lymphoma driver genes contribute to autoimmune disease.

DNA methylation is required to maintain both DNA replication timing precision and 3D genome organization integrity.

DNA replication timing and three-dimensional (3D) genome organization are associated with distinct epigenome patterns across large domains. However, whether alterations in the epigenome, in particular cancer-related DNA hypomethylation, affects higher-order levels of genome architecture is still unclear. Here, using Repli-Seq, single-cell Repli-Seq, and Hi-C, we show that genome-wide methylation loss is associated with both concordant loss of replication timing precision and deregulation of 3D genome organization. Notably, we find distinct disruption in 3D genome compartmentalization, striking gains in cell-to-cell replication timing heterogeneity and loss of allelic replication timing in cancer hypomethylation models, potentially through the gene deregulation of DNA replication and genome organization pathways. Finally, we identify ectopic H3K4me3-H3K9me3 domains from across large hypomethylated domains, where late replication is maintained, which we purport serves to protect against catastrophic genome reorganization and aberrant gene transcription. Our results highlight a potential role for the methylome in the maintenance of 3D genome regulation.

TAB-seq and ACE-seq Data Processing for Genome-Wide DNA hydroxymethylation Profiling.

5-Methylcytosine (5mC) is one of the most abundant and well-studied chemical DNA modifications of vertebrate genomes. 5mC plays an essential role in genome regulation including: silencing of retroelements, X chromosome inactivation, and heterochromatin stability. Furthermore, 5mC shapes the activity of cis-regulatory elements crucial for cell fate determination. TET enzymes can oxidize 5mC to form 5-hydroxymethylcytosine (5hmC), thereby adding an additional layer of complexity to the DNA methylation landscape dynamics. The advent of techniques enabling genome-wide 5hmC profiling provided critical insights into its genomic distribution, scope, and function. These methods include immunoprecipitation, chemical labeling and capture-based approaches, as well as single-nucleotide 5hmC profiling techniques such as TET-assisted bisulfite sequencing (TAB-seq) and APOBEC-coupled epigenetic sequencing (ACE-seq). Here we provide a detailed protocol for computational analysis required for the genomic alignment of TAB-seq and ACE-seq data, 5hmC calling, and statistical analysis.

Canfam_GSD: De novo chromosome-length genome assembly of the German Shepherd Dog (Canis lupus familiaris) using a combination of long reads, optical mapping, and Hi-C.

BACKGROUND: The German Shepherd Dog (GSD) is one of the most common breeds on earth and has been bred for its utility and intelligence. It is often first choice for police and military work, as well as protection, disability assistance, and search-and-rescue. Yet, GSDs are well known to be susceptible to a range of genetic diseases that can interfere with their training. Such diseases are of particular concern when they occur later in life, and fully trained animals are not able to continue their duties. FINDINGS: Here, we provide the draft genome sequence of a healthy German Shepherd female as a reference for future disease and evolutionary studies. We generated this improved canid reference genome (CanFam_GSD) utilizing a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. The GSD assembly is ∼80 times as contiguous as the current canid reference genome (20.9 vs 0.267 Mb contig N50), containing far fewer gaps (306 vs 23,876) and fewer scaffolds (429 vs 3,310) than the current canid reference genome CanFamv3.1. Two chromosomes (4 and 35) are assembled into single scaffolds with no gaps. BUSCO analyses of the genome assembly results show that 93.0% of the conserved single-copy genes are complete in the GSD assembly compared with 92.2% for CanFam v3.1. Homology-based gene annotation increases this value to ∼99%. Detailed examination of the evolutionarily important pancreatic amylase region reveals that there are most likely 7 copies of the gene, indicative of a duplication of 4 ancestral copies and the disruption of 1 copy. CONCLUSIONS: GSD genome assembly and annotation were produced with major improvement in completeness, continuity, and quality over the existing canid reference. This resource will enable further research related to canine diseases, the evolutionary relationships of canids, and other aspects of canid biology.

The DNA methylation landscape in cancer.

As one of the most abundant and well-studied epigenetic modifications, DNA methylation plays an essential role in normal development and cellular biology. Global alterations to the DNA methylation landscape contribute to alterations in the transcriptome and deregulation of cellular pathways. Indeed, improved methods to study DNA methylation patterning and dynamics at base pair resolution and across individual DNA molecules on a genome-wide scale has highlighted the scope of change to the DNA methylation landscape in disease states, particularly during tumorigenesis. More recently has been the development of DNA hydroxymethylation profiling techniques, which allows differentiation between 5mC and 5hmC profiles and provides further insights into DNA methylation dynamics and remodeling in tumorigenesis. In this review, we describe the distribution of DNA methylation and DNA hydroxymethylation in different genomic contexts, first in normal cells, and how this is altered in cancer. Finally, we discuss DNA methylation profiling technologies and the most recent advances in single-cell methods, bisulfite-free approaches and ultra-long read sequencing techniques.

Retention of paternal DNA methylome in the developing zebrafish germline.

Two waves of DNA methylation reprogramming occur during mammalian embryogenesis; during preimplantation development and during primordial germ cell (PGC) formation. However, it is currently unclear how evolutionarily conserved these processes are. Here we characterise the DNA methylomes of zebrafish PGCs at four developmental stages and identify retention of paternal epigenetic memory, in stark contrast to the findings in mammals. Gene expression profiling of zebrafish PGCs at the same developmental stages revealed that the embryonic germline is defined by a small number of markers that display strong developmental stage-specificity and that are independent of DNA methylation-mediated regulation. We identified promoters that are specifically targeted by DNA methylation in somatic and germline tissues during vertebrate embryogenesis and that are frequently misregulated in human cancers. Together, these detailed methylome and transcriptome maps of the zebrafish germline provide insight into vertebrate DNA methylation reprogramming and enhance our understanding of the relationships between germline fate acquisition and oncogenesis.

DNA Hypermethylation Encroachment at CpG Island Borders in Cancer Is Predisposed by H3K4 Monomethylation Patterns.

Promoter CpG islands are typically unmethylated in normal cells, but in cancer a proportion are subject to hypermethylation. Using methylome sequencing we identified CpG islands that display partial methylation encroachment across the 5' or 3' CpG island borders. CpG island methylation encroachment is widespread in prostate and breast cancer and commonly associates with gene suppression. We show that the pattern of H3K4me1 at CpG island borders in normal cells predicts the different modes of cancer CpG island hypermethylation. Notably, genetic manipulation of Kmt2d results in concordant alterations in H3K4me1 levels and CpG island border DNA methylation encroachment. Our findings suggest a role for H3K4me1 in the demarcation of CpG island methylation borders in normal cells, which become eroded in cancer.

Integrated epigenomic analysis stratifies chromatin remodellers into distinct functional groups.

BACKGROUND: ATP-dependent chromatin remodelling complexes are responsible for establishing and maintaining the positions of nucleosomes. Chromatin remodellers are targeted to chromatin by transcription factors and non-coding RNA to remodel the chromatin into functional states. However, the influence of chromatin remodelling on shaping the functional epigenome is not well understood. Moreover, chromatin remodellers have not been extensively explored as a collective group across two-dimensional and three-dimensional epigenomic layers. RESULTS: Here, we have integrated the genome-wide binding profiles of eight chromatin remodellers together with DNA methylation, nucleosome positioning, histone modification and Hi-C chromosomal contacts to reveal that chromatin remodellers can be stratified into two functional groups. Group 1 (BRG1, SNF2H, CHD3 and CHD4) has a clear preference for binding at 'actively marked' chromatin and Group 2 (BRM, INO80, SNF2L and CHD1) for 'repressively marked' chromatin. We find that histone modifications and chromatin architectural features, but not DNA methylation, stratify the remodellers into these functional groups. CONCLUSIONS: Our findings suggest that chromatin remodelling events are synchronous and that chromatin remodellers themselves should be considered simultaneously and not as individual entities in isolation or necessarily by structural similarity, as they are traditionally classified. Their coordinated function should be considered by preference for chromatin features in order to gain a more accurate and comprehensive picture of chromatin regulation.

Comprehensive evaluation of genome-wide 5-hydroxymethylcytosine profiling approaches in human DNA.

BACKGROUND: The discovery that 5-methylcytosine (5mC) can be oxidized to 5-hydroxymethylcytosine (5hmC) by the ten-eleven translocation (TET) proteins has prompted wide interest in the potential role of 5hmC in reshaping the mammalian DNA methylation landscape. The gold-standard bisulphite conversion technologies to study DNA methylation do not distinguish between 5mC and 5hmC. However, new approaches to mapping 5hmC genome-wide have advanced rapidly, although it is unclear how the different methods compare in accurately calling 5hmC. In this study, we provide a comparative analysis on brain DNA using three 5hmC genome-wide approaches, namely whole-genome bisulphite/oxidative bisulphite sequencing (WG Bis/OxBis-seq), Infinium HumanMethylation450 BeadChip arrays coupled with oxidative bisulphite (HM450K Bis/OxBis) and antibody-based immunoprecipitation and sequencing of hydroxymethylated DNA (hMeDIP-seq). We also perform loci-specific TET-assisted bisulphite sequencing (TAB-seq) for validation of candidate regions. RESULTS: We show that whole-genome single-base resolution approaches are advantaged in providing precise 5hmC values but require high sequencing depth to accurately measure 5hmC, as this modification is commonly in low abundance in mammalian cells. HM450K arrays coupled with oxidative bisulphite provide a cost-effective representation of 5hmC distribution, at CpG sites with 5hmC levels >~10%. However, 5hmC analysis is restricted to the genomic location of the probes, which is an important consideration as 5hmC modification is commonly enriched at enhancer elements. Finally, we show that the widely used hMeDIP-seq method provides an efficient genome-wide profile of 5hmC and shows high correlation with WG Bis/OxBis-seq 5hmC distribution in brain DNA. However, in cell line DNA with low levels of 5hmC, hMeDIP-seq-enriched regions are not detected by WG Bis/OxBis or HM450K, either suggesting misinterpretation of 5hmC calls by hMeDIP or lack of sensitivity of the latter methods. CONCLUSIONS: We highlight both the advantages and caveats of three commonly used genome-wide 5hmC profiling technologies and show that interpretation of 5hmC data can be significantly influenced by the sensitivity of methods used, especially as the levels of 5hmC are low and vary in different cell types and different genomic locations.