Allele-specific DNA methylation of disease susceptibility genes in Japanese patients with inflammatory bowel disease

Haplotype-resolved genome assembly of Ficus carica L. reveals allele-specific expression in the fruit

In this study, we produced a haplotype-phased genome sequence of the fig tree (Ficus carica L.), a non-Rosaceae fruit tree model species, providing a systematic overview of the organization of a heterozygous diploid genome and, for the first time in a fruit tree, evidence of allelic expression direction-shifting among haplotypes. The genome was used for whole genome analysis of heterozygosis, allelic cytosine methylation, and expression profiles in peel and pulp fruit tissues. The two pseudo-haplotypes spanned approximately 355 and 346 Mbp, respectively, and 97% of the sequences were associated with 13 chromosome pairs of the fig tree. Overall, the methylation profile in peel and pulp tissues showed no variations between the homologous chromosomes. However, we detected differential DNA methylation within defined heterozygous allelic gene regions, particularly in upstream regions. Among 6768 heterozygous coding sequences identified, 4024 exhibited allele-specific expression, with approximately 18% specific to the peel and 14% to the pulp. Specifically, 2715 genes were consistent, with one allele always more expressed than the other in both peel and pulp. On the contrary, 22 allele-specific expressed genes switched allele expression among the fig fruit peel and pulp tissues, indicating evidence of overdominance and suggesting that the genome can express one of the two alleles higher or lower depending on developmental or environmental triggers. Notably, these genes were associated with key biological processes, including fruit maturation regulation, seed maturation, and stress response, highlighting their potential role in the plant's developmental and adaptive functions in view of gene editing-based breeding.

Perfluorooctanesulfonic acid (PFOS) induced cancer related DNA methylation alterations in human breast cells: A whole genome methylome study

DNA methylation plays a pivotal role in cancer. The ubiquitous contaminant perfluorooctanesulfonic acid (PFOS) has been epidemiologically associated with breast cancer, and can induce proliferation and malignant transformation of normal human breast epithelial cells (MCF-10A), but the information about its effect on DNA methylation is sparse. The aim of this study was to characterize the whole-genome methylome effects of PFOS in our breast cell model and compare the findings with previously demonstrated DNA methylation alterations in breast tumor tissues. The DNA methylation profile was assessed at single CpG resolution in MCF-10A cells treated with 1 μM PFOS for 72 h by using Enzymatic Methyl sequencing (EM-seq). We found 12,591 differentially methylated CpG-sites and 13,360 differentially methylated 100 bp tiles in the PFOS exposed breast cells. These differentially methylated regions (DMRs) overlapped with 2406 genes of which 494 were long non-coding RNA and 1841 protein coding genes. We identified 339 affected genes that have been shown to display altered DNA methylation in breast cancer tissue and several other genes related to cancer development. This includes hypermethylation of GACAT3, DELEC1, CASC2, LCIIAR, MUC16, SYNE1 and hypomethylation of TTN and KMT2C. DMRs were also found in estrogen receptor genes (ESR1, ESR2, ESRRG, ESRRB, GREB1) and estrogen responsive genes (GPER1, EEIG1, RERG). The gene ontology analysis revealed pathways related to cancer phenotypes such as cell adhesion and growth. These findings improve the understanding of PFOS's potential role in breast cancer and illustrate the value of whole-genome methylome analysis in uncovering mechanisms of chemical effects, identifying biomarker candidates, and strengthening epidemiological associations, potentially impacting risk assessment.

SMAD3 rs36221701 T>C polymorphism impacts COPD susceptibility in the Kashi population

Small mother against decapentaplegic (SMAD) family member 3 (SMAD3) is well correlated with the inflammatory response of chronic obstructive pulmonary disease (COPD). A previous study indicated that the single nucleotide polymorphism (SNP) rs36221701 of SMAD3 was related to the risk of inflammatory disease. Hence, given the pathogenesis of COPD is intently associated with smoking and gene polymorphism, this study aims to analyze the relationship between SMAD3 rs36221701 and COPD susceptibility, and to explore whether the interaction is related to smoking status. We studied the association between the rs36221701 and rs34307601 of SMAD3 and COPD susceptibility, a total of 541 COPD patients and 534 controls of the Uyghur population were recruited at the First People's Hospital and the village of Kashi. The interrelation of the two SNPs with the risk of COPD was determined by calculating odds ratio (OR) and 95% confidence interval (95% CI). We found a significant association between the rs36221701 and COPD risk in the non-smoking population. TC genotype showed a significant decreased association with COPD risk (OR = 0.59, 95% CI = 0.41-0.83, P < 0.05), but CC genotype can increased the COPD risk (OR > 1, P < 0.05). In addition, COPD susceptibility was not related to the genetic variations in the rs34307601 (P > 0.05). In conclusion, we confirmed that the SMAD3 rs36221701 may be associated with COPD susceptibility in the Chinese Uyghur population, especially among non-smokers. Our data provide new light for the relationship between SMAD3 gene polymorphisms and COPD susceptibility in the Chinese Uyghur population.

Perspectives on Allele-Specific Expression

Diploidy has profound implications for population genetics and susceptibility to genetic diseases. Although two copies are present for most genes in the human genome, they are not necessarily both active or active at the same level in a given individual. Genomic imprinting, resulting in exclusive or biased expression in favor of the allele of paternal or maternal origin, is now believed to affect hundreds of human genes. A far greater number of genes display unequal expression of gene copies due to cis-acting genetic variants that perturb gene expression. The availability of data generated by RNA sequencing applied to large numbers of individuals and tissue types has generated unprecedented opportunities to assess the contribution of genetic variation to allelic imbalance in gene expression. Here we review the insights gained through the analysis of these data about the extent of the genetic contribution to allelic expression imbalance, the tools and statistical models for gene expression imbalance, and what the results obtained reveal about the contribution of genetic variants that alter gene expression to complex human diseases and phenotypes.

MethHaplo: combining allele-specific DNA methylation and SNPs for haplotype region identification

Background

DNA methylation is an important epigenetic modification that plays a critical role in most eukaryotic organisms. Parental alleles in haploid genomes may exhibit different methylation patterns, which can lead to different phenotypes and even different therapeutic and drug responses to diseases. However, to our knowledge, no software is available for the identification of DNA methylation haplotype regions with combined allele-specific DNA methylation, single nucleotide polymorphisms (SNPs) and high-throughput chromosome conformation capture (Hi-C) data.

Results

In this paper, we developed a new method, MethHaplo, that identify DNA methylation haplotype regions with allele-specific DNA methylation and SNPs from whole-genome bisulfite sequencing (WGBS) data. Our results showed that methylation haplotype regions were ten times longer than haplotypes with SNPs only. When we integrate WGBS and Hi-C data, MethHaplo could call even longer haplotypes.

Conclusions

This study illustrates the usefulness of methylation haplotypes. By constructing methylation haplotypes for various cell lines, we provide a clearer picture of the effect of DNA methylation on gene expression, histone modification and three-dimensional chromosome structure at the haplotype level. Our method could benefit the study of parental inheritance-related disease and hybrid vigor in agriculture.

SMAD3 Hypomethylation as a Biomarker for Early Prediction of Colorectal Cancer

The incidence and mortality rates of colorectal cancer (CRC) have been high in recent years. Prevention and early detection are crucial for decreasing the death rate. Therefore, this study aims to characterize the alteration patterns of mothers against decapentaplegic homolog 3 (SMAD3) in patients with CRC and its applications in early detection by using a genome-wide methylation array to identify an aberrant hypomethylation site in the intron position of the SMAD3 gene. Quantitative methylation-specific polymerase chain reaction showed that hypomethylated SMAD3 occurred in 91.4% (501/548) of Taiwanese CRC tissues and 66.6% of benign tubular adenoma polyps. In addition, SMAD3 hypomethylation was observed in 94.7% of patients with CRC from The Cancer Genome Atlas dataset. A decrease in circulating cell-free methylation SMAD3 was detected in 70% of CRC patients but in only 20% of healthy individuals. SMAD3 mRNA expression was low in 42.9% of Taiwanese CRC tumor tissues but high in 29.4% of tumors compared with paired adjacent normal tissues. Hypomethylated SMAD3 was found in cancers of the digestive system, such as liver cancer, gastric cancer, and colorectal cancer, but not in breast cancer, endometrial cancer, and lung cancer. In conclusion, SMAD3 hypomethylation is a potential diagnostic marker for CRC in Western and Asian populations.

The DNA methylome of inflammatory bowel disease (IBD) reflects intrinsic and extrinsic factors in intestinal mucosal cells

Abnormal DNA methylation has been described in human inflammatory conditions of the gastrointestinal tract, such as inflammatory bowel disease (IBD). As other complex diseases, IBD results from the balance between genetic predisposition and environmental exposures. As such, DNA methylation may be the consequence (and potential effector) of both, genetic susceptibility variants and/or environmental signals such as cytokine exposure. We attempted to discern between these two non-excluding possibilities by performing a combined analysis of published DNA methylation data in intestinal mucosal cells of IBD and control samples. We identified abnormal DNA methylation at different levels: deviation from mean methylation signals at site and region levels, and differential variability. A fraction of such changes is associated with genetic polymorphisms linked to IBD susceptibility. In addition, by comparing with another intestinal inflammatory condition (i.e., coeliac disease) we propose that aberrant DNA methylation can also be the result of unspecific processes such as chronic inflammation. Our characterization suggests that IBD methylomes combine intrinsic and extrinsic responses in intestinal mucosal cells, and could point to knowledge-based biomarkers of IBD detection and progression.

Non-coding DNA in IBD: from sequence variation in DNA regulatory elements to novel therapeutic potential

Genome-wide association studies have identified over 200 loci associated with IBD. We and others have recently shown that, in addition to variants in protein-coding genes, the majority of the associated loci are related to DNA regulatory elements (DREs). These findings add a dimension to the already complex genetic background of IBD. In this review we summarise the existing evidence on the role of DREs in IBD. We discuss how epigenetic research can be used in candidate gene approaches that take non-coding variants into account and can help to pinpoint the essential pathways and cell types in the pathogenesis of IBD. Despite the increased level of genetic complexity, these findings can contribute to novel therapeutic options that target transcription factor binding and enhancer activity. Finally, we summarise the future directions and challenges of this emerging field.

Correction: Allele-specific DNA methylation of disease susceptibility genes in Japanese patients with inflammatory bowel disease

[This corrects the article DOI: 10.1371/journal.pone.0194036.].

Crosstalk of Genetic Variants, Allele-Specific DNA Methylation, and Environmental Factors for Complex Disease Risk

Over the past decades, genome-wide association studies (GWAS) have identified thousands of phenotype-associated DNA sequence variants for potential explanations of inter-individual phenotypic differences and disease susceptibility. However, it remains a challenge for translating the associations into causative mechanisms for complex diseases, partially due to the involved variants in the noncoding regions and the inconvenience of functional studies in human population samples. So far, accumulating evidence has suggested a complex crosstalk among genetic variants, allele-specific binding of transcription factors (ABTF), and allele-specific DNA methylation patterns (ASM), as well as environmental factors for disease risk. This review aims to summarize the current studies regarding the interactions of the aforementioned factors with a focus on epigenetic insights. We present two scenarios of single nucleotide polymorphisms (SNPs) in coding regions and non-coding regions for disease risk, via potentially impacting epigenetic patterns. While a SNP in a coding region may confer disease risk via altering protein functions, a SNP in non-coding region may cause diseases, via SNP-altering ABTF, ASM, and allele-specific gene expression (ASE). The allelic increases or decreases of gene expression are key for disease risk during development. Such ASE can be achieved via either a "SNP-introduced ABTF to ASM" or a "SNP-introduced ASM to ABTF." Together with our additional in-depth review on insulator CTCF, we are convinced to propose a working model that the small effect of a SNP acts through altered ABTF and/or ASM, for ASE and eventual disease outcome (named as a "SNP intensifier" model). In summary, the significance of complex crosstalk among genetic factors, epigenetic patterns, and environmental factors requires further investigations for disease susceptibility.

Methylenetetrahydrofolate Dehydrogenase 1 Silencing Expedites the Apoptosis of Non-Small Cell Lung Cancer Cells via Modulating DNA Methylation

BACKGROUND Non-small cell lung cancer (NSCLC) accounts for about 85% of all types of lung cancer. Methylenetetrahydrofolate dehydrogenase 1 (MTHFD1) is involved in DNA methylation, and DNA methylation is related to tumorigenesis. The role of MTHFD1 in NSCLC was examined in our study. MATERIAL AND METHODS The correlation between the expression of MTHFD1 and the clinicopathological features of patients diagnosed with lung cancer was investigated using the chi-square test. The viability and apoptosis of NCI-H1299 cells was respectively detected using cell counting kit-8 and flow cytometry assays. The expression levels of MTHFD1, apoptosis-related factors and DNA methyltransferase-related factors were assessed by quantitative real-time PCR (qRT-PCR) and western blot assays. RESULTS We found that MTHFD1 expression in the tumor tissues and cells was higher than that of adjacent normal tissues and cells. The survival time of patients with high MTHFD1 expression was shorter than those with low MTHFD1 expression. The expression level of MTHFD1 was related to tumor size, TNM stage, histologic grade, and metastasis, but not linked to gender and age. Besides, si-MTHFD1 significantly decreased the viability of cells in a time-dependent manner, and increased cell apoptosis. When cells were transfected with MTHFD1-siRNA, the levels of surviving and B-cell lymphoma-2 (Bcl-2) were attenuated, while p53 and Bcl-2 associated X protein (Bax) levels were enhanced. Moreover, si-MTHFD1 markedly downregulated the expression levels of DNA methyltransferase 1 (DNMT1), DNMT3a, and DNMT3b. CONCLUSIONS Collectively, our results proved that MTHFD1 silencing obviously reduced the proliferation and enhanced the apoptosis of NSCLC via suppressing DNA methylation.