EWAS Open Platform: integrated data, knowledge and toolkit for epigenome-wide association study

A multi-trait epigenome-wide association study identified DNA methylation signature of inflammation among people with HIV

Inflammation underlies many conditions causing excess morbidity and mortality among people with HIV (PWH). A handful of single-trait epigenome-wide association studies (EWAS) have suggested that inflammation is associated with DNA methylation (DNAm) among PWH. Multi-trait EWAS may further improve statistical power and reveal pathways in common between different inflammatory markers. We conducted single-trait EWAS of three inflammatory markers (soluble CD14, D-dimers, and interleukin 6) in the Veteran Aging Cohort Study (n = 920). The study population was all male PWH with an average age of 51 years, and 82.3% self-reported as Black. We then applied two multi-trait EWAS methods-CPASSOC and OmniTest-to combine single-trait EWAS results. CPASSOC and OmniTest identified 189 and 157 inflammation-associated DNAm sites respectively, of which 112 overlapped. Among the identified sites, 56% were not significant in any single-trait EWAS. Top sites were mapped to inflammation-related genes including IFITM1 , PARP9 and STAT1 . These genes were significantly enriched in pathways such as "type I interferon signaling" and "immune response to virus". We demonstrate that multi-trait EWAS can improve the discovery of inflammation-associated DNAm sites, genes, and pathways. These DNAm sites suggest molecular mechanisms in response to inflammation associated with HIV and might hold the key to addressing persistent inflammation in PWH.

Epigenomic insights into common human disease pathology

The epigenome-the chemical modifications and chromatin-related packaging of the genome-enables the same genetic template to be activated or repressed in different cellular settings. This multi-layered mechanism facilitates cell-type specific function by setting the local sequence and 3D interactive activity level. Gene transcription is further modulated through the interplay with transcription factors and co-regulators. The human body requires this epigenomic apparatus to be precisely installed throughout development and then adequately maintained during the lifespan. The causal role of the epigenome in human pathology, beyond imprinting disorders and specific tumour suppressor genes, was further brought into the spotlight by large-scale sequencing projects identifying that mutations in epigenomic machinery genes could be critical drivers in both cancer and developmental disorders. Abrogation of this cellular mechanism is providing new molecular insights into pathogenesis. However, deciphering the full breadth and implications of these epigenomic changes remains challenging. Knowledge is accruing regarding disease mechanisms and clinical biomarkers, through pathogenically relevant and surrogate tissue analyses, respectively. Advances include consortia generated cell-type specific reference epigenomes, high-throughput DNA methylome association studies, as well as insights into ageing-related diseases from biological 'clocks' constructed by machine learning algorithms. Also, 3rd-generation sequencing is beginning to disentangle the complexity of genetic and DNA modification haplotypes. Cell-free DNA methylation as a cancer biomarker has clear clinical utility and further potential to assess organ damage across many disorders. Finally, molecular understanding of disease aetiology brings with it the opportunity for exact therapeutic alteration of the epigenome through CRISPR-activation or inhibition.

Epigenetic modifications in the development of bronchopulmonary dysplasia: a review

While perinatal medicine advancements have bolstered survival outcomes for premature infants, bronchopulmonary dysplasia (BPD) continues to threaten their long-term health. Gene-environment interactions, mediated by epigenetic modifications such as DNA methylation, histone modification, and non-coding RNA regulation, take center stage in BPD pathogenesis. Recent discoveries link methylation variations across biological pathways with BPD. Also, the potential reversibility of histone modifications fuels new treatment avenues. The review also highlights the promise of utilizing mesenchymal stem cells and their exosomes as BPD therapies, given their ability to modulate non-coding RNA, opening novel research and intervention possibilities. IMPACT: The complexity and universality of epigenetic modifications in the occurrence and development of bronchopulmonary dysplasia were thoroughly discussed. Both molecular and cellular mechanisms contribute to the diverse nature of epigenetic changes, suggesting the need for deeper biochemical techniques to explore these molecular alterations. The utilization of innovative cell-specific drug delivery methods like exosomes and extracellular vesicles holds promise in achieving precise epigenetic regulation.

Methylome-wide and meQTL analysis helps to distinguish treatment response from non-response and pathogenesis markers in schizophrenia

Schizophrenia is a complex condition with entwined genetic and epigenetic risk factors, posing a challenge to disentangle the intermixed pathological and therapeutic epigenetic signatures. To resolve this, we performed 850K methylome-wide and 700K genome-wide studies on the same set of schizophrenia patients by stratifying them into responders, non-responders, and drug-naïve patients. The key genes that signified the response were followed up using real-time gene expression studies to understand the effect of antipsychotics at the gene transcription level. The study primarily implicates hypermethylation in therapeutic response and hypomethylation in the drug-non-responsive state. Several differentially methylated sites and regions colocalized with the schizophrenia genome-wide association study (GWAS) risk genes and variants, supporting the convoluted gene-environment association. Gene ontology and protein-protein interaction (PPI) network analyses revealed distinct patterns that differentiated the treatment response from drug resistance. The study highlights the strong involvement of several processes related to nervous system development, cell adhesion, and signaling in the antipsychotic response. The ability of antipsychotic medications to alter the pathology by modulating gene expression or methylation patterns is evident from the general increase in the gene expression of response markers and histone modifiers and the decrease in class II human leukocyte antigen (HLA) genes following treatment with varying concentrations of medications like clozapine, olanzapine, risperidone, and haloperidol. The study indicates a directional overlap of methylation markers between pathogenesis and therapeutic response, thereby suggesting a careful distinction of methylation markers of pathogenesis from treatment response. In addition, there is a need to understand the trade-off between genetic and epigenetic observations. It is suggested that methylomic changes brought about by drugs need careful evaluation for their positive effects on pathogenesis, course of disease progression, symptom severity, side effects, and refractoriness.

Gbdmr: identifying differentially methylated CpG regions in the human genome via generalized beta regressions

BACKGROUND

DNA methylation is a biochemical process in which a methyl group is added to the cytosine-phosphate-guanine (CpG) site on DNA molecules without altering the DNA sequence. Multiple CpG sites in a certain genome region can be differentially methylated across phenotypes. Identifying these differentially methylated CpG regions (DMRs) associated with the phenotypes contributes to disease prediction and precision medicine development.

RESULTS

We propose a novel DMR detection algorithm, gbdmr. In contrast to existing methods under a linear regression framework, gbdmr assumes that DNA methylation levels follow a generalized beta distribution. We compare gbdmr to alternative approaches via simulations and real data analyses, including dmrff, a new DMR detection approach that shows promising performance among competitors, and the traditional EWAS that focuses on single CpG sites. Our simulations demonstrate that gbdmr is superior to the other two when the correlation between neighboring CpG sites is strong, while dmrff shows a higher power when the correlation is weak. We provide an explanation of these phenomena from a theoretical perspective. We further applied the three methods to multiple real DNA methylation datasets. One is from a birth cohort study undertaken on the Isle of Wight, United Kingdom, and the other two are from the Gene Expression Omnibus database repository. Overall, gbdmr identifies more DMR CpGs linked to phenotypes than dmrff, and the simulated results support the findings.

CONCLUSIONS

Gbdmr is an innovative method for detecting DMRs based on generalized beta regression. It demonstrated notable advantages over dmrff and traditional EWAS, particularly when adjacent CpGs exhibited moderate to strong correlations. Our real data analyses and simulated findings highlight the reliability of gbdmr as a robust DMR detection tool. The gbdmr approach is accessible and implemented by R on GitHub: https://github.com/chengzhouwu/gbdmr .

The association between DNA methylation and human height and a prospective model of DNA methylation-based height prediction

As a vital anthropometric characteristic, human height information not only helps to understand overall developmental status and genetic risk factors, but is also important for forensic DNA phenotyping. We utilized linear regression analysis to test the association between each CpG probe and the height phenotype. Next, we designed a methylation sequencing panel targeting 959 CpGs and subsequent height inference models were constructed for the Chinese population. A total of 11,730 height-associated sites were identified. By employing KPCA and deep neural networks, a prediction model was developed, of which the cross-validation RMSE, MAE and R2 were 5.62 cm, 4.45 cm and 0.64, respectively. Genetic factors could explain 39.4% of the methylation level variance of sites used in the height inference models. Collectively, we demonstrated an association between height and DNA methylation status through an EWAS analysis. Targeted methylation sequencing of only 959 CpGs combined with deep learning techniques could provide a model to estimate human height with higher accuracy than SNP-based prediction models.

Integrative multi-omics analysis of genomic, epigenomic, and metabolomics data leads to new insights for Attention-Deficit/Hyperactivity Disorder

The evolving field of multi-omics combines data and provides methods for simultaneous analysis across several omics levels. Here, we integrated genomics (transmitted and non-transmitted polygenic scores [PGSs]), epigenomics, and metabolomics data in a multi-omics framework to identify biomarkers for Attention-Deficit/Hyperactivity Disorder (ADHD) and investigated the connections among the three omics levels. We first trained single- and next multi-omics models to differentiate between cases and controls in 596 twins (cases = 14.8%) from the Netherlands Twin Register (NTR) demonstrating reasonable in-sample prediction through cross-validation. The multi-omics model selected 30 PGSs, 143 CpGs, and 90 metabolites. We confirmed previous associations of ADHD with glucocorticoid exposure and the transmembrane protein family TMEM, show that the DNA methylation of the MAD1L1 gene associated with ADHD has a relation with parental smoking behavior, and present novel findings including associations between indirect genetic effects and CpGs of the STAP2 gene. However, out-of-sample prediction in NTR participants (N = 258, cases = 14.3%) and in a clinical sample (N = 145, cases = 51%) did not perform well (range misclassification was [0.40, 0.57]). The results highlighted connections between omics levels, with the strongest connections between non-transmitted PGSs, CpGs, and amino acid levels and show that multi-omics designs considering interrelated omics levels can help unravel the complex biology underlying ADHD.

Improvements in lung function following vitamin C supplementation to pregnant smokers are associated with buccal DNA methylation at 5 years of age

BACKGROUND

We previously reported in the "Vitamin C to Decrease the Effects of Smoking in Pregnancy on Infant Lung Function" randomized clinical trial (RCT) that vitamin C (500 mg/day) supplementation to pregnant smokers is associated with improved respiratory outcomes that persist through 5 years of age. The objective of this study was to assess whether buccal cell DNA methylation (DNAm), as a surrogate for airway epithelium, is associated with vitamin C supplementation, improved lung function, and decreased occurrence of wheeze.

METHODS

We conducted epigenome-wide association studies (EWAS) using Infinium MethylationEPIC arrays and buccal DNAm from 158 subjects (80 placebo; 78 vitamin C) with pulmonary function testing (PFT) performed at the 5-year visit. EWAS were performed on (1) vitamin C treatment, (2) forced expiratory flow between 25 and 75% of expired volume (FEF25-75), and (3) offspring wheeze. Models were adjusted for sex, race, study site, gestational age at randomization (≤ OR > 18 weeks), proportion of epithelial cells, and latent covariates in addition to child length at PFT in EWAS for FEF25-75. We considered FDR p < 0.05 as genome-wide significant and nominal p < 0.001 as candidates for downstream analyses. Buccal DNAm measured in a subset of subjects at birth and near 1 year of age was used to determine whether DNAm signatures originated in utero, or emerged with age.

RESULTS

Vitamin C treatment was associated with 457 FDR significant (q < 0.05) differentially methylated CpGs (DMCs; 236 hypermethylated; 221 hypomethylated) and 53 differentially methylated regions (DMRs; 26 hyper; 27 hypo) at 5 years of age. FEF25-75 was associated with one FDR significant DMC (cg05814800), 1,468 candidate DMCs (p < 0.001), and 44 DMRs. Current wheeze was associated with 0 FDR-DMCs, 782 candidate DMCs, and 19 DMRs (p < 0.001). In 365/457 vitamin C FDR significant DMCs at 5 years of age, there was no significant interaction between time and treatment.

CONCLUSIONS

Vitamin C supplementation to pregnant smokers is associated with buccal DNA methylation in offspring at 5 years of age, and most methylation signatures appear to be persistent from the prenatal period. Buccal methylation at 5 years was also associated with current lung function and occurrence of wheeze, and these functionally associated loci are enriched for vitamin C associated loci. Clinical trial registration ClinicalTrials.gov, NCT01723696 and NCT03203603.

DNA-methylome-derived epigenetic fingerprint as an immunophenotype indicator of durable clinical immunotherapeutic benefits in head and neck squamous cell carcinoma

BACKGROUND

Cancer immunotherapy provides durable response and improves survival in a subset of head and neck squamous cell carcinoma (HNSC) patients, which may due to discriminative tumor microenvironment (TME). Epigenetic regulations play critical roles in HNSC tumorigenesis, progression, and activation of functional immune cells. This study aims to identify an epigenetic signature as an immunophenotype indicator of durable clinical immunotherapeutic benefits in HNSC patients.

METHODS

Unsupervised consensus clustering approach was applied to distinguish immunophenotypes based on five immune signatures in The Cancer Genome Atlas (TCGA) HNSC cohort. Two immunophenotypes (immune 'Hot' and immune 'Cold') that had different TME features, diverse prognosis, and distinct DNA methylation patterns were recognized. Immunophenotype-related methylated signatures (IPMS) were identified by the least absolute shrinkage and selector operation algorithm. Additionally, the IPMS score by deconvolution algorithm was constructed as an immunophenotype classifier to predict clinical outcomes and immunotherapeutic response.

RESULTS

The 'Hot' HNSC immunophenotype had higher immunoactivity and better overall survival (p = 0.00055) compared to the 'Cold' tumors. The immunophenotypes had distinct DNA methylation patterns, which was closely associated with HNSC tumorigenesis and functional immune cell infiltration. 311 immunophenotype-related methylated CpG sites (IRMCs) was identified from TCGA-HNSC dataset. IPMS score model achieved a strong clinical predictive performance for classifying immunophenotypes. The area under the curve value (AUC) of the IPMS score model reached 85.9% and 89.8% in TCGA train and test datasets, respectively, and robustness was verified in five HNSC validation datasets. It was also validated as an immunophenotype classifier for predicting durable clinical benefits (DCB) in lung cancer patients who received anti-PD-1/PD-L1 immunotherapy (p = 0.017) and TCGA-SKCM patients who received distinct immunotherapy (p = 0.033).

CONCLUSIONS

This study systematically analyzed DNA methylation patterns in distinct immunophenotypes to identify IPMS with clinical prognostic potential for personalized epigenetic anticancer approaches in HNSC patients. The IPMS score model may serve as a reliable epigenome prognostic tool for clinical immunophenotyping to guide immunotherapeutic strategies in HNSC.

Novel markers and networks related to restored skeletal muscle transcriptome after bariatric surgery

OBJECTIVE

The aim of this study was to discover novel markers underlying the improvement of skeletal muscle metabolism after bariatric surgery.

METHODS

Skeletal muscle transcriptome data of lean people and people with obesity, before and 1 year after bariatric surgery, were subjected to weighted gene co-expression network analysis (WGCNA) and least absolute shrinkage and selection operator (LASSO) regression. Results of LASSO were confirmed in a replication cohort.

RESULTS

The expression levels of 440 genes differing between individuals with and without obesity were no longer different 1 year after surgery, indicating restoration. WGCNA clustered 116 genes with normalized expression in one major module, particularly correlating to weight loss and decreased plasma free fatty acids (FFA), 44 of which showed an obesity-related phenotype upon deletion in mice. Among the genes of the major module, 105 represented prominent markers for reduced FFA concentration, including 55 marker genes for decreased BMI in both the discovery and replication cohorts.

CONCLUSIONS

Previously unknown gene networks and marker genes underlined the important role of FFA in restoring muscle gene expression after bariatric surgery and further suggest novel therapeutic targets for obesity.

Comprehensive evaluation of smoking exposures and their interactions on DNA methylation

BACKGROUND

Smoking impacts DNA methylation, but data are lacking on smoking-related differential methylation by sex or dietary intake, recent smoking cessation (<1 year), persistence of differential methylation from in utero smoking exposure, and effects of environmental tobacco smoke (ETS).

METHODS

We meta-analysed data from up to 15,014 adults across 5 cohorts with DNA methylation measured in blood using Illumina's EPIC array for current smoking (2560 exposed), quit < 1 year (500 exposed), in utero (286 exposed), and ETS exposure (676 exposed). We also evaluated the interaction of current smoking with sex or diet (fibre, folate, and vitamin C).

FINDINGS

Using false discovery rate (FDR < 0.05), 65,857 CpGs were differentially methylated in relation to current smoking, 4025 with recent quitting, 594 with in utero exposure, and 6 with ETS. Most current smoking CpGs attenuated within a year of quitting. CpGs related to in utero exposure in adults were enriched for those previously observed in newborns. Differential methylation by current smoking at 4-71 CpGs may be modified by sex or dietary intake. Nearly half (35-50%) of differentially methylated CpGs on the 450 K array were associated with blood gene expression. Current smoking and in utero smoking CpGs implicated 3049 and 1067 druggable targets, including chemotherapy drugs.

INTERPRETATION

Many smoking-related methylation sites were identified with Illumina's EPIC array. Most signals revert to levels observed in never smokers within a year of cessation. Many in utero smoking CpGs persist into adulthood. Smoking-related druggable targets may provide insights into cancer treatment response and shared mechanisms across smoking-related diseases.

FUNDING

Intramural Research Program of the National Institutes of Health, Norwegian Ministry of Health and Care Services and the Ministry of Education and Research, Chief Scientist Office of the Scottish Government Health Directorates and the Scottish Funding Council, Medical Research Council UK and the Wellcome Trust.

Blood epigenome-wide association studies of suicide attempt in adults with bipolar disorder

Suicide attempt (SA) risk is elevated in individuals with bipolar disorder (BD), and DNA methylation patterns may serve as possible biomarkers of SA. We conducted epigenome-wide association studies (EWAS) of blood DNA methylation associated with BD and SA. DNA methylation was measured at >700,000 positions in a discovery cohort of n = 84 adults with BD with a history of SA (BD/SA), n = 79 adults with BD without history of SA (BD/non-SA), and n = 76 non-psychiatric controls (CON). EWAS revealed six differentially methylated positions (DMPs) and seven differentially methylated regions (DMRs) between BD/SA and BD/non-SA, with multiple immune-related genes implicated. There were no epigenome-wide significant differences when BD/SA and BD/non-SA were each compared to CON, and patterns suggested that epigenetics differentiating BD/SA from BD/non-SA do not differentiate BD/non-SA from CON. Weighted gene co-methylation network analysis and trait enrichment analysis of the BD/SA vs. BD/non-SA contrast further corroborated immune system involvement, while gene ontology analysis implicated calcium signalling. In an independent replication cohort of n = 48 BD/SA and n = 47 BD/non-SA, fold changes at the discovery cohort's significant sites showed moderate correlation across cohorts and agreement on direction. In both cohorts, classification accuracy for SA history among individuals with BD was highest when methylation at the significant CpG sites as well as information from clinical interviews were combined, with an AUC of 88.8% (CI = 83.8-93.8%) and 82.1% (CI = 73.6-90.5%) for the combined epigenetic-clinical classifier in the discovery and replication cohorts, respectively. Our results provide novel insight to the role of immune system functioning in SA and BD and also suggest that integrating information from multiple levels of analysis holds promise to improve risk assessment for SA in adults with BD.

A Pregnancy and Childhood Epigenetics Consortium (PACE) meta-analysis highlights potential relationships between birth order and neonatal blood DNA methylation

Higher birth order is associated with altered risk of many disease states. Changes in placentation and exposures to in utero growth factors with successive pregnancies may impact later life disease risk via persistent DNA methylation alterations. We investigated birth order with Illumina DNA methylation array data in each of 16 birth cohorts (8164 newborns) with European, African, and Latino ancestries from the Pregnancy and Childhood Epigenetics Consortium. Meta-analyzed data demonstrated systematic DNA methylation variation in 341 CpGs (FDR adjusted P < 0.05) and 1107 regions. Forty CpGs were located within known quantitative trait loci for gene expression traits in blood, and trait enrichment analysis suggested a strong association with immune-related, transcriptional control, and blood pressure regulation phenotypes. Decreasing fertility rates worldwide with the concomitant increased proportion of first-born children highlights a potential reflection of birth order-related epigenomic states on changing disease incidence trends.

Database Resources of the National Genomics Data Center, China National Center for Bioinformation in 2024

The National Genomics Data Center (NGDC), which is a part of the China National Center for Bioinformation (CNCB), provides a family of database resources to support the global academic and industrial communities. With the rapid accumulation of multi-omics data at an unprecedented pace, CNCB-NGDC continuously expands and updates core database resources through big data archiving, integrative analysis and value-added curation. Importantly, NGDC collaborates closely with major international databases and initiatives to ensure seamless data exchange and interoperability. Over the past year, significant efforts have been dedicated to integrating diverse omics data, synthesizing expanding knowledge, developing new resources, and upgrading major existing resources. Particularly, several database resources are newly developed for the biodiversity of protists (P10K), bacteria (NTM-DB, MPA) as well as plant (PPGR, SoyOmics, PlantPan) and disease/trait association (CROST, HervD Atlas, HALL, MACdb, BioKA, BioKA, RePoS, PGG.SV, NAFLDkb). All the resources and services are publicly accessible at https://ngdc.cncb.ac.cn.

HALL: a comprehensive database for human aging and longevity studies

Diverse individuals age at different rates and display variable susceptibilities to tissue aging, functional decline and aging-related diseases. Centenarians, exemplifying extreme longevity, serve as models for healthy aging. The field of human aging and longevity research is rapidly advancing, garnering significant attention and accumulating substantial data in recent years. Omics technologies, encompassing phenomics, genomics, transcriptomics, proteomics, metabolomics and microbiomics, have provided multidimensional insights and revolutionized cohort-based investigations into human aging and longevity. Accumulated data, covering diverse cells, tissues and cohorts across the lifespan necessitates the establishment of an open and integrated database. Addressing this, we established the Human Aging and Longevity Landscape (HALL), a comprehensive multi-omics repository encompassing a diverse spectrum of human cohorts, spanning from young adults to centenarians. The core objective of HALL is to foster healthy aging by offering an extensive repository of information on biomarkers that gauge the trajectory of human aging. Moreover, the database facilitates the development of diagnostic tools for aging-related conditions and empowers targeted interventions to enhance longevity. HALL is publicly available at https://ngdc.cncb.ac.cn/hall/index.

A multi-omics data analysis workflow packaged as a FAIR Digital Object

BACKGROUND

Applying good data management and FAIR (Findable, Accessible, Interoperable, and Reusable) data principles in research projects can help disentangle knowledge discovery, study result reproducibility, and data reuse in future studies. Based on the concepts of the original FAIR principles for research data, FAIR principles for research software were recently proposed. FAIR Digital Objects enable discovery and reuse of Research Objects, including computational workflows for both humans and machines. Practical examples can help promote the adoption of FAIR practices for computational workflows in the research community. We developed a multi-omics data analysis workflow implementing FAIR practices to share it as a FAIR Digital Object.

FINDINGS

We conducted a case study investigating shared patterns between multi-omics data and childhood externalizing behavior. The analysis workflow was implemented as a modular pipeline in the workflow manager Nextflow, including containers with software dependencies. We adhered to software development practices like version control, documentation, and licensing. Finally, the workflow was described with rich semantic metadata, packaged as a Research Object Crate, and shared via WorkflowHub.

CONCLUSIONS

Along with the packaged multi-omics data analysis workflow, we share our experiences adopting various FAIR practices and creating a FAIR Digital Object. We hope our experiences can help other researchers who develop omics data analysis workflows to turn FAIR principles into practice.

Triggering pyroptosis enhances the antitumor efficacy of PARP inhibitors in prostate cancer

PURPOSE

PARP inhibitors have revolutionized the treatment landscape for advanced prostate cancer (PCa) patients who harboring mutations in homologous recombination repair (HRR) genes. However, the molecular mechanisms underlying PARP inhibitors function beyond DNA damage repair pathways remain elusive, and identifying novel predictive targets that favorably respond to PARP inhibitors in PCa is an active area of research.

METHODS

The expression of GSDME in PCa cell lines and human PCa samples was determined by western blotting. Targeted bisulfite sequencing, gene enrichment analysis (GSEA), clone formation, construction of the stably transfected cell lines, lactate dehydrogenase (LDH) assay, western blotting as well as a mouse model of subcutaneous xenografts were used to investigate the role of GSDME in PCa. The combinational therapeutic effect of olaparib and decitabine was determined using both in vitro and in vivo experiments.

RESULTS

We have found low expression of GSDME in PCa. Interestingly, we demonstrated that GSDME activity is robustly induced in olaparib-treated cells undergoing pyroptosis, and that high methylation of the GSDME promoter dampens its activity in PCa cells. Intriguingly, genetically overexpressing GSDME does not inhibit tumor cell proliferation but instead confers sensitivity to olaparib. Furthermore, pharmacological treatment with the combination of olaparib and decitabine synergistically induces GSDME expression and cleavage through caspase-3 activation, thus promoting pyroptosis and enhancing anti-tumor response, ultimately resulting in tumor remission.

CONCLUSION

Our findings highlight a novel therapeutic strategy for enhancing the long-term response to olaparib beyond HRR-deficient tumors in PCa, underscoring the critical role of GSDME in regulating tumorigenesis.

Vimentin promotes glioma progression and maintains glioma cell resistance to oxidative phosphorylation inhibition

PURPOSE

Glioma has been demonstrated as one of the most malignant intracranial tumors and currently there is no effective treatment. Based on our previous RNA-sequencing data for oxidative phosphorylation (OXPHOS)-inhibition resistant and OXPHOS-inhibition sensitive cancer cells, we found that vimentin (VIM) is highly expressed in the OXPHOS-inhibition resistant cancer cells, especially in glioma cancer cells. Further study of VIM in the literature indicates that it plays important roles in cancer progression, immunotherapy suppression, cancer stemness and drug resistance. However, its role in glioma remains elusive. This study aims to decipher the role of VIM in glioma, especially its role in OXPHOS-inhibition sensitivity, which may provide a promising therapeutic target for glioma treatment.

METHODS

The expression of VIM in glioma and the normal tissue has been obtained from The Cancer Genome Atlas (TCGA) database, and further validated in Human Protein Atlas (HPA) and Chinese Glioma Genome Atlas (CGGA). And the single-cell sequencing data was obtained from TISCH2. The immune infiltration was calculated via Tumor Immune Estimation Resource (TIMER), Estimation of Stromal and Immune Cells in Malignant Tumors using Expression Data (ESTIMATE) and ssGSEA, and the Immunophenoscore (IPS) was calculated via R package. The differentiated expressed genes were analyzed including GO/KEGG and Gene Set Enrichment Analysis (GSEA) between the VIM-high and -low groups. The methylation of VIM was checked at the EWAS and Methsurv. The correlation between VIM expression and cancer stemness was obtained from SangerBox. We also employed DepMap data and verified the role of VIM by knocking down it in VIM-high glioma cell and over-expressing it in VIM-low glioma cells to check the cell viability.

RESULTS

Vim is highly expressed in the glioma patients compared to normal samples and its high expression negatively correlates with patients' survival. The DNA methylation in VIM promoters in glioma patients is lower than that in the normal samples. High VIM expression positively correlates with the immune infiltration and tumor progression. Furthermore, Vim is expressed high in the OXPHOS-inhibition glioma cancer cells and low in the OXPHOS-inhibition sensitive ones and its expression maintains the OXPHOS-inhibition resistance.

CONCLUSIONS

In conclusion, we comprehensively deciphered the role of VIM in the progression of glioma and its clinical outcomes. Thus provide new insights into targeting VIM in glioma cancer immunotherapy in combination with the current treatment.

A consistent pattern of slide effects in Illumina DNA methylation BeadChip array data

Background: Recent studies have identified thousands of associations between DNA methylation CpGs and complex diseases/traits, emphasizing the critical role of epigenetics in understanding disease aetiology and identifying biomarkers. However, association analyses based on methylation array data are susceptible to batch/slide effects, which can lead to inflated false positive rates or reduced statistical powerResults: We use multiple DNA methylation datasets based on the popular Illumina Infinium MethylationEPIC BeadChip array to describe consistent patterns and the joint distribution of slide effects across CpGs, confirming and extending previous results. The susceptible CpGs overlap with the Illumina Infinium HumanMethylation450 BeadChip array content.Conclusions: Our findings reveal systematic patterns in slide effects. The observations provide further insights into the characteristics of these effects and can improve existing adjustment approaches.

Deciphering the DNA methylome in women with PCOS diagnosed using the new international evidence-based guidelines

STUDY QUESTION

Is there any methylome alteration in women with PCOS who were diagnosed using the new international evidence-based guidelines?

SUMMARY ANSWER

A total of 264 differentially methylated probes (DMPs) and 53 differentially methylated regions (DMRs) were identified in patients with PCOS and healthy controls.

WHAT IS KNOWN ALREADY

PCOS is a common endocrine disorder among women of reproductive age and polycystic ovarian morphology (PCOM) is one of the main features of the disease. Owing to the availability of more sensitive ultrasound machines, the traditional diagnosis of PCOM according to the Rotterdam criteria (≥12 antral follicles per ovary) is currently debated as there is a risk of overdiagnosis. The new international evidence-based guidelines set the threshold for PCOM as ≥20 antral follicles per ovary when using endovaginal ultrasound transducers with a frequency bandwidth that includes 8 MHz. However, current DNA methylation studies in PCOS are still based on the Rotterdam criteria. This study aimed to explore aberrant DNA methylation in patients diagnosed with PCOS according to the new evidence-based guidelines.

STUDY DESIGN, SIZE, DURATION

This cross-sectional case-control study included 34 PCOS cases diagnosed using new international evidence-based guidelines and 36 controls.

PARTICIPANTS/MATERIALS, SETTING, METHODS

A total of 70 women, including 34 PCOS cases and 36 controls, were recruited. DNA extracted from whole blood samples of participants were profiled using array technology. Data quality control, preprocessing, annotation, and statistical analyses were performed. Least absolute shrinkage and selection operator (LASSO) regression were used to build a PCOS diagnosis model with DNA methylation sites.

MAIN RESULTS AND THE ROLE OF CHANCE

We identified 264 DMPs between PCOS cases and controls, which were mainly located in intergenic regions or gene bodies of the genome, CpG open sea sites, and heterochromatin of functional elements. Pathway enrichment analysis showed that DMPs were significantly enriched in biological processes involved in triglyceride regulation. Three of these DMPs overlapped with the PCOS susceptibility genes thyroid adenoma-associated protein (THADA), aminopeptidase O (AOPEP), and tripartite motif family-like protein 2 (TRIML2). Fifty-three DMRs were identified and their annotated genes were largely enriched in allograft rejection, thyroid hormone production, and peripheral downstream signaling effects. Two DMRs were closely related to the PCOS susceptibility genes, potassium voltage-gated channel subfamily A member 4 (KCNA4) and farnesyl-diphosphate farnesyltransferase 1 (FDFT1). Finally, based on LASSO regression, we built a methylation marker model with high accuracy for PCOS diagnosis (AUC=0.952).

LIMITATIONS, REASONS FOR CAUTION

The study cohort was single-center and the sample size was relatively limited. Further analyses with a larger number of participants are required.

WIDER IMPLICATIONS OF THE FINDINGS

This is the first study to identify DNA methylation alterations in women with PCOS diagnosed using the new international evidence-based guideline, and it provided new molecular insight into the application of the new guidelines.

STUDY FUNDING/COMPETING INTEREST(S)

This study was supported by the National Key Research and Development Program of China (2021YFC2700400), Basic Science Center Program of NSFC (31988101), CAMS Innovation Fund for Medical Sciences (2021-I2M-5-001), National Natural Science Foundation of China (32370916, 82071606, 82101707, 82192874, and 31871509), Shandong Provincial Key Research and Development Program (2020ZLYS02), Taishan Scholars Program of Shandong Province (ts20190988), and Fundamental Research Funds of Shandong University. The authors declare no conflicts of interest.

TRIAL REGISTRATION NUMBER

N/A.

Epigenome-wide methylation and progression to high-grade cervical intraepithelial neoplasia (CIN2+): a prospective cohort study in the United States

BACKGROUND

Methylation levels may be associated with and serve as markers to predict risk of progression of precancerous cervical lesions. We conducted an epigenome-wide association study (EWAS) of CpG methylation and progression to high-grade cervical intraepithelial neoplasia (CIN2 +) following an abnormal screening test.

METHODS

A prospective US cohort of 289 colposcopy patients with normal or CIN1 enrollment histology was assessed. Baseline cervical sample DNA was analyzed using Illumina HumanMethylation 450K (n = 76) or EPIC 850K (n = 213) arrays. Participants returned at provider-recommended intervals and were followed up to 5 years via medical records. We assessed continuous CpG M values for 9 cervical cancer-associated genes and time-to-progression to CIN2+. We estimated CpG-specific time-to-event ratios (TTER) and hazard ratios using adjusted, interval-censored Weibull accelerated failure time models. We also conducted an exploratory EWAS to identify novel CpGs with false discovery rate (FDR) < 0.05.

RESULTS

At enrollment, median age was 29.2 years; 64.0% were high-risk HPV-positive, and 54.3% were non-white. During follow-up (median 24.4 months), 15 participants progressed to CIN2+. Greater methylation levels were associated with a shorter time-to-CIN2+ for CADM1 cg03505501 (TTER = 0.28; 95%CI 0.12, 0.63; FDR = 0.03) and RARB Cluster 1 (TTER = 0.46; 95% CI 0.29, 0.71; FDR = 0.01). There was evidence of similar trends for DAPK1 cg14286732, PAX1 cg07213060, and PAX1 Cluster 1. The EWAS detected 336 novel progression-associated CpGs, including those located in CpG islands associated with genes FGF22, TOX, COL18A1, GPM6A, XAB2, TIMP2, GSPT1, NR4A2, and APBB1IP.

CONCLUSIONS

Using prospective time-to-event data, we detected associations between CADM1-, DAPK1-, PAX1-, and RARB-related CpGs and cervical disease progression, and we identified novel progression-associated CpGs.

IMPACT

Methylation levels at novel CpG sites may help identify individuals with ≤CIN1 histology at higher risk of progression to CIN2+ and inform risk-based cervical cancer screening guidelines.

Genomic gain/methylation modification/hsa-miR-132-3p increases RRS1 overexpression in liver hepatocellular carcinoma

This study aimed to determine the upstream regulatory factors affecting ribosome biogenesis regulator 1 homolog (RRS1) expression and the development and prognosis of liver hepatocellular carcinoma (LIHC). The expression profiles of RRS1 were evaluated in pan-cancer tissues and liver tumor cell lines. The associations of RRS1 with pan-cancer survival, immune infiltrations, immune checkpoints, and drug sensitivity were identified. We explored the potential upstream regulatory mechanisms of RRS1 expression. Hsa-miR-132-3p knockdown, CCK-8 assays, transwell, and wound healing assays were performed to validate the regulatory effect of hsa-miR-132-3p on RRS1 expression and the development of LIHC. Our findings demonstrated that RRS1 was significantly elevated in 27 types of cancers. RRS1 predicts a poor outcome of LIHC, lung adenocarcinoma, head and neck cancer, and kidney papillary cell carcinoma. RRS1 expression showed a significant association with immune cell infiltrates and the expression of immune checkpoints-related genes in LIHC tissues. Increased RRS1 expression may have a negative effect on these anticancer drugs of LIHC. Low methylation of the RRS1 promoter and its genomic gain may elevate RRS1 expression and predict poor prognosis for LIHC. Increased hsa-miR-132-3p expression may elevate RRS1 expression and result in poor prognosis for LIHC. Hsa-miR-132-3p inhibition can decrease RRS1 expression and the development of liver tumor cell lines. Low methylation of the RRS1 promoter, RRS1 genomic gain, and hsa-miR-132-3p upregulation in LIHC may promote RRS1 upregulation and thus lead to the development and poor prognosis for LIHC. RRS1 is a promising therapeutic target for LIHC.

Exploring the Functional Basis of Epigenetic Aging in Relation to Body Fat Phenotypes in the Norfolk Island Cohort

DNA methylation is an epigenetic factor that is modifiable and can change over a lifespan. While many studies have identified methylation sites (CpGs) related to aging, the relationship of these to gene function and age-related disease phenotypes remains unclear. This research explores this question by testing for the conjoint association of age-related CpGs with gene expression and the relation of these to body fat phenotypes. The study included blood-based gene transcripts and intragenic CpG methylation data from Illumina 450 K arrays in 74 healthy adults from the Norfolk Island population. First, a series of regression analyses were performed to detect associations between gene transcript level and intragenic CpGs and their conjoint relationship with age. Second, we explored how these age-related expression CpGs (eCpGs) correlated with obesity-related phenotypes, including body fat percentage, body mass index, and waist-to-hip ratio. We identified 35 age-related eCpGs associated with age. Of these, ten eCpGs were associated with at least one body fat phenotype. Collagen Type XI Alpha 2 Chain (COL11A2), Complement C1s (C1s), and four and a half LIM domains 2 (FHL2) genes were among the most significant genes with multiple eCpGs associated with both age and multiple body fat phenotypes. The COL11A2 gene contributes to the correct assembly of the extracellular matrix in maintaining the healthy structural arrangement of various components, with the C1s gene part of complement systems functioning in inflammation. Moreover, FHL2 expression was upregulated under hypermethylation in both blood and adipose tissue with aging. These results suggest new targets for future studies and require further validation to confirm the specific function of these genes on body fat regulation.

Caveolin-1 promotes glioma progression and maintains its mitochondrial inhibition resistance

BACKGROUND

Glioma is a lethal brain cancer and lacking effective therapies. Challenges include no effective therapeutic target, intra- and intertumoral heterogeneity, inadequate effective drugs, and an immunosuppressive microenvironment, etc. Deciphering the pathogenesis of gliomas and finding out the working mechanisms are urgent and necessary for glioma treatment. Identification of prognostic biomarkers and targeting the biomarker genes will be a promising therapy.

METHODS

From our RNA-sequencing data of the oxidative phosphorylation (OXPHOS)-inhibition sensitive and OXPHOS-resistant cell lines, we found that the scaffolding protein caveolin 1 (CAV1) is highly expressed in the resistant group but not in the sensitive group. By comprehensive analysis of our RNA sequencing data, Whole Genome Bisulfite Sequencing (WGBS) data and public databases, we found that CAV1 is highly expressed in gliomas and its expression is positively related with pathological processes, higher CAV1 predicts shorter overall survival.

RESULTS

Further analysis indicated that (1) the differentiated genes in CAV1-high groups are enriched in immune infiltration and immune response; (2) CAV1 is positively correlated with tumor metastasis markers; (3) the methylation level of CAV1 promoters in glioma group is lower in higher stage than that in lower stage; (4) CAV1 is positively correlated with glioma stemness; (5) higher expression of CAV1 renders the glioma cells' resistant to oxidative phosphorylation inhibitors.

CONCLUSION

Therefore, we identified a key gene CAV1 and deciphered its function in glioma progression and prognosis, proposing that CAV1 may be a therapeutic target for gliomas.

Integrative multi-omics and systems bioinformatics in translational neuroscience: A data mining perspective

Bioinformatic analysis of large and complex omics datasets has become increasingly useful in modern day biology by providing a great depth of information, with its application to neuroscience termed neuroinformatics. Data mining of omics datasets has enabled the generation of new hypotheses based on differentially regulated biological molecules associated with disease mechanisms, which can be tested experimentally for improved diagnostic and therapeutic targeting of neurodegenerative diseases. Importantly, integrating multi-omics data using a systems bioinformatics approach will advance the understanding of the layered and interactive network of biological regulation that exchanges systemic knowledge to facilitate the development of a comprehensive human brain profile. In this review, we first summarize data mining studies utilizing datasets from the individual type of omics analysis, including epigenetics/epigenomics, transcriptomics, proteomics, metabolomics, lipidomics, and spatial omics, pertaining to Alzheimer's disease, Parkinson's disease, and multiple sclerosis. We then discuss multi-omics integration approaches, including independent biological integration and unsupervised integration methods, for more intuitive and informative interpretation of the biological data obtained across different omics layers. We further assess studies that integrate multi-omics in data mining which provide convoluted biological insights and offer proof-of-concept proposition towards systems bioinformatics in the reconstruction of brain networks. Finally, we recommend a combination of high dimensional bioinformatics analysis with experimental validation to achieve translational neuroscience applications including biomarker discovery, therapeutic development, and elucidation of disease mechanisms. We conclude by providing future perspectives and opportunities in applying integrative multi-omics and systems bioinformatics to achieve precision phenotyping of neurodegenerative diseases and towards personalized medicine.

Blood epigenome-wide association studies of suicide attempt in adults with bipolar disorder

Suicide attempt (SA) risk is elevated in individuals with bipolar disorder (BD), and DNA methylation patterns may serve as possible biomarkers of SA. We conducted epigenome-wide association studies (EWAS) of blood DNA methylation associated with BD and SA. DNA methylation was measured at > 700,000 positions in a discovery cohort of n = 84 adults with BD with a history of SA (BD/SA), n = 79 adults with BD without history of SA (BD/non-SA), and n = 76 non-psychiatric controls (CON). EWAS revealed six differentially methylated positions (DMPs) and seven differentially methylated regions (DMRs) between BD/SA and BD/non-SA, with multiple immune-related genes implicated. There were no epigenome-wide significant differences when BD/SA and BD/non-SA were each compared to CON, and patterns suggested that epigenetics differentiating BD/SA from BD/non-SA do not differentiate BD/non-SA from CON. Weighted gene co-methylation network analysis and trait enrichment analysis of the BD/SA vs. BD/non-SA contrast further corroborated immune system involvement, while gene ontology analysis implicated calcium signalling. In an independent replication cohort of n = 48 BD/SA and n = 47 BD/non-SA, fold-changes at the discovery cohort's significant sites showed moderate correlation across cohorts and agreement on direction. In both cohorts, classification accuracy for SA history among individuals with BD was highest when methylation at the significant CpG sites as well as information from clinical interviews were combined, with an AUC of 88.8% (CI = 83.8-93.8%) and 82.1% (CI = 73.6-90.5%) for the combined epigenetic-clinical predictor in the discovery and replication cohorts, respectively. Our results provide novel insight to the role of immune system functioning in SA and BD and also suggest that integrating information from multiple levels of analysis holds promise to improve risk assessment for SA in adults with BD.

Epigenomic profiling of isolated blood cell types reveals highly specific B cell smoking signatures and links to disease risk

BACKGROUND

Tobacco smoking alters the DNA methylation profiles of immune cells which may underpin some of the pathogenesis of smoking-associated diseases. To link smoking-driven epigenetic effects in specific immune cell types with disease risk, we isolated six leukocyte subtypes, CD14+ monocytes, CD15+ granulocytes, CD19+ B cells, CD4+ T cells, CD8+ T cells, and CD56+ natural killer cells, from whole blood of 67 healthy adult smokers and 74 nonsmokers for epigenome-wide association study (EWAS) using Illumina 450k and EPIC methylation arrays.

RESULTS

Numbers of smoking-associated differentially methylated sites (smCpGs) at genome-wide significance (p < 1.2 × 10-7) varied widely across cell types, from 5 smCpGs in CD8+ T cells to 111 smCpGs in CD19+ B cells. We found unique smoking effects in each cell type, some of which were not apparent in whole blood. Methylation-based deconvolution to estimate B cell subtypes revealed that smokers had 7.2% (p = 0.033) less naïve B cells. Adjusting for naïve and memory B cell proportions in EWAS and RNA-seq allowed the identification of genes enriched for B cell activation-related cytokine signaling pathways, Th1/Th2 responses, and hematopoietic cancers. Integrating with large-scale public datasets, 62 smCpGs were among CpGs associated with health-relevant EWASs. Furthermore, 74 smCpGs had reproducible methylation quantitative trait loci single nucleotide polymorphisms (SNPs) that were in complete linkage disequilibrium with genome-wide association study SNPs, associating with lung function, disease risks, and other traits.

CONCLUSIONS

We observed blood cell-type-specific smCpGs, a naïve-to-memory shift among B cells, and by integrating genome-wide datasets, we identified their potential links to disease risks and health traits.

New genetic and epigenetic insights into the chemokine system: the latest discoveries aiding progression toward precision medicine

Over the past thirty years, the importance of chemokines and their seven-transmembrane G protein-coupled receptors (GPCRs) has been increasingly recognized. Chemokine interactions with receptors trigger signaling pathway activity to form a network fundamental to diverse immune processes, including host homeostasis and responses to disease. Genetic and nongenetic regulation of both the expression and structure of chemokines and receptors conveys chemokine functional heterogeneity. Imbalances and defects in the system contribute to the pathogenesis of a variety of diseases, including cancer, immune and inflammatory diseases, and metabolic and neurological disorders, which render the system a focus of studies aiming to discover therapies and important biomarkers. The integrated view of chemokine biology underpinning divergence and plasticity has provided insights into immune dysfunction in disease states, including, among others, coronavirus disease 2019 (COVID-19). In this review, by reporting the latest advances in chemokine biology and results from analyses of a plethora of sequencing-based datasets, we outline recent advances in the understanding of the genetic variations and nongenetic heterogeneity of chemokines and receptors and provide an updated view of their contribution to the pathophysiological network, focusing on chemokine-mediated inflammation and cancer. Clarification of the molecular basis of dynamic chemokine-receptor interactions will help advance the understanding of chemokine biology to achieve precision medicine application in the clinic.

The X-factor in ART: does the use of assisted reproductive technologies influence DNA methylation on the X chromosome?

BACKGROUND

Assisted reproductive technologies (ART) may perturb DNA methylation (DNAm) in early embryonic development. Although a handful of epigenome-wide association studies of ART have been published, none have investigated CpGs on the X chromosome. To bridge this knowledge gap, we leveraged one of the largest collections of mother-father-newborn trios of ART and non-ART (natural) conceptions to date to investigate sex-specific DNAm differences on the X chromosome. The discovery cohort consisted of 982 ART and 963 non-ART trios from the Norwegian Mother, Father, and Child Cohort Study (MoBa). To verify our results from the MoBa cohort, we used an external cohort of 149 ART and 58 non-ART neonates from the Australian 'Clinical review of the Health of adults conceived following Assisted Reproductive Technologies' (CHART) study. The Illumina EPIC array was used to measure DNAm in both datasets. In the MoBa cohort, we performed a set of X-chromosome-wide association studies ('XWASs' hereafter) to search for sex-specific DNAm differences between ART and non-ART newborns. We tested several models to investigate the influence of various confounders, including parental DNAm. We also searched for differentially methylated regions (DMRs) and regions of co-methylation flanking the most significant CpGs. Additionally, we ran an analogous model to our main model on the external CHART dataset.

RESULTS

In the MoBa cohort, we found more differentially methylated CpGs and DMRs in girls than boys. Most of the associations persisted after controlling for parental DNAm and other confounders. Many of the significant CpGs and DMRs were in gene-promoter regions, and several of the genes linked to these CpGs are expressed in tissues relevant for both ART and sex (testis, placenta, and fallopian tube). We found no support for parental DNAm-dependent features as an explanation for the observed associations in the newborns. The most significant CpG in the boys-only analysis was in UBE2DNL, which is expressed in testes but with unknown function. The most significant CpGs in the girls-only analysis were in EIF2S3 and AMOT. These three loci also displayed differential DNAm in the CHART cohort.

CONCLUSIONS

Genes that co-localized with the significant CpGs and DMRs associated with ART are implicated in several key biological processes (e.g., neurodevelopment) and disorders (e.g., intellectual disability and autism). These connections are particularly compelling in light of previous findings indicating that neurodevelopmental outcomes differ in ART-conceived children compared to those naturally conceived.

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).

Epigenome-Wide Association Studies of the Fractional Exhaled Nitric Oxide and Bronchodilator Drug Response in Moderate-to-Severe Pediatric Asthma

Asthma is the most prevalent pediatric chronic disease. Bronchodilator drug response (BDR) and fractional exhaled nitric oxide (FeNO) are clinical biomarkers of asthma. Although DNA methylation (DNAm) contributes to asthma pathogenesis, the influence of DNAm on BDR and FeNO is scarcely investigated. This study aims to identify DNAm markers in whole blood associated either with BDR or FeNO in pediatric asthma. We analyzed 121 samples from children with moderate-to-severe asthma. The association of genome-wide DNAm with BDR and FeNO has been assessed using regression models, adjusting for age, sex, ancestry, and tissue heterogeneity. Cross-tissue validation was assessed in 50 nasal samples. Differentially methylated regions (DMRs) and enrichment in traits and biological pathways were assessed. A false discovery rate (FDR) < 0.1 and a genome-wide significance threshold of p < 9 × 10−8 were used to control for false-positive results. The CpG cg12835256 (PLA2G12A) was genome-wide associated with FeNO in blood samples (coefficient= −0.015, p = 2.53 × 10−9) and nominally associated in nasal samples (coefficient = −0.015, p = 0.045). Additionally, three CpGs were suggestively associated with BDR (FDR < 0.1). We identified 12 and four DMRs associated with FeNO and BDR (FDR < 0.05), respectively. An enrichment in allergic and inflammatory processes, smoking, and aging was observed. We reported novel associations of DNAm markers associated with BDR and FeNO enriched in asthma-related processes.

Promises and challenges in pharmacoepigenetics

Pharmacogenetics, the study of how interindividual genetic differences affect drug response, does not explain all observed heritable variance in drug response. Epigenetic mechanisms, such as DNA methylation, and histone acetylation may account for some of the unexplained variances. Epigenetic mechanisms modulate gene expression and can be suitable drug targets and can impact the action of nonepigenetic drugs. Pharmacoepigenetics is the field that studies the relationship between epigenetic variability and drug response. Much of this research focuses on compounds targeting epigenetic mechanisms, called epigenetic drugs, which are used to treat cancers, immune disorders, and other diseases. Several studies also suggest an epigenetic role in classical drug response; however, we know little about this area. The amount of information correlating epigenetic biomarkers to molecular datasets has recently expanded due to technological advances, and novel computational approaches have emerged to better identify and predict epigenetic interactions. We propose that the relationship between epigenetics and classical drug response may be examined using data already available by (1) finding regions of epigenetic variance, (2) pinpointing key epigenetic biomarkers within these regions, and (3) mapping these biomarkers to a drug-response phenotype. This approach expands on existing knowledge to generate putative pharmacoepigenetic relationships, which can be tested experimentally. Epigenetic modifications are involved in disease and drug response. Therefore, understanding how epigenetic drivers impact the response to classical drugs is important for improving drug design and administration to better treat disease.

Database Resources of the National Genomics Data Center, China National Center for Bioinformation in 2023

The National Genomics Data Center (NGDC), part of the China National Center for Bioinformation (CNCB), provides a family of database resources to support global academic and industrial communities. With the explosive accumulation of multi-omics data generated at an unprecedented rate, CNCB-NGDC constantly expands and updates core database resources by big data archive, integrative analysis and value-added curation. In the past year, efforts have been devoted to integrating multiple omics data, synthesizing the growing knowledge, developing new resources and upgrading a set of major resources. Particularly, several database resources are newly developed for infectious diseases and microbiology (MPoxVR, KGCoV, ProPan), cancer-trait association (ASCancer Atlas, TWAS Atlas, Brain Catalog, CCAS) as well as tropical plants (TCOD). Importantly, given the global health threat caused by monkeypox virus and SARS-CoV-2, CNCB-NGDC has newly constructed the monkeypox virus resource, along with frequent updates of SARS-CoV-2 genome sequences, variants as well as haplotypes. All the resources and services are publicly accessible at https://ngdc.cncb.ac.cn.

Integrative Multi-omics Analysis of Childhood Aggressive Behavior

This study introduces and illustrates the potential of an integrated multi-omics approach in investigating the underlying biology of complex traits such as childhood aggressive behavior. In 645 twins (cases = 42%), we trained single- and integrative multi-omics models to identify biomarkers for subclinical aggression and investigated the connections among these biomarkers. Our data comprised transmitted and two non-transmitted polygenic scores (PGSs) for 15 traits, 78,772 CpGs, and 90 metabolites. The single-omics models selected 31 PGSs, 1614 CpGs, and 90 metabolites, and the multi-omics model comprised 44 PGSs, 746 CpGs, and 90 metabolites. The predictive accuracy for these models in the test (N = 277, cases = 42%) and independent clinical data (N = 142, cases = 45%) ranged from 43 to 57%. We observed strong connections between DNA methylation, amino acids, and parental non-transmitted PGSs for ADHD, Autism Spectrum Disorder, intelligence, smoking initiation, and self-reported health. Aggression-related omics traits link to known and novel risk factors, including inflammation, carcinogens, and smoking.

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.

EpiVisR: exploratory data analysis and visualization in epigenome-wide association analyses

Background

With the widespread availability of microarray technology for epigenetic research, methods for calling differentially methylated probes or differentially methylated regions have become effective tools to analyze this type of data. Furthermore, visualization is usually employed for quality check of results and for further insights. Expert knowledge is required to leverage capabilities of these methods. To overcome this limitation and make visualization in epigenetic research available to the public, we designed EpiVisR.

Results

The EpiVisR tool allows to select and visualize combinations of traits (i.e., concentrations of chemical compounds) and differentially methylated probes/regions. It supports various modes of enriched presentation to get the most knowledge out of existing data: (1) enriched Manhattan plot and enriched volcano plot for selection of probes, (2) trait-methylation plot for visualization of selected trait values against methylation values, (3) methylation profile plot for visualization of a selected range of probes against selected trait values as well as, (4) correlation profile plot for selection and visualization of further probes that are correlated to the selected probe. EpiVisR additionally allows exporting selected data to external tools for tasks such as network analysis.

Conclusion

The key advantage of EpiVisR is the annotation of data in the enriched plots (and tied tables) as well as linking to external data sources for further integrated data analysis. Using the EpiVisR approach will allow users to integrate data from traits with epigenetic analyses that are connected by belonging to the same individuals. Merging data from various data sources among the same cohort and visualizing them will enable users to gain more insights from existing data.

MethBank 4.0: an updated database of DNA methylation across a variety of species

DNA methylation, as the most intensively studied epigenetic mark, regulates gene expression in numerous biological processes including development, aging, and disease. With the rapid accumulation of whole-genome bisulfite sequencing data, integrating, archiving, analyzing, and visualizing those data becomes critical. Since its first publication in 2015, MethBank has been continuously updated to include more DNA methylomes across more diverse species. Here, we present MethBank 4.0 (https://ngdc.cncb.ac.cn/methbank/), which reports an increase of 309% in data volume, with 1449 single-base resolution methylomes of 23 species, covering 236 tissues/cell lines and 15 biological contexts. Value-added information, such as more rigorous quality evaluation, more standardized metadata, and comprehensive downstream annotations have been integrated in the new version. Moreover, expert-curated knowledge modules of featured differentially methylated genes associated with biological contexts and methylation analysis tools have been incorporated as new components of MethBank. In addition, MethBank 4.0 is equipped with a series of new web interfaces to browse, search, and visualize DNA methylation profiles and related information. With all these improvements, we believe the updated MethBank 4.0 will serve as a fundamental resource to provide a wide range of data services for the global research community.

ASCancer Atlas: a comprehensive knowledgebase of alternative splicing in human cancers

Alternative splicing (AS) is a fundamental process that governs almost all aspects of cellular functions, and dysregulation in this process has been implicated in tumor initiation, progression and treatment resistance. With accumulating studies of carcinogenic mis-splicing in cancers, there is an urgent demand to integrate cancer-associated splicing changes to better understand their internal cross-talks and functional consequences from a global view. However, a resource of key functional AS events in human cancers is still lacking. To fill the gap, we developed ASCancer Atlas (https://ngdc.cncb.ac.cn/ascancer), a comprehensive knowledgebase of aberrant splicing in human cancers. Compared to extant databases, ASCancer Atlas features a high-confidence collection of 2006 cancer-associated splicing events experimentally proved to promote tumorigenesis, a systematic splicing regulatory network, and a suit of multi-scale online analysis tools. For each event, we manually curated the functional axis including upstream splicing regulators, splicing event annotations, downstream oncogenic effects, and possible therapeutic strategies. ASCancer Atlas also houses about 2 million computationally putative splicing events. Additionally, a user-friendly web interface was built to enable users to easily browse, search, visualize, analyze, and download all splicing events. Overall, ASCancer Atlas provides a unique resource to study the functional roles of splicing dysregulation in human cancers.

Multi-Omics Analysis of the Tumor Microenvironment in Liver Metastasis of Colorectal Cancer Identified FJX1 as a Novel Biomarker

Colorectal cancer incidence and mortality have increased in recent years, with more than half of patients who died of colorectal cancer developing liver metastases. Consequently, colorectal cancer liver metastasis is the focus of clinical treatment, as well as being the most difficult. The primary target genes related to colorectal cancer liver metastasis were via bioinformatics analysis. First, five prognosis-related genes, CTAG1A, CSTL1, FJX1, IER5L, and KLHL35, were identified through screening, and the prognosis of the CSTL1, FJX1, IER5L, and KLHL35 high expression group was considerably poorer than that of the low expression group. Furthermore, the clinical correlation analysis revealed that in distinct pathological stages T, N, and M, the mRNA expression levels of CSTL1, IER5L, and KLHL35 were higher than in normal tissues. Finally, a correlation study of the above genes and clinical manifestations revealed that FJX1 was strongly linked to colorectal cancer liver metastasis. FJX1 is thought to affect chromogenic modification enzymes, the Notch signaling system, cell senescence, and other signaling pathways, according to KEGG enrichment analysis. FJX1 may be a critical target in colorectal cancer metastasis, and thus has the potential as a new biomarker to predict and treat colorectal cancer liver metastases.

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.

Identification of COVID-19-Associated DNA Methylation Variations by Integrating Methylation Array and scRNA-Seq Data at Cell-Type Resolution

Single-cell transcriptome studies have revealed immune dysfunction in COVID-19 patients, including lymphopenia, T cell exhaustion, and increased levels of pro-inflammatory cytokines, while DNA methylation plays an important role in the regulation of immune response and inflammatory response. The specific cell types of immune responses regulated by DNA methylation in COVID-19 patients will be better understood by exploring the COVID-19 DNA methylation variation at the cell-type level. Here, we developed an analytical pipeline to explore single-cell DNA methylation variations in COVID-19 patients by transferring bulk-tissue-level knowledge to the single-cell level. We discovered that the methylation variations in the whole blood of COVID-19 patients showed significant cell-type specificity with remarkable enrichment in gamma-delta T cells and presented a phenomenon of hypermethylation and low expression. Furthermore, we identified five genes whose methylation variations were associated with several cell types. Among them, S100A9, AHNAK, and CX3CR1 have been reported as potential COVID-19 biomarkers previously, and the others (TRAF3IP3 and LFNG) are closely associated with the immune and virus-related signaling pathways. We propose that they might serve as potential epigenetic biomarkers for COVID-19 and could play roles in important biological processes such as the immune response and antiviral activity.

Identifying the Potential Roles of PBX4 in Human Cancers Based on Integrative Analysis

PBX4 belongs to the pre-B-cell leukemia homeobox (PBX) transcription factors family and acts as a transcriptional cofactor of HOX proteins participating in several pathophysiological processes. Recent studies have revealed that the dysregulation of PBX4 is closely related to multiple diseases, especially cancers. However, the research on PBX4’s potential roles in 33 cancers from the Cancer Genome Atlas (TCGA) is still insufficient. Therefore, we performed a comprehensive pan-cancer analysis to explore the roles of PBX4with multiple public databases. Our results showed that PBX4 was differentially expressed in 17 types of human cancer and significantly correlated to the pathological stage, tumor grade, and immune and molecular subtypes. We used the Kaplan–Meier plotter and PrognoScan databases to find the significant associations between PBX4 expression and prognostic values of multiple cancers. It was also found that PBX4 expression was statistically related to mutation status, DNA methylation, immune infiltration, drug sensitivity, and immune checkpoint blockade (ICB) therapy. Additionally, we found that PBX4 was involved in different functional states of multiple cancers from the single-cell resolution perspective. Enrichment analysis results showed that PBX4-related genes were enriched in the cell cycle process, MAPK cascade, ncRNA metabolic process, positive regulation of GTPase activity, and regulation of lipase activity and mainly participated in the pathways of cholesterol metabolism, base excision repair, herpes simplex virus 1 infection, transcriptional misregulation in cancer, and Epstein–Barr virus infection. Altogether, our integrative analysis could help in better understanding the potential roles of PBX4 in different human cancers.

An epigenome-wide study of DNA methylation profiles and lung function among American Indians in the Strong Heart Study

BACKGROUND

Epigenetic modifications, including DNA methylation (DNAm), are often related to environmental exposures, and are increasingly recognized as key processes in the pathogenesis of chronic lung disease. American Indian communities have a high burden of lung disease compared to the national average. The objective of this study was to investigate the association of DNAm and lung function in the Strong Heart Study (SHS). We conducted a cross-sectional study of American Indian adults, 45-74 years of age who participated in the SHS. DNAm was measured using the Illumina Infinium Human MethylationEPIC platform at baseline (1989-1991). Lung function was measured via spirometry, including forced expiratory volume in 1 s (FEV1) and forced vital capacity (FVC), at visit 2 (1993-1995). Airflow limitation was defined as FEV1 < 70% predicted and FEV1/FVC < 0.7, restriction was defined as FEV1/FVC > 0.7 and FVC < 80% predicted, and normal spirometry was defined as FEV1/FVC > 0.7, FEV1 > 70% predicted, FVC > 80% predicted. We used elastic-net models to select relevant CpGs for lung function and spirometry-defined lung disease. We also conducted bioinformatic analyses to evaluate the biological plausibility of the findings.

RESULTS

Among 1677 participants, 21.2% had spirometry-defined airflow limitation and 13.6% had spirometry-defined restrictive pattern lung function. Elastic-net models selected 1118 Differentially Methylated Positions (DMPs) as predictors of airflow limitation and 1385 for restrictive pattern lung function. A total of 12 DMPs overlapped between airflow limitation and restrictive pattern. EGFR, MAPK1 and PRPF8 genes were the most connected nodes in the protein-protein interaction network. Many of the DMPs targeted genes with biological roles related to lung function such as protein kinases.

CONCLUSION

We found multiple differentially methylated CpG sites associated with chronic lung disease. These signals could contribute to better understand molecular mechanisms involved in lung disease, as assessed systemically, as well as to identify patterns that could be useful for diagnostic purposes. Further experimental and longitudinal studies are needed to assess whether DNA methylation has a causal role in lung disease.

Incorporating local ancestry improves identification of ancestry-associated methylation signatures and meQTLs in African Americans

Here we report three epigenome-wide association studies (EWAS) of DNA methylation on self-reported race, global genetic ancestry, and local genetic ancestry in admixed Americans from three sets of samples, including internal and external replications (Ntotal = 1224). Our EWAS on local ancestry (LA) identified the largest number of ancestry-associated DNA methylation sites and also featured the highest replication rate. Furthermore, by incorporating ancestry origins of genetic variations, we identified 36 methylation quantitative trait loci (meQTL) clumps for LA-associated CpGs that cannot be captured by a model that assumes identical genetic effects across ancestry origins. Lead SNPs at 152 meQTL clumps had significantly different genetic effects in the context of an African or European ancestry background. Local ancestry information enables superior capture of ancestry-associated methylation signatures and identification of ancestry-specific genetic effects on DNA methylation. These findings highlight the importance of incorporating local ancestry for EWAS in admixed samples from multi-ancestry cohorts.

GMQN: A Reference-Based Method for Correcting Batch Effects and Probe Bias in HumanMethylation BeadChip

The Illumina HumanMethylation BeadChip is one of the most cost-effective methods to quantify DNA methylation levels at single-base resolution across the human genome, which makes it a routine platform for epigenome-wide association studies. It has accumulated tens of thousands of DNA methylation array samples in public databases, providing great support for data integration and further analysis. However, the majority of public DNA methylation data are deposited as processed data without background probes which are widely used in data normalization. Here, we present Gaussian mixture quantile normalization (GMQN), a reference based method for correcting batch effects as well as probe bias in the HumanMethylation BeadChip. Availability and implementation: https://github.com/MengweiLi-project/gmqn.

Database Resources of the National Genomics Data Center, China National Center for Bioinformation in 2022

The National Genomics Data Center (NGDC), part of the China National Center for Bioinformation (CNCB), provides a family of database resources to support global research in both academia and industry. With the explosively accumulated multi-omics data at ever-faster rates, CNCB-NGDC is constantly scaling up and updating its core database resources through big data archive, curation, integration and analysis. In the past year, efforts have been made to synthesize the growing data and knowledge, particularly in single-cell omics and precision medicine research, and a series of resources have been newly developed, updated and enhanced. Moreover, CNCB-NGDC has continued to daily update SARS-CoV-2 genome sequences, variants, haplotypes and literature. Particularly, OpenLB, an open library of bioscience, has been established by providing easy and open access to a substantial number of abstract texts from PubMed, bioRxiv and medRxiv. In addition, Database Commons is significantly updated by cataloguing a full list of global databases, and BLAST tools are newly deployed to provide online sequence search services. All these resources along with their services are publicly accessible at https://ngdc.cncb.ac.cn.