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

Prediction of deleterious non-synonymous SNPs of human MDC1 gene: an in silico approach

MDC1 (Mediator of DNA damage Checkpoint protein 1) functions to facilitate the localization of numerous DNA damage response (DDR) components to DNA double-strand break sites. MDC1 is an integral component in preserving genomic stability and appropriate DDR regulation. There haven't been systematic investigations of MDC1 mutations that induce cancer and genomic instability. Variations in nsSNPs have the potential to modify the protein chemistry and their function. Describing functional SNPs in disease-associated genes presents a significant conundrum for investigators, it is possible to assess potential functional SNPs before conducting larger population examinations. Multiple sequences and structure-based bioinformatics strategies were implemented in the current in-silico investigation to discern potential nsSNPs of the MDC1 genes. The nsSNPs were identified with SIFT, SNAP2, Align GVGD, PolyPhen-2, and PANTHER, and their stability was determined with MUpro. The conservation, solvent accessibility, and structural effects of the mutations were identified with ConSurf, NetSurfP-2.0, and SAAFEC-SEQ respectively. Cancer-related analysis of the nsSNPs was conducted using cBioPortal and TCGA web servers. The present study appraised five nsSNPs (P1426T, P69S, P194R, P203L, and H131Y) as probably mutilating due to their existence in highly conserved regions and propensity to deplete protein stability. The nsSNPs P194R, P203L, and H131Y were concluded as deleterious and possibly damaging from the 5 prediction tools. The functional nsSNP P194R mutation is associated with skin cutaneous melanoma while no significant records were found for other nsSNPs. The present study concludes that the highly deleterious P194R mutations can potentially induce genomic instability and contribute to various cancers' pathogenesis. Developing drugs targeting these mutations can undoubtedly be advantageous in large population-based studies, particularly in the development of precision medicine.

AIMER: A SNP-independent software for identifying imprinting-like allelic methylated regions from DNA methylome

Genomic imprinting is essential for mammalian growth and embryogenesis. High-throughput bisulfite sequencing accompanied with parental haplotype-specific information allows analysis of imprinted genes and imprinting control regions (ICRs) on a large scale. Currently, although several allelic methylated regions (AMRs) detection software were developed, methods for detecting imprinted AMRs is still limited. Here, we developed a SNP-independent statistical approach, AIMER, to detect imprinting-like AMRs. By using the mouse frontal cortex methylome as input, we demonstrated that AIMER performs very well in detecting known germline ICRs compared with other methods. Furthermore, we found the putative parental AMRs AIMER detected could be distinguished from sequence-dependent AMRs. Finally, we found a novel germline imprinting-like AMR using WGBS data from 17 distinct mouse tissue samples. The results indicate that AIMER is a good choice for detecting imprinting-like (parent-of-origin-dependent) AMRs. We hope this method will be helpful for future genomic imprinting studies. The Python source code for our project is now publicly available on both GitHub (https://github.com/ZhaoLab-TMU/AIMER) and Gitee (https://gitee.com/zhaolab_tmu/AIMER).

Approaches based on miRNAs in Behçet's Disease: Unveiling pathogenic mechanisms, diagnostic strategies, and therapeutic applications

Behçet's Disease (BD) is an intricate medical puzzle, captivating researchers with its enigmatic pathogenesis. This complex ailment, distinguished by recurrent mouth and genital lesions, eye irritation, and skin injuries, presents a substantial obstacle to therapeutic research. This review explores the complex interaction of microRNAs (miRNAs) with BD, highlighting their crucial involvement in the disease's pathophysiology. miRNAs, recognized for regulatory influence in diverse biological processes, hold a pivotal position in the molecular mechanisms of autoimmune diseases, such as BD. The exploration begins with examining miRNA biogenic pathways and functions, establishing a foundational understanding of their regulatory mechanisms. Shifting to the molecular landscape governing BD, the review highlights miRNA-mediated impacts on critical signaling pathways like Notch, nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB), and protein kinase B (AKT)/mammalian target of rapamycin (mTOR), offering insights into intricate pathophysiological mechanisms. Dissecting the immunological landscape reveals the profound influence of miRNAs on BD, shedding light on the intricate modulation of immune responses and offering novel perspectives on disease etiology and progression. Beyond molecular intricacies, the review explores the clinical relevance of miRNAs in BD, emphasizing their potential as diagnostic and prognostic indicators. The discussion extends to the promising realm of miRNA-based therapeutic interventions, highlighting their potential in alleviating symptoms and altering disease progression. This comprehensive review, serving as a valuable resource for researchers, clinicians, and stakeholders, aims to decipher the intricate molecular tapestry of BD and explore the therapeutic potential of miRNAs.

Anemia and Its Connections to Inflammation in Older Adults: A Review

Anemia is a common hematological disorder that affects 12% of the community-dwelling population, 40% of hospitalized patients, and 47% of nursing home residents. Our understanding of the impact of inflammation on iron metabolism and erythropoiesis is still lacking. In older adults, anemia can be divided into nutritional deficiency anemia, bleeding anemia, and unexplained anemia. The last type of anemia might be caused by reduced erythropoietin (EPO) activity, progressive EPO resistance of bone marrow erythroid progenitors, and the chronic subclinical pro-inflammatory state. Overall, one-third of older patients with anemia demonstrate a nutritional deficiency, one-third have a chronic subclinical pro-inflammatory state and chronic kidney disease, and one-third suffer from anemia of unknown etiology. Understanding anemia's pathophysiology in people aged 65 and over is crucial because it contributes to frailty, falls, cognitive decline, decreased functional ability, and higher mortality risk. Inflammation produces adverse effects on the cells of the hematological system. These effects include iron deficiency (hypoferremia), reduced EPO production, and the elevated phagocytosis of erythrocytes by hepatic and splenic macrophages. Additionally, inflammation causes enhanced eryptosis due to oxidative stress in the circulation. Identifying mechanisms behind age-related inflammation is essential for a better understanding and preventing anemia in older adults.

Polymorphism in epigenetic regulating genes in relation to periodontitis, number of teeth, and levels of high-sensitivity C-reactive protein and glycated hemoglobin: The Tromsø Study 2015-2016

BACKGROUND

The aim of this study was to investigate the association between periodontitis and four single nucleotide polymorphisms (SNPs) in genes involved in epigenetic regulation of DNA, and between these same SNPs and tooth loss, high-sensitivity C-reactive protein (hs-CRP), and glycated hemoglobin (HbA1c) levels.

METHODS

We included participants with periodontal examination (n = 3633, aged: 40-93 years) from the Tromsø Study seventh survey (2015-2016), Norway. Periodontitis was defined according to the 2017 AAP/EFP classification system as no periodontitis, grades A, B, or C. Salivary DNA was extracted and genotyping was performed to investigate four SNPs (rs2288349, rs35474715, rs34023346, and rs10010325) in the sequence of the genes DNMT1, IDH2, TET1, and TET2. Association between SNPs and periodontitis was analyzed by logistic regression adjusted for age, sex, and smoking. Subgroup analyses on participants aged 40-49 years were performed.

RESULTS

In participants aged 40-49 years, homozygous carriage of minor A-allele of rs2288349 (DNMT1) was associated with decreased susceptibility to periodontitis (grade A: odds ratio [OR] 0.55; p = 0.014: grade B/C OR 0.48; p = 0.004). The minor A-allele of rs10010325 (TET2) was associated with increased susceptibility to periodontitis (grade A OR 1.69; p = 0.035: grade B/C OR 1.90; p = 0.014). In the entire sample, homozygous carriage of the G-allele of rs35474715 (IDH2) was associated with having ≤24 teeth (OR 1.31; p = 0.018). Homozygous carriage of the A-allele of TET2 was associated with hs-CRP≥3 mg/L (OR 1.37; p = 0.025) and HbA1c≥6.5% (OR 1.62; p = 0.028).

CONCLUSIONS

In this Norwegian population, there were associations between polymorphism in genes related to DNA methylation and periodontitis, tooth loss, low-grade inflammation, and hyperglycemia.

Integrated analysis of human DNA methylation, gene expression, and genomic variation in iMETHYL database using kernel tensor decomposition-based unsupervised feature extraction

Integrating gene expression, DNA methylation, and genomic variants simultaneously without location coincidence (i.e., irrespective of distance from each other) or pairwise coincidence (i.e., direct identification of triplets of gene expression, DNA methylation, and genomic variants, and not integration of pairwise coincidences) is difficult. In this study, we integrated gene expression, DNA methylation, and genome variants from the iMETHYL database using the recently proposed kernel tensor decomposition-based unsupervised feature extraction method with limited computational resources (i.e., short CPU time and small memory requirements). Our methods do not require prior knowledge of the subjects because they are fully unsupervised in that unsupervised tensor decomposition is used. The selected genes and genomic variants were significantly targeted by transcription factors that were biologically enriched in KEGG pathway terms as well as in the intra-related regulatory network. The proposed method is promising for integrated analyses of gene expression, methylation, and genomic variants with limited computational resources.

Impact of Plasmodium falciparum infection on DNA methylation of circulating immune cells

The regulation of immune cell responses to infection is a complex process that involves various molecular mechanisms, including epigenetic regulation. DNA methylation has been shown to play central roles in regulating gene expression and modulating cell response during infection. However, the nature and extent to which DNA methylation is involved in the host immune response in human malaria remains largely unknown. Here, we present a longitudinal study investigating the temporal dynamics of genome-wide in vivo DNA methylation profiles using 189 MethylationEPIC 850 K profiles from 66 children in Burkina Faso, West Africa, sampled three times: before infection, during symptomatic parasitemia, and after malaria treatment. The results revealed major changes in the DNA methylation profiles of children in response to both Plasmodium falciparum infection and malaria treatment, with widespread hypomethylation of CpGs upon infection (82% of 6.8 K differentially methylated regions). We document a remarkable reversal of CpG methylation profiles upon treatment to pre-infection states. These changes implicate divergence in core immune processes, including the regulation of lymphocyte, neutrophil, and myeloid leukocyte function. Integrative DNA methylation-mRNA analysis of a top differentially methylated region overlapping the pro-inflammatory gene TNF implicates DNA methylation of TNF cis regulatory elements in the molecular mechanisms of TNF regulation in human malaria. Our results highlight a central role of epigenetic regulation in mounting the host immune response to P. falciparum infection and in response to malaria treatment.

Widespread allele-specific topological domains in the human genome are not confined to imprinted gene clusters

BACKGROUND

There is widespread interest in the three-dimensional chromatin conformation of the genome and its impact on gene expression. However, these studies frequently do not consider parent-of-origin differences, such as genomic imprinting, which result in monoallelic expression. In addition, genome-wide allele-specific chromatin conformation associations have not been extensively explored. There are few accessible bioinformatic workflows for investigating allelic conformation differences and these require pre-phased haplotypes which are not widely available.

RESULTS

We developed a bioinformatic pipeline, "HiCFlow," that performs haplotype assembly and visualization of parental chromatin architecture. We benchmarked the pipeline using prototype haplotype phased Hi-C data from GM12878 cells at three disease-associated imprinted gene clusters. Using Region Capture Hi-C and Hi-C data from human cell lines (1-7HB2, IMR-90, and H1-hESCs), we can robustly identify the known stable allele-specific interactions at the IGF2-H19 locus. Other imprinted loci (DLK1 and SNRPN) are more variable and there is no "canonical imprinted 3D structure," but we could detect allele-specific differences in A/B compartmentalization. Genome-wide, when topologically associating domains (TADs) are unbiasedly ranked according to their allele-specific contact frequencies, a set of allele-specific TADs could be defined. These occur in genomic regions of high sequence variation. In addition to imprinted genes, allele-specific TADs are also enriched for allele-specific expressed genes. We find loci that have not previously been identified as allele-specific expressed genes such as the bitter taste receptors (TAS2Rs).

CONCLUSIONS

This study highlights the widespread differences in chromatin conformation between heterozygous loci and provides a new framework for understanding allele-specific expressed genes.

DNMT3B rs2424913 as a Risk Factor for Congenital Heart Defects in Down Syndrome

Impairments of the genes that encode enzymes that are involved in one-carbon metabolism because of the presence of gene polymorphisms can affect the methylation pattern. The altered methylation profiles of the genes involved in cardiogenesis may result in congenital heart defects (CHDs). The aim of this study was to investigate the association between the MTHFR rs1801133, MTHFR rs1801131, MTRR rs1801394, DNMT1 rs2228611, DNMT3A rs1550117, DNMT3B rs1569686, and DNMT3B rs2424913 gene polymorphisms and congenital heart defects in Down syndrome (DS) individuals. The study was conducted on 350 participants, including 134 DS individuals with CHDs (DSCHD+), 124 DS individuals without CHDs (DSCHD−), and 92 individuals with non-syndromic CHD. The genotyping was performed using the PCR–RFLP method. A statistically significant higher frequency of the DNMT3B rs2424913 TT in the DSCHD+ individuals was observed. The DNMT3B rs2424913 TT genotype, as well as the T allele, had significantly higher frequencies in the individuals with DS and atrial septal defects (ASDs) in comparison with the individuals with DS and other CHDs. Furthermore, our results indicate a statistically significant effect of the DNMT3B rs1569686 TT genotype in individuals with non-syndromic CHDs. The results of the study suggest that the DNMT3B rs2424913 TT genotypes may be a possible predisposing factor for CHDs in DS individuals, and especially those with ASDs.

The Mutagenic Consequences of DNA Methylation within and across Generations

DNA methylation is an epigenetic modification with wide-ranging consequences across the life of an organism. This modification can be stable, persisting through development despite changing environmental conditions. However, in other contexts, DNA methylation can also be flexible, underlying organismal phenotypic plasticity. One underappreciated aspect of DNA methylation is that it is a potent mutagen; methylated cytosines mutate at a much faster rate than other genetic motifs. This mutagenic property of DNA methylation has been largely ignored in eco-evolutionary literature, despite its prevalence. Here, we explore how DNA methylation induced by environmental and other factors could promote mutation and lead to evolutionary change at a more rapid rate and in a more directed manner than through stochastic genetic mutations alone. We argue for future research on the evolutionary implications of DNA methylation driven mutations both within the lifetime of organisms, as well as across timescales.

DNA methylation as a mediator of genetic and environmental influences on Parkinson's disease susceptibility: Impacts of alpha-Synuclein, physical activity, and pesticide exposure on the epigenome

Parkinson's disease (PD) is a neurodegenerative disorder with a complex etiology and increasing prevalence worldwide. As PD is influenced by a combination of genetic and environment/lifestyle factors in approximately 90% of cases, there is increasing interest in identification of the interindividual mechanisms underlying the development of PD as well as actionable lifestyle factors that can influence risk. This narrative review presents an outline of the genetic and environmental factors contributing to PD risk and explores the possible roles of cytosine methylation and hydroxymethylation in the etiology and/or as early-stage biomarkers of PD, with an emphasis on epigenome-wide association studies (EWAS) of PD conducted over the past decade. Specifically, we focused on variants in the SNCA gene, exposure to pesticides, and physical activity as key contributors to PD risk. Current research indicates that these factors individually impact the epigenome, particularly at the level of CpG methylation. There is also emerging evidence for interaction effects between genetic and environmental contributions to PD risk, possibly acting across multiple omics layers. We speculated that this may be one reason for the poor replicability of the results of EWAS for PD reported to date. Our goal is to provide direction for future epigenetics studies of PD to build upon existing foundations and leverage large datasets, new technologies, and relevant statistical approaches to further elucidate the etiology of this disease.

Epigenetic regulation of fetal brain development in pig

How fetal brain development is regulated at the molecular level is not well understood. Due to ethical challenges associated with research on the human fetus, large animals particularly pigs are increasingly used to study development and disorders of fetal brain. The pig fetal brain grows rapidly during the last ∼ 50 days before birth which is around day 60 (d60) of pig gestation. But what regulates the onset of accelerated growth of the brain is unknown. The current study tests the hypothesis that epigenetic alteration around d60 is involved in the onset of rapid growth of fetal brain of pig. To test this hypothesis, DNA methylation changes of fetal brain was assessed in a genome-wide manner by Enzymatic Methyl-seq (EM-seq) during two gestational periods (GP): d45 vs. d60 (GP1) and d60 vs. d90 (GP2). The cytosine-guanine (CpG) methylation data was analyzed in an integrative manner with the RNA-seq data generated from the same brain samples from our earlier study. A neural network based modeling approach was implemented to learn changes in methylation patterns of the differentially expressed genes, and then predict methylations of the brain in a genome-wide manner during rapid growth. This approach identified specific methylations that changed in a mutually informative manner during rapid growth of the fetal brain. These methylations were significantly overrepresented in specific genic as well as intergenic features including CpG islands, introns, and untranslated regions. In addition, sex-bias methylations of known single nucleotide polymorphic sites were also identified in the fetal brain ide during rapid growth.

Bioinformatic Prioritization and Functional Annotation of GWAS-Based Candidate Genes for Primary Open-Angle Glaucoma

Background: Primary open-angle glaucoma (POAG) is the most prevalent glaucoma subtype, but its exact etiology is still unknown. In this study, we aimed to prioritize the most likely ‘causal’ genes and identify functional characteristics and underlying biological pathways of POAG candidate genes. Methods: We used the results of a large POAG genome-wide association analysis study from GERA and UK Biobank cohorts. First, we performed systematic gene-prioritization analyses based on: (i) nearest genes; (ii) nonsynonymous single-nucleotide polymorphisms; (iii) co-regulation analysis; (iv) transcriptome-wide association studies; and (v) epigenomic data. Next, we performed functional enrichment analyses to find overrepresented functional pathways and tissues. Results: We identified 142 prioritized genes, of which 64 were novel for POAG. BICC1, AFAP1, and ABCA1 were the most highly prioritized genes based on four or more lines of evidence. The most significant pathways were related to extracellular matrix turnover, transforming growth factor-β, blood vessel development, and retinoic acid receptor signaling. Ocular tissues such as sclera and trabecular meshwork showed enrichment in prioritized gene expression (>1.5 fold). We found pleiotropy of POAG with intraocular pressure and optic-disc parameters, as well as genetic correlation with hypertension and diabetes-related eye disease. Conclusions: Our findings contribute to a better understanding of the molecular mechanisms underlying glaucoma pathogenesis and have prioritized many novel candidate genes for functional follow-up studies.

The Prospective Application of Melatonin in Treating Epigenetic Dysfunctional Diseases

In the past, different human disorders were described by scientists from the perspective of either environmental factors or just by genetically related mechanisms. The rise in epigenetic studies and its modifications, i.e., heritable alterations in gene expression without changes in DNA sequences, have now been confirmed in diseases. Modifications namely, DNA methylation, posttranslational histone modifications, and non-coding RNAs have led to a better understanding of the coaction between epigenetic alterations and human pathologies. Melatonin is a widely-produced indoleamine regulator molecule that influences numerous biological functions within many cell types. Concerning its broad spectrum of actions, melatonin should be investigated much more for its contribution to the upstream and downstream mechanistic regulation of epigenetic modifications in diseases. It is, therefore, necessary to fill the existing gaps concerning corresponding processes associated with melatonin with the physiological abnormalities brought by epigenetic modifications. This review outlines the findings on melatonin’s action on epigenetic regulation in human diseases including neurodegenerative diseases, diabetes, cancer, and cardiovascular diseases. It summarizes the ability of melatonin to act on molecules such as proteins and RNAs which affect the development and progression of diseases.

Alpha-synuclein overexpression induces epigenomic dysregulation of glutamate signaling and locomotor pathways

Parkinson’s disease (PD) is a neurological disorder with complex interindividual etiology that is becoming increasingly prevalent worldwide. Elevated alpha-synuclein levels can increase risk of PD and may influence epigenetic regulation of PD pathways. Here, we report genome-wide DNA methylation and hydroxymethylation alterations associated with overexpression of two PD-linked alpha-synuclein variants (wild-type and A30P) in LUHMES cells differentiated to dopaminergic neurons. Alpha-synuclein altered DNA methylation at thousands of CpGs and DNA hydroxymethylation at hundreds of CpGs in both genotypes, primarily in locomotor behavior and glutamate signaling pathway genes. In some cases, epigenetic changes were associated with transcription. SMITE network analysis incorporating H3K4me1 ChIP-seq to score DNA methylation and hydroxymethylation changes across promoters, enhancers, and gene bodies confirmed epigenetic and transcriptional deregulation of glutamate signaling modules in both genotypes. Our results identify distinct and shared impacts of alpha-synuclein variants on the epigenome, and associate alpha-synuclein with the epigenetic etiology of PD.

Maternal LINE-1 DNA Methylation in Early Spontaneous Preterm Birth

Despite considerable effort aimed at decreasing the incidence of spontaneous preterm birth (SPTB), it remains the leading cause of infant mortality and morbidity. The aim of this study was to evaluate maternal LINE-1 DNA methylation (DNAm), along with DNMT polymorphisms and factors proposed to modulate DNAm, in patients who delivered early preterm. This case-control study included women who delivered spontaneously early preterm (23-33^6/7 weeks of gestation), and control women. DNAm was analyzed in peripheral blood lymphocytes by quantification of LINE-1 DNAm using the MethyLight method. There was no significant difference in LINE-1 DNAm between patients with early PTB and controls. Among the investigated predictors, only the history of previous PTB was significantly associated with LINE-1 DNAm in PTB patients (β = -0.407; R² = 0.131; p = 0.011). The regression analysis showed the effect of DNMT3B rs1569686 TT+TG genotypes on LINE-1 DNAm in patients with familial PTB (β = -0.524; R² = 0.275; p = 0.037). Our findings suggest novel associations of maternal LINE-1 DNA hypomethylation with DNMT3B rs1569686 T allele. These results also contribute to the understanding of a complex (epi)genetic and environmental relationship underlying the early PTB.

Allele-specific DNA methylation maps in monozygotic twins discordant for psychiatric disorders reveal that disease-associated switching at the EIPR1 regulatory loci modulates neural function

The non-Mendelian features of phenotypic variations within monozygotic twins are likely complicated by environmental modifiers of genetic effects that have yet to be elucidated. Here, we performed methylome and genome analyses of blood DNA from psychiatric disorder-discordant monozygotic twins to study how allele-specific methylation (ASM) mediates phenotypic variations. We identified that thousands of genetic variants with ASM imbalances exhibit phenotypic variation-associated switching at regulatory loci. These ASMs have plausible causal associations with psychiatric disorders through effects on interactions between transcription factors, DNA methylations, and other epigenomic markers and then contribute to dysregulated gene expression, which eventually increases disease susceptibility. Moreover, we also experimentally validated the model that the rs4854158 alternative C allele at an ASM switching regulatory locus of EIPR1 encoding endosome-associated recycling protein-interacting protein 1, is associated with demethylation and higher RNA expression and shows lower TF binding affinities in unaffected controls. An epigenetic ASM switching induces C allele hypermethylation and then recruits repressive Polycomb repressive complex 2 (PRC2), reinforces trimethylation of lysine 27 on histone 3 and inhibits its transcriptional activity, thus leading to downregulation of EIPR1 in schizophrenia. Moreover, disruption of rs4854158 induces gain of EIPR1 function and promotes neural development and vesicle trafficking. Our study provides a powerful framework for identifying regulatory risk variants and contributes to our understanding of the interplay between genetic and epigenetic variants in mediating psychiatric disorder susceptibility.

Genetic and epigenetic associations of ANRIL with coronary artery disease and risk factors

BACKGROUND

Both DNA genotype and methylation of antisense non-coding RNA in the INK4 locus (ANRIL) have been robustly associated with coronary artery disease (CAD), but the interdependent mechanisms of genotype and methylation remain unclear.

METHODS

Eighteen tag single nucleotide polymorphisms (SNPs) of ANRIL were genotyped in a matched case-control study (cases 503 and controls 503). DNA methylation of ANRIL and the INK4/ARF locus (p14^ARF, p15^INK4b and p16^INK4a) was measured using pyrosequencing in the same set of samples (cases 100 and controls 100).

RESULTS

Polymorphisms of ANRIL (rs1004638, rs1333048 and rs1333050) were significantly associated with CAD (p < 0.05). The incidence of CAD, multi-vessel disease, and modified Gensini scores demonstrated a strong, direct association with ANRIL gene dosage (p < 0.05). There was no significant association between ANRIL polymorphisms and myocardial infarction/acute coronary syndrome (MI/ACS) (p > 0.05). Methylation levels of ANRIL were similar between the two studied groups (p > 0.05), but were different in the rs1004638 genotype, with AA and AT genotype having a higher level of ANRIL methylation (pos4, p = 0.006; pos8, p = 0.019). Further Spearman analyses indicated that methylation levels of ANRIL were positively associated with systolic blood pressure (pos6, r = 0.248, p = 0.013), diastolic blood pressure (pos3, r = 0.213, p = 0.034; pos6, r = 0.220, p = 0.028), and triglyceride (pos4, r = 0.253, p = 0.013), and negatively associated with high-density lipoprotein cholesterol (pos2, r = - 0.243, p = 0.017). Additionally, we identified 12 transcription factor binding sites (TFBS) within the methylated ANRIL region, and functional annotation indicated these TFBS were associated with basal transcription. Methylation at the INK4/ARF locus was not associated with ANRIL genotype.

CONCLUSIONS

These results indicate that ANRIL genotype (tag SNPs rs1004638, rs1333048 and rs1333050) mainly affects coronary atherosclerosis, but not MI/ACS. There may be allele-related DNA methylation and allele-related binding of transcription factors within the ANRIL promoter.

Epigenetic marks in suicide: a review

Suicide is a complex phenomenon and a global public health problem that involves several biological factors that could contribute to the pathophysiology of suicide. There is evidence that epigenetic factors influence some psychiatric disorders, suggesting a predisposition to suicide or suicidal behavior. Here, we review studies of molecular mechanisms of suicide in an epigenetic perspective in the postmortem brain of suicide completers and peripheral blood cells of suicide attempters. Besides, we include studies of gene-specific DNA methylation, epigenome-wide association, histone modification, and interfering RNAs as epigenetic factors. This review provides an overview of the epigenetic mechanisms described in different biological systems related to suicide, contributing to an understanding of the genetic regulation in suicide. We conclude that epigenetic marks are potential biomarkers in suicide, and they could become attractive therapeutic targets due to their reversibility and importance in regulating gene expression.

Disease-associated genetic variants in the regulatory regions of human genes: mechanisms of action on transcription and genomic resources for dissecting these mechanisms

Whole genome and whole exome sequencing technologies play a very important role in the studies of the genetic aspects of the pathogenesis of various diseases. The ample use of genome-wide and exome-wide association study methodology (GWAS and EWAS) made it possible to identify a large number of genetic variants associated with diseases. This information is accumulated in the databases like GWAS central, GWAS catalog, OMIM, ClinVar, etc. Most of the variants identified by the GWAS technique are located in the noncoding regions of the human genome. According to the ENCODE project, the fraction of regions in the human genome potentially involved in transcriptional control is many times greater than the fraction of coding regions. Thus, genetic variation in noncoding regions of the genome can increase the susceptibility to diseases by disrupting various regulatory elements (promoters, enhancers, silencers, insulator regions, etc.). However, identification of the mechanisms of influence of pathogenic genetic variants on the diseases risk is difficult due to a wide variety of regulatory elements. The present review focuses on the molecular genetic mechanisms by which pathogenic genetic variants affect gene expression. At the same time, attention is concentrated on the transcriptional level of regulation as an initial step in the expression of any gene. A triggering event mediating the effect of a pathogenic genetic variant on the level of gene expression can be, for example, a change in the functional activity of transcription factor binding sites (TFBSs) or DNA methylation change, which, in turn, affects the functional activity of promoters or enhancers. Dissecting the regulatory roles of polymorphic loci have been impossible without close integration of modern experimental approaches with computer analysis of a growing wealth of genetic and biological data obtained using omics technologies. The review provides a brief description of a number of the most well-known public genomic information resources containing data obtained using omics technologies, including (1) resources that accumulate data on the chromatin states and the regions of transcription factor binding derived from ChIP-seq experiments; (2) resources containing data on genomic loci, for which allele-specific transcription factor binding was revealed based on ChIP-seq technology; (3) resources containing in silico predicted data on the potential impact of genetic variants on the transcription factor binding sites

Non-Coding Variants in Cancer: Mechanistic Insights and Clinical Potential for Personalized Medicine

The cancer genome is characterized by extensive variability, in the form of Single Nucleotide Polymorphisms (SNPs) or structural variations such as Copy Number Alterations (CNAs) across wider genomic areas. At the molecular level, most SNPs and/or CNAs reside in non-coding sequences, ultimately affecting the regulation of oncogenes and/or tumor-suppressors in a cancer-specific manner. Notably, inherited non-coding variants can predispose for cancer decades prior to disease onset. Furthermore, accumulation of additional non-coding driver mutations during progression of the disease, gives rise to genomic instability, acting as the driving force of neoplastic development and malignant evolution. Therefore, detection and characterization of such mutations can improve risk assessment for healthy carriers and expand the diagnostic and therapeutic toolbox for the patient. This review focuses on functional variants that reside in transcribed or not transcribed non-coding regions of the cancer genome and presents a collection of appropriate state-of-the-art methodologies to study them.

Whole-Genome Differentially Hydroxymethylated DNA Regions among Twins Discordant for Cardiovascular Death

Epigenetics is a mechanism underlying cardiovascular disease. It is unknown whether DNA hydroxymethylation is prospectively associated with the risk for cardiovascular death independent of germline and common environment. Male twin pairs middle-aged in 1969–1973 and discordant for cardiovascular death through December 31, 2014, were included. Hydroxymethylation was quantified in buffy coat DNA collected in 1986–1987. The 1893 differentially hydroxymethylated regions (DhMRs) were identified after controlling for blood leukocyte subtypes and age among 12 monozygotic (MZ) pairs (Benjamini–Hochberg False Discovery Rate < 0.01), of which the 102 DhMRs were confirmed with directionally consistent log₂-fold changes and p < 0.01 among additional 7 MZ pairs. These signature 102 DhMRs, independent of the germline, were located on all chromosomes except for chromosome 21 and the Y chromosome, mainly within/overlapped with intergenic regions and introns, and predominantly hyper-hydroxymethylated. A binary linear classifier predicting cardiovascular death among 19 dizygotic pairs was identified and equivalent to that generated from MZ via the 2D transformation. Computational bioinformatics discovered pathways, phenotypes, and DNA motifs for these DhMRs or their subtypes, suggesting that hydroxymethylation was a pathophysiological mechanism underlying cardiovascular death that might be influenced by genetic factors and warranted further investigations of mechanisms of these signature regions in vivo and in vitro.

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.

Elevation of DNA Methylation in the Promoter Regions of the Brain-Derived Neurotrophic Factor Gene is Associated with Heroin Addiction

To study the potential role of brain-derived neurotrophic factor (BDNF) methylation in heroin addiction, we first detected the methylation level of seven CpG islands that included 106 CpG sites in the promoter regions of BDNF from 120 people addicted to heroin and 113 controls. Methylation quantitative trait locus (mQTL) analysis was then employed to determine the association between the single-nucleotide polymorphism rs6265, a well-known locus shown to be correlated with heroin addiction, and the methylation levels of these CpG sites. Finally, we used the JASPAR database to predict whether transcription factors could bind to these CpG sites. We found that the methylation levels of CpG islands 6 and 7 and the methylation levels of BDNF_45 and BDNF_80 were significantly higher in the heroin addiction group than in the control group. We also found that rs6265 was an mQTL and was associated with the methylation level of BDNF_58. Using the JASPAR database, we found that ALX homeobox 3 (ALX3), achaete-scute family bHLH transcription factor 1 (ASCL1) and aryl hydrocarbon receptor nuclear translocator 2 (ARNT2) could bind to CpG island 6, and ALX3 could bind to CpG island 7. In summary, we showed that increased DNA methylation in the promoter regions of the BDNF gene was associated with heroin addiction in Han Chinese.

miR-122-5p Mediates Fluoride-Induced Osteoblast Activation by Targeting CDK4

Chronic intake of fluoride, existing in the environment, may cause endemic fluorosis, which is characterized by the occurrence of skeletal and dental fluorosis. However, the pathogenesis of fluorosis has not yet been elucidated. Abnormal osteoblast proliferation and activation have a pivotal role in bone turnover disorders which are linked to skeletal fluorosis. MicroRNAs are involved in fundamental cellular processes, including cell proliferation. Based on our previous study, population study and in vitro experiments were designed to understand the effect of miR-122-5p on osteoblast activation in skeletal fluorosis through targeting cyclin-dependent kinase 4 (CDK4). In human populations with coal-burning type fluoride exposure, the results showed that miR-122-5p was downregulated but CDK4 expression was upregulated and miR-122-5p was negatively correlated with CDK4 expression. Furthermore, in human osteoblasts treated with sodium fluoride, we demonstrated that miR-122-5p mediated osteoblast activation of skeletal fluorosis via upregulation of the CDK4 protein. In support of this, dual-luciferase reporter assay showed that miR-122-5p modulated CDK4 protein levels by targeting its 3'-untranslated region. These findings show, for the first time, that miR-122-5p may be involved in the cause and development of skeletal fluorosis by targeting CDK4.

Exploring Common Therapeutic Targets for Neurodegenerative Disorders Using Transcriptome Study

Alzheimer's disease (AD) and Parkinson's disease (PD) are well-known neuronal degenerative disorders that share common pathological events. Approved medications alleviate symptoms but do not address the root cause of the disease. Energy dysfunction in the neuronal population leads to various pathological events and ultimately results in neuronal death. Identifying common therapeutic targets for these disorders may help in the drug discovery process. The Brodmann area 9 (BA9) region is affected in both the disease conditions and plays an essential role in cognitive, motor, and memory-related functions. Analyzing transcriptome data of BA9 provides deep insights related to common pathological pathways involved in AD and PD. In this work, we map the preprocessed BA9 fastq files generated by RNA-seq for disease and control samples with reference hg38 genomic assembly and identify common variants and differentially expressed genes (DEG). These variants are predominantly located in the 3' UTR (non-promoter) region, affecting the conserved transcription factor (TF) binding motifs involved in the methylation and acetylation process. We have constructed BA9-specific functional interaction networks, which show the relationship between TFs and DEGs. Based on expression signature analysis, we propose that MAPK1, VEGFR1/FLT1, and FGFR1 are promising drug targets to restore blood-brain barrier functionality by reducing neuroinflammation and may save neurons.

Epigenome-wide association study on asthma and chronic obstructive pulmonary disease overlap reveals aberrant DNA methylations related to clinical phenotypes

We hypothesized that epigenetics is a link between smoking/allergen exposures and the development of Asthma and chronic obstructive pulmonary disease (ACO). A total of 75 of 228 COPD patients were identified as ACO, which was independently associated with increased exacerbations. Microarray analysis identified 404 differentially methylated loci (DML) in ACO patients, and 6575 DML in those with rapid lung function decline in a discovery cohort. In the validation cohort, ACO patients had hypermethylated PDE9A (+ 30,088)/ZNF323 (− 296), and hypomethylated SEPT8 (− 47) genes as compared with either pure COPD patients or healthy non-smokers. Hypermethylated TIGIT (− 173) gene and hypomethylated CYSLTR1 (+ 348)/CCDC88C (+ 125,722)/ADORA2B (+ 1339) were associated with severe airflow limitation, while hypomethylated IFRD1 (− 515) gene with frequent exacerbation in all the COPD patients. Hypermethylated ZNF323 (− 296) / MPV17L (+ 194) and hypomethylated PTPRN2 (+ 10,000) genes were associated with rapid lung function decline. In vitro cigarette smoke extract and ovalbumin concurrent exposure resulted in specific DNA methylation changes of the MPV17L / ZNF323 genes, while 5-aza-2′-deoxycytidine treatment reversed promoter hypermethylation-mediated MPV17L under-expression accompanied with reduced apoptosis and decreased generation of reactive oxygen species. Aberrant DNA methylations may constitute a determinant for ACO, and provide a biomarker of airflow limitation, exacerbation, and lung function decline.

Differentially methylated regions within lung cancer risk loci are enriched in deregulated enhancers

Genome-wide association studies (GWAS) have identified SNPs linked with lung cancer risk. Our aim was to discover the genes, non-coding RNAs, and regulatory elements within GWAS-identified risk regions that are deregulated in non-small cell lung carcinoma (NSCLC) to identify novel, clinically targetable genes and mechanisms in carcinogenesis. A targeted bisulphite-sequencing approach was used to comprehensively investigate DNA methylation changes occurring within lung cancer risk regions in 17 NSCLC and adjacent normal tissue pairs. We report differences in differentially methylated regions between adenocarcinoma and squamous cell carcinoma. Among the minimal regions found to be differentially methylated in at least 50% of the patients, 7 candidates were replicated in 2 independent cohorts (n = 27 and n = 87) and the potential of 6 as methylation-dependent regulatory elements was confirmed by functional assays. This study contributes to understanding the pathways implicated in lung cancer initiation and progression, and provides new potential targets for cancer treatment.

DNA methylation differences associated with social anxiety disorder and early life adversity

Social anxiety disorder (SAD) is a psychiatric disorder characterized by extensive fear in social situations. Multiple genetic and environmental factors are known to contribute to its pathogenesis. One of the main environmental risk factors is early life adversity (ELA). Evidence is emerging that epigenetic mechanisms such as DNA methylation might play an important role in the biological mechanisms underlying SAD and ELA. To investigate the relationship between ELA, DNA methylation, and SAD, we performed an epigenome-wide association study for SAD and ELA examining DNA from whole blood of a cohort of 143 individuals using DNA methylation arrays. We identified two differentially methylated regions (DMRs) associated with SAD located within the genes SLC43A2 and TNXB. As this was the first epigenome-wide association study for SAD, it is worth noting that both genes have previously been associated with panic disorder. Further, we identified two DMRs associated with ELA within the SLC17A3 promoter region and the SIAH3 gene and several DMRs that were associated with the interaction of SAD and ELA. Of these, the regions within C2CD2L and MRPL28 showed the largest difference in DNA methylation. Lastly, we found that two DMRs were associated with both the severity of social anxiety and ELA, however, neither of them was found to mediate the contribution of ELA to SAD later in life. Future studies are needed to replicate our findings in independent cohorts and to investigate the biological pathways underlying these effects.

The correlation of salivary telomere length and single nucleotide polymorphisms of the ADIPOQ, SIRT1 and FOXO3A genes with lifestyle-related diseases in a Japanese population

Background

It has been reported that genetic factors are associated with risk factors and onset of lifestyle-related diseases, but this finding is still the subject of much debate.

Objective

The aim of the present study was to investigate the correlation of genetic factors, including salivary telomere length and three single nucleotide polymorphisms (SNPs) that may influence lifestyle-related diseases, with lifestyle-related diseases themselves.

Methods

In one year at a single facility, relative telomere length and SNPs were determined by using monochrome multiplex quantitative polymerase chain reaction and TaqMan SNP Genotyping Assays, respectively, and were compared with lifestyle-related diseases in 120 Japanese individuals near our university.

Results

In men and all participants, age was inversely correlated with relative telomere length with respective p values of 0.049 and 0.034. In men, the frequency of hypertension was significantly higher in the short relative telomere length group than in the long group with unadjusted p value of 0.039, and the difference in the frequency of hypertension between the two groups was of borderline statistical significance after adjustment for age (p = 0.057). Furthermore, in men and all participants, the sum of the number of affected lifestyle-related diseases, including hypertension, was significantly higher in the short relative telomere length group than in the long group, with p values of 0.004 and 0.029, respectively. For ADIPOQ rs1501299, men’s ankle brachial index was higher in the T/T genotype than in the G/G and G/T genotypes, with p values of 0.001 and 0.000, respectively. For SIRT1 rs7895833, men’s body mass index and waist circumference and all participants’ brachial-ankle pulse wave velocity were higher in the A/G genotype than in the G/G genotype, with respective p values of 0.048, 0.032 and 0.035. For FOXO3A rs2802292, women’s body temperature and all participants’ saturation of peripheral oxygen were lower in the G/T genotype than in the T/T genotype, with respective p values of 0.039 and 0.032. However, relative telomere length was not associated with physiological or anthropometric measurements except for height in men (p = 0.016). ADIPOQ rs1501299 in men, but not the other two SNPs, was significantly associated with the sum of the number of affected lifestyle-related diseases (p = 0.013), by genotype. For each SNPs, there was no significant difference in the frequency of hypertension or relative telomere length by genotype.

Conclusion

Relative telomere length and the three types of SNPs determined using saliva have been shown to be differentially associated with onset of and measured risk factors for lifestyle-related diseases consisting mainly of cardiovascular diseases and cancer.

A Primer on DNA Methylation and Its Potential to Impact Maternal Depression Risk and Assessment During Pregnancy and the Postpartum

Depression onset during and after pregnancy is prevalent and associated with significant implications for maternal, child, and family health. Although environmental risk factors important to the expression of pregnancy-related depression are well known, knowledge of the genetic underpinning is limited. Given the joint contribution of environmental and genetic factors to depression risk liability, DNA methylation presents itself as an ideal biomarker to investigate basic mechanisms and opportunities for translational research to care for pregnancy-related depression health outcomes. This article is an introduction to DNA methylation and its potential to serve as a marker of depression risk during pregnancy and the postpartum. This commentary discusses current clinical uses of DNA methylation-based testing and how it may be applied to perinatal depression clinical care and management.

A review of linc00673 as a novel lncRNA for tumor regulation

Long non-coding RNAs (LncRNAs) act as regulators and play important roles in a variety of biological processes. These regulators constitute a huge information network among genes and participate in the pathophysiological process of human diseases. Increasing evidence has demonstrated that LncRNA, as an oncogene or tumor suppressor gene, is closely related to the occurrence and development of tumors. Linc00673 is a recently discovered LncRNA molecule that is dysregulated in several solid tumors. Moreover, its genetic polymorphism is believed to affect the susceptibility of a population to the corresponding cancer species. This article summarizes the role of Linc00673 in different human cancers and its molecular mechanisms with a focus on the characteristics of Linc00673 and the existing literature on it while highlighting the future research directions for Linc00673. Linc00673 has the potential to become a feasible clinical diagnostic and prognostic marker toward providing a new molecular therapeutic target for cancer patients.

DNA methylation analysis of candidate genes associated with dementia in peripheral blood

Aim: To investigate whether genes implicated in dementia pathogenesis are differently methylated in peripheral blood. Materials & methods: Participants included 160 cognitively healthy individuals aged 70+ years: 73 who were subsequently diagnosed with dementia and 87 controls matched on age, gender, education, smoking and baseline cognition. A total of 49 participants also provided blood samples at diagnosis. Blood DNA methylation of APOE, APP, BDNF, PIN1, SNCA and TOMM40 was examined. Results: A total of 56 of 299 probes were differentially methylated in dementia compared with controls and 39 probes prior to diagnosis. The greatest effect size was in APP (cg19423170, Δ-8.32%, adjusted p = 0.009 at diagnosis; cg19933173, Δ-4.18%, adjusted p < 0.0001 prediagnosis). Conclusion: Genes implicated in dementia pathogenesis show differential blood methylation in dementia, even prior to diagnosis.

Genetic Aspects of Inflammation and Immune Response in Stroke

Genetic determinants play important role in the complex processes of inflammation and immune response in stroke and could be studied in different ways. Inflammation and immunomodulation are associated with repair processes in ischemic stroke, and together with the concept of preconditioning are promising modes of stroke treatment. One of the important aspects to be considered in the recovery of patients after the stroke is a genetic predisposition, which has been studied extensively. Polymorphisms in a number of candidate genes, such as IL-6, BDNF, COX2, CYPC19, and GPIIIa could be associated with stroke outcome and recovery. Recent GWAS studies pointed to the variant in genesPATJ and LOC as new genetic markers of long term outcome. Epigenetic regulation of immune response in stroke is also important, with mechanisms of histone modifications, DNA methylation, and activity of non-coding RNAs. These complex processes are changing from acute phase over the repair to establishing homeostasis or to provoke exaggerated reaction and death. Pharmacogenetics and pharmacogenomics of stroke cures might also be evaluated in the context of immuno-inflammation and brain plasticity. Potential novel genetic treatment modalities are challenged but still in the early phase of the investigation.

SNPs in Sites for DNA Methylation, Transcription Factor Binding, and miRNA Targets Leading to Allele-Specific Gene Expression and Contributing to Complex Disease Risk: A Systematic Review

INTRODUCTION

The complex genetic diversity among human populations results from an assortment of factors acting at various sequential levels, including mutations, population migrations, genetic drift, and selection. Although there are a plethora of DNA sequence variations identified through genome-wide association studies (GWAS), the challenge remains to explain the mechanisms underlying interindividual phenotypic disparity accounting for disease susceptibility. Single nucleotide polymorphisms (SNPs) present in the sites for DNA methylation, transcription factor (TF) binding, or miRNA targets can alter the gene expression. The systematic review aimed to evaluate the complex crosstalk among SNPs, miRNAs, DNA methylation, and TFs for complex multifactorial disease risk.

METHODS

PubMed and Scopus databases were used from inception until May 15, 2019. Initially, screening of articles involved studies assessing the interaction of SNPs with TFs, DNA methylation, or miRNAs resulting in allele-specific gene expression in complex multifactorial diseases. We also included the studies which provided experimental validation of the interaction of SNPs with each of these factors. The results from various studies on multifactorial diseases were assessed.

RESULTS

A total of 11 articles for SNPs interacting with DNA methylation, 30 articles for SNPs interacting with TFs, and 11 articles for SNPs in miRNA binding sites were selected. The interactions of SNPs with epigenetic factors were found to be implicated in different types of cancers, autoimmune diseases, cardiovascular diseases, diabetes, and asthma.

CONCLUSION

The systematic review provides evidence for the interplay between genetic and epigenetic risk factors through allele-specific gene expression in various complex multifactorial diseases.

CpG-SNP site methylation regulates allele-specific expression of MTHFD1 gene in type 2 diabetes

The interaction of genetic and epigenetic mechanisms is one of the underlying causes of phenotypic variability in complex diseases such as type 2 diabetes (T2D). To explore the influence of genetic and epigenetic changes in T2D, we examined the effect of methylation of CpG-SNP sites on allele-specific expression (ASE) in one-carbon metabolism pathway genes in T2D. Case-control study was conducted on 860 individuals (430 T2D and 430 controls). CpG-SNPs shortlisted through in silico analysis were genotyped using tetra ARMS PCR and validated using Sanger DNA sequencing. Global DNA methylation was carried out using RP-HPLC. Promoter DNA methylation and CpG site-specific methylation were carried out using bisulfite sequencing. mRNA expression and ASE were examined by SYBR green and TaqMan assay, respectively. Four exonic CpG-SNPs of MTHFD1, MTRR, and GGH genes were identified in folate pathway genes. Among these, MTHFD1 rs2236225 showed significant association with T2D independent of obesity, displayed ASE, and correlated with CpG-SNP site-specific methylation when compared with controls. Our results demonstrate that SNP rs2236225 in the CpG site of MTHFD1, which regulates allele-specific gene expression in PBMCs is methylation dependent and may perturb one-carbon metabolism pathway in T2D subjects.

Epigenetic Vulnerability of Insulator CTCF Motifs at Parkinson's Disease-Associated Genes in Response to Neurotoxicant Rotenone

CCCTC-binding factor (CTCF) is a regulatory protein that binds DNA to control spatial organization and transcription. The sequence-specific binding of CTCF is variable and is impacted by nearby epigenetic patterns. It has been demonstrated that non-coding genetic variants cluster with CTCF sites in topological associating domains and thus can affect CTCF activity on gene expression. Therefore, environmental factors that alter epigenetic patterns at CTCF binding sites may dictate the interaction of non-coding genetic variants with regulatory proteins. To test this mechanism, we treated human cell line HEK293 with rotenone for 24 h and characterized its effect on global epigenetic patterns specifically at regulatory regions of Parkinson's disease (PD) risk loci. We used RNA sequencing to examine changes in global transcription and identified over 2000 differentially expressed genes (DEGs, >1.5-fold change, FDR < 0.05). Among these DEGs, 13 were identified as PD-associated genes according to Genome-wide association studies meta-data. We focused on eight genes that have non-coding risk variants and a prominent CTCF binding site. We analyzed methylation of a total of 165 CGs surrounding CTCF binding sites and detected differential methylation (|>1%|, q < 0.05) in 45 CGs at 7 PD-associated genes. Of these 45 CGs, 47% were hypomethylated and 53% were hypermethylated. Interestingly, 5 out of the 7 genes had correlated gene upregulation with CG hypermethylation at CTCF and gene downregulation with CG hypomethylation at CTCF. We also investigated active H3K27ac surrounding the same CTCF binding sites within these seven genes. We observed a significant increase in H3K27ac in four genes (FDR < 0.05). Three genes (PARK2, GPRIN3, FER) showed increased CTCF binding in response to rotenone. Our data indicate that rotenone alters regulatory regions of PD-associated genes through changes in epigenetic patterns, and these changes impact high-order chromatin organization to increase the influence of non-coding variants on genome integrity and cellular survival.

Vitamin D receptor gene methylation in patients with hand, foot, and mouth disease caused by enterovirus 71

To evaluate the epigenetic regulation of the VDR gene in enterovirus 71 (EV71)-associated severe hand, foot, and mouth disease (HFMD), a total of 116 patients with EV71-HFMD, including 58 with mild EV71-HFMD and 58 with severe EV71-HFMD, as well as 60 healthy controls, were enrolled in this study. Quantitative real-time PCR was used to measure the relative levels of VDR mRNA expression, and the methylation status of the VDR promoter was assessed using a MethylTarget™ assay. The DNA methylation levels of the VDR promoter in children with EV71-associated severe HFMD were lower than those in the healthy controls and in children with mild HFMD (P < 0.05). Hypomethylation at CpG site 133 and hypermethylation at the CpG 42 sites and 68 downregulated VDR expression. Moreover, the methylation level of VDR could be used for differential diagnosis of mild and severe EV71-associated HFMD (AUC₅₆, 0.73; AUC₆₈, 0.699; AUC₄₂, 0.694; AUC₆₆, 0.693). VDR expression and promoter methylation were associated with the progression of EV71 infection. Determining the VDR promoter status might help clinicians initiate the appropriate strategy for treatment of EV71-associated HFMD.

DAMEfinder: a method to detect differential allele-specific methylation

Background

DNA methylation is a highly studied epigenetic signature that is associated with regulation of gene expression, whereby genes with high levels of promoter methylation are generally repressed. Genomic imprinting occurs when one of the parental alleles is methylated, i.e., when there is inherited allele-specific methylation (ASM). A special case of imprinting occurs during X chromosome inactivation in females, where one of the two X chromosomes is silenced, to achieve dosage compensation between the sexes. Another more widespread form of ASM is sequence dependent (SD-ASM), where ASM is linked to a nearby heterozygous single nucleotide polymorphism (SNP).

Results

We developed a method to screen for genomic regions that exhibit loss or gain of ASM in samples from two conditions (treatments, diseases, etc.). The method relies on the availability of bisulfite sequencing data from multiple samples of the two conditions. We leverage other established computational methods to screen for these regions within a new R package called DAMEfinder. It calculates an ASM score for all CpG sites or pairs in the genome of each sample, and then quantifies the change in ASM between conditions. It then clusters nearby CpG sites with consistent change into regions. In the absence of SNP information, our method relies only on reads to quantify ASM. This novel ASM score compares favorably to current methods that also screen for ASM. Not only does it easily discern between imprinted and non-imprinted regions, but also females from males based on X chromosome inactivation. We also applied DAMEfinder to a colorectal cancer dataset and observed that colorectal cancer subtypes are distinguishable according to their ASM signature. We also re-discover known cases of loss of imprinting.

Conclusion

We have designed DAMEfinder to detect regions of differential ASM (DAMEs), which is a more refined definition of differential methylation, and can therefore help in breaking down the complexity of DNA methylation and its influence in development and disease.

Genomic resources for dissecting the role of non-protein coding variation in gene-environment interactions

The majority of single nucleotide variants (SNVs) identified in Genome Wide Association Studies (GWAS) fall within non-protein coding DNA and have the potential to alter gene expression. Non-protein coding DNA can control gene expression by acting as transcription factor (TF) binding sites or by regulating the organization of DNA into chromatin. SNVs in non-coding DNA sequences can disrupt TF binding and chromatin structure and this can result in pathology. Further, environmental health studies have shown that exposure to xenobiotics can disrupt the ability of TFs to regulate entire gene networks and result in pathology. However, there is a large amount of interindividual variability in exposure-linked health outcomes. One explanation for this heterogeneity is that genetic variation and exposure combine to disrupt gene regulation, and this eventually manifests in disease. Many resources exist that annotate common variants from GWAS and combine them with conservation, functional genomics, and TF binding data. These annotation tools provide clues regarding the biological implications of an SNV, as well as lead to the generation of hypotheses regarding potentially disrupted target genes, epigenetic markers, pathways, and cell types. Collectively this information can be used to predict how SNVs can alter an individual's response to exposure and disease risk. A basic understanding of the regulatory information contained within non-protein coding DNA is needed to predict the biological consequences of SNVs, and to determine how these SNVs impact exposure-related disease. We hope that this review will aid in the characterization of disease-associated genetic variation in the non-protein coding genome.

Remote modulation of lncRNA GCLET by risk variant at 16p13 underlying genetic susceptibility to gastric cancer

The biological effects of susceptibility loci are rarely reported in gastric tumorigenesis. We conducted a large-scale cross-ancestry genetic study in 18,852 individuals and identified the potential causal variant rs3850997 T>G at 16p13 significantly associated with a decreased risk of gastric cancer [odds ratio (OR) = 0.87, 95% confidence interval (CI) = 0.83 to 0.91, P = 2.13 × 10-9]. This risk effect was mediated through the mapped long noncoding RNA GCLET (Gastric Cancer Low-Expressed Transcript; ORindirect = 0.987, 95% CI = 0.975 to 0.999, P = 0.018). Mechanistically, rs3850997 exerted an allele-specific long-range regulatory effect on GCLET by affecting the binding affinity of CTCF. Furthermore, GCLET increased FOXP2 expression by competing with miR-27a-3p, and this regulation remarkably affected in vitro, in vivo, and clinical gastric cancer phenotypes. The findings highlight the genetic functions and implications for the etiology and pathology of cancers.

Missing heritability in Parkinson's disease: the emerging role of non-coding genetic variation

Parkinson's disease (PD) is a neurodegenerative disorder caused by a complex interplay of genetic and environmental factors. For the stratification of PD patients and the development of advanced clinical trials, including causative treatments, a better understanding of the underlying genetic architecture of PD is required. Despite substantial efforts, genome-wide association studies have not been able to explain most of the observed heritability. The majority of PD-associated genetic variants are located in non-coding regions of the genome. A systematic assessment of their functional role is hampered by our incomplete understanding of genotype-phenotype correlations, for example through differential regulation of gene expression. Here, the recent progress and remaining challenges for the elucidation of the role of non-coding genetic variants is reviewed with a focus on PD as a complex disease with multifactorial origins. The function of gene regulatory elements and the impact of non-coding variants on them, and the means to map these elements on a genome-wide level, will be delineated. Moreover, examples of how the integration of functional genomic annotations can serve to identify disease-associated pathways and to prioritize disease- and cell type-specific regulatory variants will be given. Finally, strategies for functional validation and considerations for suitable model systems are outlined. Together this emphasizes the contribution of rare and common genetic variants to the complex pathogenesis of PD and points to remaining challenges for the dissection of genetic complexity that may allow for better stratification, improved diagnostics and more targeted treatments for PD in the future.

Targeted lipid nanoparticles for RNA therapeutics and immunomodulation in leukocytes

Abnormalities in leukocytes' function are associated with many immune related disorders, such as cancer, autoimmunity and susceptibility to infectious diseases. Recent developments in Genome-wide-association-studies give rise to new opportunities for novel therapeutics. RNA-based modalities, that allow a selective genetic manipulation in vivo, are powerful tools for personalized medicine, enabling downregulation or expression of relevant proteins. Yet, RNA-based therapeutics requires a delivery modality to facilitate the stability, uptake and intracellular release of the RNA molecules. The use of lipid nanoparticles as a drug delivery approach improves the payloads' stability, pharmacokinetics, bio-distribution and therapeutic benefit while reducing side effects. Moreover, a wide variety of targeting moieties allow a precise and modular manipulation of gene expression, together with the ability to identify and selectively affect disease-relevant leukocytes-subsets. Altogether, RNA-based therapeutics, targeting leukocytes subsets, is believed to be one of the most promising therapeutic concepts of the near future, addressing pressing issues in cancer and inflammation heterogeneity.

hOGG1 promoter methylation, hOGG1 genetic variants and their interactions for risk of coal-borne arsenicosis: A case-control study

To identify the effect of hOGG1 methylation, Ser326Cys polymorphism and their interactions on the risk of coal-borne arsenicosis, 113 coal-borne arsenicosis subjects and 55 reference subjects were recruited. Urinary arsenic contents were analyzed with ICP-MS. hOGG1 methylation and Ser326Cys polymorphism was measured by mehtylation-specific PCR and restriction fragment length polymorphism PCR in PBLCs, respectively. The results showed that the prevalence of methylated hOGG1 and variation genotype (326 ^Ser/Cys & 326 ^Cys/Cys) were increased with raised levels of urinary arsenic in arsenicosis subjects. Increased prevalence of methylated hOGG1 and variation genotype were associated with raised risk of arsenicosis. Moreover, the results revealed that variant genotype might increase the susceptibility to hOGG1 methylation. The interactions of methylated hOGG1 and variation genotype were also found to contribute to increased risk of arsenicosis. Taken together, hOGG1 hypermethylation, hOGG1 variants and their interactions might be potential biomarkers for evaluating risk of coal-borne arsenicosis.

The conserved DNMT1-dependent methylation regions in human cells are vulnerable to neurotoxicant rotenone exposure

BACKGROUND

Allele-specific DNA methylation (ASM) describes genomic loci that maintain CpG methylation at only one inherited allele rather than having coordinated methylation across both alleles. The most prominent of these regions are germline ASMs (gASMs) that control the expression of imprinted genes in a parent of origin-dependent manner and are associated with disease. However, our recent report reveals numerous ASMs at non-imprinted genes. These non-germline ASMs are dependent on DNA methyltransferase 1 (DNMT1) and strikingly show the feature of random, switchable monoallelic methylation patterns in the mouse genome. The significance of these ASMs to human health has not been explored. Due to their shared allelicity with gASMs, herein, we propose that non-traditional ASMs are sensitive to exposures in association with human disease.

RESULTS

We first explore their conservancy in the human genome. Our data show that our putative non-germline ASMs were in conserved regions of the human genome and located adjacent to genes vital for neuronal development and maturation. We next tested the hypothesized vulnerability of these regions by exposing human embryonic kidney cell HEK293 with the neurotoxicant rotenone for 24 h. Indeed,14 genes adjacent to our identified regions were differentially expressed from RNA-sequencing. We analyzed the base-resolution methylation patterns of the predicted non-germline ASMs at two neurological genes, HCN2 and NEFM, with potential to increase the risk of neurodegeneration. Both regions were significantly hypomethylated in response to rotenone.

CONCLUSIONS

Our data indicate that non-germline ASMs seem conserved between mouse and human genomes, overlap important regulatory factor binding motifs, and regulate the expression of genes vital to neuronal function. These results support the notion that ASMs are sensitive to environmental factors such as rotenone and may alter the risk of neurological disease later in life by disrupting neuronal development.

Analyzing the Relationship between Cohort and Case-Control Study Results Based on Model for Multiple Pathogenic Factors

Objective

Although the relative risk from a prospective cohort study is numerically approximate to the odds ratio from a case-control study for a low-probability event, a definite relationship between case-control and cohort studies cannot be confirmed. In this study, we established a different model to determine the relationship between case-control and cohort studies.

Methods

Two analysis models (the cross-sectional model and multiple pathogenic factor model) were established. Incidences in both the exposure group and the nonexposure group in a cohort study were compared with the frequency of the observed factor in each group (diseased and nondiseased) in a case-control study.

Results

The relationship between the results of a case-control study and a cohort study is as follows: Pe=(Pd∗m)/(Pc∗m)/(Pd∗m)/(Pn=(m)/(∗PdPc∗m)/(Pd∗m)/(Pe and Pn represent the incidence in the exposed group and nonexposed group, respectively, from the cohort study, while Pd and Pc represent the observed frequencies in the disease group and the control group, respectively, for the case-control study; finally, m)/(

Conclusions

There is a definite relationship between the results of case-control and cohort studies assessing the same exposure. The outcomes of case-control studies can be translated into cohort study data.

Epigenetic biomarkers and preterm birth

Preterm birth (PTB) is a major public health challenge, and novel, sensitive approaches to predict PTB are still evolving. Epigenomic markers are being explored as biomarkers of PTB because of their molecular stability compared to gene expression. This approach is also relatively new compared to gene-based diagnostics, which relies on mutations or single nucleotide polymorphisms. The fundamental principle of epigenome diagnostics is that epigenetic reprogramming in the target tissue (e.g. placental tissue) might be captured by more accessible surrogate tissue (e.g. blood) using biochemical epigenome assays on circulating DNA that incorporate methylation, histone modifications, nucleosome positioning, and/or chromatin accessibility. Epigenomic-based biomarkers may hold great potential for early identification of the majority of PTBs that are not associated with genetic variants or mutations. In this review, we discuss recent advances made in the development of epigenome assays focusing on its potential exploration for association and prediction of PTB. We also summarize population-level cohort studies conducted in the USA and globally that provide opportunities for genetic and epigenetic marker development for PTB. In addition, we summarize publicly available epigenome resources and published PTB studies. We particularly focus on ongoing genome-wide DNA methylation and epigenome-wide association studies. Finally, we review the limitations of current research, the importance of establishing a comprehensive biobank, and possible directions for future studies in identifying effective epigenome biomarkers to enhance health outcomes for pregnant women at risk of PTB and their infants.

Novel pharmacological targets in drug development for the treatment of anxiety and anxiety-related disorders

Current medication for anxiety disorders is suboptimal in terms of efficiency and tolerability, highlighting the need for improved drug treatments. In this review an overview of drugs being studied in different phases of clinical trials for their potential in the treatment of fear-, anxiety- and trauma-related disorders is presented. One strategy followed in drug development is refining and improving compounds interacting with existing anxiolytic drug targets, such as serotonergic and prototypical GABAergic benzodiazepines. A more innovative approach involves the search for compounds with novel mechanisms of anxiolytic action using the growing knowledge base concerning the relevant neurocircuitries and neurobiological mechanisms underlying pathological fear and anxiety. The target systems evaluated in clinical trials include glutamate, endocannabinoid and neuropeptide systems, as well as ion channels and targets derived from phytochemicals. Examples of promising novel candidates currently in clinical development for generalised anxiety disorder, social anxiety disorder, panic disorder, obsessive compulsive disorder or post-traumatic stress disorder include ketamine, riluzole, xenon with one common pharmacological action of modulation of glutamatergic neurotransmission, as well as the neurosteroid aloradine. Finally, compounds such as D-cycloserine, MDMA, L-DOPA and cannabinoids have shown efficacy in enhancing fear-extinction learning in humans. They are thus investigated in clinical trials as an augmentative strategy for speeding up and enhancing the long-term effectiveness of exposure-based psychotherapy, which could render chronic anxiolytic drug treatment dispensable for many patients. These efforts are indicative of a rekindled interest and renewed optimism in the anxiety drug discovery field, after decades of relative stagnation.