5-Hydroxymethylcytosine (5hmC) is an important mammalian DNA epigenetic modification that has been linked to gene regulation and cancer pathogenesis. Here we explored the diagnostic potential of 5hmC in circulating cell-free DNA (cfDNA) using a sensitive chemical labeling-based low-input shotgun sequencing approach. We sequenced cell-free 5hmC from 49 patients of seven different cancer types and found distinct features that could be used to predict cancer types and stages with high accuracy. Specifically, we discovered that lung cancer leads to a progressive global loss of 5hmC in cfDNA, whereas hepatocellular carcinoma and pancreatic cancer lead to disease-specific changes in the cell-free hydroxymethylome. Our proof-of-principle results suggest that cell-free 5hmC signatures may potentially be used not only to identify cancer types but also to track tumor stage in some cancers.
Histone modifications and DNA methylation are two major epigenetic factors regulating gene expression. However, the mechanism in which DNA methylation and histone modifications co-regulate gene expression was little studied. In our study, classifications of DNA methylation and gene expression showed the complicated relationship between gene expression and epigenetic factors. A Bayesian network was constructed by using the high-resolution maps of histone modifications, DNA methylation and gene expression in human CD4+ T cells to deduce causal and combinatorial relationships among them. PolII was found as the only direct regulator to gene expression, which was not found in prior studies. Our Bayesian network showed that epigenetic factors such as H3K4me3, H3K27me3 and DNA methylation are key regulators of gene expression, though indirectly. However they were considered to combinatorially stabilize the state and structure of chromatin.
Gene expression regulation by synergistic epigenetic modifications and genetic elements has gain substantial attention for ten years. However, epigenetic and genetic alterations have also been considered to participate in carcinogenesis. The recent development of next-generation sequencing to explore various human cancers has uncovered many mutations in genes that control the epigenome and many epigenetic patterns that control the genome. The genetic alterations tend to change the epigenetic patterns including DNA methylation, histone modifications and ultimately to change gene expression. The epigenetic alterations contribute to cancer by inducing genetic mutations and impairing DNA repair functions. The development of next-generation high-throughput technologies makes the discovery of tumor-suppressor genes and oncogenes more easily. With the aid of Bioinformatics, a large number of studies have suggested that cancer is a multilayer and complex disease involving both genetic mutations and epigenetic changes. We illustrate the recent progress of cancer pathogenesis and explain why cancer genome and epigenome complement each other in cancer progression. The crosstalk between genome and epigenome is a new viewpoint that may be useful in cancer therapy.
Wiedemann-Rautenstrauch syndrome (WDRTS) is an extremely rare autosomal recessive neonatal disorder. Currently, over 50 cases with variable phenotypes of WDRTS have been reported. In our cohort of prenatal and postnatal growth retardation, a female proband was found to have general growth retardation, neurocutaneous syndrome, and anemia. Karyotype test and array-CGH detected no obvious chromosomal aberrations. Trio-based whole-exome sequencing (Trio-WES) identified bi-allelic compound mutations in the coding sequence (CDS) of POLR3A gene (c.3342C > T, p.Ser1114 = and c.3718G > A, p.Gly1240Ser). For the mild anemia phenotype, the underlying causal genetic factors could be attributed to the compound heterozygous mutations in FANCA gene (c.2832dup, p.Ala945CysfsTer6 and c.1902 T > G, p.Asp634Glu). Mini-gene reporter assays revealed that the synonymous variant of POLR3A and the missense variant of FANCA could affect pre-mRNA splicing of each gene. For POLR3A, the synonymous mutation (c.3342C > T, p.Ser1114=) generated three types of aberrant isoforms. Therefore, the female patient was finally diagnosed as WDRTS caused by POLR3A. For FANCA, the missense variant (c.1902 T > G, p.Asp634Glu) disrupted the normal splicing between exon 21 and 22, and produced two types of abnormal isoforms, one carrying the 1902G and the other spliced between exon 21 and 23 to exclude exon 22. Network analysis showed that POLR3A and FANCA could be STRINGed, indicating both proteins might collaborate for some unknown functions. Current investigation would broaden the knowledge for clinicians and genetic counselors and remind them to interpret those synonymous or predicted "benign" variants more carefully.
To detect potential mutation of the ASPM gene in a Chinese pedigree affected with autosomal recessive primary microcephaly 5 (MCPH5).Peripheral venous blood samples were collected from the proband and her parents. Amniotic fluid sample was also collected upon her mother' s subsequent pregnancy. Following extraction of genomic DNA, PCR and Sanger sequencing were carried out to identify potential variants of the ASPM gene.The proband was found to harbor compound heterozygous variants of the ASPM gene, namely c.8214dupT (p.Q2739fs) in exon 18 and c.9541C>T (p.R3181X) in exon 23, which were respectively inherited from her father and mother. The fetus has found to have inherited the c.9541C>T (p.R3181X) variant only.The c.8214dupT (p.Q2739fs) and c.9541C>T (p.R3181X) compound heterozygous variants of the ASPM gene probably underlay the pathogenesis of MCPH5 in this patient. Above finding has enabled genetic counseling and prenatal diagnosis for her family.
Abstract 5-Hydroxymethylcytosine (5hmC) is an important mammalian DNA epigenetic modification that has been linked to gene regulation and cancer pathogenesis. Here we explored the diagnostic potential of 5hmC in circulating cell-free DNA (cfDNA) using a sensitive chemical labeling-based low-input shotgun sequencing approach. We sequenced cell-free 5hmC from 49 patients of seven different cancer types and found distinct features that could be used to predict cancer types and stages with high accuracy. Specifically, we discovered that lung cancer leads to a progressive global loss of 5hmC in cfDNA, whereas hepatocellular carcinoma and pancreatic cancer lead to disease-specific changes in the cell-free hydroxymethylome. Our proof of principle results suggest that cell-free 5hmC signatures may potentially be used not only to identify cancer types but also to track tumor stage in some cancers. One Sentence Summary Analyzing the epigenetic modification 5-hydroxymethylcysoine in circulating cell-free DNA reveals tumor tissue of origin and stages for cancer diagnostics.
Abstract Exploring the epigenetic regulation mechanism of colorectal cancer (CRC) from the perspective of N6-methyladenosine (m6A) modification may provide a new target for tumor therapy. Analysis using high-throughput RNA-seq profile from TCGA found that the gene expression of Methyltransferase-like 3 (METTL3) was significantly upregulated among 20 m6A binding proteins in CRC, which was also validated in CRC cancer tissues and cell lines. Moreover, transcriptome sequencing in METTL3 knockdown cells using CRISPR/Cas9 editing suggested that EphA2 and VEGFA were differential expression, which were enriched in the vasculature development, PI3K/AKT and ERK1/2 signal pathway through the functional enrichment analysis. The results in vitro revealed that METTL3 as the m6A “writers” participates the methylation of EphA2 and VEGFA, which were recognized by the m6A “readers”, insulin-like growth factor 2 mRNA binding protein 2/3 (IGF2BP2/3), to prevent their mRNA degradation. In addition, EphA2 and VEGFA targeted by METTL3 via different IGF2BP-dependent mechanisms were found to promote vasculogenic mimicry (VM) formation via PI3K/AKT/mTOR and ERK1/2 signaling in CRC. The study suggests that intervention with m6A-binding proteins (METTL3 and IGF2BP2/3) may provide a potential diagnostic or prognostic target of VM-based anti-metastasis drugs for CRC.
Methylation quantitative trait loci (mQTLs) are the genetic variants that may affect the DNA methylation patterns of CpG sites. However, their roles in influencing the disturbances of smoking-related epigenetic changes have not been well established. This study was conducted to address whether mQTLs exist in the vicinity of smoking-related CpG sites (± 50 kb) and to examine their associations with smoking exposure and all-cause mortality in older adults. We obtained DNA methylation profiles in whole blood samples by Illumina Infinium Human Methylation 450 BeadChip array of two independent subsamples of the ESTHER study (discovery set, n = 581; validation set, n = 368) and their corresponding genotyping data using the Illumina Infinium OncoArray BeadChip. After correction for multiple testing (FDR), we successfully identified that 70 out of 151 previously reported smoking-related CpG sites were significantly associated with 192 SNPs within the 50 kb search window of each locus. The 192 mQTLs significantly influenced the active smoking-related DNA methylation changes, with percentage changes ranging from 0.01 to 18.96%, especially for the weakly/moderately smoking-related CpG sites. However, these identified mQTLs were not directly associated with active smoking exposure or all-cause mortality. Our findings clearly demonstrated that if not dealt with properly, the mQTLs might impair the power of epigenetic-based models of smoking exposure to a certain extent. In addition, such genetic variants could be the key factor to distinguish between the heritable and smoking-induced impact on epigenome disparities. These mQTLs are of special importance when DNA methylation markers measured by Illumina Infinium assay are used for any comparative population studies related to smoking-related cancers and chronic diseases.
Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder. Aberrant DNA methylation has been observed in ASD but the mechanisms remain largely unknown. Here, we employed discordant monozygotic twins to investigate the contribution of DNA methylation to ASD etiology. Genome-wide DNA methylation analysis was performed using samples obtained from five pairs of ASD-discordant monozygotic twins, which revealed a total of 2,397 differentially methylated genes. Further, such gene list was annotated with Kyoto Encyclopedia of Genes and Genomes and demonstrated predominant activation of neurotrophin signaling pathway in ASD-discordant monozygotic twins. The methylation of SH2B1 gene was further confirmed in the ASD-discordant, ASD-concordant monozygotic twins, and a set of 30 pairs of sporadic case-control by bisulfite-pyrosequencing. The results showed that there was a greater DNA methylation difference in ASD-discordant monozygotic twins than ASD-concordant monozygotic twins. Further, verification of the Chr.16:28856743 of SH2B1 showed significant differences in DNA methylation between case and control. These results suggest abnormal methylation of SH2B1 is associated with ASD etiology. Our data suggest that it might be worthwhile to further explore the functions of SH2B1 and related genes of neurotrophin signaling pathway in ASD.
PTC is not generally considered a lethal disease, but prone to recurrence as the prognosis. Hashimoto's thyroiditis (HT) is an important factor that affects the prognosis of papillary thyroid carcinoma (PTC). It is crucial to find biomarkers to identify the combination of HT with PTC and to predict the prognosis.RNASeq data from the Cancer Genome Atlas (TCGA) database was used to screen differentially expressed genes (DEGs) of PTC with HT via the edgeR package of R software version 3.3.0. Also, the DEGs were applied to the DAVID web-based tool to determine the enrichment of gene functions via Gene Ontology (GO) analysis and to identify associated pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. By constructing protein interaction networks within Cytoscape software, we screened candidate genes and explored possible relationships with the clinical phenotype of PTC. Finally, additional thyroid tissue samples were collected to verify the results above.After analyzing the RNA-Seq data of PTC patients from the Cancer Genomic Atlas, 497 differentially expressed PTC genes were found to be associated with HT, of which protein tyrosine phosphatase receptor type C (PTPRC), KIT, and COL1A1 were associated with tumor size and lymph node metastasis (p < 0.05). Verification of these results with another 30 thyroid tissues of clinical PTC patients revealed that the expression level of PTPRC in the PTC with HT group was higher than that in the PTC without HT group (p < 0.05) and the ROC curve showed a good discrimination (area under the curve = 0.846). However, the correlation with the clinical phenotype was not statistically significant (p > 0.05).These data suggest that upregulation of PTPRC enhances the incidence of HT associated with PTC and is also predictive of a poor prognosis.