DNA methylation dynamics and dysregulation delineated by high-throughput profiling in the mouse

We have developed a mouse DNA methylation array that contains 296,070 probes representing the diversity of mouse DNA methylation biology. We present a mouse methylation atlas as a rich reference resource of 1,239 DNA samples encompassing distinct tissues, strains, ages, sexes, and pathologies. We describe applications for comparative epigenomics, genomic imprinting, epigenetic inhibitors, patient-derived xenograft assessment, backcross tracing, and epigenetic clocks. We dissect DNA methylation processes associated with differentiation, aging, and tumorigenesis. Notably, we find that tissue-specific methylation signatures localize to binding sites for transcription factors controlling the corresponding tissue development. Age-associated hypermethylation is enriched at regions of Polycomb repression, while hypomethylation is enhanced at regions bound by cohesin complex members. Apc ^Min/+ polyp-associated hypermethylation affects enhancers regulating intestinal differentiation, while hypomethylation targets AP-1 binding sites. This Infinium Mouse Methylation BeadChip (version MM285) is widely accessible to the research community and will accelerate high-sample-throughput studies in this important model organism.

Evaluation of whole-genome DNA methylation sequencing library preparation protocols

Background

With rapidly dropping sequencing cost, the popularity of whole-genome DNA methylation sequencing has been on the rise. Multiple library preparation protocols currently exist. We have performed 22 whole-genome DNA methylation sequencing experiments on snap frozen human samples, and extensively benchmarked common library preparation protocols for whole-genome DNA methylation sequencing, including three traditional bisulfite-based protocols and a new enzyme-based protocol. In addition, different input DNA quantities were compared for two kits compatible with a reduced starting quantity. In addition, we also present bioinformatic analysis pipelines for sequencing data from each of these library types.

Results

An assortment of metrics were collected for each kit, including raw read statistics, library quality and uniformity metrics, cytosine retention, and CpG beta value consistency between technical replicates. Overall, the NEBNext Enzymatic Methyl-seq and Swift Accel-NGS Methyl-Seq kits performed quantitatively better than the other two protocols. In addition, the NEB and Swift kits performed well at low-input amounts, validating their utility in applications where DNA is the limiting factor.

Results

The NEBNext Enzymatic Methyl-seq kit appeared to be the best option for whole-genome DNA methylation sequencing of high-quality DNA, closely followed by the Swift kit, which potentially works better for degraded samples. Further, a general bioinformatic pipeline is applicable across the four protocols, with the exception of extra trimming needed for the Swift Biosciences’s Accel-NGS Methyl-Seq protocol to remove the Adaptase sequence.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00401-y.

Partial deletion of the sulfate transporter SLC13A1 is associated with an osteochondrodysplasia in the Miniature Poodle breed

A crippling dwarfism was first described in the Miniature Poodle in Great Britain in 1956. Here, we resolve the genetic basis of this recessively inherited disorder. A case-control analysis (8:8) of genotype data from 173 k SNPs revealed a single associated locus on CFA14 (P(raw) <10(-8)). All affected dogs were homozygous for an ancestral haplotype consistent with a founder effect and an identical-by-descent mutation. Systematic failure of nine, nearly contiguous SNPs, was observed solely in affected dogs, suggesting a deletion was the causal mutation. A 130-kb deletion was confirmed both by fluorescence in situ hybridization (FISH) analysis and by cloning the physical breakpoints. The mutation was perfectly associated in all cases and obligate heterozygotes. The deletion ablated all but the first exon of SLC13A1, a sodium/sulfate symporter responsible for regulating serum levels of inorganic sulfate. Our results corroborate earlier findings from an Slc13a1 mouse knockout, which resulted in hyposulfatemia and syndromic defects. Interestingly, the metabolic disorder in Miniature Poodles appears to share more clinical signs with a spectrum of human disorders caused by SLC26A2 than with the mouse Slc13a1 model. SLC26A2 is the primary sodium-independent sulfate transporter in cartilage and bone and is important for the sulfation of proteoglycans such as aggregan. We propose that disruption of SLC13A1 in the dog similarly causes undersulfation of proteoglycans in the extracellular matrix (ECM), which impacts the conversion of cartilage to bone. A co-dominant DNA test of the deletion was developed to enable breeders to avoid producing affected dogs and to selectively eliminate the mutation from the gene pool.

Spectral karyotyping of sarcomas and fibroblasts derived from Ink4a/Arf-deficient mice reveals chromosomal instability in vitro

The Ink4a/Arf locus is functionally linked to the Rb and p53 pathways through the action of its two gene products. Mouse models null for this locus show rapid onset of cancer with a preponderance of lymphomas and sarcomas. We report on a study of cell lines derived from sarcomas arising in Ink4a/Arf null mice. The cytogenetics of these lines was monitored over the course of serial passage. Results indicate that early passage cells are relatively normal. However, after multiple passages chromosomal instability becomes apparent as evidenced by increasing tetraploidy and aneuploidy, and the concomitant loss of clonality. To further evaluate the effect of Ink4a/Arf-deficiency on chromosomal stability in vitro, we isolated Ink4a/Arf deficient primary murine embryonic fibroblasts (MEFs), serially passaged them, and analyzed their chromosomal stability by spectral karyotyping (a 24-color chromosome paint-FISH technique). We found that chromosomal instability in Ink4a/Arf deficient MEFs developed with the same timing as seen in cell lines derived from Ink4a/Arf deficient sarcomas. Thus, chromosomal instability seen in Ink4a/Arf deficient tumors in vitro may be unrelated to the original phenotype of the tumor in vivo. Therefore, interpretation of cytogenetic data from cell lines derived from Ink4a/Arf deficient tumors should be done on early passage cells.