Genetic variation influencing DNA methylation provides insights into molecular mechanisms regulating genomic function

Analysis of the interlink between glucose-6-phosphate dehydrogenase (G6PD) and lung cancer through multi-omics databases

Glucose-6-Phosphate Dehydrogenase (G6PD) is a crucial enzyme that executes the pentose phosphate pathway. Due to its critical nodal position in the metabolic network, it is associated with different forms of cancer tumorigeneses and progression. Nonetheless, its functional role and molecular mechanism in lung cancer remain unknown. The present study provides intricate information associated with G6PD and Lung Cancer. Varieties of public datasets were retrieved by us, including UALCAN, TCGA, cBioPortal, and the UCSC Xena browser. The data obtained were used to assess the expression of G6PD, its clinical features, epigenetic regulation, relationship with tumour infiltration, tumour mutation burden, microsatellite instability, tumour microenvironment, immune checkpoint genes, genomic alteration, and patient's overall survival rate. The present study revealed that the G6PD expression was correlated with the clinical features of lung cancer including disease stage, race, sex, age, smoking habits, and lymph node metastasis. Moreover, the expression profile of G6PD also imparts epigenetic changes by modulating the DNA promoter methylation activity. Methylation of promoters changes the expression of various transcription factors, genes leading to an influence on the immune system. These events linked with G6PD-related mutational gene alterations (FAM3A, LAG3, p53, KRAS). The entire circumstance influences the patient's overall survival rate and poor prognosis. Functional investigation using STRING, GO, and KEGG found that G6PD primarily engages in hallmark functions (metabolism, immunological responses, proliferation, apoptosis, p53, HIF-1, FOXO, PI3K-AKT signaling). This work provides a wide knowledge of G6PD's function in lung cancer, as well as a theoretical foundation for possible prognostic therapeutic markers.

Genetic effects on the skin methylome in healthy older twins

Whole-skin DNA methylation variation has been implicated in several diseases, including melanoma, but its genetic basis has not yet been fully characterized. Using bulk skin tissue samples from 414 healthy female UK twins, we performed twin-based heritability and methylation quantitative trait loci (meQTL) analyses for >400,000 DNA methylation sites. We find that the human skin DNA methylome is on average less heritable than previously estimated in blood and other tissues (mean heritability: 10.02%). meQTL analysis identified local genetic effects influencing DNA methylation at 18.8% (76,442) of tested CpG sites, as well as 1,775 CpG sites associated with at least one distal genetic variant. As a functional follow-up, we performed skin expression QTL (eQTL) analyses in a partially overlapping sample of 604 female twins. Colocalization analysis identified over 3,500 shared genetic effects affecting thousands of CpG sites (10,067) and genes (4,475). Mediation analysis of putative colocalized gene-CpG pairs identified 114 genes with evidence for eQTL effects being mediated by DNA methylation in skin, including in genes implicating skin disease such as ALOX12 and CSPG4. We further explored the relevance of skin meQTLs to skin disease and found that skin meQTLs and CpGs under genetic influence were enriched for multiple skin-related genome-wide and epigenome-wide association signals, including for melanoma and psoriasis. Our findings give insights into the regulatory landscape of epigenomic variation in skin.

Population Diversity at the Single-Cell Level

Population-scale single-cell genomics is a transformative approach for unraveling the intricate links between genetic and cellular variation. This approach is facilitated by cutting-edge experimental methodologies, including the development of high-throughput single-cell multiomics and advances in multiplexed environmental and genetic perturbations. Examining the effects of natural or synthetic genetic variants across cellular contexts provides insights into the mutual influence of genetics and the environment in shaping cellular heterogeneity. The development of computational methodologies further enables detailed quantitative analysis of molecular variation, offering an opportunity to examine the respective roles of stochastic, intercellular, and interindividual variation. Future opportunities lie in leveraging long-read sequencing, refining disease-relevant cellular models, and embracing predictive and generative machine learning models. These advancements hold the potential for a deeper understanding of the genetic architecture of human molecular traits, which in turn has important implications for understanding the genetic causes of human disease.

Genome-Wide Methylation Analysis Reveals a KCNK3-Prominent Causal Cascade on Hypertension

BACKGROUND

Despite advances in understanding hypertension's genetic structure, how noncoding genetic variants influence it remains unclear. Studying their interaction with DNA methylation is crucial to deciphering this complex disease's genetic mechanisms.

METHODS

We investigated the genetic and epigenetic interplay in hypertension using whole-genome bisulfite sequencing. Methylation profiling in 918 males revealed allele-specific methylation and methylation quantitative trait loci. We engineered rs1275988T/C mutant mice using CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 (CRISPR-associated protein 9), bred them for homozygosity, and subjected them to a high-salt diet. Telemetry captured their cardiovascular metrics. Protein-DNA interactions were elucidated using DNA pull-downs, mass spectrometry, and Western blots. A wire myograph assessed vascular function, and analysis of the Kcnk3 gene methylation highlighted the mutation's role in hypertension.

RESULTS

We discovered that DNA methylation-associated genetic effects, especially in non-cytosine-phosphate-guanine (non-CpG) island and noncoding distal regulatory regions, significantly contribute to hypertension predisposition. We identified distinct methylation quantitative trait locus patterns in the hypertensive population and observed that the onset of hypertension is influenced by the transmission of genetic effects through the demethylation process. By evidence-driven prioritization and in vivo experiments, we unearthed rs1275988 in a cell type-specific enhancer as a notable hypertension causal variant, intensifying hypertension through the modulation of local DNA methylation and consequential alterations in Kcnk3 gene expression and vascular remodeling. When exposed to a high-salt diet, mice with the rs1275988C/C genotype exhibited exacerbated hypertension and significant vascular remodeling, underscored by increased aortic wall thickness. The C allele of rs1275988 was associated with elevated DNA methylation levels, driving down the expression of the Kcnk3 gene by attenuating Nr2f2 (nuclear receptor subfamily 2 group F member 2) binding at the enhancer locus.

CONCLUSIONS

Our research reveals new insights into the complex interplay between genetic variations and DNA methylation in hypertension. We underscore hypomethylation's potential in hypertension onset and identify rs1275988 as a causal variant in vascular remodeling. This work advances our understanding of hypertension's molecular mechanisms and encourages personalized health care strategies.

A systematic review on the contribution of DNA methylation to hearing loss

BACKGROUND

DNA methylation may have a regulatory role in monogenic sensorineural hearing loss and complex, polygenic phenotypic forms of hearing loss, including age-related hearing impairment or Meniere disease. The purpose of this systematic review is to critically assess the evidence supporting a functional role of DNA methylation in phenotypes associated with hearing loss.

RESULTS

The search strategy yielded a total of 661 articles. After quality assessment, 25 records were selected (12 human DNA methylation studies, 5 experimental animal studies and 8 studies reporting mutations in the DNMT1 gene). Although some methylation studies reported significant differences in CpG methylation in diverse gene promoters associated with complex hearing loss phenotypes (ARHI, otosclerosis, MD), only one study included a replication cohort that supported a regulatory role for CpG methylation in the genes TCF25 and POLE in ARHI. Conversely, several studies have independently confirmed pathogenic mutations within exon 21 of the DNMT1 gene, which encodes the DNA (cytosine-5)-methyltransferase 1 enzyme. This methylation enzyme is strongly associated with a rare disease defined by autosomal dominant cerebellar ataxia, deafness and narcolepsy (ADCA-DN). Of note, rare variants in DNMT1 and DNMT3A genes have also been reported in noise-induced hearing loss.

CONCLUSIONS

Evidence supporting a functional role for DNA methylation in hearing loss is limited to few genes in complex disorders such as ARHI. Mutations in the DNMT1 gene are associated with ADCA-DN, suggesting the CpG methylation in hearing loss genes deserves further attention in hearing research.

mLiftOver: harmonizing data across Infinium DNA methylation platforms

MOTIVATION

Infinium DNA methylation BeadChips are widely used for genome-wide DNA methylation profiling at the population scale. Recent updates to probe content and naming conventions in the EPIC version 2 (EPICv2) arrays have complicated integrating new data with previous Infinium array platforms, such as the MethylationEPIC (EPIC) and the HumanMethylation450 (HM450) BeadChip.

RESULTS

We present mLiftOver, a user-friendly tool that harmonizes probe ID, methylation level, and signal intensity data across different Infinium platforms. It manages probe replicates, missing data imputation, and platform-specific bias for accurate data conversion. We validated the tool by applying HM450-based cancer classifiers to EPICv2 cancer data, achieving high accuracy. Additionally, we successfully integrated EPICv2 healthy tissue data with legacy HM450 data for tissue identity analysis and produced consistent copy number profiles in cancer cells.

AVAILABILITY AND IMPLEMENTATION

mLiftOver is implemented R and available in the Bioconductor package SeSAMe (version 1.21.13+): https://bioconductor.org/packages/release/bioc/html/sesame.html. Analysis of EPIC and EPICv2 platform-specific bias and high-confidence mapping is available at https://github.com/zhou-lab/InfiniumAnnotationV1/raw/main/Anno/EPICv2/EPICv2ToEPIC_conversion.tsv.gz. The source code is available at https://github.com/zwdzwd/sesame/blob/devel/R/mLiftOver.R under the MIT license.

bsgenova: an accurate, robust, and fast genotype caller for bisulfite-sequencing data

BACKGROUND

Bisulfite sequencing (BS-Seq) is a fundamental technique for characterizing DNA methylation profiles. Genotype calling from bisulfite-converted BS-Seq data allows allele-specific methylation analysis and the concurrent exploration of genetic and epigenetic profiles. Despite various methods have been proposed, single nucleotide polymorphisms (SNPs) calling from BS-Seq data, particularly for SNPs on chromosome X and in the presence of contaminative data, poses ongoing challenges.

RESULTS

We introduce bsgenova, a novel SNP caller tailored for bisulfite sequencing data, employing a Bayesian multinomial model. The performance of bsgenova is assessed by comparing SNPs called from real-world BS-Seq data with those from corresponding whole-genome sequencing (WGS) data across three human cell lines. bsgenova is both sensitive and precise, especially for chromosome X, compared with three existing methods. Moreover, in the presence of low-quality reads, bsgenova outperforms other methods notably. In addition, bsgenova is meticulously implemented, leveraging matrix imputation and multi-process parallelization. Compared to existing methods, bsgenova stands out for its speed and efficiency in memory and disk usage. Furthermore, bsgenova integrates bsextractor, a methylation extractor, enhancing its flexibility and expanding its utility.

CONCLUSIONS

We introduce bsgenova for SNP calling from bisulfite-sequencing data. The source code is available at https://github.com/hippo-yf/bsgenova under license GPL-3.0.

Cascade testing for hereditary cancer in Singapore: how population genomics help guide clinical policy

Hereditary Cancer makes up around 5-10% of all cancers. It is important to diagnose hereditary cancer in a timely fashion, as not only do patients require long-term care from a young age, but their relatives also require management. The main approach to capture at-risk relatives is cascade testing. It involves genetic testing of relatives of the first detected carrier of a pathogenic variant in a family i.e. the proband. The current standard of care for cascade testing is a patient-mediated approach. Probands are then advised to inform and encourage family members to undergo genetic testing. In Singapore, cascade testing is inefficient, around 10-15%, lower than the 30% global average. Here, we describe the cascade testing process and its effort to increase testing in Singapore. Precision Health Research, Singapore (PRECISE), was set up to coordinate Singapore's National Precision Medicine strategy and has awarded five clinical implementation pilots, with one of them seeking to identify strategies for how cascade testing for hereditary cancer can be increased in a safe and cost-efficient manner. Achieving this will be done through addressing barriers such as cost, manpower shortages, exploring a digital channel for contacting at-risk relatives, and getting a deeper insight into why genetic testing gets declined. If successful, it will likely result in care pathways that are a cost-effective public health intervention for identifying individuals at risk. Surveillance and management of those unaffected at-risk individuals, if caught early, will result in improved patient outcomes, and further reduce the healthcare burden for the economy.

A review of the role of epigenetic studies for intramuscular fat deposition in beef cattle

Intramuscular fat (IMF) deposition profoundly influences meat quality and economic value in beef cattle production. Meanwhile, contemporary developments in epigenetics have opened new outlooks for understanding the molecular basics of IMF regulation, and it has become a key area of research for world scholars. Therefore, the aim of this paper was to provide insight and synthesis into the intricate relationship between epigenetic mechanisms and IMF deposition in beef cattle. The methodology involves a thorough analysis of existing literature, including pertinent books, academic journals, and online resources, to provide a comprehensive overview of the role of epigenetic studies in IMF deposition in beef cattle. This review summarizes the contemporary studies in epigenetic mechanisms in IMF regulation, high-resolution epigenomic mapping, single-cell epigenomics, multi-omics integration, epigenome editing approaches, longitudinal studies in cattle growth, environmental epigenetics, machine learning in epigenetics, ethical and regulatory considerations, and translation to industry practices from perspectives of IMF deposition in beef cattle. Moreover, this paper highlights DNA methylation, histone modifications, acetylation, phosphorylation, ubiquitylation, non-coding RNAs, DNA hydroxymethylation, epigenetic readers, writers, and erasers, chromatin immunoprecipitation followed by sequencing, whole genome bisulfite sequencing, epigenome-wide association studies, and their profound impact on the expression of crucial genes governing adipogenesis and lipid metabolism. Nutrition and stress also have significant influences on epigenetic modifications and IMF deposition. The key findings underscore the pivotal role of epigenetic studies in understanding and enhancing IMF deposition in beef cattle, with implications for precision livestock farming and ethical livestock management. In conclusion, this review highlights the crucial significance of epigenetic pathways and environmental factors in affecting IMF deposition in beef cattle, providing insightful information for improving the economics and meat quality of cattle production.

MSA: scalable DNA methylation screening BeadChip for high-throughput trait association studies

The Infinium DNA Methylation BeadChips have significantly contributed to population-scale epigenetics research by enabling epigenome-wide trait association discoveries. Here, we design, describe, and experimentally verify a new iteration of this technology, the Methylation Screening Array (MSA), to focus on human trait screening and discovery. This array utilizes extensive data from previous Infinium platform-based epigenome-wide association studies (EWAS). It incorporates knowledge from the latest single-cell and cell type-resolution whole genome methylome profiles. The MSA is engineered to achieve scalable screening of epigenetics-trait association in an ultra-high sample throughput. Our design encompassed diverse human trait associations, including those with genetic, cellular, environmental, and demographical variables and human diseases such as genetic, neurodegenerative, cardiovascular, infectious, and immune diseases. We comprehensively evaluated this array's reproducibility, accuracy, and capacity for cell-type deconvolution and supporting 5-hydroxymethylation profiling in diverse human tissues. Our first atlas data using this platform uncovered the complex chromatin and tissue contexts of DNA modification variations and genetic variants linked to human phenotypes.

Identification and characterization of whole blood gene expression and splicing quantitative trait loci during early to mid-lactation of dairy cattle

BACKGROUND

Characterization of regulatory variants (e.g., gene expression quantitative trait loci, eQTL; gene splicing QTL, sQTL) is crucial for biologically interpreting molecular mechanisms underlying loci associated with complex traits. However, regulatory variants in dairy cattle, particularly in specific biological contexts (e.g., distinct lactation stages), remain largely unknown. In this study, we explored regulatory variants in whole blood samples collected during early to mid-lactation (22-150 days after calving) of 101 Holstein cows and analyzed them to decipher the regulatory mechanisms underlying complex traits in dairy cattle.

RESULTS

We identified 14,303 genes and 227,705 intron clusters expressed in the white blood cells of 101 cattle. The average heritability of gene expression and intron excision ratio explained by cis-SNPs is 0.28 ± 0.13 and 0.25 ± 0.13, respectively. We identified 23,485 SNP-gene expression pairs and 18,166 SNP-intron cluster pairs in dairy cattle during early to mid-lactation. Compared with the 2,380,457 cis-eQTLs reported to be present in blood in the Cattle Genotype-Tissue Expression atlas (CattleGTEx), only 6,114 cis-eQTLs (P < 0.05) were detected in the present study. By conducting colocalization analysis between cis-e/sQTL and the results of genome-wide association studies (GWAS) from four traits, we identified a cis-e/sQTL (rs109421300) of the DGAT1 gene that might be a key marker in early to mid-lactation for milk yield, fat yield, protein yield, and somatic cell score (PP4 > 0.6). Finally, transcriptome-wide association studies (TWAS) revealed certain genes (e.g., FAM83H and TBC1D17) whose expression in white blood cells was significantly (P < 0.05) associated with complex traits.

CONCLUSIONS

This study investigated the genetic regulation of gene expression and alternative splicing in dairy cows during early to mid-lactation and provided new insights into the regulatory mechanisms underlying complex traits of economic importance.

Analysis of blood methylation quantitative trait loci in East Asians reveals ancestry-specific impacts on complex traits

Methylation quantitative trait loci (mQTLs) are essential for understanding the role of DNA methylation changes in genetic predisposition, yet they have not been fully characterized in East Asians (EAs). Here we identified mQTLs in whole blood from 3,523 Chinese individuals and replicated them in additional 1,858 Chinese individuals from two cohorts. Over 9% of mQTLs displayed specificity to EAs, facilitating the fine-mapping of EA-specific genetic associations, as shown for variants associated with height. Trans-mQTL hotspots revealed biological pathways contributing to EA-specific genetic associations, including an ERG-mediated 233 trans-mCpG network, implicated in hematopoietic cell differentiation, which likely reflects binding efficiency modulation of the ERG protein complex. More than 90% of mQTLs were shared between different blood cell lineages, with a smaller fraction of lineage-specific mQTLs displaying preferential hypomethylation in the respective lineages. Our study provides new insights into the mQTL landscape across genetic ancestries and their downstream effects on cellular processes and diseases/traits.

Cellular senescence gene TACC3 associated with colorectal cancer risk via genetic and DNA methylated alteration

Cell senescence genes play a vital role in the pathogenesis of colorectal cancer, a process that may involve the triggering of genetic variations and reversible phenotypes caused by epigenetic modifications. However, the specific regulatory mechanisms remain unclear. Using CellAge and The Cancer Genome Atlas databases and in-house RNA-seq data, DNA methylation-modified cellular senescence genes (DMCSGs) were validated by Support Vector Machine and correlation analyses. In 1150 cases and 1342 controls, we identified colorectal cancer risk variants in DMCSGs. The regulatory effects of gene, variant, and DNA methylation were explored through dual-luciferase and 5-azacytidine treatment experiments, complemented by multiple database analyses. Biological functions of key gene were evaluated via cell proliferation assays, SA-β-gal staining, senescence marker detection, and immune infiltration analyses. The genetic variant rs4558926 in the downstream of TACC3 was significantly associated with colorectal cancer risk (OR = 1.35, P = 3.22 × 10-4). TACC3 mRNA expression increased due to rs4558926 C > G and decreased DNA methylation levels. The CpG sites in the TACC3 promoter region were regulated by rs4558926. TACC3 knockdown decreased proliferation and senescence in colorectal cancer cells. In addition, subjects with high-TACC3 expression presented an immunosuppressive microenvironment. These findings provide insights into the involvement of genetic variants of cellular senescence genes in the development and progression of colorectal cancer.

Low-input and single-cell methods for Infinium DNA methylation BeadChips

The Infinium BeadChip is the most widely used DNA methylome assay technology for population-scale epigenome profiling. However, the standard workflow requires over 200 ng of input DNA, hindering its application to small cell-number samples, such as primordial germ cells. We developed experimental and analysis workflows to extend this technology to suboptimal input DNA conditions, including ultra-low input down to single cells. DNA preamplification significantly enhanced detection rates to over 50% in five-cell samples and ∼25% in single cells. Enzymatic conversion also substantially improved data quality. Computationally, we developed a method to model the background signal's influence on the DNA methylation level readings. The modified detection P-value calculation achieved higher sensitivities for low-input datasets and was validated in over 100 000 public diverse methylome profiles. We employed the optimized workflow to query the demethylation dynamics in mouse primordial germ cells available at low cell numbers. Our data revealed nuanced chromatin states, sex disparities, and the role of DNA methylation in transposable element regulation during germ cell development. Collectively, we present comprehensive experimental and computational solutions to extend this widely used methylation assay technology to applications with limited DNA.

Conserved methylation signatures associate with the tumor immune microenvironment and immunotherapy response

BACKGROUND

Aberrant DNA methylation is a major characteristic of cancer genomes. It remains unclear which biological processes determine epigenetic reprogramming and how these processes influence the variants in the cancer methylome, which can further impact cancer phenotypes.

METHODS

We performed pairwise permutations of 381,900 loci in 569 paired DNA methylation profiles of cancer tissue and matched normal tissue from The Cancer Genome Atlas (TCGA) and defined conserved differentially methylated positions (DMPs) based on the resulting null distribution. Then, we derived independent methylation signatures from 2,465 cancer-only methylation profiles from the TCGA and 241 cell line-based methylation profiles from the Genomics of Drug Sensitivity in Cancer (GDSC) cohort using nonnegative matrix factorization (NMF). We correlated DNA methylation signatures with various clinical and biological features, including age, survival, cancer stage, tumor immune microenvironment factors, and immunotherapy response. We inferred the determinant genes of these methylation signatures by integrating genomic and transcriptomic data and evaluated the impact of these signatures on cancer phenotypes in independent bulk and single-cell RNA/methylome cohorts.

RESULTS

We identified 7,364 differentially methylated positions (2,969 Hyper-DMPs and 4,395 Hypo-DMPs) in nine cancer types from the TCGA. We subsequently retrieved three highly conserved, independent methylation signatures (Hyper-MS1, Hypo-MS1, and Hypo-MS4) from cancer tissues and cell lines based on these Hyper and Hypo-DMPs. Our data suggested that Hypo-MS4 activity predicts poor survival and is associated with immunotherapy response and distant tumor metastasis, and Hypo-MS4 activity is related to TP53 mutation and FOXA1 binding specificity. In addition, we demonstrated a correlation between the activities of Hypo-MS4 in cancer cells and the fractions of regulatory CD4 + T cells with the expression levels of immunological genes in the tumor immune microenvironment.

CONCLUSIONS

Our findings demonstrated that the methylation signatures of distinct biological processes are associated with immune activity in the cancer microenvironment and predict immunotherapy response.

Genetic control of DNA methylation is largely shared across European and East Asian populations

DNA methylation is an ideal trait to study the extent of the shared genetic control across ancestries, effectively providing hundreds of thousands of model molecular traits with large QTL effect sizes. We investigate cis DNAm QTLs in three European (n = 3701) and two East Asian (n = 2099) cohorts to quantify the similarities and differences in the genetic architecture across populations. We observe 80,394 associated mQTLs (62.2% of DNAm probes with significant mQTL) to be significant in both ancestries, while 28,925 mQTLs (22.4%) are identified in only a single ancestry. mQTL effect sizes are highly conserved across populations, with differences in mQTL discovery likely due to differences in allele frequency of associated variants and differing linkage disequilibrium between causal variants and assayed SNPs. This study highlights the overall similarity of genetic control across ancestries and the value of ancestral diversity in increasing the power to detect associations and enhancing fine mapping resolution.

Alternative polyadenylation quantitative trait methylation mapping in human cancers provides clues into the molecular mechanisms of APA

Genetic variants are involved in the orchestration of alternative polyadenylation (APA) events, while the role of DNA methylation in regulating APA remains unclear. We generated a comprehensive atlas of APA quantitative trait methylation sites (apaQTMs) across 21 different types of cancer (1,612 to 60,219 acting in cis and 4,448 to 142,349 in trans). Potential causal apaQTMs in non-cancer samples were also identified. Mechanistically, we observed a strong enrichment of cis-apaQTMs near polyadenylation sites (PASs) and both cis- and trans-apaQTMs in proximity to transcription factor (TF) binding regions. Through the integration of ChIP-signals and RNA-seq data from cell lines, we have identified several regulators of APA events, acting either directly or indirectly, implicating novel functions of some important genes, such as TCF7L2, which is known for its involvement in type 2 diabetes and cancers. Furthermore, we have identified a vast number of QTMs that share the same putative causal CpG sites with five different cancer types, underscoring the roles of QTMs, including apaQTMs, in the process of tumorigenesis. DNA methylation is extensively involved in the regulation of APA events in human cancers. In an attempt to elucidate the potential underlying molecular mechanisms of APA by DNA methylation, our study paves the way for subsequent experimental validations into the intricate biological functions of DNA methylation in APA regulation and the pathogenesis of human cancers. To present a comprehensive catalog of apaQTM patterns, we introduce the Pancan-apaQTM database, available at https://pancan-apaqtm-zju.shinyapps.io/pancanaQTM/.

QTL mapping of human retina DNA methylation identifies 87 gene-epigenome interactions in age-related macular degeneration

DNA methylation provides a crucial epigenetic mark linking genetic variations to environmental influence. We have analyzed array-based DNA methylation profiles of 160 human retinas with co-measured RNA-seq and >8 million genetic variants, uncovering sites of genetic regulation in cis (37,453 methylation quantitative trait loci and 12,505 expression quantitative trait loci) and 13,747 DNA methylation loci affecting gene expression, with over one-third specific to the retina. Methylation and expression quantitative trait loci show non-random distribution and enrichment of biological processes related to synapse, mitochondria, and catabolism. Summary data-based Mendelian randomization and colocalization analyses identify 87 target genes where methylation and gene-expression changes likely mediate the genotype effect on age-related macular degeneration. Integrated pathway analysis reveals epigenetic regulation of immune response and metabolism including the glutathione pathway and glycolysis. Our study thus defines key roles of genetic variations driving methylation changes, prioritizes epigenetic control of gene expression, and suggests frameworks for regulation of macular degeneration pathology by genotype-environment interaction in retina.

New insights into the role of immunity and inflammation in diabetic kidney disease in the omics era

Diabetic kidney disease (DKD) is becoming the leading cause of chronic kidney disease, especially in the industrialized world. Despite mounting evidence has demonstrated that immunity and inflammation are highly involved in the pathogenesis and progression of DKD, the underlying mechanisms remain incompletely understood. Substantial molecules, signaling pathways, and cell types participate in DKD inflammation, by integrating into a complex regulatory network. Most of the studies have focused on individual components, without presenting their importance in the global or system-based processes, which largely hinders clinical translation. Besides, conventional technologies failed to monitor the different behaviors of resident renal cells and immune cells, making it difficult to understand their contributions to inflammation in DKD. Recently, the advancement of omics technologies including genomics, epigenomics, transcriptomics, proteomics, and metabolomics has revolutionized biomedical research, which allows an unbiased global analysis of changes in DNA, RNA, proteins, and metabolites in disease settings, even at single-cell and spatial resolutions. They help us to identify critical regulators of inflammation processes and provide an overview of cell heterogeneity in DKD. This review aims to summarize the application of multiple omics in the field of DKD and emphasize the latest evidence on the interplay of inflammation and DKD revealed by these technologies, which will provide new insights into the role of inflammation in the pathogenesis of DKD and lead to the development of novel therapeutic approaches and diagnostic biomarkers.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects

Bulk-tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, and context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from the blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell-type proportions, we demonstrate that cell-type iQTLs could be considered as proxies for cell-type-specific QTL effects, particularly for the most abundant cell type in the tissue. The interpretation of age iQTLs, however, warrants caution because the moderation effect of age on the genotype and molecular phenotype association could be mediated by changes in cell-type composition. Finally, we show that cell-type iQTLs contribute to cell-type-specific enrichment of diseases that, in combination with additional functional data, could guide future functional studies. Overall, this study highlights the use of iQTLs to gain insights into the context specificity of regulatory effects.

Genome-wide QTL mapping across three tissues highlights several Alzheimer's and Parkinson's disease loci potentially acting via DNA methylation

DNA methylation (DNAm) is an epigenetic mark with essential roles in disease development and predisposition. Here, we created genome-wide maps of methylation quantitative trait loci (meQTL) in three peripheral tissues and used Mendelian randomization (MR) analyses to assess the potential causal relationships between DNAm and risk for two common neurodegenerative disorders, i.e. Alzheimer's disease (AD) and Parkinson's disease (PD). Genome-wide single nucleotide polymorphism (SNP; ~5.5M sites) and DNAm (~850K CpG sites) data were generated from whole blood (n=1,058), buccal (n=1,527) and saliva (n=837) specimens. We identified between 11 and 15 million genome-wide significant (p<10-14) SNP-CpG associations in each tissue. Combining these meQTL GWAS results with recent AD/PD GWAS summary statistics by MR strongly suggests that the previously described associations between PSMC3, PICALM, and TSPAN14 and AD may be founded on differential DNAm in or near these genes. In addition, there is strong, albeit less unequivocal, support for causal links between DNAm at PRDM7 in AD as well as at KANSL1/MAPT in AD and PD. Our study adds valuable insights on AD/PD pathogenesis by combining two high-resolution "omics" domains, and the meQTL data shared along with this publication will allow like-minded analyses in other diseases.

Pleiotropic influence of DNA methylation QTLs on physiological and ageing traits

DNA methylation is influenced by genetic and non-genetic factors. Here, we chart quantitative trait loci (QTLs) that modulate levels of methylation at highly conserved CpGs using liver methylome data from mouse strains belonging to the BXD family. A regulatory hotspot on chromosome 5 had the highest density of trans-acting methylation QTLs (trans-meQTLs) associated with multiple distant CpGs. We refer to this locus as meQTL.5a. Trans-modulated CpGs showed age-dependent changes and were enriched in developmental genes, including several members of the MODY pathway (maturity onset diabetes of the young). The joint modulation by genotype and ageing resulted in a more 'aged methylome' for BXD strains that inherited the DBA/2J parental allele at meQTL.5a. Further, several gene expression traits, body weight, and lipid levels mapped to meQTL.5a, and there was a modest linkage with lifespan. DNA binding motif and protein-protein interaction enrichment analyses identified the hepatic nuclear factor, Hnf1a (MODY3 gene in humans), as a strong candidate. The pleiotropic effects of meQTL.5a could contribute to variations in body size and metabolic traits, and influence CpG methylation and epigenetic ageing that could have an impact on lifespan.

Circulating tumor DNA as liquid biopsy in lung cancer: Biological characteristics and clinical integration

Lung cancer maintains high morbidity and mortality rate globally despite significant advancements in diagnosis and treatment in the era of precision medicine. Pathological analysis of tumor tissue, the current gold standard for lung cancer diagnosis, is intrusive and intrinsically confined to evaluating the limited amount of tissues that could be physically extracted. However, tissue biopsy has several limitations, including the invasiveness of the procedure and difficulty in obtaining samples for patients at advanced stages., there Additionally,has been no major breakthrough in tumor biomarkers with high specificity and sensitivity, particularly for early-stage lung cancer. Liquid biopsy has been considered a feasible auxiliary tool for tearly dianosis, evaluating treatment responses and monitoring prognosis of lung cancer. Circulating tumor DNA (ctDNA), an ideal biomarker of liquid biopsy, has emerged as one of the most reliable tools for monitoring tumor processes at molecular levels. Herein, this review focuses on tumor heterogeneity to elucidate the superiority of liquid biopsy and retrospectively discussdeciphersolution. We systematically elaborate ctDNA biological characteristics, introduce methods for ctDNA detection, and discuss the current role of plasma ctDNA in lung cancer management. Finally, we summarize the drawbacks of ctDNA analysis and highlight its potential clinical application in lung cancer.

Epigenetics insights from perceived facial aging

Facial aging is the most visible manifestation of aging. People desire to look younger than others of the same chronological age. Hence, perceived age is often used as a visible marker of aging, while biological age, often estimated by methylation markers, is used as an objective measure of age. Multiple epigenetics-based clocks have been developed for accurate estimation of general biological age and the age of specific organs, including the skin. However, it is not clear whether the epigenetic biomarkers (CpGs) used in these clocks are drivers of aging processes or consequences of aging. In this proof-of-concept study, we integrate data from GWAS on perceived facial aging and EWAS on CpGs measured in blood. By running EW Mendelian randomization, we identify hundreds of putative CpGs that are potentially causal to perceived facial aging with similar numbers of damaging markers that causally drive or accelerate facial aging and protective methylation markers that causally slow down or protect from aging. We further demonstrate that while candidate causal CpGs have little overlap with known epigenetics-based clocks, they affect genes or proteins with known functions in skin aging, such as skin pigmentation, elastin, and collagen levels. Overall, our results suggest that blood methylation markers reflect facial aging processes, and thus can be used to quantify skin aging and develop anti-aging solutions that target the root causes of aging.

Investigation of the impact of biogeographic ancestry on DNA methylation based age predictions comparing a Middle East and a Central European population

DNA methylation based age prediction is a new method in the toolbox of forensic genetics. Typically, the method is applied in the course of police investigation e.g. to predict the age of an unknown person that has left a biological trace at a crime scene. The method can also be used to answer other forensic questions, for example to estimate the age of unknown human bodies in the course of the identification process. In the present study, we tested for a potential impact of biogeographic ancestry (BGA) on age predictions using five age dependent methylated CpG sites within the genetic regions of ELOVL2, MIR29B2CHG, FHL2, KLF14 and TRIM59. We collected 102 blood samples each from donors living in Iraq, Middle East (ME) and Germany, Central Europe (EU). Both sample sets were matched in sex and age ranging from 18 to 68 years with exactly one male and female sample per year of age. All samples were analyzed by bisulfite pyrosequencing applying a multiplex pre-amplification strategy based on a single input of 35 ng converted DNA in the PCR. For the CpGs in MIR29B2CHG, FHL2 and KLF14, we observed significantly different methylation levels between the two populations. While we were able to train two highly accurate prediction models for the respective population with mean absolute deviations between predicted and actual ages (MAD) of 3.34 years for the ME model, and 2.72 years for the EU model, we found an absolute prediction difference between the two population specific models of more than 4 years. A combined model for both populations compensated the methylation difference between the two populations, providing MADs of prediction of only 3.81 years for ME and 3.31 years for EU samples. In total, the results of the present study strongly support the benefit of BGA information for more reliable methylation based age predictions.

Polymorphic variants involved in methylation regulation: a strategy to discover risk loci for pancreatic ductal adenocarcinoma

INTRODUCTION

Only a small number of risk factors for pancreatic ductal adenocarcinoma (PDAC) has been established. Several studies identified a role of epigenetics and of deregulation of DNA methylation. DNA methylation is variable across a lifetime and in different tissues; nevertheless, its levels can be regulated by genetic variants like methylation quantitative trait loci (mQTLs), which can be used as a surrogate.

MATERIALS AND METHODS

We scanned the whole genome for mQTLs and performed an association study in 14 705 PDAC cases and 246 921 controls. The methylation data were obtained from whole blood and pancreatic cancer tissue through online databases. We used the Pancreatic Cancer Cohort Consortium and the Pancreatic Cancer Case-Control Consortium genome-wide association study (GWAS) data as discovery phase and the Pancreatic Disease Research consortium, the FinnGen project and the Japan Pancreatic Cancer Research consortium GWAS as replication phase.

RESULTS

The C allele of 15q26.1-rs12905855 showed an association with a decreased risk of PDAC (OR=0.90, 95% CI 0.87 to 0.94, p=4.93×10-8 in the overall meta-analysis), reaching genome-level statistical significance. 15q26.1-rs12905855 decreases the methylation of a 'C-phosphate-G' (CpG) site located in the promoter region of the RCCD1 antisense (RCCD1-AS1) gene which, when expressed, decreases the expression of the RCC1 domain-containing (RCCD1) gene (part of a histone demethylase complex). Thus, it is possible that the rs12905855 C-allele has a protective role in PDAC development through an increase of RCCD1 gene expression, made possible by the inactivity of RCCD1-AS1.

CONCLUSION

We identified a novel PDAC risk locus which modulates cancer risk by controlling gene expression through DNA methylation.

Low-input and single-cell methods for Infinium DNA methylation BeadChips

The Infinium BeadChip is the most widely used DNA methylome assay technology for population-scale epigenome profiling. However, the standard workflow requires over 200 ng of input DNA, hindering its application to small cell-number samples, such as primordial germ cells. We developed experimental and analysis workflows to extend this technology to suboptimal input DNA conditions, including ultra-low input down to single cells. DNA preamplification significantly enhanced detection rates to over 50% in five-cell samples and ∼25% in single cells. Enzymatic conversion also substantially improved data quality. Computationally, we developed a method to model the background signal's influence on the DNA methylation level readings. The modified detection p -values calculation achieved higher sensitivities for low-input datasets and was validated in over 100,000 public datasets with diverse methylation profiles. We employed the optimized workflow to query the demethylation dynamics in mouse primordial germ cells available at low cell numbers. Our data revealed nuanced chromatin states, sex disparities, and the role of DNA methylation in transposable element regulation during germ cell development. Collectively, we present comprehensive experimental and computational solutions to extend this widely used methylation assay technology to applications with limited DNA.

Blood-based genome-wide DNA methylation correlations across body-fat- and adiposity-related biochemical traits

The recent increase in obesity levels across many countries is likely to be driven by nongenetic factors. The epigenetic modification DNA methylation (DNAm) may help to explore this, as it is sensitive to both genetic and environmental exposures. While the relationship between DNAm and body-fat traits has been extensively studied, there is limited literature on the shared associations of DNAm variation across such traits. Akin to genetic correlation estimates, here, we introduce an approach to evaluate the similarities in DNAm associations between traits: DNAm correlations. As DNAm can be both a cause and consequence of complex traits, DNAm correlations have the potential to provide insights into trait relationships above that currently obtained from genetic and phenotypic correlations. Utilizing 7,519 unrelated individuals from Generation Scotland with DNAm from the EPIC array, we calculated DNAm correlations between body-fat- and adiposity-related traits by using the bivariate OREML framework in the OSCA software. For each trait, we also estimated the shared contribution of DNAm between sexes. We identified strong, positive DNAm correlations between each of the body-fat traits (BMI, body-fat percentage, and waist-to-hip ratio, ranging from 0.96 to 1.00), finding larger associations than those identified by genetic and phenotypic correlations. We identified a significant deviation from 1 in the DNAm correlations for BMI between males and females, with sex-specific DNAm changes associated with BMI identified at eight DNAm probes. Employing genome-wide DNAm correlations to evaluate the similarities in the associations of DNAm with complex traits has provided insight into obesity-related traits beyond that provided by genetic correlations.

Pan-cancer analysis of G6PD carcinogenesis in human tumors

Glucose-6-phosphate dehydrogenase (G6PD) is involved in the catalytic pentose phosphate pathway (PPP), which is closely related to energy metabolism. G6PD plays a crucial role in many types of cancer, but the specific molecular mechanisms of G6PD in cancer remain unclear. Therefore, we investigated the potential oncogenic role of G6PD in various tumors based on The Cancer Genome Atlas (TCGA), the cBioPortal datasets, the University of California Santa Cruz (UCSC) Xena browser, and the UALCAN-based online tool. G6PD was highly expressed in several cancer tissues (hepatocellular carcinoma, glioma, and breast cancer) compared with normal tissues and was significantly associated with poor prognosis of hepatocellular carcinoma, clear cell renal cell carcinoma, and breast cancer. Promoter methylation levels of G6PD were lower in Bladder Urothelial Carcinoma (BLCA) (P = 2.77e-02), breast invasive carcinoma (BRCA) (P = 1.62e-12), kidney renal clear cell carcinoma (KIRC) (P = 4.23e-02), kidney renal papillary cell carcinoma (KIRP) (P = 2.64e-03), liver hepatocellular carcinoma (LIHC) (P = 1.76e-02), stomach adenocarcinoma (STAD) (P = 3.50e-02), testicular germ cell tumors (TGCT) (P = 1.62e-12), higher in prostate adenocarcinoma (PRAD) (P = 1.81e-09), and uterine corpus endometrial carcinoma (UCEC) (P = 2.96e-04) compared with corresponding normal tissue samples. G6PD expression was positively correlated with the infiltration level of immune cells in most tumors, suggesting that G6PD may be involved in tumor immune infiltration. In addition, the functional mechanism of G6PD also involves 'Carbon metabolism', 'Glycolysis/Gluconeogenesis', 'Pentose phosphate pathway', and 'Central carbon pathway metabolism in cancer signaling pathway'. This pan-cancer study provides a relatively broad understanding of the oncogenic role of G6PD in various tumors and presents a theoretical basis for the development of G6PD inhibitors as therapeutic drugs for multiple cancers.

MicroRNAs as biomarkers for early diagnosis, targeting and prognosis of prostate cancer

Globally, prostate cancer (PC) is leading cause of cancer-related mortality in men worldwide. Despite significant advances in the treatment and management of this disease, the cure rates for PC remains low, largely due to late detection. PC detection is mostly reliant on prostate-specific antigen (PSA) and digital rectal examination (DRE); however, due to the low positive predictive value of current diagnostics, there is an urgent need to identify new accurate biomarkers. Recent studies support the biological role of microRNAs (miRNAs) in the initiation and progression of PC, as well as their potential as novel biomarkers for patients' diagnosis, prognosis, and disease relapse. In the advanced stages, cancer-cell-derived small extracellular vesicles (SEVs) may constitute a significant part of circulating vesicles and cause detectable changes in the plasma vesicular miRNA profile. Recent computational model for the identification of miRNA biomarkers discussed. In addition, accumulating evidence indicates that miRNAs can be utilized to target PC cells. In this article, the current understanding of the role of microRNAs and exosomes in the pathogenesis and their significance in PC prognosis, early diagnosis, chemoresistance, and treatment are reviewed.

Genetic impacts on DNA methylation help elucidate regulatory genomic processes

BACKGROUND

Pinpointing genetic impacts on DNA methylation can improve our understanding of pathways that underlie gene regulation and disease risk.

RESULTS

We report heritability and methylation quantitative trait locus (meQTL) analysis at 724,499 CpGs profiled with the Illumina Infinium MethylationEPIC array in 2358 blood samples from three UK cohorts. Methylation levels at 34.2% of CpGs are affected by SNPs, and 98% of effects are cis-acting or within 1 Mbp of the tested CpG. Our results are consistent with meQTL analyses based on the former Illumina Infinium HumanMethylation450 array. Both SNPs and CpGs with meQTLs are overrepresented in enhancers, which have improved coverage on this platform compared to previous approaches. Co-localisation analyses across genetic effects on DNA methylation and 56 human traits identify 1520 co-localisations across 1325 unique CpGs and 34 phenotypes, including in disease-relevant genes, such as USP1 and DOCK7 (total cholesterol levels), and ICOSLG (inflammatory bowel disease). Enrichment analysis of meQTLs and integration with expression QTLs give insights into mechanisms underlying cis-meQTLs (e.g. through disruption of transcription factor binding sites for CTCF and SMC3) and trans-meQTLs (e.g. through regulating the expression of ACD and SENP7 which can modulate DNA methylation at distal sites).

CONCLUSIONS

Our findings improve the characterisation of the mechanisms underlying DNA methylation variability and are informative for prioritisation of GWAS variants for functional follow-ups. The MeQTL EPIC Database and viewer are available online at https://epicmeqtl.kcl.ac.uk .

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects

Bulk tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, while context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell type proportions, we demonstrate that cell type iQTLs could be considered as proxies for cell type-specific QTL effects. The interpretation of age iQTLs, however, warrants caution as the moderation effect of age on the genotype and molecular phenotype association may be mediated by changes in cell type composition. Finally, we show that cell type iQTLs contribute to cell type-specific enrichment of diseases that, in combination with additional functional data, may guide future functional studies. Overall, this study highlights iQTLs to gain insights into the context-specificity of regulatory effects.

QTL mapping of human retina DNA methylation identifies 87 gene-epigenome interactions in age-related macular degeneration

DNA methylation (DNAm) provides a crucial epigenetic mark linking genetic variations to environmental influence. We analyzed array-based DNAm profiles of 160 human retinas with co-measured RNA-seq and > 8 million genetic variants, uncovering sites of genetic regulation in cis (37,453 mQTLs and 12,505 eQTLs) and 13,747 eQTMs (DNAm loci affecting gene expression), with over one-third specific to the retina. mQTLs and eQTMs show non-random distribution and enrichment of biological processes related to synapse, mitochondria, and catabolism. Summary data-based Mendelian randomization and colocalization analyses identify 87 target genes where methylation and gene-expression changes likely mediate the genotype effect on age-related macular degeneration (AMD). Integrated pathway analysis reveals epigenetic regulation of immune response and metabolism including the glutathione pathway and glycolysis. Our study thus defines key roles of genetic variations driving methylation changes, prioritizes epigenetic control of gene expression, and suggests frameworks for regulation of AMD pathology by genotype-environment interaction in retina.

Variation in histone configurations correlates with gene expression across nine inbred strains of mice

The Diversity Outbred (DO) mice and their inbred founders are widely used models of human disease. However, although the genetic diversity of these mice has been well documented, their epigenetic diversity has not. Epigenetic modifications, such as histone modifications and DNA methylation, are important regulators of gene expression and, as such, are a critical mechanistic link between genotype and phenotype. Therefore, creating a map of epigenetic modifications in the DO mice and their founders is an important step toward understanding mechanisms of gene regulation and the link to disease in this widely used resource. To this end, we performed a strain survey of epigenetic modifications in hepatocytes of the DO founders. We surveyed four histone modifications (H3K4me1, H3K4me3, H3K27me3, and H3K27ac), as well as DNA methylation. We used ChromHMM to identify 14 chromatin states, each of which represents a distinct combination of the four histone modifications. We found that the epigenetic landscape is highly variable across the DO founders and is associated with variation in gene expression across strains. We found that epigenetic state imputed into a population of DO mice recapitulated the association with gene expression seen in the founders, suggesting that both histone modifications and DNA methylation are highly heritable mechanisms of gene expression regulation. We illustrate how DO gene expression can be aligned with inbred epigenetic states to identify putative cis-regulatory regions. Finally, we provide a data resource that documents strain-specific variation in the chromatin state and DNA methylation in hepatocytes across nine widely used strains of laboratory mice.

meQTL mapping in the GENOA study reveals genetic determinants of DNA methylation in African Americans

Identifying genetic variants that are associated with variation in DNA methylation, an analysis commonly referred to as methylation quantitative trait locus (meQTL) mapping, is an important first step towards understanding the genetic architecture underlying epigenetic variation. Most existing meQTL mapping studies have focused on individuals of European ancestry and are underrepresented in other populations, with a particular absence of large studies in populations with African ancestry. We fill this critical knowledge gap by performing a large-scale cis-meQTL mapping study in 961 African Americans from the Genetic Epidemiology Network of Arteriopathy (GENOA) study. We identify a total of 4,565,687 cis-acting meQTLs in 320,965 meCpGs. We find that 45% of meCpGs harbor multiple independent meQTLs, suggesting potential polygenic genetic architecture underlying methylation variation. A large percentage of the cis-meQTLs also colocalize with cis-expression QTLs (eQTLs) in the same population. Importantly, the identified cis-meQTLs explain a substantial proportion (median = 24.6%) of methylation variation. In addition, the cis-meQTL associated CpG sites mediate a substantial proportion (median = 24.9%) of SNP effects underlying gene expression. Overall, our results represent an important step toward revealing the co-regulation of methylation and gene expression, facilitating the functional interpretation of epigenetic and gene regulation underlying common diseases in African Americans.

Comprehensive analysis of ZNF family genes in prognosis, immunity, and treatment of esophageal cancer

BACKGROUND

As a common malignant tumor, esophageal carcinoma (ESCA) has a low early diagnosis rate and poor prognosis. This study aimed to construct the prognostic features composed of ZNF family genes to effectively predict the prognosis of ESCA patients.

METHODS

The mRNA expression matrix and clinical data were downloaded from TCGA and GEO database. Using univariate Cox analysis, lasso regression and multivariate Cox analysis, we screened six prognosis-related ZNF family genes to construct the prognostic model. We then used Kaplan-Meier plot, time-dependent receiver operating characteristic (ROC), multivariable Cox regression analysis of clinical information, and nomogram to evaluate the prognostic value within and across sets, separately and combined. We also validated the prognostic value of the six-gene signature using GSE53624 dataset. The different immune status was observed in the single sample Gene Set Enrichment Analysis (ssGSEA). Finally, real-time quantitative PCR was used to detect the expression of six prognostic ZNF genes in twelve pairs of ESCA and adjacent normal tissues.

RESULTS

A six prognosis-related ZNF family genes model consisted of ZNF91, ZNF586, ZNF502, ZNF865, ZNF106 and ZNF225 was identified. Multivariable Cox regression analysis revealed that six prognosis-related ZNF family genes were independent prognostic factors for overall survival of ESCA patients in TCGA and GSE53624. Further, a prognostic nomogram including the riskScore, age, gender, T, stage was constructed, and TCGA/GSE53624-based calibration plots indicated its excellent predictive performance. Drug Sensitivity and ssGSEA analysis showed that the six genes model was closely related to immune cells infiltration and could be used as a potential predictor of chemotherapy sensitivity.

CONCLUSION

We identified six prognosis-related ZNF family genes model of ESCA, which provide evidence for individualized prevention and treatment.

The Biology of Stress Intolerance in Patients with Chronic Pain—State of the Art and Future Directions

Stress has been consistently linked to negative impacts on physical and mental health. More specifically, patients with chronic pain experience stress intolerance, which is an exacerbation or occurrence of symptoms in response to any type of stress. The pathophysiological mechanisms underlying this phenomenon remain unsolved. In this state-of-the-art paper, we summarised the role of the autonomic nervous system (ANS) and hypothalamus-pituitary-adrenal (HPA) axis, the two major stress response systems in stress intolerance. We provided insights into such mechanisms based on evidence from clinical studies in both patients with chronic pain, showing dysregulated stress systems, and healthy controls supported by preclinical studies, highlighting the link between these systems and symptoms of stress intolerance. Furthermore, we explored the possible regulating role for (epi)genetic mechanisms influencing the ANS and HPA axis. The link between stress and chronic pain has become an important area of research as it has the potential to inform the development of interventions to improve the quality of life for individuals living with chronic pain. As stress has become a prevalent concern in modern society, understanding the connection between stress, HPA axis, ANS, and chronic health conditions such as chronic pain is crucial to improve public health and well-being.

3D genomics and its applications in precision medicine

Three-dimensional (3D) genomics is an emerging discipline that studies the three-dimensional structure of chromatin and the three-dimensional and functions of genomes. It mainly studies the three-dimensional conformation and functional regulation of intranuclear genomes, such as DNA replication, DNA recombination, genome folding, gene expression regulation, transcription factor regulation mechanism, and the maintenance of three-dimensional conformation of genomes. Self-chromosomal conformation capture (3C) technology has been developed, and 3D genomics and related fields have developed rapidly. In addition, chromatin interaction analysis techniques developed by 3C technologies, such as paired-end tag sequencing (ChIA-PET) and whole-genome chromosome conformation capture (Hi-C), enable scientists to further study the relationship between chromatin conformation and gene regulation in different species. Thus, the spatial conformation of plant, animal, and microbial genomes, transcriptional regulation mechanisms, interaction patterns of chromosomes, and the formation mechanism of spatiotemporal specificity of genomes are revealed. With the help of new experimental technologies, the identification of key genes and signal pathways related to life activities and diseases is sustaining the rapid development of life science, agriculture, and medicine. In this paper, the concept and development of 3D genomics and its application in agricultural science, life science, and medicine are introduced, which provides a theoretical basis for the study of biological life processes.

Comparative analysis of genome-scale, base-resolution DNA methylation profiles across 580 animal species

Methylation of cytosines is a prototypic epigenetic modification of the DNA. It has been implicated in various regulatory mechanisms across the animal kingdom and particularly in vertebrates. We mapped DNA methylation in 580 animal species (535 vertebrates, 45 invertebrates), resulting in 2443 genome-scale DNA methylation profiles of multiple organs. Bioinformatic analysis of this large dataset quantified the association of DNA methylation with the underlying genomic DNA sequence throughout vertebrate evolution. We observed a broadly conserved link with two major transitions-once in the first vertebrates and again with the emergence of reptiles. Cross-species comparisons focusing on individual organs supported a deeply conserved association of DNA methylation with tissue type, and cross-mapping analysis of DNA methylation at gene promoters revealed evolutionary changes for orthologous genes. In summary, this study establishes a large resource of vertebrate and invertebrate DNA methylomes, it showcases the power of reference-free epigenome analysis in species for which no reference genomes are available, and it contributes an epigenetic perspective to the study of vertebrate evolution.

DNA methylation QTL mapping across diverse human tissues provides molecular links between genetic variation and complex traits

Studies of DNA methylation (DNAm) in solid human tissues are relatively scarce; tissue-specific characterization of DNAm is needed to understand its role in gene regulation and its relevance to complex traits. We generated array-based DNAm profiles for 987 human samples from the Genotype-Tissue Expression (GTEx) project, representing 9 tissue types and 424 subjects. We characterized methylome and transcriptome correlations (eQTMs), genetic regulation in cis (mQTLs and eQTLs) across tissues and e/mQTLs links to complex traits. We identified mQTLs for 286,152 CpG sites, many of which (>5%) show tissue specificity, and mQTL colocalizations with 2,254 distinct GWAS hits across 83 traits. For 91% of these loci, a candidate gene link was identified by integration of functional maps, including eQTMs, and/or eQTL colocalization, but only 33% of loci involved an eQTL and mQTL present in the same tissue type. With this DNAm-focused integrative analysis, we contribute to the understanding of molecular regulatory mechanisms in human tissues and their impact on complex traits.

Using epigenomics to understand cellular responses to environmental influences in diseases

It is a generally accepted model that environmental influences can exert their effects, at least in part, by changing the molecular regulators of transcription that are described as epigenetic. As there is biochemical evidence that some epigenetic regulators of transcription can maintain their states long term and through cell division, an epigenetic model encompasses the idea of maintenance of the effect of an exposure long after it is no longer present. The evidence supporting this model is mostly from the observation of alterations of molecular regulators of transcription following exposures. With the understanding that the interpretation of these associations is more complex than originally recognised, this model may be oversimplistic; therefore, adopting novel perspectives and experimental approaches when examining how environmental exposures are linked to phenotypes may prove worthwhile. In this review, we have chosen to use the example of nonalcoholic fatty liver disease (NAFLD), a common, complex human disease with strong environmental and genetic influences. We describe how epigenomic approaches combined with emerging functional genetic and single-cell genomic techniques are poised to generate new insights into the pathogenesis of environmentally influenced human disease phenotypes exemplified by NAFLD.

Multi-omics revolution to promote plant breeding efficiency

Crop production is the primary goal of agricultural activities, which is always taken into consideration. However, global agricultural systems are coming under increasing pressure from the rising food demand of the rapidly growing world population and changing climate. To address these issues, improving high-yield and climate-resilient related-traits in crop breeding is an effective strategy. In recent years, advances in omics techniques, including genomics, transcriptomics, proteomics, and metabolomics, paved the way for accelerating plant/crop breeding to cope with the changing climate and enhance food production. Optimized omics and phenotypic plasticity platform integration, exploited by evolving machine learning algorithms will aid in the development of biological interpretations for complex crop traits. The precise and progressive assembly of desire alleles using precise genome editing approaches and enhanced breeding strategies would enable future crops to excel in combating the changing climates. Furthermore, plant breeding and genetic engineering ensures an exclusive approach to developing nutrient sufficient and climate-resilient crops, the productivity of which can sustainably and adequately meet the world's food, nutrition, and energy needs. This review provides an overview of how the integration of omics approaches could be exploited to select crop varieties with desired traits.

DNA methylation and aeroallergen sensitization: The chicken or the egg?

Background

DNA methylation (DNAm) is considered a plausible pathway through which genetic and environmental factors may influence the development of allergies. However, causality has yet to be determined as it is unknown whether DNAm is rather a cause or consequence of allergic sensitization. Here, we investigated the direction of the observed associations between well-known environmental and genetic determinants of allergy, DNAm, and aeroallergen sensitization using a combination of high-dimensional and causal mediation analyses.

Methods

Using prospectively collected data from the German LISA birth cohort from two time windows (6–10 years: N = 234; 10–15 years: N = 167), we tested whether DNAm is a cause or a consequence of aeroallergen sensitization (specific immunoglobulin E > 0.35kU/l) by conducting mediation analyses for both effect directions using maternal smoking during pregnancy, family history of allergies, and a polygenic risk score (PRS) for any allergic disease as exposure variables. We evaluated individual CpG sites (EPIC BeadChip) and allergy-related methylation risk scores (MRS) as potential mediators in the mediation analyses. We applied three high-dimensional mediation approaches (HIMA, DACT, gHMA) and validated results using causal mediation analyses. A replication of results was attempted in the Swedish BAMSE cohort.

Results

Using high-dimensional methods, we identified five CpGs as mediators of prenatal exposures to sensitization with significant (adjusted p < 0.05) indirect effects in the causal mediation analysis (maternal smoking: two CpGs, family history: one, PRS: two). None of these CpGs could be replicated in BAMSE. The effect of family history on allergy-related MRS was significantly mediated by aeroallergen sensitization (proportions mediated: 33.7–49.6%), suggesting changes in DNAm occurred post-sensitization.

Conclusion

The results indicate that DNAm may be a cause or consequence of aeroallergen sensitization depending on genomic location. Allergy-related MRS, identified as a potential cause of sensitization, can be considered as a cross-sectional biomarker of disease. Differential DNAm in individual CpGs, identified as mediators of the development of sensitization, could be used as clinical predictors of disease development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-022-01332-5.

Incorporation of DNA methylation quantitative trait loci (mQTLs) in epigenome-wide association analysis: application to birthweight effects in neonatal whole blood

BACKGROUND

Epigenome-wide association studies (EWAS) have helped to define the associations between DNA methylation and many clinicopathologic and developmental traits. Since DNA methylation is affected by genetic variation at certain loci, EWAS associations may be potentially influenced by genetic effects. However, a formal assessment of the value of incorporating genetic variation in EWAS evaluations is lacking especially for multiethnic populations.

METHODS

Using single nucleotide polymorphism (SNP) from Illumina Omni Express or Affymetrix PMDA arrays and DNA methylation data from the Illumina 450 K or EPIC array from 1638 newborns of diverse genetic ancestries, we generated DNA methylation quantitative trait loci (mQTL) databases for both array types. We then investigated associations between neonatal DNA methylation and birthweight (incorporating gestational age) using EWAS modeling, and reported how EWAS results were influenced by controlling for mQTLs.

RESULTS

For CpGs on the 450 K array, an average of 15.4% CpGs were assigned as mQTLs, while on the EPIC array, 23.0% CpGs were matched to mQTLs (adjusted P value < 0.05). The CpGs associated with SNPs were enriched in the CpG island shore regions. Correcting for mQTLs in the EWAS model for birthweight helped to increase significance levels for top hits. For CpGs overlapping genes associated with birthweight-related pathways (nutrition metabolism, biosynthesis, for example), accounting for mQTLs changed their regression coefficients more dramatically (> 20%) than for other random CpGs.

CONCLUSION

DNA methylation levels at circa 20% CpGs in the genome were affected by common SNP genotypes. EWAS model fit significantly improved when taking these genetic effects into consideration. Genetic effects were stronger on CpGs overlapping genetic elements associated with control of gene expression.

Identification of novel meQTLs strongly associated with rheumatoid arthritis by large-scale epigenome-wide analysis

Rheumatoid arthritis (RA) is highly heritable, and previous studies have suggested that genetic variation may affect susceptibility to RA by altering epigenetic modifications (e.g. DNA methylation). Here we examined how genetic variation influences DNA methylation (DNAm) in RA by integrating individual genetic variation and DNAm data. Epigenome-wide meQTL (methylation quantitative trait loci) analysis was performed on 354 RA patients and 335 controls, scanning 30,101,744 relationships between 62 SNPs and 485,512 DNA methylation sites. Two regulatory relationship pairs (FDR < 0.05) showed very strong associations with RA risk. One was rs10796216-cg00475509, and the DNAm decreased by 0.0168 per addition of allele rs10796216-A. The other was rs6546473-cg13358873, for which a 0.0365 reduction of DNAm at cg13358873 was observed for each addition of allele rs6546473-A, and lower DNAm was found to be significantly associated with RA risk (P = 2.0407e-28). Moreover, both pairs of meQTL showed a strong regulatory relationship only in RA samples, so they can be subsequently considered as risk markers for RA. In conclusion, our integrated analysis of genetic and epigenetic variation suggests that genetic variation may affect the risk of RA by regulating DNA methylation levels. Alterations of DNAm at cg00475509 and cg13358873 loci conferred by rs10796216-A and rs6546473-A allele may suggest a potential risk for RA. Our results deepen our understanding of the genetic and epigenetic mechanisms of RA and provide novel associations that can be further investigated in future studies.

Elucidating the genetic architecture of DNA methylation to identify promising molecular mechanisms of disease

DNA methylation commonly occurs at cytosine-phosphate-guanine sites (CpGs) that can serve as biomarkers for many diseases. We analyzed whole genome sequencing data to identify DNA methylation quantitative trait loci (mQTLs) in 4126 Framingham Heart Study participants. Our mQTL mapping identified 94,362,817 cis-mQTLvariant-CpG pairs (for 210,156 unique autosomal CpGs) at P < 1e-7 and 33,572,145 trans-mQTL variant-CpG pairs (for 213,606 unique autosomal CpGs) at P < 1e-14. Using cis-mQTL variants for 1258 CpGs associated with seven cardiovascular disease (CVD) risk factors, we found 104 unique CpGs that colocalized with at least one CVD trait. For example, cg11554650 (PPP1R18) colocalized with type 2 diabetes, and was driven by a single nucleotide polymorphism (rs2516396). We performed Mendelian randomization (MR) analysis and demonstrated 58 putatively causal relations of CVD risk factor-associated CpGs to one or more risk factors (e.g., cg05337441 [APOB] with LDL; MR P = 1.2e-99, and 17 causal associations with coronary artery disease (e.g. cg08129017 [SREBF1] with coronary artery disease; MR P = 5e-13). We also showed that three CpGs, e.g., cg14893161 (PM20D1), are putatively causally associated with COVID-19 severity. To assist in future analyses of the role of DNA methylation in disease pathogenesis, we have posted a comprehensive summary data set in the National Heart, Lung, and Blood Institute's BioData Catalyst.

Network reconstruction for trans acting genetic loci using multi-omics data and prior information

BACKGROUND

Molecular measurements of the genome, the transcriptome, and the epigenome, often termed multi-omics data, provide an in-depth view on biological systems and their integration is crucial for gaining insights in complex regulatory processes. These data can be used to explain disease related genetic variants by linking them to intermediate molecular traits (quantitative trait loci, QTL). Molecular networks regulating cellular processes leave footprints in QTL results as so-called trans-QTL hotspots. Reconstructing these networks is a complex endeavor and use of biological prior information can improve network inference. However, previous efforts were limited in the types of priors used or have only been applied to model systems. In this study, we reconstruct the regulatory networks underlying trans-QTL hotspots using human cohort data and data-driven prior information.

METHODS

We devised a new strategy to integrate QTL with human population scale multi-omics data. State-of-the art network inference methods including BDgraph and glasso were applied to these data. Comprehensive prior information to guide network inference was manually curated from large-scale biological databases. The inference approach was extensively benchmarked using simulated data and cross-cohort replication analyses. Best performing methods were subsequently applied to real-world human cohort data.

RESULTS

Our benchmarks showed that prior-based strategies outperform methods without prior information in simulated data and show better replication across datasets. Application of our approach to human cohort data highlighted two novel regulatory networks related to schizophrenia and lean body mass for which we generated novel functional hypotheses.

CONCLUSIONS

We demonstrate that existing biological knowledge can improve the integrative analysis of networks underlying trans associations and generate novel hypotheses about regulatory mechanisms.

The immune factors driving DNA methylation variation in human blood

Epigenetic changes are required for normal development, yet the nature and respective contribution of factors that drive epigenetic variation in humans remain to be fully characterized. Here, we assessed how the blood DNA methylome of 884 adults is affected by DNA sequence variation, age, sex and 139 factors relating to life habits and immunity. Furthermore, we investigated whether these effects are mediated or not by changes in cellular composition, measured by deep immunophenotyping. We show that DNA methylation differs substantially between naïve and memory T cells, supporting the need for adjustment on these cell-types. By doing so, we find that latent cytomegalovirus infection drives DNA methylation variation and provide further support that the increased dispersion of DNA methylation with aging is due to epigenetic drift. Finally, our results indicate that cellular composition and DNA sequence variation are the strongest predictors of DNA methylation, highlighting critical factors for medical epigenomics studies.

A Methylation Diagnostic Model Based on Random Forests and Neural Networks for Asthma Identification

Background

Asthma significantly impacts human life and health as a chronic disease. Traditional treatments for asthma have several limitations. Artificial intelligence aids in cancer treatment and may also accelerate our understanding of asthma mechanisms. We aimed to develop a new clinical diagnosis model for asthma using artificial neural networks (ANN).

Methods

Datasets (GSE85566, GSE40576, and GSE13716) were downloaded from Gene Expression Omnibus (GEO) and identified differentially expressed CpGs (DECs) enriched by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Random forest (RF) and ANN algorithms further identified gene characteristics and built clinical models. In addition, two external validation datasets (GSE40576 and GSE137716) were used to validate the diagnostic ability of the model.

Results

The methylation analysis tool (ChAMP) considered DECs that were up-regulated (n =121) and down-regulated (n =20). GO results showed enrichment of actin cytoskeleton organization and cell-substrate adhesion, shigellosis, and serotonergic synapses. RF (random forest) analysis identified 10 crucial DECs (cg05075579, cg20434422, cg03907390, cg00712106, cg05696969, cg22862094, cg11733958, cg00328720, and cg13570822). ANN constructed the clinical model according to 10 DECs. In two external validation datasets (GSE40576 and GSE137716), the Area Under Curve (AUC) for GSE137716 was 1.000, and AUC for GSE40576 was 0.950, confirming the reliability of the model.

Conclusion

Our findings provide new methylation markers and clinical diagnostic models for asthma diagnosis and treatment.

Digital methylation assessments of alcohol and cigarette consumption account for common variance in accelerated epigenetic ageing

Smoking and Heavy Alcohol Consumption (HAC) are established risk factors for myriad complex disorders of ageing. Yet many prior studies of Epigenetic Ageing (EA) have shown only modest effects of smoking and drinking on accelerated ageing. One potential reason for this conundrum might be the reliance of some prior EA studies on self-reported substance use, which may be unreliable in many samples. To test whether novel, non-self-reported indices would show a stronger association of smoking and HAC to EA, we used methylation sensitive digital PCR (MSdPCR) and data from 437 African American subjects from Wave 7 of the Family and Community Health Study Offspring Cohort to examine the effects of subjective and objective measures of smoking and HAC on 7 indices of EA. Because of limited overall correlations between the various EA indices, we examined patterns of association separately for each index. Consistent with expectations, MSdPCR assessments of smoking and HAC, but not self-reported alcohol consumption, were strongly correlated with accelerated EA. MSdPCR assessments of smoking and HAC accounted for 57% of GrimAge acceleration and the shared variance in GrimAge and DunedinPOAM accelerated EA. We conclude that MSdPCR assessments of smoking and HAC are valuable tools for understanding EA, represent directly targetable conditions for the prevention of premature ageing, and substantially improve upon self-reported assessment of smoking and HAC.