Association between gut microbiota and CpG island methylator phenotype in colorectal cancer.

The intestinal microbiota is an important environmental factor implicated in CRC development. Intriguingly, modulation of DNA methylation by gut microbiota has been reported in preclinical models, although the relationship between tumor-infiltrating bacteria and CIMP status is currently unexplored. In this study, we investigated tumor-associated bacteria in 203 CRC tumor cases and validated the findings using The Cancer Genome Atlas datasets. We assessed the abundance of Bacteroides fragilis, Escherichia coli, Fusobacterium nucleatum, and Klebsiella pneumoniae through qPCR analysis and observed enrichment of all four bacterial species in CRC samples. Notably, except for E. coli, all exhibited significant enrichment in cases of CIMP. This enrichment was primarily driven by a subset of cases distinguished by high levels of these bacteria, which we labeled as "Superhigh". The bacterial Superhigh status showed a significant association with CIMP (odds ratio 3.1, p-value = 0.013) and with MLH1 methylation (odds ratio 4.2, p-value = 0.0025). In TCGA CRC cases (393 tumor and 45 adj. normal), bacterial taxa information was extracted from non-human whole exome sequencing reads, and the bacterial Superhigh status was similarly associated with CIMP (odds ratio 2.9, p < 0.001) and MLH1 methylation (odds ratio 3.5, p < 0.001). Finally, 16S ribosomal RNA gene sequencing revealed high enrichment of Bergeyella spp. C. concisus, and F. canifelinum in CIMP-Positive tumor cases. Our findings highlight that specific bacterial taxa may influence DNA methylation, particularly in CpG islands, and contribute to the development and progression of CIMP in colorectal cancer.

Genome-wide screening and functional validation of methylation barriers near promoters.

CpG islands near promoters are normally unmethylated despite being surrounded by densely methylated regions. Aberrant hypermethylation of these CpG islands has been associated with the development of various human diseases. Although local genetic elements have been speculated to play a role in protecting promoters from methylation, only a limited number of methylation barriers have been identified. In this study, we conducted an integrated computational and experimental investigation of colorectal cancer methylomes. Our study revealed 610 genes with disrupted methylation barriers. Genomic sequences of these barriers shared a common 41-bp sequence motif (MB-41) that displayed homology to the chicken HS4 methylation barrier. Using the CDKN2A (P16) tumor suppressor gene promoter, we validated the protective function of MB-41 and showed that loss of such protection led to aberrant hypermethylation. Our findings highlight a novel sequence signature of cis-acting methylation barriers in the human genome that safeguard promoters from silencing.

MC180295 is a highly potent and selective CDK9 inhibitor with preclinical in vitro and in vivo efficacy in cancer.

BACKGROUND: Inhibition of cyclin-dependent kinase 9 (CDK9), a novel epigenetic target in cancer, can reactivate epigenetically silenced genes in cancer by dephosphorylating the SWI/SNF chromatin remodeler BRG1. Here, we characterized the anti-tumor efficacy of MC180295, a newly developed CDK9 inhibitor. METHODS: In this study, we explored the pharmacokinetics of MC180295 in mice and rats, and tested the anti-tumor efficacy of MC180295, and its enantiomers, in multiple cancer cell lines and mouse models. We also combined CDK9 inhibition with a DNA methyltransferase (DNMT) inhibitor, decitabine, in multiple mouse models, and tested MC180295 dependence on T cells. Drug toxicity was measured by checking body weights and complete blood counts. RESULTS: MC180295 had high specificity for CDK9 and high potency against multiple neoplastic cell lines (median IC50 of 171 nM in 46 cell lines representing 6 different malignancies), with the highest potency seen in AML cell lines derived from patients with MLL translocations. MC180295 is a racemic mixture of two enantiomers, MC180379 and MC180380, with MC180380 showing higher potency in a live-cell epigenetic assay. Both MC180295 and MC180380 showed efficacy in in vivo AML and colon cancer xenograft models, and significant synergy with decitabine in both cancer models. Lastly, we found that CDK9 inhibition-mediated anti-tumoral effects were partially dependent on CD8 + T cells in vivo, indicating a significant immune component to the response. CONCLUSIONS: MC180380, an inhibitor of cyclin-dependent kinase 9 (CDK9), is an efficacious anti-cancer agent worth advancing further toward clinical use.

Use of genome-wide DNA methylation analysis to identify prognostic CpG site markers associated with longer survival time in dogs with multicentric high-grade B-cell lymphoma.

BACKGROUND: DNA methylation analysis might identify prognostic CpG sites in CHOP-treated dogs with multicentric high-grade B-cell lymphoma (MHGL) with heterogenous prognosis. OBJECTIVE: To identify prognostic CpG sites of MHGL through genome-wide DNA methylation analysis with pyrosequencing validation. ANIMALS: Test group: 24 dogs. Validation group: 100 dogs. All client-owned dogs were diagnosed with MHGL and treated with CHOP chemotherapy. METHODS: Cohort study. DNA was extracted from lymph node samples obtained via FNA. Genome-wide DNA methylation analysis using Digital Restriction Enzyme Analysis of Methylation (DREAM) was performed on the test group to identify differentially methylated CpG sites (DMCs). Bisulfite pyrosequencing was used to measure methylation status of candidate DMCs in the validation group. Median survival times (MST) were analyzed using Kaplan-Meier (log-rank) product limit method. RESULTS: DREAM analyzed 101 576 CpG sites. Hierarchical clustering of 16 262 CpG sites in test group identified group with better prognosis (MST = 55-477 days vs 10-301 days, P = .007). Volcano plot identified 1371 differentially methylated CpG sites (DMCs). DMC near the genes of FAM213A (DMC-F) and PHLPP1 (DMC-P) were selected as candidates. Bisulfite-pyrosequencing performed on validation group showed group with methylation level of DMC-F < 40% had favorable prognosis (MST = 11-1072 days vs 8-1792 days, P = .01), whereas group with the methylation level combination of DMC-F < 40% plus DMC-P < 10% had excellent prognosis (MST = 18-1072 days vs 8-1792 days, P = .009). CONCLUSION AND CLINICAL IMPORTANCE: Methylation status of prognostic CpG sites delineate canine MGHL cases with longer MST, providing owners with information on expectations of potential improved treatment outcomes.

Postmortem toxicology findings from the Camden Opioid Research Initiative.

The United States continues to be impacted by decades of an opioid misuse epidemic, worsened by the COVID-19 pandemic and by the growing prevalence of highly potent synthetic opioids (HPSO) such as fentanyl. In instances of a toxicity event, first-response administration of reversal medications such as naloxone can be insufficient to fully counteract the effects of HPSO, particularly when there is co-occurring substance use. In an effort to characterize and study this multi-faceted problem, the Camden Opioid Research Initiative (CORI) has been formed. The CORI study has collected and analyzed post-mortem toxicology data from 42 cases of decedents who expired from opioid-related toxicity in the South New Jersey region to characterize substance use profiles. Co-occurring substance use, whether by intent or through possible contamination of the illicit opioid supply, is pervasive among deaths due to opioid toxicity, and evidence of medication-assisted treatment is scarce. Nearly all (98%) of the toxicology cases show the presence of the HPSO, fentanyl, and very few (7%) results detected evidence of medication-assisted treatment for opioid use disorder, such as buprenorphine or methadone, at the time of death. The opioid toxicity reversal drug, naloxone, was detected in 19% of cases, but 100% of cases expressed one or more stimulants, and sedatives including xylazine were detected in 48% of cases. These results showing complex substance use profiles indicate that efforts at mitigating the opioid misuse epidemic must address the complications presented by co-occurring stimulant and other substance use, and reduce barriers to and stigmas of seeking effective medication-assisted treatments.

Diverse genome-wide DNA methylation alterations in canine hepatocellular tumours.

BACKGROUND: Canine hepatocellular tumours (HCTs) are common primary liver tumours. However, the exact mechanisms of tumourigenesis remain unclear. Although some genetic mutations have been reported, DNA methylation alterations in canine HCT have not been well studied. OBJECTIVES: In this study, we aimed to analyse the DNA methylation status of canine HCT. METHODS: Tissues from 33 hepatocellular carcinomas, 3 hepatocellular adenomas, 1 nodular hyperplasia, 21 non-tumour livers from the patients and normal livers from 5 healthy dogs were used. We analysed the DNA methylation levels of 72,367 cytosine-guanine dinucleotides (CpG sites) in all 63 samples. RESULTS AND CONCLUSIONS: Although a large fraction of CpG sites that were highly methylated in the normal liver became hypomethylated in tumours from most patients, we also found some patients with less remarkable change or no change in DNA methylation. Hierarchical clustering analysis revealed that 32 of 37 tumour samples differed from normal livers, although the remaining 5 tumour livers fell into the same cluster as normal livers. In addition, the number of hypermethylated genes in tumour livers varied among tumour cases, suggesting various DNA methylation patterns in different tumour groups. However, patient and clinical parameters, such as age, were not associated with DNA methylation status. In conclusion, we found that HCTs undergo aberrant and diverse patterns of genome-wide DNA methylation compared with normal liver tissue, suggesting a complex epigenetic mechanism in canine HCT.

Folic Acid Supplementation Promotes Hypomethylation in Both the Inflamed Colonic Mucosa and Colitis-Associated Dysplasia.

PURPOSE: The purpose of this study was to assess the effect of folic acid (FA) supplementation on colitis-associated colorectal cancer (CRC) using the azoxymethane/dextran sulfate sodium (AOM/DSS) model. METHODS: Mice were fed a chow containing 2 mg/kg FA at baseline and randomized after the first DSS treatment to receive 0, 2, or 8 mg/kg FA chow for 16 weeks. Colon tissue was collected for histopathological evaluation, genome-wide methylation analyses (Digital Restriction Enzyme Assay of Methylation), and gene expression profiling (RNA-Seq). RESULTS: A dose-dependent increase in the multiplicity of colonic dysplasias was observed, with the multiplicity of total and polypoid dysplasias higher (64% and 225%, respectively) in the 8 mg FA vs. the 0 mg FA group (p < 0.001). Polypoid dysplasias were hypomethylated, as compared to the non-neoplastic colonic mucosa (p < 0.05), irrespective of FA treatment. The colonic mucosa of the 8 mg FA group was markedly hypomethylated as compared to the 0 mg FA group. Differential methylation of genes involved in Wnt/ß-catenin and MAPK signaling resulted in corresponding alterations in gene expression within the colonic mucosa. CONCLUSIONS: High-dose FA created an altered epigenetic field effect within the non-neoplastic colonic mucosa. The observed decrease in site-specific DNA methylation altered oncogenic pathways and promoted colitis-associated CRC.

DNA methylation entropy as a measure of stem cell replication and aging.

BACKGROUND: Epigenetic marks are encoded by DNA methylation and accumulate errors as organisms age. This drift correlates with lifespan, but the biology of how this occurs is still unexplained. We analyze DNA methylation with age in mouse intestinal stem cells and compare them to nonstem cells. RESULTS: Age-related changes in DNA methylation are identical in stem and nonstem cells, affect most prominently CpG islands and correlate weakly with gene expression. Age-related DNA methylation entropy, measured by the Jensen-Shannon Distribution, affects up to 25% of the detectable CpG sites and is a better measure of aging than individual CpG methylation. We analyze this entropy as a function of age in seven other tissues (heart, kidney, skeletal muscle, lung, liver, spleen, and blood) and it correlates strikingly with tissue-specific stem cell division rates. Thus, DNA methylation drift and increased entropy with age are primarily caused by and are sensors for, stem cell replication in adult tissues. CONCLUSIONS: These data have implications for the mechanisms of tissue-specific functional declines with aging and for the development of DNA-methylation-based biological clocks.

TET1 and TDG Suppress Inflammatory Response in Intestinal Tumorigenesis: Implications for Colorectal Tumors With the CpG Island Methylator Phenotype.

BACKGROUND & AIMS: Aberrant DNA methylation is frequent in colorectal cancer (CRC), but underlying mechanisms and pathologic consequences are poorly understood. METHODS: We disrupted active DNA demethylation genes Tet1 and/or Tdg from ApcMin mice and characterized the methylome and transcriptome of colonic adenomas. Data were compared to human colonic adenocarcinomas (COAD) in The Cancer Genome Atlas. RESULTS: There were increased numbers of small intestinal adenomas in ApcMin mice expressing the TdgN151A allele, whereas Tet1-deficient and Tet1/TdgN151A-double heterozygous ApcMin colonic adenomas were larger with features of erosion and invasion. We detected reduction in global DNA hypomethylation in colonic adenomas from Tet1- and Tdg-mutant ApcMin mice and hypermethylation of CpG islands in Tet1-mutant ApcMin adenomas. Up-regulation of inflammatory, immune, and interferon response genes was present in Tet1- and Tdg-mutant colonic adenomas compared to control ApcMin adenomas. This up-regulation was also seen in murine colonic organoids and human CRC lines infected with lentiviruses expressing TET1 or TDG short hairpin RNA. A 127-gene inflammatory signature separated colonic adenocarcinomas into 4 groups, closely aligned with their microsatellite or chromosomal instability and characterized by different levels of DNA methylation and DNMT1 expression that anticorrelated with TET1 expression. Tumors with the CpG island methylator phenotype (CIMP) had concerted high DNMT1/low TET1 expression. TET1 or TDG knockdown in CRC lines enhanced killing by natural killer cells. CONCLUSIONS: Our findings reveal a novel epigenetic regulation, linked to the type of genomic instability, by which TET1/TDG-mediated DNA demethylation decreases methylation levels and inflammatory/interferon/immune responses. CIMP in CRC is triggered by an imbalance of methylating activities over demethylating activities. These mice represent a model of CIMP CRC.

Novel treatment strategy of targeting epigenetic dysregulation in pancreatic neuroendocrine tumors.

BACKGROUND: Epigenetic dysregulation is an integral step in the progression of pancreatic neuroendocrine tumors. We hypothesized that tumor suppressor repression by DNA methyltransferase 1 in pancreatic neuroendocrine tumors could be targeted with epigenetic treatment. METHODS: Resected pancreatic neuroendocrine tumors from 32 patients were stained for DNA methyltransferase 1 and scored. Human (BON1) and murine (STC) pancreatic neuroendocrine tumor cells were treated with DNA methyltransferase 1 inhibitor 5-azacytidine and chemotherapeutic agents 5-fluorouracil and temozolomide. Cell proliferation assay and tumor suppressor gene analysis were performed with qRT-PCR and Clarion S microarray. Tumor measurements were compared in a murine treatment model. RESULTS: High DNA methyltransferase scores were associated with high Ki-67 (6.7% vs 70.6% P < .01), mitotic rate (0.0% vs 31.3%), and grade (20.0% vs 80.4%, P < .01). Treatment with 5-azacytidine and chemotherapy resulted in a reduction of cell proliferation compared to chemotherapy alone in BON1 (77.3% vs 53.1%, P < .001) and STC (73.4% vs 34.2%, P < .001). Treatment with 5-azacytidine and chemotherapy resulted in upregulation of tumor suppressors CDKN1A (7.6 rel. fold, P < .001), BRCA2 (4.3 rel. fold, P < .001), and CDH1 (6.0 rel. fold, P = .026) in BON1 and CDKN1a (14.5 rel. fold, P < .001) and CDH (17.5 rel. fold, P < .001) in STC. In microarray, 5-azacytidine drove global genetic changes in combination treatment. In vivo tumors treated with chemotherapy measured 88.6 ± 19.54 mm3 vs 52.89 ± 10.51 mm3 in those treated with combination therapy (P = .009). CONCLUSION: Epigenetic dysregulation with DNA methyltransferase 1 is associated with pancreatic neuroendocrine tumors and is a potential targetable strategy. 5-azacytidine and chemotherapy in combination can reduce cell proliferation, upregulate silenced tumor suppressor genes, and decrease in vivo tumors in pancreatic neuroendocrine tumors.

Hypomethylation in MTNR1B: a novel epigenetic marker for atherosclerosis profiling using stenosis radiophenotype and blood inflammatory cells.

BACKGROUND: Changes in gene-specific promoter methylation may result from aging and environmental influences. Atherosclerosis is associated with aging and environmental effects. Thus, promoter methylation profiling may be used as an epigenetic tool to evaluate the impact of aging and the environment on atherosclerosis development. However, gene-specific methylation changes are currently inadequate epigenetic markers for predicting atherosclerosis and cardiovascular disease pathogenesis. RESULTS: We profiled and validated changes in gene-specific promoter methylation associated with atherosclerosis using stenosis radiophenotypes of cranial vessels and blood inflammatory cells rather than direct sampling of atherosclerotic plaques. First, we profiled gene-specific promoter methylation changes using digital restriction enzyme analysis of methylation (DREAM) sequencing in peripheral blood mononuclear cells from eight samples each of cranial vessels with and without severe-stenosis radiophenotypes. Using DREAM sequencing profiling, 11 tags were detected in the promoter regions of the ACVR1C, ADCK5, EFNA2, ENOSF1, GLS2, KNDC1, MTNR1B, PACSIN3, PAX8-AS1, TLDC1, and ZNF7 genes. Using methylation evaluation, we found that EFNA2, ENOSF1, GLS2, KNDC1, MTNR1B, PAX8-AS1, and TLDC1 showed > 5% promoter methylation in non-plaque intima, atherosclerotic vascular tissues, and buffy coats. Using logistic regression analysis, we identified hypomethylation of MTNR1B as an independent variable for the stenosis radiophenotype prediction model by combining it with traditional atherosclerosis risk factors including age, hypertension history, and increases in creatinine, lipoprotein (a), and homocysteine. We performed fivefold cross-validation of the prediction model using 384 patients with ischemic stroke (50 [13%] no-stenosis and 334 [87%] > 1 stenosis radiophenotype). For the cross-validation, the training dataset included 70% of the dataset. The prediction model showed an accuracy of 0.887, specificity to predict stenosis radiophenotype of 0.940, sensitivity to predict no-stenosis radiophenotype of 0.533, and area under receiver operating characteristic curve of 0.877 to predict stenosis radiophenotype from the test dataset including 30% of the dataset. CONCLUSIONS: We identified and validated MTNR1B hypomethylation as an epigenetic marker to predict cranial vessel atherosclerosis using stenosis radiophenotypes and blood inflammatory cells rather than direct atherosclerotic plaque sampling.

Non-pathogenic microbiota accelerate age-related CpG Island methylation in colonic mucosa.

DNA methylation is an epigenetic process altered in cancer and ageing. Age-related methylation drift can be used to estimate lifespan and can be influenced by extrinsic factors such as diet. Here, we report that non-pathogenic microbiota accelerate age-related methylation drift in the colon when compared with germ-free mice. DNA methylation analyses showed that microbiota and IL10KO were associated with changes in 5% and 4.1% of CpG sites, while mice with both factors had 18% alterations. Microbiota, IL10KO, and their combination altered 0.4%, 0.4%, and 4% of CpG island methylation, respectively. These are comparable to what is seen in colon cancer. Ageing changes were accelerated in the IL10KO mice with microbiota, and the affected genes were more likely to be altered in colon cancer. Thus, the microbiota affect DNA methylation of the colon in patterns reminiscent of what is observed in ageing and colorectal cancer.

DNA polymerase θ protects leukemia cells from metabolically induced DNA damage.

Leukemia cells accumulate DNA damage, but altered DNA repair mechanisms protect them from apoptosis. We showed here that formaldehyde generated by serine/1-carbon cycle metabolism contributed to the accumulation of toxic DNA-protein crosslinks (DPCs) in leukemia cells, especially in driver clones harboring oncogenic tyrosine kinases (OTKs: FLT3(internal tandem duplication [ITD]), JAK2(V617F), BCR-ABL1). To counteract this effect, OTKs enhanced the expression of DNA polymerase theta (POLθ) via ERK1/2 serine/threonine kinase-dependent inhibition of c-CBL E3 ligase-mediated ubiquitination of POLθ and its proteasomal degradation. Overexpression of POLθ in OTK-positive cells resulted in the efficient repair of DPC-containing DNA double-strand breaks by POLθ-mediated end-joining. The transforming activities of OTKs and other leukemia-inducing oncogenes, especially of those causing the inhibition of BRCA1/2-mediated homologous recombination with and without concomitant inhibition of DNA-PK-dependent nonhomologous end-joining, was abrogated in Polq-/- murine bone marrow cells. Genetic and pharmacological targeting of POLθ polymerase and helicase activities revealed that both activities are promising targets in leukemia cells. Moreover, OTK inhibitors or DPC-inducing drug etoposide enhanced the antileukemia effect of POLθ inhibitor in vitro and in vivo. In conclusion, we demonstrated that POLθ plays an essential role in protecting leukemia cells from metabolically induced toxic DNA lesions triggered by formaldehyde, and it can be targeted to achieve a therapeutic effect.

Single cell transcriptomic analysis reveals cellular diversity of murine esophageal epithelium.

Although morphologic progression coupled with expression of specific molecular markers has been characterized along the esophageal squamous differentiation gradient, the molecular heterogeneity within cell types along this trajectory has yet to be classified at the single cell level. To address this knowledge gap, we perform single cell RNA-sequencing of 44,679 murine esophageal epithelial, to identify 11 distinct cell populations as well as pathways alterations along the basal-superficial axis and in each individual population. We evaluate the impact of aging upon esophageal epithelial cell populations and demonstrate age-associated mitochondrial dysfunction. We compare single cell transcriptomic profiles in 3D murine organoids and human esophageal biopsies with that of murine esophageal epithelium. Finally, we employ pseudotemporal trajectory analysis to develop a working model of cell fate determination in murine esophageal epithelium. These studies provide comprehensive molecular perspective on the cellular heterogeneity of murine esophageal epithelium in the context of homeostasis and aging.

Genome-wide DNA methylation profile in feline haematological tumours: A preliminary study.

Although DNA methylation has been analysed in few studies for a limited number of loci in cats with diseases, genome-wide profile of DNA methylation has never been addressed. The hypothesis for this study is that next-generation sequencing with sequential digestion of genomic DNA with SmaI and XmaI enzymes could provide highly quantitative information on methylation levels in cats. Using blood from four healthy control cats and two disease cats as well as three feline lymphoma/leukemia cell lines, approximately 74-94 thousand CpG sites across the cat genome could be analysed. CpG sites in CpG island (CGI) were broadly either methylated or unmethylated in normal blood, while CpG sites in non-CpG islands (NCGI) are largely methylated. Lymphoma cell lines showed thousands of CpG sites with gain of methylation at normally unmethylated CGI sites and loss of methylation at normally methylated NCGI sites. Hypermethylated CpG sites located at promoter regions included genes annotated with 'developmental process' and 'anatomical structure morphogenesis' such as HOXD10. This highly quantitative method would be suitable for studies of DNA methylation changes not only in cancer but also in other common diseases in cats.

Obese status is associated with accelerated DNA methylation change in peripheral blood of senior dogs.

Obesity and its associated comorbidities constitute a major and growing health problem worldwide not only involved with people but also dogs and cats. Although few genetic mutations have been associated with obesity in dogs, molecular mechanism remains to be clearly understood. Given the fact that DNA methylation leads to gene expression variability and has plasticity affected by metabolic phenotypes such as obesity in human, the objective of this study is to identify obesity-associated differentially methylated cytosine-phosphate-guanine (CpG) dinucleotide sites in dogs. With genome-wide DNA methylation analysis using next-generation sequencing for blood samples from fourteen Miniature dachshunds with body condition score (BCS) 4-5 and BCS ≥6, over 100,000 sites could be analysed to identify genomic locations of differentially methylated CpG sites. As a result, 191 differentially methylated CpG sites (89 CpG sites were hypermethylated in BCS ≥6 and 102 were hypermethylated in BCS 4-5) were identified. These sites included promoter regions of Kisspeptin receptor (KISS1R) and Calcyphosine 2 (CAPS2) genes which were subsequently validated by bisulfite-pyrosequencing for another set of 157 dog blood samples. KISS1R methylation levels were found to be higher in BCS ≥6 group than BCS 4-5 in senior (>84 months) dogs. Especially male dogs but not female dogs as well as uncastrated male dogs but not castrated male dogs showed this trend. DNA methylation of KISS1R gene will be useful for understanding of comprehensive epigenetic change in obese dogs.

DNA methylation landscape of 16 canine somatic tissues by methylation-sensitive restriction enzyme-based next generation sequencing.

DNA methylation plays important functions in gene expression regulation that is involved in individual development and various diseases. DNA methylation has been well studied in human and model organisms, but only limited data exist in companion animals like dog. Using methylation-sensitive restriction enzyme-based next generation sequencing (Canine DREAM), we obtained canine DNA methylation maps of 16 somatic tissues from two dogs. In total, we evaluated 130,861 CpG sites. The majority of CpG sites were either highly methylated (> 70%, 52.5-64.6% of all CpG sites analyzed) or unmethylated (< 30%, 22.5-28.0% of all CpG sites analyzed) which are methylation patterns similar to other species. The overall methylation status of CpG sites across the 32 methylomes were remarkably similar. However, the tissue types were clearly defined by principle component analysis and hierarchical clustering analysis with DNA methylome. We found 6416 CpG sites located closely at promoter region of genes and inverse correlation between DNA methylation and gene expression of these genes. Our study provides basic dataset for DNA methylation profiles in dogs.

Genome-wide DNA methylation analysis of dogs with high lead exposure living near a lead mining area in Kabwe, Zambia.

Lead (Pb) is a heavy metal that has been proven to be toxic to both animals and humans. Genom-wide DNA methylation in domestic dogs exposed to high levels of Pb in Kabwe, Zambia was analyzed in this study. Using next-generation sequencing on samples from 20 domestic dogs (mean blood Pb concentration: 43.6 µg/dL and 7.2 µg/dL in the high and low exposure groups), a digital restriction enzyme analysis of methylation was performed to identify the genomic locations of differentially methylated CpG sites. A validation study on an additional 20 dogs followed (blood Pb concentration: 4.9-29.7 µg/dL). The cluster analysis resolved two broad clusters indicating high and low Pb exposure. The study identified 827 (1.2%) CpG sites with differences in methylation (101 CpG sites were hypermethylated in the low exposure group and 726 were hypermethylated in the high exposure group). The sites corresponded to 26 genes with differentially methylated CpG sites at their promoter regions, including the NGF gene. The methylation of four CpG sites was validated using bisulfite pyrosequencing. The results indicate that aberrant hypermethylation is prevalent in dogs exposed to Pb. The altered DNA methylation of the genes identified in this study contributes to a greater understanding of the epigenetic changes caused by Pb exposure and highlights novel biomarker discoveries across species.