iMETHYL: an integrative database of human DNA methylation, gene expression, and genomic variation

Knowledge-guided learning methods for integrative analysis of multi-omics data

Integrative analysis of multi-omics data has the potential to yield valuable and comprehensive insights into the molecular mechanisms underlying complex diseases such as cancer and Alzheimer's disease. However, a number of analytical challenges complicate multi-omics data integration. For instance, -omics data are usually high-dimensional, and sample sizes in multi-omics studies tend to be modest. Furthermore, when genes in an important pathway have relatively weak signal, it can be difficult to detect them individually. There is a growing body of literature on knowledge-guided learning methods that can address these challenges by incorporating biological knowledge such as functional genomics and functional proteomics into multi-omics data analysis. These methods have been shown to outperform their counterparts that do not utilize biological knowledge in tasks including prediction, feature selection, clustering, and dimension reduction. In this review, we survey recently developed methods and applications of knowledge-guided multi-omics data integration methods and discuss future research directions.

Torch-eCpG: a fast and scalable eQTM mapper for thousands of molecular phenotypes with graphical processing units

BACKGROUND

Gene expression may be regulated by the DNA methylation of regulatory elements in cis, distal, and trans regions. One method to evaluate the relationship between DNA methylation and gene expression is the mapping of expression quantitative trait methylation (eQTM) loci (also called expression associated CpG loci, eCpG). However, no open-source tools are available to provide eQTM mapping. In addition, eQTM mapping can involve a large number of comparisons which may prevent the analyses due to limitations of computational resources. Here, we describe Torch-eCpG, an open-source tool to perform eQTM mapping that includes an optimized implementation that can use the graphical processing unit (GPU) to reduce runtime.

RESULTS

We demonstrate the analyses using the tool are reproducible, up to 18 × faster using the GPU, and scale linearly with increasing methylation loci.

CONCLUSIONS

Torch-eCpG is a fast, reliable, and scalable tool to perform eQTM mapping. Source code for Torch-eCpG is available at https://github.com/kordk/torch-ecpg .

jMorp: Japanese Multi-Omics Reference Panel update report 2023

Modern medicine is increasingly focused on personalized medicine, and multi-omics data is crucial in understanding biological phenomena and disease mechanisms. Each ethnic group has its unique genetic background with specific genomic variations influencing disease risk and drug response. Therefore, multi-omics data from specific ethnic populations are essential for the effective implementation of personalized medicine. Various prospective cohort studies, such as the UK Biobank, All of Us and Lifelines, have been conducted worldwide. The Tohoku Medical Megabank project was initiated after the Great East Japan Earthquake in 2011. It collects biological specimens and conducts genome and omics analyses to build a basis for personalized medicine. Summary statistical data from these analyses are available in the jMorp web database (https://jmorp.megabank.tohoku.ac.jp), which provides a multidimensional approach to the diversity of the Japanese population. jMorp was launched in 2015 as a public database for plasma metabolome and proteome analyses and has been continuously updated. The current update will significantly expand the scale of the data (metabolome, genome, transcriptome, and metagenome). In addition, the user interface and backend server implementations were rewritten to improve the connectivity between the items stored in jMorp. This paper provides an overview of the new version of the jMorp.

Unlocking the secrets: the power of methylation-based cfDNA detection of tissue damage in organ systems

BACKGROUND

Detecting organ and tissue damage is essential for early diagnosis, treatment decisions, and monitoring disease progression. Methylation-based assays offer a promising approach, as DNA methylation patterns can change in response to tissue damage. These assays have potential applications in early detection, monitoring disease progression, evaluating treatment efficacy, and assessing organ viability for transplantation. cfDNA released into the bloodstream upon tissue or organ injury can serve as a biomarker for damage. The epigenetic state of cfDNA, including DNA methylation patterns, can provide insights into the extent of tissue and organ damage.

CONTENT

Firstly, this review highlights DNA methylation as an extensively studied epigenetic modification that plays a pivotal role in processes such as cell growth, differentiation, and disease development. It then presents a variety of highly precise 5-mC methylation detection techniques that serve as powerful tools for gaining profound insights into epigenetic alterations linked with tissue damage. Subsequently, the review delves into the mechanisms underlying DNA methylation changes in organ and tissue damage, encompassing inflammation, oxidative stress, and DNA damage repair mechanisms. Next, it addresses the current research status of cfDNA methylation in the detection of specific organ tissues and organ damage. Finally, it provides an overview of the multiple steps involved in identifying specific methylation markers associated with tissue and organ damage for clinical trials. This review will explore the mechanisms and current state of research on cfDNA methylation-based assay detecting organ and tissue damage, the underlying mechanisms, and potential applications in clinical practice.

Lack of Methylation Changes in GJB2 and RB1 Non-coding Regions of Cochlear Implant Patients with Sensorineural Hearing Loss

Objective

Recent advances in epigenetic studies continue to reveal novel mechanisms of gene regulation and control, however little is known on the role of epigenetics in sensorineural hearing loss (SNHL) in humans. We aimed to investigate the methylation patterns of two regions, one in RB1 and another in GJB2 in Filipino patients with SNHL compared to hearing control individuals.

Methods

We investigated an RB1 promoter region that was previously identified as differentially methylated in children with SNHL and lead exposure. Additionally, we investigated a sequence in an enhancer-like region within GJB2 that contains four CpGs in close proximity. Bisulfite conversion was performed on salivary DNA samples from 15 children with SNHL and 45 unrelated ethnically-matched individuals. We then performed methylation-specific real-time PCR analysis (qMSP) using TaqMan® probes to determine percentage methylation of the two regions.

Results

Using qMSP, both our cases and controls had zero methylation at the targeted GJB2 and RB1 regions.

Conclusion

Our study showed no changes in methylation at the selected CpG regions in RB1 and GJB2 in the two comparison groups with or without SNHL. This may be due to a lack of environmental exposures to these target regions. Other epigenetic marks may be present around these regions as well as those of other HL-associated genes.

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.

Association between CORIN methylation and hypertension in Chinese adults

Background

Corin, a physical activator of atrial natriuretic peptide, has been associated with hypertension with unclear mechanisms. Here, we aimed to examine whether CORIN gene methylation was involved in the underlying molecular mechanisms.

Methods

DNA methylation levels of CORIN were measured by target bisulfite sequencing using genomic DNA isolated from peripheral blood mononuclear cells in 2498 participants in the Gusu cohort (discovery sample) and 1771 independent participants (replication sample). We constructed a mediation model with DNA methylation as the predictor, serum corin as the mediator, and hypertension as the outcome, adjusting for covariates. Multiple testing was controlled by false discovery rate (FDR) approach.

Results

Of the 9 CpGs assayed, hypermethylation at all CpGs were significantly associated with a lower level of blood pressure in the discovery sample and eight associations were also significant in the replication sample (all FDR-adjusted p<0.05). Serum corin mediated approximately 3.07% (p=0.004), 6.25% (p=0.002) and 10.11% (p=0.034) of the associations of hypermethylation at one CpG (Chr4:47840096) with systolic and diastolic blood pressure, and hypertension, respectively. All these mediations passed the causal inference test.

Conclusions

These results suggest that hypermethylation in the CORIN gene is associated with a lower odds of prevalent hypertension and may be involved in the role of corin in blood pressure regulation.

Epigenome-wide association study of plasma lipids in West Africans: the RODAM study

BACKGROUND

DNA-methylation has been associated with plasma lipid concentration in populations of diverse ethnic backgrounds, but epigenome-wide association studies (EWAS) in West-Africans are lacking. The aim of this study was to identify DNA-methylation loci associated with plasma lipids in Ghanaians.

METHODS

We conducted an EWAS using Illumina 450k DNA-methylation array profiles of extracted DNA from 663 Ghanaian participants. Differentially methylated positions (DMPs) were examined for association with plasma total cholesterol (TC), LDL-cholesterol, HDL-cholesterol, and triglycerides concentrations using linear regression models adjusted for age, sex, body mass index, diabetes mellitus, and technical covariates. Findings were replicated in independent cohorts of different ethnicities.

FINDINGS

We identified one significantly associated DMP with triglycerides (cg19693031 annotated to TXNIP, regression coefficient beta -0.26, false discovery rate adjusted p-value 0.001), which replicated in-silico in South African Batswana, African American, and European populations. From the top five DMPs with the lowest nominal p-values, two additional DMPs for triglycerides (CPT1A, ABCG1), two DMPs for LDL-cholesterol (EPSTI1, cg13781819), and one for TC (TXNIP) replicated. With the exception of EPSTI1, these loci are involved in lipid transport/metabolism or are known GWAS-associated loci. The top 5 DMPs per lipid trait explained 9.5% in the variance of TC, 8.3% in LDL-cholesterol, 6.1% in HDL-cholesterol, and 11.0% in triglycerides.

INTERPRETATION

The top DMPs identified in this study are in loci that play a role in lipid metabolism across populations, including West-Africans. Future studies including larger sample size, longitudinal study design and translational research is needed to increase our understanding on the epigenetic regulation of lipid metabolism among West-African populations.

FUNDING

European Commission under the Framework Programme (grant number: 278901).

CRISPR/Cas9 genome editing demonstrates functionality of the autoimmunity-associated SNP rs12946510

Genome-wide association studies (GWAS) map genetic associations of complex traits with precision limited to a linkage disequilibrium group. To translate GWAS results into new understanding of disease mechanisms, individual causative polymorphisms and their target genes should be identified. CRISPR/Cas9 genome editing can be used to create isogenic cell lines bearing alternative genotypes of candidate single-nucleotide polymorphisms to test their causality and to reveal gene targets. An intergenic polymorphism rs12946510 is associated with multiple sclerosis, inflammatory bowel disease and asthma. We created sublines of the T-helper cell line bearing alternative genotypes of rs12946510 and showed that its risk ("T") allele is associated with lower expression of IKZF3 and ORMDL3 genes and reduced cell activation. Our editing procedure can become an effective tool for discovering new genes involved in pathogenesis of complex diseases.

Epigenome-wide association study of BMI in Black populations from InterGEN and GENOA

OBJECTIVE

Obesity is a significant public health concern across the globe. Research investigating epigenetic mechanisms related to obesity and obesity-associated conditions has identified differences that may contribute to cellular dysregulation that accelerates the development of disease. However, few studies include Black women, who experience the highest incidence of obesity and early onset of cardiometabolic disorders.

METHODS

The association of BMI with epigenome-wide DNA methylation (DNAm) was examined using the 850K Illumina EPIC BeadChip in two Black populations (Intergenerational Impact of Genetic and Psychological Factors on Blood Pressure [InterGEN], n = 239; and The Genetic Epidemiology Network of Arteriopathy [GENOA] study, n = 961) using linear mixed-effects regression models adjusted for batch effects, cell type heterogeneity, population stratification, and confounding factors.

RESULTS

Cross-sectional analysis of the InterGEN discovery cohort identified 28 DNAm sites significantly associated with BMI, 24 of which had not been previously reported. Of these, 17 were replicated using the GENOA study. In addition, a meta-analysis, including both the InterGEN and GENOA cohorts, identified 658 DNAm sites associated with BMI with false discovery rate < 0.05. In a meta-analysis of Black women, we identified 628 DNAm sites significantly associated with BMI. Using a more stringent significance threshold of Bonferroni-corrected p value 0.05, 65 and 61 DNAm sites associated with BMI were identified from the combined sex and female-only meta-analyses, respectively.

CONCLUSIONS

This study suggests that BMI is associated with differences in DNAm among women that can be identified with DNA extracted from salivary (discovery) and peripheral blood (replication) samples among Black populations across two cohorts.

An explorative epigenome-wide association study of plasma renin and aldosterone concentration in a Ghanaian population: the RODAM study

BACKGROUND

The epigenetic regulation of the renin-angiotensin-aldosterone system (RAAS) potentially plays a role in the pathophysiology underlying the high burden of hypertension in sub-Saharan Africans (SSA). Here we report the first epigenome-wide association study (EWAS) of plasma renin and aldosterone concentrations and the aldosterone-to-renin ratio (ARR).

METHODS

Epigenome-wide DNA methylation was measured using the Illumina 450K array on whole blood samples of 68 Ghanaians. Differentially methylated positions (DMPs) were assessed for plasma renin concentration, aldosterone, and ARR using linear regression models adjusted for age, sex, body mass index, diabetes mellitus, hypertension, and technical covariates. Additionally, we extracted methylation loci previously associated with hypertension, kidney function, or that were annotated to RAAS-related genes and associated these with renin and aldosterone concentration.

RESULTS

We identified one DMP for renin, ten DMPs for aldosterone, and one DMP associated with ARR. Top DMPs were annotated to the PTPRN2, SKIL, and KCNT1 genes, which have been reported in relation to cardiometabolic risk factors, atherosclerosis, and sodium-potassium handling. Moreover, EWAS loci previously associated with hypertension, kidney function, or RAAS-related genes were also associated with renin, aldosterone, and ARR.

CONCLUSION

In this first EWAS on RAAS hormones, we identified DMPs associated with renin, aldosterone, and ARR in a SSA population. These findings are a first step in understanding the role of DNA methylation in regulation of the RAAS in general and in a SSA population specifically. Replication and translational studies are needed to establish the role of these DMPs in the hypertension burden in SSA populations.

SNP-by-CpG Site Interactions in ABCA7 Are Associated with Cognition in Older African Americans

SNPs in ABCA7 confer the largest genetic risk for Alzheimer's Disease (AD) in African Americans (AA) after APOE ε4. However, the relationship between ABCA7 and cognitive function has not been thoroughly examined. We investigated the effects of five known AD risk SNPs and 72 CpGs in ABCA7, as well as their interactions, on general cognitive function (cognition) in 634 older AA without dementia from Genetic Epidemiology Network of Arteriopathy (GENOA). Using linear mixed models, no SNP or CpG was associated with cognition after multiple testing correction, but five CpGs were nominally associated (p < 0.05). Four SNP-by-CpG interactions were associated with cognition (FDR q < 0.1). Contrast tests show that methylation is associated with cognition in some genotype groups (p < 0.05): a 1% increase at cg00135882 and cg22271697 is associated with a 0.68 SD decrease and 0.14 SD increase in cognition for those with the rs3764647 GG/AG (p = 0.004) and AA (p = 2 × 10-4) genotypes, respectively. In addition, a 1% increase at cg06169110 and cg17316918 is associated with a 0.37 SD decrease (p = 2 × 10-4) and 0.33 SD increase (p = 0.004), respectively, in cognition for those with the rs115550680 GG/AG genotype. While AD risk SNPs in ABCA7 were not associated with cognition in this sample, some have interactions with proximal methylation on cognition.

A critical review of datasets and computational suites for improving cancer theranostics and biomarker discovery

Cancer has been constantly evolving and so is the research pertaining to cancer diagnosis and therapeutic regimens. Early detection and specific therapeutics are the key features of modern cancer therapy. These requirements can only be fulfilled with the integration of diverse high-throughput technologies. Integration of advanced omics methodology involving genomics, epigenomics, proteomics, and transcriptomics provide a clear understanding of multi-faceted cancer. In the past few years, tremendous high-throughput data have been generated from cancer genomics and epigenomic analyses, which on further methodological analyses can yield better biological insights. The major epigenetic alterations reported in cancer are DNA methylation levels, histone post-translational modifications, and epi-miRNA regulating the oncogenes and tumor suppressor genes. While the genomic analyses like gene expression profiling, cancer gene prediction, and genome annotation divulge the genetic alterations in oncogenes or tumor suppressor genes. Also, systems biology approach using biological networks is being extensively used to identify novel cancer biomarkers. Therefore, integration of these multi-dimensional approaches will help to identify potential diagnostic and therapeutic biomarkers. Here, we reviewed the critical databases and tools dedicated to various epigenomic and genomic alterations in cancer. The review further focuses on the multi-omics resources available for further validating the identified cancer biomarkers. We also highlighted the tools for cancer biomarker discovery using a systems biology approach utilizing genomic and epigenomic data. Biomarkers predicted using such integrative approaches are shown to be more clinically relevant.

Using empirical biological knowledge to infer regulatory networks from multi-omics data

Background

Integration of multi-omics data can provide a more complex view of the biological system consisting of different interconnected molecular components, the crucial aspect for developing novel personalised therapeutic strategies for complex diseases. Various tools have been developed to integrate multi-omics data. However, an efficient multi-omics framework for regulatory network inference at the genome level that incorporates prior knowledge is still to emerge.

Results

We present IntOMICS, an efficient integrative framework based on Bayesian networks. IntOMICS systematically analyses gene expression, DNA methylation, copy number variation and biological prior knowledge to infer regulatory networks. IntOMICS complements the missing biological prior knowledge by so-called empirical biological knowledge, estimated from the available experimental data. Regulatory networks derived from IntOMICS provide deeper insights into the complex flow of genetic information on top of the increasing accuracy trend compared to a published algorithm designed exclusively for gene expression data. The ability to capture relevant crosstalks between multi-omics modalities is verified using known associations in microsatellite stable/instable colon cancer samples. Additionally, IntOMICS performance is compared with two algorithms for multi-omics regulatory network inference that can also incorporate prior knowledge in the inference framework. IntOMICS is also applied to detect potential predictive biomarkers in microsatellite stable stage III colon cancer samples.

Conclusions

We provide IntOMICS, a framework for multi-omics data integration using a novel approach to biological knowledge discovery. IntOMICS is a powerful resource for exploratory systems biology and can provide valuable insights into the complex mechanisms of biological processes that have a vital role in personalised medicine.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04891-9.

Association between CORIN methylation and hypertension in Chinese adults

BACKGROUND

Corin, a physical activator of atrial natriuretic peptide, has been associated with hypertension with unclear mechanisms. Here, we aimed to examine whether CORIN gene methylation was involved in the underlying molecular mechanisms.

METHODS

DNA methylation levels of CORIN were measured by target bisulfite sequencing using genomic DNA isolated from peripheral blood mononuclear cells in 2498 participants in the Gusu cohort (discovery sample) and 1771 independent participants (replication sample). We constructed a mediation model with DNA methylation as the predictor, serum corin as the mediator, and hypertension as the outcome, adjusting for covariates. Multiple testing was controlled by false discovery rate (FDR) approach.

RESULTS

Of the 9 CpGs assayed, hypermethylation at all CpGs were significantly associated with a lower level of blood pressure in the discovery sample and eight associations were also significant in the replication sample (all FDR-adjusted p<0.05). Serum corin mediated approximately 3.07% (p=0.004), 6.25% (p=0.002) and 10.11% (p=0.034) of the associations of hypermethylation at one CpG (Chr4:47840096) with systolic and diastolic blood pressure, and hypertension, respectively. All these mediations passed the causal inference test.

CONCLUSIONS

These results suggest that hypermethylation in the CORIN gene is associated with a lower odds of prevalent hypertension and may be involved in the role of corin in blood pressure regulation.

Genome wide exploration of the methylome in aggressive B-cell lymphoma in Golden Retrievers reveals a conserved hypermethylome

Few recurrent DNA mutations are seen in aggressive canine B cell lymphomas (cBCL), suggesting other frequent drivers. The methylated island recovery assay (MIRA-seq) or methylated CpG-binding domain sequencing (MBD-seq) was used to define the genome-wide methylation profiles in aggressive cBCL in Golden Retrievers to determine if cBCL can be better defined by epigenetic changes than by DNA mutations. DNA hypermethylation patterns were relatively homogenous within cBCL samples in Golden Retrievers, in different breeds and in geographical regions. Aberrant hypermethylation is thus suspected to be a central and early event in cBCL lymphomagenesis. Distinct subgroups within cBCL in Golden Retrievers were not identified with DNA methylation profiles. In comparison, the methylome profile of human DLBCL (hDLBCL) is relatively heterogeneous. Only moderate similarity between hDLBCL and cBCL was seen and cBCL likely cannot be accurately classified into the subtypes seen in hDLBCL. Genes with hypermethylated regions in the promoter-TSS-first exon of cBCL compared to normal B cells often also had additional hyper- and hypomethylated regions distributed throughout the gene suggesting non-randomized repeat targeting of key genes by epigenetic mechanisms. The prevalence of hypermethylation in transcription factor families in aggressive cBCL may represent a fundamental step in lymphomagenesis.

An epigenome-wide association study of insulin resistance in African Americans

Background

African Americans have a high risk for type 2 diabetes (T2D) and insulin resistance. Studies among other population groups have identified DNA methylation loci associated with insulin resistance, but data in African Americans are lacking. Using DNA methylation profiles of blood samples obtained from the Illumina Infinium® HumanMethylation450 BeadChip, we performed an epigenome-wide association study to identify DNA methylation loci associated with insulin resistance among 136 non-diabetic, unrelated African American men (mean age 41.6 years) from the Howard University Family Study.

Results

We identified three differentially methylated positions (DMPs) for homeostatic model assessment of insulin resistance (HOMA-IR) at 5% FDR. One DMP (cg14013695, HOXA5) is a known locus among Mexican Americans, while the other two DMPs are novel—cg00456326 (OSR1; beta = 0.027) and cg20259981 (ST18; beta = 0.010). Although the cg00456326 DMP is novel, the OSR1 gene has previously been found associated with both insulin resistance and T2D in Europeans. The genes HOXA5 and ST18 have been implicated in biological processes relevant to insulin resistance. Differential methylation at the significant HOXA5 and OSR1 DMPs is associated with differences in gene expression in the iMETHYL database. Analysis of differentially methylated regions (DMRs) did not identify any epigenome-wide DMRs for HOMA-IR. We tested transferability of HOMA-IR associated DMPs from five previous EWAS in Mexican Americans, Indian Asians, Europeans, and European ancestry Americans. Out of the 730 previously reported HOMA-IR DMPs, 47 (6.4%) were associated with HOMA-IR in this cohort of African Americans.

Conclusions

The findings from our study suggest substantial differences in DNA methylation patterns associated with insulin resistance across populations. Two of the DMPs we identified in African Americans have not been reported in other populations, and we found low transferability of HOMA-IR DMPs reported in other populations in African Americans. More work in African-ancestry populations is needed to confirm our findings as well as functional analyses to understand how such DNA methylation alterations contribute to T2D pathology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-022-01309-4.

DNA Methylation of the Natriuretic Peptide System Genes and Ischemic Stroke: Gene-Based and Gene Set Analyses

Background and Objectives

The natriuretic peptide (NP) system has been considered an important regulator for ischemic stroke (IS) with a limited clinical implication. A better understanding of the underlying molecular mechanisms is urgent. Here, we aimed to examine the role of DNA methylation of NP system genes in IS.

Methods

DNA methylation at promoter regions of 4 core NP system genes, e.g., CORIN, FURIN, NPPA, and NPPB, was measured by targeted bisulfite sequencing in 853 patients with IS and 918 controls. We first examined the association between DNA methylation at each single CpG and IS, followed by gene-based and gene set analyses to examine the joint associations of DNA methylation at multiple CpGs in a gene or all 4 genes as a pathway with IS.

Results

After control of covariates and multiple testing, DNA methylation at 19 of the 36 assayed CpGs was individually associated with IS at q < 0.05. Higher average methylation levels at the targeted regions of CORIN (odds ratio [OR] = 0.64, 95% confidence interval [CI]: 0.56–0.73), FURIN (OR = 0.78, 95% CI: 0.69–0.88), and NPPA (OR = 0.78, 95% CI: 0.69–0.88) were associated with a lower odds of IS (all q < 0.05). The truncated product method revealed the same gene-based associations (all q < 0.05) and found that DNA methylation at all 4 NP system genes together was jointly associated with IS (p = 0.0001).

Discussion

DNA methylation at NP system genes was downregulated in patients with IS. Our results may unravel a molecular mechanism underlying the regulating effect of the NP system on IS and highlight the relevance of testing the joint effect of multiple CpGs in the epigenetic analysis.

FURIN Promoter Methylation Predicts the Risk of Incident Diabetes: A Prospective Analysis in the Gusu Cohort

Background

Furin has been associated with diabetes but the underlying mechanisms are unclear. As a mediator linking fixed genome and dynamic environment, DNA methylation of its coding gene FURIN may be involved. Here, we aimed to examine the prospective association between DNA methylation in FURIN promoter and incident diabetes during 4 years of follow-up in Chinese adults.

Methods

DNA methylation levels in FURIN promoter were quantified by target bisulfite sequencing using peripheral blood from 1836 participants in the Gusu cohort who were free of diabetes at baseline. To examine the association between DNA methylation levels in FURIN promoter and incident diabetes, we constructed a logistic regression model adjusting for the conventional factors. Multiple testing was controlled by adjusting for the total number of CpG sites assayed using the false-discovery rate approach.

Results

Among the 1836 participants free of diabetes at baseline, 109 (5.94%) participants developed diabetes during the average of 4 years of follow-up. Hypermethylation at two of the eight CpG sites assayed in the FURIN promoter was associated with an increased risk of diabetes, after multivariable adjustment and multiple testing correction. Every 5% increment in methylation levels at CpG1 and CpG2 were associated with a 22% (OR=1.22, 95%CI: 1.05-1.43, P=0.009, q=0.038) and 39% (OR=1.39, 95%CI: 1.08-1.77, P=0.009, q=0.038) higher risk of incident diabetes, respectively. The gene-based association analysis revealed that DNA methylation at multiple CpG loci was jointly associated with incident diabetes (P<0.001). Using the average methylation level of the 8 CpG loci in FURIN promoter revealed a similar association (OR=1.28, 95% CI: 1.02-1.62, P=0.037).

Conclusions

These results suggested that the hypermethylation levels in FURIN promoter were associated with an increased risk for incident diabetes in Chinese adults.

A single-nucleotide-polymorphism in the 5′-flanking region of MSX1 gene as a predictive marker candidate for platinum-based therapy of esophageal carcinoma

Background:

Platinum derivatives are important treatment options for patients with esophageal carcinoma (EC), and a predictive marker for platinum-based therapy is needed for precision medicine.

Patients and methods:

This study contained two cohorts consisting of EC patients treated using platinum-based chemoradiation therapy (CRT) as the first-line and another external cohort of nationwide clinicogenomic data from the BioBank Japan (BBJ).

Results:

Genome-wide association study (GWAS) of therapeutic outcomes, refractory disease or not, following platinum-based CRT as first-line in 94 patients in the first cohort suggested the association of 89 SNPs using p < 0.0001. The top 10 SNPs selected from each chromosomal region by odds ratio were evaluated for progression-free survival (PFS) and overall survival (OS) hazard ratios in the first cohort, resulting in four candidates (p < 0.0025). The four selected candidates were re-evaluated in another cohort of 24 EC patients, which included patients prospectively enrolled in this study to fulfill the sample size statistically suggested by the results of the first cohort, and of the four, only rs3815544 was replicated (p < 0.0125). Furthermore, this candidate genotype of rs3815544 proceeded to the re-evaluation study in an external cohort consisting of EC patients treated with platinum derivatives and/or by radiation therapy as the first-line treatment in BBJ, which confirmed that the alternative allele (G) of rs3815544 was statistically associated with non-response (SD or PD) to platinum-based therapy in EC patients (odds ratio = 1.801, p = 0.048). The methylation QTL database as well as online clinicogenomic databases suggested that the region including rs3815544 may regulate MSX1 expression through CpG methylation, and this down-regulation was statistically associated with poor prognosis after platinum-based therapies for EC.

Conclusion:

rs3815544 is a novel candidate predictive marker for platinum-based EC therapy.

Low MICA Gene Expression Confers an Increased Risk of Graves' Disease: A Mendelian Randomization Study

Background: Expression of natural killer group 2 member D (NKG2D) ligand (NKG2DL) plays a major role as a "danger signal" on stressed cells to promote removal of the latter by NKG2D-expressing cytotoxic lymphocytes. NKG2DL expression has been found in peripheral immune cells as well, such as in macrophages; however, the effect of this expression is yet to be determined. Methods: We determined instrumental variables (IVs; R² <0.01 in linkage disequilibrium), explaining the major variance in major histocompatibility complex class I chain-related protein A (MICA) and B (MICB) gene expression levels from the expression-quantitative trait locus (eQTL) of NKG2DLs based on the RNA-seq analysis of peripheral blood mononuclear cells (PBMCs) from 381 Japanese. Simultaneously, the target outcomes were filtered by PheWAS from 58 health risks, using a community-based cohort study composed of 44,739 Japanese residents. Finally, we estimated the causal effect of gene expression levels on the outcomes using the Mendelian randomization approach. Results: We determined nine and four IVs, explaining 87.6% and 33.0% of MICA and MICB gene expression levels, respectively. In the association test, we identified 10 or 13 significant outcomes associated with the MICA or MICB eQTLs, respectively, as well as the causal effect of MICA expression on Graves' disease (GD) (p = 4.2 × 10^-3; odds ratio per 1 S.D. difference in the expression: 0.983 [confidence interval: 0.971-0.995]), using the weighted median estimator, without significant pleiotropy (p > 0.05), and the results were consistent across the sensitivity analyses. Conclusions: Our study provide novel evidence associating NKG2DL expression with GD, an autoimmune thyroiditis; direction of the effect indicated the immunoregulatory role of MICA expression in PBMCs, suggesting the importance of further functional assays in inflammatory diseases.

Parallel analysis of transcription, integration, and sequence of single HIV-1 proviruses

HIV-1-infected cells that persist despite antiretroviral therapy (ART) are frequently considered "transcriptionally silent," but active viral gene expression may occur in some cells, challenging the concept of viral latency. Applying an assay for profiling the transcriptional activity and the chromosomal locations of individual proviruses, we describe a global genomic and epigenetic map of transcriptionally active and silent proviral species and evaluate their longitudinal evolution in persons receiving suppressive ART. Using genome-wide epigenetic reference data, we show that proviral transcriptional activity is associated with activating epigenetic chromatin features in linear proximity of integration sites and in their inter- and intrachromosomal contact regions. Transcriptionally active proviruses were actively selected against during prolonged ART; however, this pattern was violated by large clones of virally infected cells that may outcompete negative selection forces through elevated intrinsic proliferative activity. Our results suggest that transcriptionally active proviruses are dynamically evolving under selection pressure by host factors.

Analysis of the DNA methylation pattern of the promoter region of calcitonin gene-related peptide 1 gene in patients with episodic migraine: An exploratory case-control study

Recent studies suggested that epigenetic mechanisms, including DNA methylation, may be involved in migraine pathogenesis. The calcitonin gene-related peptide (CGRP), encoded by calcitonin gene-related peptide 1 (CALCA) gene, plays a key role in the disease. The aim of the study was to evaluate DNA methylation of CALCA gene in patients with episodic migraine. 22 patients with episodic migraine (F/M 15/7, mean age 39.7 ± 13.4 years) and 20 controls (F/M 12/8, mean age 40.5 ± 14.8 years) were recruited. Genomic DNA was extracted from peripheral blood. Cytosine-to-thymine conversion was obtained with sodium bisulfite. The methylation pattern of two CpG islands in the promoter region of CALCA gene was analyzed. No difference of methylation of the 30 CpG sites at the distal region of CALCA promoter was observed between migraineurs and controls. Interestingly, in patients with episodic migraine the methylation level was lower in 2 CpG sites at the proximal promoter region (CpG -1461, p = 0.037, and -1415, p = 0.035, respectively). Furthermore, DNA methylation level at different CpG sites correlates with several clinical characteristics of the disease, as age at onset, presence of nausea/vomiting, depression and anxiety (p < 0.05). In conclusion, we found that DNA methylation profile in two CpG sites at the proximal promoter region of CALCA is lower in migraineurs when compared to controls. Intriguingly, the -1415 hypomethylated unit is located at the CREB binding site, a nuclear transcription factor. In addition, we found a correlation between the level of CALCA methylation and several clinical features of migraine. Further studies with larger sample size are needed to confirm these results.

DNA Methylation Abnormalities and Altered Whole Transcriptome Profiles after Switching from Combustible Tobacco Smoking to Heated Tobacco Products

BACKGROUND

The use of heated tobacco products (HTP) has increased exponentially in Japan since 2016; however, their effects on health remain a major concern.

METHODS

Tsuruoka Metabolome Cohort Study participants (n = 11,002) were grouped on the basis of their smoking habits as never smokers (NS), past smokers (PS), combustible tobacco smokers (CS), and HTP users for <2 years. Peripheral blood mononuclear cells were collected from 52 participants per group matched to HTP users using propensity scores, and DNA and RNA were purified from the samples. DNA methylation (DNAm) analysis of the 17 smoking-associated DNAm biomarker genes (such as AHRR, F2RL3, LRRN3, and GPR15), as well as whole transcriptome analysis, was performed.

RESULTS

Ten of the 17 genes were significantly hypomethylated in CS and HTP users compared with NS, among which AHRR, F2RL3, and RARA showed intermediate characteristics between CS and NS; nonetheless, AHRR expression was significantly higher in CS than in the other three groups. Conversely, LRRN3 and GPR15 were more hypomethylated in HTP users than in NS, and GPR15 expression was markedly upregulated in all the groups when compared with that in NS.

CONCLUSIONS

HTP users (switched from CS <2 years) display abnormal DNAm and transcriptome profiles, albeit to a lesser extent than the CS. However, because the molecular genetic effects of long-term HTP use are still unknown, long-term molecular epidemiologic studies are needed.

IMPACT

This study provides new insights into the molecular genetic effects on DNAm and transcriptome profiles in HTP users who switched from CS.

Longitudinal DNA methylation dynamics as a practical indicator in clinical epigenetics

BACKGROUND

One of the fundamental assumptions of DNA methylation in clinical epigenetics is that DNA methylation status can change over time with or without interplay with environmental and clinical conditions. However, little is known about how DNA methylation status changes over time under ordinary environmental and clinical conditions. In this study, we revisited the high frequency longitudinal DNA methylation data of two Japanese males (24 time-points within three months) and characterized the longitudinal dynamics.

RESULTS

The results showed that the majority of CpGs on Illumina HumanMethylation450 BeadChip probe set were longitudinally stable over the time period of three months. Focusing on dynamic and stable CpGs extracted from datasets, dynamic CpGs were more likely to be reported as epigenome-wide association study (EWAS) markers of various traits, especially those of immune- and inflammatory-related traits; meanwhile, the stable CpGs were enriched in metabolism-related genes and were less likely to be EWAS markers, indicating that the stable CpGs are stable both in the short-term within individuals and under various environmental and clinical conditions.

CONCLUSIONS

This study indicates that CpGs with different stabilities are involved in different functions and traits, and thus, they are potential indicators that can be applied for clinical epigenetic studies to outline underlying mechanisms.

Integration features of intact latent HIV-1 in CD4+ T cell clones contribute to viral persistence

Latent intact HIV-1 proviruses persist in a small subset of long-lived CD4+ T cells that can undergo clonal expansion in vivo. Expanded clones of CD4+ T cells dominate latent reservoirs in individuals on long-term antiretroviral therapy (ART) and represent a major barrier to HIV-1 cure. To determine how integration landscape might contribute to latency, we analyzed integration sites of near full length HIV-1 genomes from individuals on long-term ART, focusing on individuals whose reservoirs are highly clonal. We find that intact proviruses in expanded CD4+ T cell clones are preferentially integrated within Krüppel-associated box (KRAB) domain-containing zinc finger (ZNF) genes. ZNF genes are associated with heterochromatin in memory CD4+ T cells; nevertheless, they are expressed in these cells under steady-state conditions. In contrast to genes carrying unique integrations, ZNF genes carrying clonal intact integrations are down-regulated upon cellular activation. Together, the data suggest selected genomic sites, including ZNF genes, can be especially permissive for maintaining HIV-1 latency during memory CD4+ T cell expansion.

Toward a Framework for Assessing Privacy Risks in Multi-Omic Research and Databases

While the accumulation and increased circulation of genomic data have captured much attention over the past decade, privacy risks raised by the diversification and integration of omics have been largely overlooked. In this paper, we propose the outline of a framework for assessing privacy risks in multi-omic research and databases. Following a comparison of privacy risks associated with genomic and epigenomic data, we dissect ten privacy risk-impacting omic data properties that affect either the risk of re-identification of research participants, or the sensitivity of the information potentially conveyed by biological data. We then propose a three-step approach for the assessment of privacy risks in the multi-omic era. Thus, we lay grounds for a data property-based, 'pan-omic' approach that moves away from genetic exceptionalism. We conclude by inviting our peers to refine these theoretical foundations, put them to the test in their respective fields, and translate our approach into practical guidance.

Dynamic DNA Methylation Changes in the COMT Gene Promoter Region in Response to Mental Stress and Its Modulation by Transcranial Direct Current Stimulation

Changes in epigenetic modifications present a mechanism how environmental factors, such as the experience of stress, can alter gene regulation. While stress-related disorders have consistently been associated with differential DNA methylation, little is known about the time scale in which these alterations emerge. We investigated dynamic DNA methylation changes in whole blood of 42 healthy male individuals in response to a stressful cognitive task, its association with concentration changes in cortisol, and its modulation by transcranial direct current stimulation (tDCS). We observed a continuous increase in COMT promotor DNA methylation which correlated with higher saliva cortisol levels and was still detectable one week later. However, this lasting effect was suppressed by concurrent activity-enhancing anodal tDCS to the dorsolateral prefrontal cortex. Our findings support the significance of gene-specific DNA methylation in whole blood as potential biomarkers for stress-related effects. Moreover, they suggest alternative molecular mechanisms possibly involved in lasting behavioral effects of tDCS.

Transcriptomic Changes in Hot Spring Frog Tadpoles (Buergeria otai) in Response to Heat Stress

Buergeria frog tadpoles exhibit high thermal tolerance and are occasionally found in water pools that temporarily exceed 40°C. With the aim of understanding how they can cope with the severe heat stress, we performed RNA-seq of three heat-treated (38°C) and three control (25°C) tadpoles and compared their transcriptomic profiles. We identified 382 differentially expressed transcripts. A protein-protein interaction (PPI) network analysis of these transcripts further identified hub proteins involved in protein degradation, stress granule assembly, and global suppression of DNA transcription and mRNA translation. Along with the avoidance behavior against high water temperature, these endurance mechanisms potentially support tadpoles to survive in high temperatures for short periods of time. Similar mechanisms may exist in many other amphibian species whose habitats are prone to high temperatures.

Hypomethylation of AHRR (cg05575921) Is Related to Smoking Status in the Mexican Mestizo Population

Tobacco smoking results in a multifactorial disease involving environmental and genetic factors; epigenome-wide association studies (EWAS) show changes in DNA methylation levels due to cigarette consumption, partially reversible upon tobacco smoking cessation. Therefore, methylation levels could predict smoking status. This study aimed to evaluate the DNA methylation level of cg05575921 (AHRR) and cg23771366 (PRSS23) and their correlation with lung function variables, cigarette consumption, and nicotine addiction in the Mexican smoking population. We included 114 non-smokers (NS) and 102 current tobacco smokers (TS); we then further subclassified them as heavy smokers (HS) (n = 53) and light smokers (LS) (n = 49). We used restriction enzymes (MspI/HpaII) and qPCR to determine the DNA methylation level. We observed significant hypomethylation of cg05575921 in smokers compared to NS (p = 0.003); further analysis found a difference between HS and NS (p = 0.02). We did not observe differences between other groups or a positive correlation between methylation levels and age, BMI, cigarette consumption, nicotine addiction, or lung function. In conclusion, the cg05575921 site of AHRR is significantly hypomethylated in Mexican smokers, especially in HS (≥20 cigarettes per day).

Genome-wide DNA methylation analysis on C-reactive protein among Ghanaians suggests molecular links to the emerging risk of cardiovascular diseases

Molecular mechanisms at the intersection of inflammation and cardiovascular diseases (CVD) among Africans are still unknown. We performed an epigenome-wide association study to identify loci associated with serum C-reactive protein (marker of inflammation) among Ghanaians and further assessed whether differentially methylated positions (DMPs) were linked to CVD in previous reports, or to estimated CVD risk in the same population. We used the Illumina Infinium® HumanMethylation450 BeadChip to obtain DNAm profiles of blood samples in 589 Ghanaians from the RODAM study (without acute infections, not taking anti-inflammatory medications, CRP levels < 40 mg/L). We then used linear models to identify DMPs associated with CRP concentrations. Post-hoc, we evaluated associations of identified DMPs with elevated CVD risk estimated via ASCVD risk score. We also performed subset analyses at CRP levels ≤10 mg/L and replication analyses on candidate probes. Finally, we assessed for biological relevance of our findings in public databases. We subsequently identified 14 novel DMPs associated with CRP. In post-hoc evaluations, we found that DMPs in PC, BTG4 and PADI1 showed trends of associations with estimated CVD risk, we identified a separate DMP in MORC2 that was associated with CRP levels ≤10 mg/L, and we successfully replicated 65 (24%) of previously reported DMPs. All DMPs with gene annotations (13) were biologically linked to inflammation or CVD traits. We have identified epigenetic loci that may play a role in the intersection between inflammation and CVD among Ghanaians. Further studies among other Africans are needed to confirm our findings.

The association between RAPSN methylation in peripheral blood and breast cancer in the Chinese population

DNA methylation in peripheral blood is associated with breast cancer (BC) but has mainly been studied in Caucasian populations. We investigated the association between blood-based methylation of receptor-associated protein of the synapse (RAPSN) and BC in Chinese population. The methylation levels of 12 RAPSN CpG sites were quantitatively evaluated by mass spectrometry in two case-control studies with 283 sporadic BC cases and 331 controls totally. The association was analyzed by logistic regression adjusted for covariants. The RAPSN methylation levels in patients with variant clinical characteristics were investigated by non-parametric tests. We found a significant association between BC and altered RAPSN methylation in blood in women at premenopausal and perimenopausal (age < 50 years old), but not in the elder women. This was approved by two independent case-control studies as well as by combining the subjects of the two studies (taken all subjects together, age < 50 years old, per 5% of methylation, odds ratio (OR) range from 1.17 to 1.30 for two CpG sites; OR = 0.75 for one CpG site; all p values < 0.02). This age-related RAPSN methylation was further modified by human epidermal growth factor receptor 2 (HER2) status (age < 50 years old, HER2 negative, per 5% of methylation, OR range from 1.27 to 1.48 for two CpG sites; OR = 0.76 for one CpG site; all p values < 0.02). We elucidated an association between BC and blood-based RAPSN methylation influenced by age and the status of HER2 in Chinese population.

DNA methylation as the link between migration and the major noncommunicable diseases: the RODAM study

Aim: We assessed epigenome-wide DNA methylation (DNAm) differences between migrant and non-migrant Ghanaians. Materials & methods: We used the Illumina Infinium® HumanMethylation450 BeadChip to profile DNAm of 712 Ghanaians in whole blood. We used linear models to detect differentially methylated positions (DMPs) associated with migration. We performed multiple post hoc analyses to validate our findings. Results: We identified 13 DMPs associated with migration (delta-beta values: 0.2-4.5%). Seven DMPs in CPLX2, EIF4E3, MEF2D, TLX3, ST8SIA1, ANG and CHRM3 were independent of extrinsic genomic influences in public databases. Two DMPs in NLRC5 were associated with duration of stay in Europe among migrants. All DMPs were biologically linked to migration-related factors. Conclusion: Our findings provide the first insights into DNAm differences between migrants and non-migrants.

jMorp updates in 2020: large enhancement of multi-omics data resources on the general Japanese population

In the Tohoku Medical Megabank project, genome and omics analyses of participants in two cohort studies were performed. A part of the data is available at the Japanese Multi Omics Reference Panel (jMorp; https://jmorp.megabank.tohoku.ac.jp) as a web-based database, as reported in our previous manuscript published in Nucleic Acid Research in 2018. At that time, jMorp mainly consisted of metabolome data; however, now genome, methylome, and transcriptome data have been integrated in addition to the enhancement of the number of samples for the metabolome data. For genomic data, jMorp provides a Japanese reference sequence obtained using de novo assembly of sequences from three Japanese individuals and allele frequencies obtained using whole-genome sequencing of 8,380 Japanese individuals. In addition, the omics data include methylome and transcriptome data from ∼300 samples and distribution of concentrations of more than 755 metabolites obtained using high-throughput nuclear magnetic resonance and high-sensitivity mass spectrometry. In summary, jMorp now provides four different kinds of omics data (genome, methylome, transcriptome, and metabolome), with a user-friendly web interface. This will be a useful scientific data resource on the general population for the discovery of disease biomarkers and personalized disease prevention and early diagnosis.

MethHC 2.0: information repository of DNA methylation and gene expression in human cancer

DNA methylation is an important epigenetic regulator in gene expression and has several roles in cancer and disease progression. MethHC version 2.0 (MethHC 2.0) is an integrated and web-based resource focusing on the aberrant methylomes of human diseases, specifically cancer. This paper presents an updated implementation of MethHC 2.0 by incorporating additional DNA methylomes and transcriptomes from several public repositories, including 33 human cancers, over 50 118 microarray and RNA sequencing data from TCGA and GEO, and accumulating up to 3586 manually curated data from >7000 collected published literature with experimental evidence. MethHC 2.0 has also been equipped with enhanced data annotation functionality and a user-friendly web interface for data presentation, search, and visualization. Provided features include clinical-pathological data, mutation and copy number variation, multiplicity of information (gene regions, enhancer regions, and CGI regions), and circulating tumor DNA methylation profiles, available for research such as biomarker panel design, cancer comparison, diagnosis, prognosis, therapy study and identifying potential epigenetic biomarkers. MethHC 2.0 is now available at http://awi.cuhk.edu.cn/∼MethHC.

Evaluation of clinical formalin-fixed paraffin-embedded tissue quality for targeted-bisulfite sequencing

Formalin-fixed paraffin-embedded (FFPE) tissues are promising biological resources for genetic research. Recent improvements in DNA extraction from FFPE samples allowed the use of these tissues for multiple sequencing methods. However, fundamental research addressing the application of FFPE-derived DNA for targeted-bisulfite sequencing (TB-seq) is lacking. Here, we evaluated the suitability of FFPE-derived DNA for TB-seq. We conducted TB-seq using FFPE-derived DNA and corresponding fresh frozen (FF) tissues of patients with kidney cancer and compared the quality of DNA, libraries, and TB-seq statistics between the two preservation methods. The approximately 600-bp average fragment size of the FFPE-derived DNA was significantly shorter than that of the FF-derived DNA. The sequencing libraries constructed using FFPE-derived DNA and the mapping ratio were approximately 10 times and 10% lower, respectively, than those constructed using FF-derived DNA. In the mapped data of FFPE-derived DNA, duplicated reads accounted for > 60% of the obtained sequence reads, with lower mean on-target coverage. Therefore, the standard TB-seq protocol is inadequate for obtaining high-quality data for epigenetic analysis from FFPE-derived DNA, and technical improvements are necessary for enabling the use of archived FFPE resources.

Cell-cycle-gated feedback control mediates desensitization to interferon stimulation

Cells use molecular circuits to interpret and respond to extracellular cues, such as hormones and cytokines, which are often released in a temporally varying fashion. In this study, we combine microfluidics, time-lapse microscopy, and computational modeling to investigate how the type I interferon (IFN)-responsive regulatory network operates in single human cells to process repetitive IFN stimulation. We found that IFN-α pretreatments lead to opposite effects, priming versus desensitization, depending on input durations. These effects are governed by a regulatory network composed of a fast-acting positive feedback loop and a delayed negative feedback loop, mediated by upregulation of ubiquitin-specific peptidase 18 (USP18). We further revealed that USP18 upregulation can only be initiated at the G1/early S phases of cell cycle upon the treatment onset, resulting in heterogeneous and delayed induction kinetics in single cells. This cell cycle gating provides a temporal compartmentalization of feedback loops, enabling duration-dependent desensitization to repetitive stimulations.

Distinct viral reservoirs in individuals with spontaneous control of HIV-1

Sustained, drug-free control of HIV-1 replication is naturally achieved in less than 0.5% of infected individuals (here termed 'elite controllers'), despite the presence of a replication-competent viral reservoir1. Inducing such an ability to spontaneously maintain undetectable plasma viraemia is a major objective of HIV-1 cure research, but the characteristics of proviral reservoirs in elite controllers remain to be determined. Here, using next-generation sequencing of near-full-length single HIV-1 genomes and corresponding chromosomal integration sites, we show that the proviral reservoirs of elite controllers frequently consist of oligoclonal to near-monoclonal clusters of intact proviral sequences. In contrast to individuals treated with long-term antiretroviral therapy, intact proviral sequences from elite controllers were integrated at highly distinct sites in the human genome and were preferentially located in centromeric satellite DNA or in Krüppel-associated box domain-containing zinc finger genes on chromosome 19, both of which are associated with heterochromatin features. Moreover, the integration sites of intact proviral sequences from elite controllers showed an increased distance to transcriptional start sites and accessible chromatin of the host genome and were enriched in repressive chromatin marks. These data suggest that a distinct configuration of the proviral reservoir represents a structural correlate of natural viral control, and that the quality, rather than the quantity, of viral reservoirs can be an important distinguishing feature for a functional cure of HIV-1 infection. Moreover, in one elite controller, we were unable to detect intact proviral sequences despite analysing more than 1.5 billion peripheral blood mononuclear cells, which raises the possibility that a sterilizing cure of HIV-1 infection, which has previously been observed only following allogeneic haematopoietic stem cell transplantation2,3, may be feasible in rare instances.

Changes in DNA Methylation and Gene Expression of Insulin and Obesity-Related Gene PIK3R1 after Roux-en-Y Gastric Bypass

Weight regulation and the magnitude of weight loss after a Roux-en-Y gastric bypass (RYGB) can be genetically determined. DNA methylation patterns and the expression of some genes can be altered after weight loss interventions, including RYGB. The present study aimed to evaluate how the gene expression and DNA methylation of PIK3R1, an obesity and insulin-related gene, change after RYGB. Blood samples were obtained from 13 women (35.9 ± 9.2 years) with severe obesity before and six months after surgical procedure. Whole blood transcriptome and epigenomic patterns were assessed by microarray-based, genome-wide technologies. A total of 1966 differentially expressed genes were identified in the pre- and postoperative periods of RYGB. From these, we observed that genes involved in obesity and insulin pathways were upregulated after surgery. Then, the PIK3R1 gene was selected for further RT-qPCR analysis and cytosine-guanine nucleotide (CpG) sites methylation evaluation. We observed that the PI3KR1 gene was upregulated, and six DNA methylation CpG sites were differently methylated after bariatric surgery. In conclusion, we found that RYGB upregulates genes involved in obesity and insulin pathways.

Establishment of diagnostic criteria for upper urinary tract urothelial carcinoma based on genome-wide DNA methylation analysis

The aim of this study was to develop a less invasive and accurate diagnostic system for upper urinary tract urothelial carcinoma (UTUC) based on genome-wide DNA methylation profiling. Genome-wide DNA methylation screening was performed using the Infinium HumanMethylation450 BeadChip, and DNA methylation quantification was verified using pyrosequencing. We analysed 26 samples of normal control urothelial tissue (C), an initial cohort of 62 samples (31 samples of non-cancerous urothelium [N] from UTUC patients and 31 samples of the corresponding UTUCs), a validation cohort of 82 samples (41 N and 41 UTUC samples), and 14 samples of urinary bladder urothelial carcinoma (BUC). In the initial cohort, we identified 2,448 CpG sites showing significant differences in DNA methylation levels between both C and UTUC and N and UTUC, but not showing differences between C and N. Among these CpG sites, 10 were located within CpG islands or their shores and shelves included in genomic domains where DNA methylation levels are stably controlled, allowing discrimination of UTUC even from BUC. Receiver operating characteristic curve analysis for discrimination of UTUC from N in these 10 CpG and neighbouring sites (37 diagnostic panels in total) yielded area under the curve values of 0.959-1.000, with a sensitivity and specificity of 86.6-100% and 93.5-100%, respectively. The diagnostic impact was successfully confirmed in the validation cohort. Our criteria were useful for diagnosis of UTUC, regardless of its clinicopathological features. Application of our criteria to voided urine samples will ultimately allow non-invasive DNA methylation diagnosis of UTUC.

MethGET: web-based bioinformatics software for correlating genome-wide DNA methylation and gene expression

Background

DNA methylation is a major epigenetic modification involved in regulating gene expression. The effects of DNA methylation on gene expression differ by genomic location and vary across kingdoms, species and environmental conditions. To identify the functional regulatory roles of DNA methylation, the correlation between DNA methylation changes and alterations in gene expression is crucial. With the advance of next-generation sequencing, genome-wide methylation and gene expression profiling have become feasible. Current bioinformatics tools for investigating such correlation are designed to the assessment of DNA methylation at CG sites. The correlation of non-CG methylation and gene expression is very limited. Some bioinformatics databases allow correlation analysis, but they are limited to specific genomes such as that of humans and do not allow user-provided data.

Results

Here, we developed a bioinformatics web tool, MethGET (Methylation and Gene Expression Teller), that is specialized to analyse the association between genome-wide DNA methylation and gene expression. MethGET is the first web tool to which users can supply their own data from any genome. It is also the tool that correlates gene expression with CG, CHG, and CHH methylation based on whole-genome bisulfite sequencing data. MethGET not only reveals the correlation within an individual sample (single-methylome) but also performs comparisons between two groups of samples (multiple-methylomes). For single-methylome analyses, MethGET provides Pearson correlations and ordinal associations to investigate the relationship between DNA methylation and gene expression. It also groups genes with different gene expression levels to view the methylation distribution at specific genomic regions. Multiple-methylome analyses include comparative analyses and heatmap representations between two groups. These functions enable the detailed investigation of the role of DNA methylation in gene regulation. Additionally, we applied MethGET to rice regeneration data and discovered that CHH methylation in the gene body region may play a role in the tissue culture process, which demonstrates the capability of MethGET for use in epigenomic research.

Conclusions

MethGET is a Python software that correlates DNA methylation and gene expression. Its web interface is publicly available at https://paoyang.ipmb.sinica.edu.tw/Software.html. The stand-alone version and source codes are available on GitHub at https://github.com/Jason-Teng/MethGET.

An epigenome-wide association study of posttraumatic stress disorder in US veterans implicates several new DNA methylation loci

BACKGROUND

Previous studies using candidate gene and genome-wide approaches have identified epigenetic changes in DNA methylation (DNAm) associated with posttraumatic stress disorder (PTSD).

METHODS

In this study, we performed an EWAS of PTSD in a cohort of Veterans (n = 378 lifetime PTSD cases and 135 controls) from the Translational Research Center for TBI and Stress Disorders (TRACTS) cohort assessed using the Illumina EPIC Methylation BeadChip which assesses DNAm at more than 850,000 sites throughout the genome. Our model included covariates for ancestry, cell heterogeneity, sex, age, and a smoking score based on DNAm at 39 smoking-associated CpGs. We also examined in EPIC-based DNAm data generated from pre-frontal cortex (PFC) tissue from the National PTSD Brain Bank (n = 72).

RESULTS

The analysis of blood samples yielded one genome-wide significant association with PTSD at cg19534438 in the gene G0S2 (p = 1.19 × 10-7, padj = 0.048). This association was replicated in an independent PGC-PTSD-EWAS consortium meta-analysis of military cohorts (p = 0.0024). We also observed association with the smoking-related locus cg05575921 in AHRR despite inclusion of a methylation-based smoking score covariate (p = 9.16 × 10-6), which replicates a previously observed PGC-PTSD-EWAS association (Smith et al. 2019), and yields evidence consistent with a smoking-independent effect. The top 100 EWAS loci were then examined in the PFC data. One of the blood-based PTSD loci, cg04130728 in CHST11, which was in the top 10 loci in blood, but which was not genome-wide significant, was significantly associated with PTSD in brain tissue (in blood p = 1.19 × 10-5, padj = 0.60, in brain, p = 0.00032 with the same direction of effect). Gene set enrichment analysis of the top 500 EWAS loci yielded several significant overlapping GO terms involved in pathogen response, including "Response to lipopolysaccharide" (p = 6.97 × 10-6, padj = 0.042).

CONCLUSIONS

The cross replication observed in independent cohorts is evidence that DNA methylation in peripheral tissue can yield consistent and replicable PTSD associations, and our results also suggest that that some PTSD associations observed in peripheral tissue may mirror associations in the brain.

Mechanism of long noncoding RNAs as transcriptional regulators in cancer

Dysregulation of gene expression, often interpreted by gene transcription as an endpoint response, is tightly associated with human cancer. Long noncoding RNAs (lncRNAs), derived from the noncoding elements in the genome and appeared no less than 200nt in length, have emerged as a novel class of pivotal regulatory component. Recently, great attention has been paid to the cancer-related lncRNAs and growing evidence have shown that lncRNAs act as key transcriptional regulators in cancer cells through diverse mechanisms. Here, we focus on the nucleus-expressed lncRNAs and summarize their molecular mechanisms in transcriptional control during tumorigenesis and cancer metastasis. Six major mechanisms will be discussed in this review: association with transcriptional factor, modulating DNA methylation or histone modification enzyme, influencing on chromatin remodelling complex, facilitating chromosomal looping, interaction with RNA polymerase and direct association with promoter.

The association of integration patterns of human papilloma virus and single nucleotide polymorphisms on immune- or DNA repair-related genes in cervical cancer patients

The present study investigated the association between single nucleotide polymorphisms (SNPs) in immune- or DNA repair-related genes and the integration pattern of human papillomavirus (HPV), a promising prognostic marker in cervical cancer. The HPV integration patterns of cervical cancer patients were determined by polymerase chain reaction and in situ hybridization, and categorized as episomal (group A), single-copy or multi-copy tandem repetition integrated (group B), and undetectable HPV types (group C). After sample and SNP quality control, 166,505 SNPs in 161 samples (38, 111, and 12 patients in groups A, B, and C, respectively) were examined. None of the SNPs reached genome-wide significance, and several candidate SNPs for future study were selected, including rs10999435 on chromosome 10q22, rs1322054 on chromosome 9q32-33, and rs10902171 on chromosome 11p15. Luciferase assay identified rs1322054 as the primary functional variant to regulate gene expression in immune cell. Further studies are needed to determine the genetic background of different integration patterns of HPV in cervical cancer patients.

Genome analyses for the Tohoku Medical Megabank Project towards establishment of personalized healthcare

Personalized healthcare (PHC) based on an individual's genetic make-up is one of the most advanced, yet feasible, forms of medical care. The Tohoku Medical Megabank (TMM) Project aims to combine population genomics, medical genetics and prospective cohort studies to develop a critical infrastructure for the establishment of PHC. To date, a TMM CommCohort (adult general population) and a TMM BirThree Cohort (birth+three-generation families) have conducted recruitments and baseline surveys. Genome analyses as part of the TMM Project will aid in the development of a high-fidelity whole-genome Japanese reference panel, in designing custom single-nucleotide polymorphism (SNP) arrays specific to Japanese, and in estimation of the biological significance of genetic variations through linked investigations of the cohorts. Whole-genome sequencing from >3,500 unrelated Japanese and establishment of a Japanese reference genome sequence from long-read data have been done. We next aim to obtain genotype data for all TMM cohort participants (>150,000) using our custom SNP arrays. These data will help identify disease-associated genomic signatures in the Japanese population, while genomic data from TMM BirThree Cohort participants will be used to improve the reference genome panel. Follow-up of the cohort participants will allow us to test the genetic markers and, consequently, contribute to the realization of PHC.

Genome-wide analysis of polymorphism × sodium interaction effect on blood pressure identifies a novel 3'-BCL11B gene desert locus

Excessive sodium intake is a global risk factor for hypertension. Sodium effects on blood pressure vary from person to person; hence, high-risk group targeting based on personal genetic information can play a complementary role to ongoing population preventive approaches to reduce sodium consumption. To identify genetic factors that modulate sodium effects on blood pressure, we conducted a population-based genome-wide interaction analysis in 8,768 Japanese subjects, which was >3 times larger than a similar previous study. We tested 7,135,436 polymorphisms in the discovery cohort, and loci that met suggestive significance were further examined in an independent replication cohort. We found that an interaction between a novel 3'-BCL11B gene desert locus and daily sodium consumption was significantly associated with systolic blood pressure in both discovery and replication cohorts under the recessive model. Further statistical analysis of rs8022678, the sentinel variant of the 3'-BCL11B gene desert locus, showed that differences in mean systolic blood pressure between high and low sodium consumption subgroups were 5.9 mm Hg (P = 8.8 × 10-12) in rs8022678 A carriers and -0.3 mm Hg (P = 0.27) in rs8022678 A non-carriers, suggesting that the rs8022678 genotype can classify persons into sodium-sensitive (A carriers) and sodium-insensitive (A non-carriers) subgroups. Our results implied that rs8022678 A carriers may receive a greater benefit from sodium-lowering interventions than non-carriers.