CGmapTools improves the precision of heterozygous SNV calls and supports allele-specific methylation detection and visualization in bisulfite-sequencing data

Distinct epigenetic features of tumor-reactive CD8+ T cells in colorectal cancer patients revealed by genome-wide DNA methylation analysis

BACKGROUND

Tumor-reactive CD8+ tumor-infiltrating lymphocytes (TILs) represent a subtype of T cells that can recognize and destroy tumor specifically. Understanding the regulatory mechanism of tumor-reactive CD8+ T cells has important therapeutic implications. Yet the DNA methylation status of this T cell subtype has not been elucidated.

RESULTS

In this study, we segregate tumor-reactive and bystander CD8+ TILs, as well as naïve and effector memory CD8+ T cell subtypes as controls from colorectal cancer patients, to compare their transcriptome and methylome characteristics. Transcriptome profiling confirms previous conclusions that tumor-reactive TILs have an exhausted tissue-resident memory signature. Whole-genome methylation profiling identifies a distinct methylome pattern of tumor-reactive CD8+ T cells, with tumor-reactive markers CD39 and CD103 being specifically demethylated. In addition, dynamic changes are observed during the transition of naïve T cells into tumor-reactive CD8+ T cells. Transcription factor binding motif enrichment analysis identifies several immune-related transcription factors, including three exhaustion-related genes (NR4A1, BATF, and EGR2) and VDR, which potentially play an important regulatory role in tumor-reactive CD8+ T cells.

CONCLUSION

Our study supports the involvement of DNA methylation in shaping tumor-reactive and bystander CD8+ TILs, and provides a valuable resource for the development of novel DNA methylation markers and future therapeutics.

IMAGE: high-powered detection of genetic effects on DNA methylation using integrated methylation QTL mapping and allele-specific analysis

Identifying genetic variants that are associated with methylation variation-an analysis commonly referred to as methylation quantitative trait locus (mQTL) mapping-is important for understanding the epigenetic mechanisms underlying genotype-trait associations. Here, we develop a statistical method, IMAGE, for mQTL mapping in sequencing-based methylation studies. IMAGE properly accounts for the count nature of bisulfite sequencing data and incorporates allele-specific methylation patterns from heterozygous individuals to enable more powerful mQTL discovery. We compare IMAGE with existing approaches through extensive simulation. We also apply IMAGE to analyze two bisulfite sequencing studies, in which IMAGE identifies more mQTL than existing approaches.

Genome-wide allele-specific methylation is enriched at gene regulatory regions in a multi-generation pedigree from the Norfolk Island isolate

Background

Allele-specific methylation (ASM) occurs when DNA methylation patterns exhibit asymmetry among alleles. ASM occurs at imprinted loci, but its presence elsewhere across the human genome is indicative of wider importance in terms of gene regulation and disease risk. Here, we studied ASM by focusing on blood-based DNA collected from 24 subjects comprising a 3-generation pedigree from the Norfolk Island genetic isolate. We applied a genome-wide bisulphite sequencing approach with a genotype-independent ASM calling method to map ASM across the genome. Regions of ASM were then tested for enrichment at gene regulatory regions using Genomic Association Test (GAT) tool.

Results

In total, we identified 1.12 M CpGs of which 147,170 (13%) exhibited ASM (P ≤ 0.05). When including contiguous ASM signal spanning ≥ 2 CpGs, this condensed to 12,761 ASM regions (AMRs). These AMRs tagged 79% of known imprinting regions and most (98.1%) co-localised with known single nucleotide variants. Notably, miRNA and lncRNA showed a 3.3- and 1.8-fold enrichment of AMRs, respectively (P < 0.005). Also, the 5′ UTR and start codons each showed a 3.5-fold enrichment of AMRs (P < 0.005). There was also enrichment of AMRs observed at subtelomeric regions of many chromosomes. Five out of 11 large AMRs localised to the protocadherin cluster on chromosome 5.

Conclusions

This study shows ASM extends far beyond genomic imprinting in humans and that gene regulatory regions are hotspots for ASM. Future studies of ASM in pedigrees should help to clarify transgenerational inheritance patterns in relation to genotype and disease phenotypes.

An in vitro attempt at precision toxicology reveals the involvement of DNA methylation alteration in ochratoxin A-induced G0/G1 phase arrest

Precision toxicology evaluates the toxicity of certain substances by isolating a small group of cells with a typical phenotype of interest followed by a single cell sequencing-based analysis. In this in vitro attempt, ochratoxin A (OTA), a typical mycotoxin and food contaminant, is found to induce G0/G1 phase cell cycle arrest in human renal proximal tubular HKC cells at a concentration of 20 μM after a 24h-treatment. A small number of G0/G1 phase HKC cells are evaluated in both the presence and absence of OTA. These cells are sorted with a flow cytometer and subjected to mRNA and DNA methylation sequencing using Smart-Seq2 and single-cell reduced-representation bisulfite sequencing (scRRBS) technology, respectively. Integrated analysis of the transcriptome and methylome profiles reveals that OTA causes abnormal expression of the essential genes that regulate G1/S phase transition, act as signal transductors in G1 DNA damage checkpoints, and associate with the anaphase-promoting complex/cyclosome. The alteration of their DNA methylation status is a significant underlying epigenetic mechanism. Furthermore, Notch signaling and Ras/MAPK/CREB pathways are found to be suppressed by OTA. This attempt at precision toxicology paves the way for a deeper understanding of OTA toxicity and provides an innovative strategy to researchers in the toxicology and pharmacology field.

Dynamic DNA methylation of ovaries during pubertal transition in gilts

Background

In female mammals, the initiation of puberty, coupling with the dramatically morphological changes in ovaries, indicates the sexual and follicular maturation. Previous studies have suggested that the disrupted DNA methylation results in the delayed puberty. However, to date, the changes in ovarian methylomes during pubertal transition have not been investigated. In this study, using gilts as a pubertal model, the genome-wide DNA methylation were profiled to explore their dynamics during pubertal transition across Pre-, In- and Post-puberty.

Results

During pubertal transition, the follicles underwent maturation and luteinization, coupled with the significant changes in the mRNA expression of DNMT1 and DNMT3a. DNA methylation levels of In-puberty were higher than that of Pre- and Post-puberty at the locations of genes and CpG islands (CGIs). Analysis of the DNA methylation changes identified 12,313, 20,960 and 17,694 differentially methylated CpGs (DMCs) for the comparisons of Pre- vs. In-, In vs. Post-, and Pre- vs. Post-puberty, respectively. Moreover, the CGIs, upstream and exonic regions showed a significant underrepresentation of DMCs, but the CGI shores, CGI shelves, intronic, downstream and intergenic regions showed a significant overrepresentation of DMCs. Furthermore, biological functions of these methylation changes enriched in PI3K-Akt signaling pathway, GnRH signaling pathway, and Insulin secretion, and the mRNA expressions of several genes of these signaling pathway, including MMP2, ESR1, GSK3B, FGF21, IGF1R, and TAC3, were significantly changed across Pre-, In- and Post-puberty in ovaries.

Conclusions

During pubertal transition in gilts, the DNA methylation changes of ovaries were likely to affect the transcription of genes related to PI3K-Akt signaling pathway, GnRH signaling pathway, and Insulin secretion. These observations can provide new insight into the epigenetic mechanism of follicular and sexual maturation during pubertal transition in mammals.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5884-x) contains supplementary material, which is available to authorized users.

High-Resolution Bisulfite-Sequencing of Peripheral Blood DNA Methylation in Early-Onset and Familial Risk Breast Cancer Patients

PURPOSE

Understanding and explaining hereditary predisposition to cancer has focused on the genetic etiology of the disease. However, mutations in known genes associated with breast cancer, such as BRCA1 and BRCA2, account for less than 25% of familial cases of breast cancer. Recently, specific epigenetic modifications at BRCA1 have been shown to promote hereditary breast cancer, but the broader potential for epigenetic contribution to hereditary breast cancer is not yet well understood.

EXPERIMENTAL DESIGN

We examined DNA methylation through deep bisulfite sequencing of CpG islands and known promoter or regulatory regions in peripheral blood DNA from 99 patients with familial or early-onset breast or ovarian cancer, 6 unaffected BRCA mutation carriers, and 49 unaffected controls.

RESULTS

In 9% of patients, we observed altered methylation in the promoter regions of genes known to be involved in cancer, including hypermethylation at the tumor suppressor PTEN and hypomethylation at the proto-oncogene TEX14. These alterations occur in the form of allelic methylation that span up to hundreds of base pairs in length.

CONCLUSIONS

Our observations suggest a broader role for DNA methylation in early-onset, familial risk breast cancer. Further studies are warranted to clarify these mechanisms and the benefits of DNA methylation screening for early risk prediction of familial cancers.

Genome-Wide DNA Methylation Analysis of Hypothalamus During the Onset of Puberty in Gilts

Although selection of the early age at puberty in gilts will make for a favorable effect on the reproductivity of sow, a large proportion of phenotypic variation in age at puberty of gilts cannot be explained by genetics. Previous studies have implicated hypothalamic DNA methylation in the onset of puberty in mammals. However, the underlying molecular mechanism regarding the regulation of the onset of puberty has remained largely unexplored in gilts. Herein, the genome-scale DNA methylation of hypothalamus was acquired, using the reduced representation bisulfite sequencing, to compare and describe the changes of DNA methylation across Pre-, In- and Post-pubertal gilts. In this study, the average methylation levels of CpGs and CpHs (where H = C, T, or A) in CpG islands- and gene-related regions were gradually decreased in hypothalamic methylomes during the pubertal transition. Comparisons of Pre- vs. In-, In- vs. Post-, and Pre- vs. Post-pubertal stage revealed that there were 85726, 92914, and 100421 differentially methylated CpGs and 5940, 14804, and 16893 differentially methylated CpHs (where H = C, T, or A) in the hypothalamic methylomes. The methylation changes of CpHs were more dynamic than that of CpGs, and methylation changes of CpGs and CpHs were likely to be, respectively, involved in the developmental processes of reproduction and the molecular processes of cellular communications in the hypothalamus. Moreover, methylation changes of CpHs were observed to overrepresent in the quantitative trait loci of age at puberty, and the biological function of these CpH methylation changes was enriched in the pancreas development in gilts. Furthermore, the mRNA levels of several differentially CpG or CpH methylated genes related to the transcription of RNA II polymerase, GnRH signaling pathway, Estrogen signaling pathway, PI3K-AKt signaling pathway, and Insulin signaling pathway, including MAX, MMP2, FGF11, IGF1R, FGF21, and GSK3B, were significantly changed across these pubertal stages in the hypothalamus. These results will help our understanding of how DNA methylation contributes to phenotypic variation of age at puberty.

Genome-wide DNA methylation analysis of pituitaries during the initiation of puberty in gilts

It has been widely recognized that the early or delayed puberty appears to display harmful effects on adult health outcomes. During the timing of puberty, pituitaries responds to the hypothalamus and then introduce the following response of ovaries in hypothalamic-pituitary-gonadal axis. DNA methylation has been recently suggested to regulate the onset of puberty in female mammals. However, to date, the changes of DNA methylation in pituitaries have not been investigated during pubertal transition. In this study, using gilts as the pubertal model, the genome-scale DNA methylation of pituitaries was profiled and compared across Pre-, In- and Post-puberty by using the reduced representation bisulfite sequencing. We found that average methylation levels of each genomic feature in Post- were lower than Pre- and In-pubertal stage in CpG context, but they were higher in In- than that in Pre- and Post-pubertal stage in CpH (where H = A, T, or C) context. The methylation patterns of CpHs were more dynamic than that of CpGs at the location of high CpG content, low CpG content promoter genes, and differently genomic CGIs. Furthermore, the differently genomic CGIs were likely to show in a similar manner in CpG context but display in a stage-specific manner in the CpH context across the Pre-, In- and Post-pubertal stage. Among these pubertal stages, 5 kb upstream regions of the transcription start sites were protected from both CpG and CpH methylation changes. 12.65% of detected CpGs were identified as the differentially methylated CpGs, regarding 4301 genes which were involved in the fundamental functions of pituitaries. 0.35% of detected CpHs were identified as differentially methylated CpHs, regarding 3691 genes which were involved in the biological functions of releasing gonadotropin hormones. These observations and analyses would provide valuable insights into epigenetic mechanism of the initiation of puberty in pituitary level.

A role for MED14 and UVH6 in heterochromatin transcription upon destabilization of silencing

The TFIIH component UVH6 and the mediator subunit MED14 are differentially required for the release of heterochromatin silencing, and MED14 regulates non-CG DNA methylation in Arabidopsis.

Constitutive heterochromatin is associated with repressive epigenetic modifications of histones and DNA which silence transcription. Yet, particular mutations or environmental changes can destabilize heterochromatin-associated silencing without noticeable changes in repressive epigenetic marks. Factors allowing transcription in this nonpermissive chromatin context remain poorly known. Here, we show that the transcription factor IIH component UVH6 and the mediator subunit MED14 are both required for heat stress–induced transcriptional changes and release of heterochromatin transcriptional silencing in Arabidopsis thaliana. We find that MED14, but not UVH6, is required for transcription when heterochromatin silencing is destabilized in the absence of stress through mutating the MOM1 silencing factor. In this case, our results raise the possibility that transcription dependency over MED14 might require intact patterns of repressive epigenetic marks. We also uncover that MED14 regulates DNA methylation in non-CG contexts at a subset of RNA-directed DNA methylation target loci. These findings provide insight into the control of heterochromatin transcription upon silencing destabilization and identify MED14 as a regulator of DNA methylation.

T-ALL leukemia stem cell 'stemness' is epigenetically controlled by the master regulator SPI1

Leukemia stem cells (LSCs) are regarded as the origins and key therapeutic targets of leukemia, but limited knowledge is available on the key determinants of LSC 'stemness'. Using single-cell RNA-seq analysis, we identify a master regulator, SPI1, the LSC-specific expression of which determines the molecular signature and activity of LSCs in the murine Pten-null T-ALL model. Although initiated by PTEN-controlled β-catenin activation, Spi1 expression and LSC 'stemness' are maintained by a β-catenin-SPI1-HAVCR2 regulatory circuit independent of the leukemogenic driver mutation. Perturbing any component of this circuit either genetically or pharmacologically can prevent LSC formation or eliminate existing LSCs. LSCs lose their 'stemness' when Spi1 expression is silenced by DNA methylation, but Spi1 expression can be reactivated by 5-AZ treatment. Importantly, similar regulatory mechanisms may be also present in human T-ALL.

A temporal transcriptome and methylome in human embryonic stem cell-derived cardiomyocytes identifies novel regulators of early cardiac development

Stem cell-based cardiogenesis has become a powerful tool to enhance our understanding of cardiac development and test novel therapeutics for cardiovascular diseases. However, transcriptional and epigenetic regulation of multiple transitional stages from pluripotent cells to committed cardiomyocytes has not yet been fully characterized. To characterize how transcription factors, lincRNAs and DNA methylation change at temporal developmental stages, and identify potential novel regulators during cardiogenesis. We utilized a previously reported protocol that yields human cardiomyocytes (hCM) with more than 90% purity from human Embryonic Stem Cells (hESC). Leveraging the purity of cells resulting from this protocol, we systematically examined how gene expression and DNA methylation programs change at temporal developmental stages during cardiogenesis. Our results provide a comprehensive view of expression changes during cardiogenesis that extend previous studies, allowing us to identify key transcription factors as well as lincRNAs that are strongly associated with cardiac differentiation. Moreover, we incorporated a simple but powerful method to screen for novel regulators of cardiogenesis solely based on expression changes and found four novel cardiac-related transcription factors, i.e., SORBS2, MITF, DPF3, and ZNF436, which have no or few prior literature reports and we were able to validate using siRNA. Our strategy of identifying novel regulators of cardiogenesis can also be easily implemented in other stem cell-based systems. Our results provide a valuable resource for understanding cardiogenesis that extends previous findings by leveraging the purity of our cell lines, which allowed us to identify four novel cardiac-related regulators.

Allele-specific epigenome maps reveal sequence-dependent stochastic switching at regulatory loci

To assess the impact of genetic variation in regulatory loci on human health, we constructed a high-resolution map of allelic imbalances in DNA methylation, histone marks, and gene transcription in 71 epigenomes from 36 distinct cell and tissue types from 13 donors. Deep whole-genome bisulfite sequencing of 49 methylomes revealed sequence-dependent CpG methylation imbalances at thousands of heterozygous regulatory loci. Such loci are enriched for stochastic switching, which is defined as random transitions between fully methylated and unmethylated states of DNA. The methylation imbalances at thousands of loci are explainable by different relative frequencies of the methylated and unmethylated states for the two alleles. Further analyses provided a unifying model that links sequence-dependent allelic imbalances of the epigenome, stochastic switching at gene regulatory loci, and disease-associated genetic variation.

Epigenetic Changes in the Regulation of Nicotiana tabacum Response to Cucumber Mosaic Virus Infection and Symptom Recovery through Single-Base Resolution Methylomes

Plants have evolved multiple mechanisms to respond to viral infection. These responses have been studied in detail at the level of host immune response and antiviral RNA silencing (RNAi). However, the possibility of epigenetic reprogramming has not been thoroughly investigated. Here, we identified the role of DNA methylation during viral infection and performed reduced representation bisulfite sequencing (RRBS) on tissues of Cucumber mosaic virus (CMV)-infected Nicotiana tabacum at various developmental stages. Differential methylated regions are enriched with CHH sequence contexts, 80% of which are located on the gene body to regulate gene expression in a temporal style. The methylated genes depressed by methyltransferase inhibition largely overlapped with methylated genes in response to viral invasion. Activation in the argonaute protein and depression in methyl donor synthase revealed the important role of dynamic methylation changes in modulating viral clearance and resistance signaling. Methylation-expression relationships were found to be required for the immune response and cellular components are necessary for the proper defense response to infection and symptom recovery.

Receptor-targeted aptamer-siRNA conjugate-directed transcriptional regulation of HIV-1

Gene-based therapies represent a promising therapeutic paradigm for the treatment of HIV-1, as they have the potential to maintain sustained viral inhibition with reduced treatment interventions. Such an option may represent a long-term treatment alternative to highly active antiretroviral therapy.

Methods: We previously described a therapeutic approach, referred to as transcriptional gene silencing (TGS), whereby small noncoding RNAs directly inhibit the transcriptional activity of HIV-1 by targeting sites within the viral promoter, specifically the 5' long terminal repeat (LTR). TGS differs from traditional RNA interference (RNAi) in that it is characterized by concomitant silent-state epigenetic marks on histones and DNA. To deliver TGS-inducing RNAs, we developed functional RNA conjugates based on the previously reported dual function of the gp120 (A-1) aptamer conjugated to 27-mer Dicer-substrate anti-HIV-1 siRNA (dsiRNA), LTR-362.

Results: We demonstrate here that high levels of processed guide RNAs localize to the nucleus in infected T lymphoblastoid CEM cell line and primary human CD4+ T-cells. Treatment of the aptamer-siRNA conjugates induced TGS with an ~10-fold suppression of viral p24 levels as measured at day 12 post infection. To explore the silencing efficacy of aptamer-siRNA conjugates in vivo, HIV-1-infected humanized NOD/SCID/IL2 rγ^null mice (hu-NSG) were treated with the aptamer-siRNA conjugates. Systemic delivery of the A-1-stick-LTR-362 27-mer siRNA conjugates suppressed HIV-1 infection and protected CD4+ T cell levels in viremia hu-NSG mice.

Principle conclusions: Collectively these data suggest that the gp120 aptamer-dsiRNA conjugate design is suitable for systemic delivery of small RNAs that can be used to suppress HIV-1.