CpG islands in vertebrate genomes

CpG island turnover events predict evolutionary changes in enhancer activity

BACKGROUND

Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. CpG islands (CGIs) have recently been shown to influence enhancer activity, and here we test how their turnover across species contributes to enhancer evolution.

RESULTS

We integrate maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and find that CGI content in enhancers is strongly associated with increased histone modification levels. CGIs show widespread turnover across species and species-specific CGIs are strongly enriched for enhancers exhibiting species-specific activity across all tissues and species. Genes associated with enhancers with species-specific CGIs show concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon.

CONCLUSIONS

Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

Genome-Wide Methylation Analysis in Two Wild-Type Non-Small Cell Lung Cancer Subgroups with Negative and High PD-L1 Expression

We conducted a pilot study to analyze the differential methylation status of 20 primary acinar adenocarcinomas of the lungs. These adenocarcinomas had to be wild type in mutation analysis and had either high (TPS > 50%; n = 10) or negative (TPS < 1%; n = 10) PD-L1 status to be integrated into our study. To examine the methylation of 866,895 specific sites, we utilized the Illumina Infinium EPIC bead chip array. Both hypermethylation and hypomethylation play significant roles in tumor development, progression, and metastasis. They also impact the formation of the tumor microenvironment, which plays a decisive role in tumor differentiation, epigenetics, dissemination, and immune evasion. The gained methylation patterns were correlated with PD-L1 expression. Our analysis has identified distinct methylation patterns in lung adenocarcinomas with high and negative PD-L1 expression. After analyzing the correlation between the methylation results of genes and promoters with their pathobiology, we found that tumors with high expression of PD-L1 tend to exhibit oncogenic effects through hypermethylation. On the other hand, tumors with negative PD-L1 expression show loss of their suppressor functions through hypomethylation. The suppressor functions of hypermethylated genes and promoters are ineffective compared to simultaneously activated dominant oncogenic mechanisms. The tumor microenvironment supports tumor growth in both groups.

Genomic context- and H2AK119 ubiquitination-dependent inheritance of human Polycomb silencing

Polycomb repressive complexes 1 and 2 (PRC1 and 2) are required for heritable repression of developmental genes. The cis- and trans-acting factors that contribute to epigenetic inheritance of mammalian Polycomb repression are not fully understood. Here, we show that, in human cells, ectopically induced Polycomb silencing at initially active developmental genes, but not near ubiquitously expressed housekeeping genes, is inherited for many cell divisions. Unexpectedly, silencing is heritable in cells with mutations in the H3K27me3 binding pocket of the Embryonic Ectoderm Development (EED) subunit of PRC2, which are known to disrupt H3K27me3 recognition and lead to loss of H3K27me3. This mode of inheritance is less stable and requires intact PRC2 and recognition of H2AK119ub1 by PRC1. Our findings suggest that maintenance of Polycomb silencing is sensitive to local genomic context and can be mediated by PRC1-dependent H2AK119ub1 and PRC2 independently of H3K27me3 recognition.

Gene-specific somatic epigenetic mosaicism of FDFT1 underlies a non-hereditary localized form of porokeratosis

Porokeratosis is a clonal keratinization disorder characterized by solitary, linearly arranged, or generally distributed multiple skin lesions. Previous studies showed that genetic alterations in MVK, PMVK, MVD, or FDPS-genes in the mevalonate pathway-cause hereditary porokeratosis, with skin lesions harboring germline and lesion-specific somatic variants on opposite alleles. Here, we identified non-hereditary porokeratosis associated with epigenetic silencing of FDFT1, another gene in the mevalonate pathway. Skin lesions of the generalized form had germline and lesion-specific somatic variants on opposite alleles in FDFT1, representing FDFT1-associated hereditary porokeratosis identified in this study. Conversely, lesions of the solitary or linearly arranged localized form had somatic bi-allelic promoter hypermethylation or mono-allelic promoter hypermethylation with somatic genetic alterations on opposite alleles in FDFT1, indicating non-hereditary porokeratosis. FDFT1 localization was uniformly diminished within the lesions, and lesion-derived keratinocytes showed cholesterol dependence for cell growth and altered expression of genes related to cell-cycle and epidermal development, confirming that lesions form by clonal expansion of FDFT1-deficient keratinocytes. In some individuals with the localized form, gene-specific promoter hypermethylation of FDFT1 was detected in morphologically normal epidermis adjacent to methylation-related lesions but not distal to these lesions, suggesting that asymptomatic somatic epigenetic mosaicism of FDFT1 predisposes certain skin areas to the disease. Finally, consistent with its genetic etiology, topical statin treatment ameliorated lesions in FDFT1-deficient porokeratosis. In conclusion, we identified bi-allelic genetic and/or epigenetic alterations of FDFT1 as a cause of porokeratosis and shed light on the pathogenesis of skin mosaicism involving clonal expansion of epigenetically altered cells.

An analysis of POMC gene methylation and expression in patients with schizophrenia

Schizophrenia is a chronic mental disorder that affects millions of people and is believed to be caused by both environmental and genetic factors. Despite extensive research, the exact mechanisms underlying schizophrenia are still unclear. Studies have shown that numerous psychiatric disorders are associated with methylation of the POMC gene, which encodes adrenocorticotropic hormone, a critical player in the hypothalamic-pituitary-adrenal axis. However, the association between DNA methylation in POMC patients and schizophrenia remains unclear. In this study, we evaluated three fragments of the POMC promoter region, including 51 CpG sites, in the peripheral blood of schizophrenia patients and healthy controls. The POMC protein level was measured via enzyme-linked immunosorbent assay (ELISA). The schizophrenia group exhibited significantly greater levels of methylation of the POMC gene than those in the control group. The methylation level of the POMC-2 fragment was significantly greater in the patient group than in the control group. There were 17 significantly hypermethylated CpG sites in the patient group. After stratification by sex, POMC methylation levels were found to be significantly greater in male schizophrenia patients than in healthy controls; the methylation levels of POMC-2 fragments were greater in the male patient group; nine CpG sites were significantly hypermethylated in the male patient group; and only one CpG site was significantly hypermethylated in the female patient group. The POMC protein level in patients was significantly lower than that in healthy controls. These findings demonstrate that the DNA methylation of POMC might be associated with the pathophysiology of schizophrenia. Overall, studying the correlation between POMC methylation and schizophrenia may contribute to the diagnosis and evaluation of neuropsychiatric disorders.

Novel islands of GGC and GCC repeats coincide with human evolution

BACKGROUND

Because of high mutation rate, overrepresentation in genic regions, and link with various neurological, neurodegenerative, and movement disorders, GGC and GCC short tandem repeats (STRs) are prone to natural selection. Among a number of lacking data, the 3-repeats of these STRs remain widely unexplored.

RESULTS

In a genome-wide search in human, here we mapped GGC and GCC STRs of ≥3-repeats, and found novel islands of up to 45 of those STRs, populating spans of 1 to 2 kb of genomic DNA. RGPD4 and NOC4L harbored the densest (GGC)3 (probability 3.09061E-71) and (GCC)3 (probability 1.72376E-61) islands, respectively, and were human-specific. We also found prime instances of directional incremented density of STRs at specific loci in human versus other species, including the FOXK2 and SKI GGC islands. The genes containing those islands significantly diverged in expression in human versus other species, and the proteins encoded by those genes interact closely in a physical interaction network, consequence of which may be human-specific characteristics such as higher order brain functions.

CONCLUSION

We report novel islands of GGC and GCC STRs of evolutionary relevance to human. The density, and in some instances, periodicity of these islands support them as a novel genomic entity, which need to be further explored in evolutionary, mechanistic, and functional platforms.

Methylation in cornea and corneal diseases: a systematic review

Corneal diseases are among the primary causes of blindness and vision loss worldwide. However, the pathogenesis of corneal diseases remains elusive, and diagnostic and therapeutic tools are limited. Thus, identifying new targets for the diagnosis and treatment of corneal diseases has gained great interest. Methylation, a type of epigenetic modification, modulates various cellular processes at both nucleic acid and protein levels. Growing evidence shows that methylation is a key regulator in the pathogenesis of corneal diseases, including inflammation, fibrosis, and neovascularization, making it an attractive potential therapeutic target. In this review, we discuss the major alterations of methylation and demethylation at the DNA, RNA, and protein levels in corneal diseases and how these dynamics contribute to the pathogenesis of corneal diseases. Also, we provide insights into identifying potential biomarkers of methylation that may improve the diagnosis and treatment of corneal diseases.

Synthetic reversed sequences reveal default genomic states

Pervasive transcriptional activity is observed across diverse species. The genomes of extant organisms have undergone billions of years of evolution, making it unclear whether these genomic activities represent effects of selection or 'noise'1-4. Characterizing default genome states could help understand whether pervasive transcriptional activity has biological meaning. Here we addressed this question by introducing a synthetic 101-kb locus into the genomes of Saccharomyces cerevisiae and Mus musculus and characterizing genomic activity. The locus was designed by reversing but not complementing human HPRT1, including its flanking regions, thus retaining basic features of the natural sequence but ablating evolved coding or regulatory information. We observed widespread activity of both reversed and native HPRT1 loci in yeast, despite the lack of evolved yeast promoters. By contrast, the reversed locus displayed no activity at all in mouse embryonic stem cells, and instead exhibited repressive chromatin signatures. The repressive signature was alleviated in a locus variant lacking CpG dinucleotides; nevertheless, this variant was also transcriptionally inactive. These results show that synthetic genomic sequences that lack coding information are active in yeast, but inactive in mouse embryonic stem cells, consistent with a major difference in 'default genomic states' between these two divergent eukaryotic cell types, with implications for understanding pervasive transcription, horizontal transfer of genetic information and the birth of new genes.

Mechanisms of DNA Methylation Regulatory Function and Crosstalk with Histone Lysine Methylation

DNA methylation is a well-studied epigenetic modification that has key roles in regulating gene expression, maintaining genome integrity, and determining cell fate. Precisely how DNA methylation patterns are established and maintained in specific cell types at key developmental stages is still being elucidated. However, research over the last two decades has contributed to our understanding of DNA methylation regulation by other epigenetic processes. Specifically, lysine methylation on key residues of histone proteins has been shown to contribute to the allosteric regulation of DNA methyltransferase (DNMT) activities. In this review, we discuss the dynamic interplay between DNA methylation and histone lysine methylation as epigenetic regulators of genome function by synthesizing key recent studies in the field. With a focus on DNMT3 enzymes, we discuss mechanisms of DNA methylation and histone lysine methylation crosstalk in the regulation of gene expression and the maintenance of genome integrity. Further, we discuss how alterations to the balance of various sites of histone lysine methylation and DNA methylation contribute to human developmental disorders and cancers. Finally, we provide perspectives on the current direction of the field and highlight areas for continued research and development.

FinaleMe: Predicting DNA methylation by the fragmentation patterns of plasma cell-free DNA

Analysis of DNA methylation in cell-free DNA reveals clinically relevant biomarkers but requires specialized protocols such as whole-genome bisulfite sequencing. Meanwhile, millions of cell-free DNA samples are being profiled by whole-genome sequencing. Here, we develop FinaleMe, a non-homogeneous Hidden Markov Model, to predict DNA methylation of cell-free DNA and, therefore, tissues-of-origin, directly from plasma whole-genome sequencing. We validate the performance with 80 pairs of deep and shallow-coverage whole-genome sequencing and whole-genome bisulfite sequencing data.

Structural Features and Physiological Associations of Human 14-3-3ζ Pseudogenes

There are about 14,000 pseudogenes that are mutated or truncated sequences resembling functional parent genes. About two-thirds of pseudogenes are processed, while others are duplicated. Although initially thought dead, emerging studies indicate they have functional and regulatory roles. We study 14-3-3ζ, an adaptor protein that regulates cytokine signaling and inflammatory diseases, including rheumatoid arthritis, cancer, and neurological disorders. To understand how 14-3-3ζ (gene symbol YWHAZ) performs diverse functions, we examined the human genome and identified nine YWHAZ pseudogenes spread across many chromosomes. Unlike the 32 kb exon-to-exon sequence in YWHAZ, all pseudogenes are much shorter and lack introns. Out of six, four YWHAZ exons are highly conserved, but the untranslated region (UTR) shows significant diversity. The putative amino acid sequence of pseudogenes is 78-97% homologous, resulting in striking structural similarities with the parent protein. The OMIM and Decipher database searches revealed chromosomal loci containing pseudogenes are associated with human diseases that overlap with the parent gene. To the best of our knowledge, this is the first report on pseudogenes of the 14-3-3 family protein and their implications for human health. This bioinformatics-based study introduces a new insight into the complexity of 14-3-3ζ's functions in biology.

Transcriptional regulation of the alternative sex hormone-binding globulin promoter by KLF4

In most mammals the major site of sex hormone-binding globulin (SHBG) synthesis is the liver wherefrom it is secreted into the bloodstream and is the primary determinant of sex steroid access to target tissues. The minor site of SHBG synthesis is the testis and in lower mammals testicular SHBG has long been known to be synthesized and secreted by Sertoli cells. However, human testicular SHBG is expressed in developing germ cells from an upstream alternative promoter (altP-SHBG). Transcripts arising from this region comprise an alternative first exon (1A) with the resultant protein confined to the acrosomal compartment of the mature spermatozoa. I have dissected the regulatory components of the alternative SHBG promoter and identified motifs that are required for optimal transcriptional activity from this region. Transcriptional activity is driven by two CACCC elements that appear to be functionally redundant. The transcription factor KLF4 interacts with promoter the region spanning these elements in vivo. Knockdown of Klf4 results in decreased altP-SHBG activity, while Klf4 overexpression relieves the effects of knockdown. Based on their shared patterns of expression in vivo, I conclude that KLF4 is a transcriptional regulator of SHBG in male germ cells.

Reconstructing DNA methylation maps of ancient populations

Studying premortem DNA methylation from ancient DNA (aDNA) provides a proxy for ancient gene activity patterns, and hence valuable information on evolutionary changes in gene regulation. Due to statistical limitations, current methods to reconstruct aDNA methylation maps are constrained to high-coverage shotgun samples, which comprise a small minority of available ancient samples. Most samples are sequenced using in-situ hybridization capture sequencing which targets a predefined set of genomic positions. Here, we develop methods to reconstruct aDNA methylation maps of samples that were not sequenced using high-coverage shotgun sequencing, by way of pooling together individuals to obtain a DNA methylation map that is characteristic of a population. We show that the resulting DNA methylation maps capture meaningful biological information and allow for the detection of differential methylation across populations. We offer guidelines on how to carry out comparative studies involving ancient populations, and how to control the rate of falsely discovered differentially methylated regions. The ability to reconstruct DNA methylation maps of past populations allows for the development of a whole new frontier in paleoepigenetic research, tracing DNA methylation changes throughout human history, using data from thousands of ancient samples.

Quality assessment of enzymatic methyl-seq library constructed using crude cell lysate

Enzymatic methyl-seq (EM-seq), an enzyme-based method, identifies genome-wide DNA methylation, which enables us to obtain reliable methylome data from purified genomic DNA by avoiding bisulfite-induced DNA damage. However, the loss of DNA during purification hinders the methylome analysis of limited samples. The crude DNA extraction method is the quickest and minimal sample loss approach for obtaining useable DNA without requiring additional dissolution and purification. However, it remains unclear whether crude DNA can be used directly for EM-seq library construction. In this study, we aimed to assess the quality of EM-seq libraries prepared directly using crude DNA. The crude DNA-derived libraries provided appropriate fragment sizes and concentrations for sequencing similar to those of the purified DNA-derived libraries. However, the sequencing results of crude samples exhibited lower reference sequence mapping efficiencies than those of the purified samples. Additionally, the lower-input crude DNA-derived sample exhibited a marginally lower cytosine-to-thymine conversion efficiency and hypermethylated pattern around gene regulatory elements than the higher-input crude DNA- or purified DNA-derived samples. In contrast, the methylation profiles of the crude and purified samples exhibited a significant correlation. Our findings indicate that crude DNA can be used as a raw material for EM-seq library construction.

The practical utility of AI-assisted molecular profiling in the diagnosis and management of cancer of unknown primary: an updated review

Cancer of unknown primary (CUP) presents a complex diagnostic challenge, characterized by metastatic tumors of unknown tissue origin and a dismal prognosis. This review delves into the emerging significance of artificial intelligence (AI) and machine learning (ML) in transforming the landscape of CUP diagnosis, classification, and treatment. ML approaches, trained on extensive molecular profiling data, have shown promise in accurately predicting tissue of origin. Genomic profiling, encompassing driver mutations and copy number variations, plays a pivotal role in CUP diagnosis by providing insights into tumor type-specific oncogenic alterations. Mutational signatures (MS), reflecting somatic mutation patterns, offer further insights into CUP diagnosis. Known MS with established etiology, such as ultraviolet (UV) light-induced DNA damage and tobacco exposure, have been identified in cases of dedifferentiated/transdifferentiated melanoma and carcinoma. Deep learning models that integrate gene expression data and DNA methylation patterns offer insights into tissue lineage and tumor classification. In digital pathology, machine learning algorithms analyze whole-slide images to aid in CUP classification. Finally, precision oncology, guided by molecular profiling, offers targeted therapies independent of primary tissue identification. Clinical trials assigning CUP patients to molecularly guided therapies, including targetable alterations and tumor mutation burden as an immunotherapy biomarker, have resulted in improved overall survival in a subset of patients. In conclusion, AI- and ML-driven approaches are revolutionizing CUP management by enhancing diagnostic accuracy. Precision oncology utilizing enhanced molecular profiling facilitates the identification of targeted therapies that transcend the need to identify the tissue of origin, ultimately improving patient outcomes.

The Versatile Attributes of MGMT: Its Repair Mechanism, Crosstalk with Other DNA Repair Pathways, and Its Role in Cancer

O6-methylguanine-DNA methyltransferase (MGMT or AGT) is a DNA repair protein with the capability to remove alkyl groups from O6-AlkylG adducts. Moreover, MGMT plays a crucial role in repairing DNA damage induced by methylating agents like temozolomide and chloroethylating agents such as carmustine, and thereby contributes to chemotherapeutic resistance when these agents are used. This review delves into the structural roles and repair mechanisms of MGMT, with emphasis on the potential structural and functional roles of the N-terminal domain of MGMT. It also explores the development of cancer therapeutic strategies that target MGMT. Finally, it discusses the intriguing crosstalk between MGMT and other DNA repair pathways.

Modeling methyl-sensitive transcription factor motifs with an expanded epigenetic alphabet

BACKGROUND

Transcription factors bind DNA in specific sequence contexts. In addition to distinguishing one nucleobase from another, some transcription factors can distinguish between unmodified and modified bases. Current models of transcription factor binding tend not to take DNA modifications into account, while the recent few that do often have limitations. This makes a comprehensive and accurate profiling of transcription factor affinities difficult.

RESULTS

Here, we develop methods to identify transcription factor binding sites in modified DNA. Our models expand the standard A/C/G/T DNA alphabet to include cytosine modifications. We develop Cytomod to create modified genomic sequences and we also enhance the MEME Suite, adding the capacity to handle custom alphabets. We adapt the well-established position weight matrix (PWM) model of transcription factor binding affinity to this expanded DNA alphabet. Using these methods, we identify modification-sensitive transcription factor binding motifs. We confirm established binding preferences, such as the preference of ZFP57 and C/EBPβ for methylated motifs and the preference of c-Myc for unmethylated E-box motifs.

CONCLUSIONS

Using known binding preferences to tune model parameters, we discover novel modified motifs for a wide array of transcription factors. Finally, we validate our binding preference predictions for OCT4 using cleavage under targets and release using nuclease (CUT&RUN) experiments across conventional, methylation-, and hydroxymethylation-enriched sequences. Our approach readily extends to other DNA modifications. As more genome-wide single-base resolution modification data becomes available, we expect that our method will yield insights into altered transcription factor binding affinities across many different modifications.

A SINE-VNTR-Alu at the LRIG2 locus is associated with proximal and distal gene expression in CRISPR and population models

SINE-VNTR-Alu (SVA) retrotransposons represent mobile regulatory elements that have the potential to influence the surrounding genome when they insert into a locus. Evolutionarily recent mobilisation has resulted in loci in the human genome where a given retrotransposon might be observed to be present or absent, termed a retrotransposon insertion polymorphism (RIP). We previously observed that an SVA RIP ~ 2 kb upstream of LRIG2 on chromosome 1, the 'LRIG2 SVA', was associated with differences in local gene expression and methylation, and that the two were correlated. Here, we have used CRISPR-mediated deletion of the LRIG2 SVA in a cell line model to validate that presence of the retrotransposon is directly affecting local expression and provide evidence that is suggestive of a modest role for the SVA in modulating nearby methylation. Additionally, in leveraging an available Hi-C dataset we observed that the LRIG2 SVA was also involved in long-range chromatin interactions with a cluster of genes ~ 300 kb away, and that expression of these genes was to varying degrees associated with dosage of the SVA in both CRISPR cell line and population models. Altogether, these data support a regulatory role for SVAs in the modulation of gene expression, with the latter potentially involving chromatin looping, consistent with the model that RIPs may contribute to interpersonal differences in transcriptional networks.

Quantifying negative selection in human 3' UTRs uncovers constrained targets of RNA-binding proteins

Many non-coding variants associated with phenotypes occur in 3' untranslated regions (3' UTRs), and may affect interactions with RNA-binding proteins (RBPs) to regulate gene expression post-transcriptionally. However, identifying functional 3' UTR variants has proven difficult. We use allele frequencies from the Genome Aggregation Database (gnomAD) to identify classes of 3' UTR variants under strong negative selection in humans. We develop intergenic mutability-adjusted proportion singleton (iMAPS), a generalized measure related to MAPS, to quantify negative selection in non-coding regions. This approach, in conjunction with in vitro and in vivo binding data, identifies precise RBP binding sites, miRNA target sites, and polyadenylation signals (PASs) under strong selection. For each class of sites, we identify thousands of gnomAD variants under selection comparable to missense coding variants, and find that sites in core 3' UTR regions upstream of the most-used PAS are under strongest selection. Together, this work improves our understanding of selection on human genes and validates approaches for interpreting genetic variants in human 3' UTRs.

FinaleMe: Predicting DNA methylation by the fragmentation patterns of plasma cell-free DNA

Analysis of DNA methylation in cell-free DNA (cfDNA) reveals clinically relevant biomarkers but requires specialized protocols and sufficient input material that limits its applicability. Millions of cfDNA samples have been profiled by genomic sequencing. To maximize the gene regulation information from the existing dataset, we developed FinaleMe, a non-homogeneous Hidden Markov Model (HMM), to predict DNA methylation of cfDNA and, therefore, tissues-of-origin directly from plasma whole-genome sequencing (WGS). We validated the performance with 80 pairs of deep and shallow-coverage WGS and whole-genome bisulfite sequencing (WGBS) data.

DNA methylation and expression of LGR6 gene in ankylosing spondylitis: A case-control study

OBJECTIVE

The objectives of the present research were to ascertain the relationship of Leucine-Rich Repeat-Containing G-Protein Coupled Receptors 6 (LGR6) methylation and transcript levels with ankylosing spondylitis (AS).

METHODS

Targeted bisulfite sequencing was applied to analyze LGR6 DNA methylation in 81 AS cases and 81 controls. Besides, the LGR6 transcription level of peripheral blood mononuclear cells (PBMCs) from 70 AS cases and 64 controls was measured utilizing quantitative real-time transcription-polymerase chain reaction (qRT-PCR).

RESULTS

The study detected the methylation levels of 43 sites in two CpG (cytosine-guanine dinucleotide) islands of LGR6 and found that LGR6 were significantly hypomethylated in AS patients (LGR6_1: P = 0.002; LGR6_2: P < 0.001). LGR6 transcript level was obviously reduced in AS (P = 0.001) and was positively related to DNA methylation level (CpG-1: P = 0.010; CpG-2: P = 0.007). Besides, the Receiver operating characteristic curve (ROC) exhibited good diagnostic performance of LGR6 methylation level (AUC = 0.676, 95% CI = 0.594-0.758, P < 0.001). Further subgroup analysis revealed that gender may affect the LGR6_1 methylation pattern.

CONCLUSION

The present study revealed that LGR6 DNA methylation dysregulation may be involved in the pathogenesis of AS from an epigenetic perspective for the first time, with the aim of providing new directions for biomarker identification and treatment development for AS patients.

Transcript CD81-215 may be a long noncoding RNA of stromal origin with tumor-promoting role in colon cancer

The role of tetraspanin CD81 in malignant transformation is best studied in colorectal cancer, and it appears that other transcripts beside the fully coding mRNA may also be dysregulated in malignant cells. Recent data from a comprehensive pan-cancer transcriptome analysis demonstrated differential activity of two alternative CD81 gene promoters in malignant versus nonmalignant gut mucosa. The promoter active in gut mucosa gives rise to transcripts CD81-203 and CD81-213, while the promoter active in colon and rectal cancer gives rise to transcripts CD81-205 and CD81-215. Our study aimed to explore the biomarker potential of the transcripts from the alternative CD81 gene promoters in colon cancer, as well as to investigate their structure and potential function using in silico tools. The analysis of the transcripts' expression in several colon cell lines cultivated in 2D and 3D and a set of colon cancer and healthy gut mucosa samples by qPCR and RNA sequencing suggested their low expression and stromal origin. Expression patterns in tumor and nontumor tissue along with in silico data suppose that the transcript CD81-215 may be a noncoding RNA of stromal origin with possible involvement in signaling related to malignant transformation.

Association between nutrient intake related to the one-carbon metabolism and colorectal cancer risk: a case-control study in the Basque Country

PURPOSE

Epidemiologic evidence for the association between methyl-donor nutrient intake and colorectal cancer (CRC) risk remains inconclusive. We aimed to examine the relationship between intake of vitamins of the B group, methionine, total choline and betaine and CRC risk, in a population from the CRC screening programme in the Basque Country.

DESIGN

This observational study included 308 patients with CRC and 308 age- and sex-matched subjects as controls. During recruitment, dietary, anthropometric, lifestyle, socioeconomic, demographic, and health status information was collected. Conditional logistic regression was used to estimate the odds ratios (ORs) for CRC risk.

RESULTS

The adjusted ORs for CRC risk decreased with higher intakes of choline and betaine (p < 0.05). After further adjustment for folate, high intake of choline and betaine remained associated with a reduced CRC risk (adjusted model for choline, OR third tertile vs first tertile = 0.45, 95% CI 0.26-0.80, p = 0.006; for betaine, OR third tertile vs first tertile = 0.27, 95% CI 0.16-0.47, p < 0.001). Regarding the other nutrients, our findings indicated a non-significant decrease in CRC risk with the high level of intake.

CONCLUSIONS

Our data suggest that choline and betaine intake influence CRC risk in the studied population.

DNA Methylation Changes in Blood Cells of Fibromyalgia and Chronic Fatigue Syndrome Patients

Purpose

Fibromyalgia (FM) and Chronic Fatigue Syndrome (CFS) affect 0.4% and 1% of society, respectively, and the prevalence of these pain syndromes is increasing. To date, no strong association between these syndromes and the genetic background of affected individuals has been shown. Therefore, it is plausible that epigenetic changes might play a role in the development of these syndromes.

Patients and Methods

Three previous studies have attempted to elaborate the involvement of genome-wide methylation changes in blood cells in the development of fibromyalgia and chronic fatigue syndrome. These studies included 22 patients with fibromyalgia and 127 patients with CFS, and the results of the studies were largely discrepant. Contradicting results of those studies may be attributed to differences in the omics data analysis approaches used in each study. We reanalyzed the data collected in these studies using an updated and coherent data-analysis framework.

Results

Overall, the methylation changes that we observed overlapped with previous results only to some extent. However, the gene set enrichment analyses based on genes annotated to methylation changes identified in each of the analyzed datasets were surprisingly coherent and uniformly associated with the physiological processes that, when affected, may result in symptoms characteristic of fibromyalgia and chronic fatigue syndrome.

Conclusion

Methylomes of the blood cells of patients with FM and CFS in three independent studies have shown methylation changes that appear to be implicated in the pathogenesis of these syndromes.

Identifying unique exposure-specific transgenerational differentially DNA methylated region epimutations in the genome using hybrid deep learning prediction models

Exposure to environmental toxicants can lead to epimutations in the genome and an increase in differential DNA methylated regions (DMRs) that have been linked to increased susceptibility to various diseases. However, the unique effect of particular toxicants on the genome in terms of leading to unique DMRs for the toxicants has been less studied. One hurdle to such studies is the low number of observed DMRs per toxicants. To address this hurdle, a previously validated hybrid deep-learning cross-exposure prediction model is trained per exposure and used to predict exposure-specific DMRs in the genome. Given these predicted exposure-specific DMRs, a set of unique DMRs per exposure can be identified. Analysis of these unique DMRs through visualization, DNA sequence motif matching, and gene association reveals known and unknown links between individual exposures and their unique effects on the genome. The results indicate the potential ability to define exposure-specific epigenetic markers in the genome and the potential relative impact of different exposures. Therefore, a computational approach to predict exposure-specific transgenerational epimutations was developed, which supported the exposure specificity of ancestral toxicant actions and provided epigenome information on the DMR sites predicted.

A novel chromatin-remodeling complex variant, dcPBAF, is involved in maintaining transcription in differentiated neurons

The Polybromo-associated BAF (BRG1- or BRM-associated factors) (PBAF) chromatin-remodeling complex is essential for transcription in mammalian cells. In this study, we describe a novel variant of the PBAF complex from differentiated neuronal cells, called dcPBAF, that differs from the canonical PBAF existing in proliferating neuroblasts. We describe that in differentiated adult neurons, a specific subunit of PBAF, PHF10, is replaced by a PHF10 isoform that lacks N- and C-terminal domains (called PHF10D). In addition, dcPBAF does not contain the canonical BRD7 subunit. dcPBAF binds promoters of the actively transcribed neuron-specific and housekeeping genes in terminally differentiated neurons of adult mice. Furthermore, in differentiated human neuronal cells, PHF10D-containing dcPBAF maintains a high transcriptional level at several neuron-specific genes.

Cancer signature ensemble integrating cfDNA methylation, copy number, and fragmentation facilitates multi-cancer early detection

Cell-free DNA (cfDNA) sequencing has demonstrated great potential for early cancer detection. However, most large-scale studies have focused only on either targeted methylation sites or whole-genome sequencing, limiting comprehensive analysis that integrates both epigenetic and genetic signatures. In this study, we present a platform that enables simultaneous analysis of whole-genome methylation, copy number, and fragmentomic patterns of cfDNA in a single assay. Using a total of 950 plasma (361 healthy and 589 cancer) and 240 tissue samples, we demonstrate that a multifeature cancer signature ensemble (CSE) classifier integrating all features outperforms single-feature classifiers. At 95.2% specificity, the cancer detection sensitivity with methylation, copy number, and fragmentomic models was 77.2%, 61.4%, and 60.5%, respectively, but sensitivity was significantly increased to 88.9% with the CSE classifier (p value < 0.0001). For tissue of origin, the CSE classifier enhanced the accuracy beyond the methylation classifier, from 74.3% to 76.4%. Overall, this work proves the utility of a signature ensemble integrating epigenetic and genetic information for accurate cancer detection.

The Role of Clonal Hematopoiesis of Indeterminant Potential and DNA (Cytosine-5)-Methyltransferase Dysregulation in Pulmonary Arterial Hypertension and Other Cardiovascular Diseases

DNA methylation is an epigenetic mechanism that regulates gene expression without altering gene sequences in health and disease. DNA methyltransferases (DNMTs) are enzymes responsible for DNA methylation, and their dysregulation is both a pathogenic mechanism of disease and a therapeutic target. DNMTs change gene expression by methylating CpG islands within exonic and intergenic DNA regions, which typically reduces gene transcription. Initially, mutations in the DNMT genes and pathologic DNMT protein expression were found to cause hematologic diseases, like myeloproliferative disease and acute myeloid leukemia, but recently they have been shown to promote cardiovascular diseases, including coronary artery disease and pulmonary hypertension. We reviewed the regulation and functions of DNMTs, with an emphasis on somatic mutations in DNMT3A, a common cause of clonal hematopoiesis of indeterminant potential (CHIP) that may also be involved in the development of pulmonary arterial hypertension (PAH). Accumulation of somatic mutations in DNMT3A and other CHIP genes in hematopoietic cells and cardiovascular tissues creates an inflammatory environment that promotes cardiopulmonary diseases, even in the absence of hematologic disease. This review summarized the current understanding of the roles of DNMTs in maintenance and de novo methylation that contribute to the pathogenesis of cardiovascular diseases, including PAH.

QUFIND: tool for comparative prediction and mining of G4 quadruplexes overlapping with CpG islands

G-quadruplexes (G4s) are secondary structures in DNA that have been shown to be involved in gene regulation. They play a vital role in the cellular processes and several pathogens including bacteria, fungi, and viruses have also been shown to possess G4s that help them in their pathogenesis. Additionally, cross-talk among the CpG islands and G4s has been shown to influence biological processes. The virus-encoded G4s are affected by the mutational landscape leading to the formation/deletion of these G4s. Therefore, understanding and predicting these multivariate effects on traditional and non-traditional quadruplexes forms an important area of research, that is, yet to be investigated. We have designed a user-friendly webserver QUFIND (http://soodlab.com/qufinder/) that can predict traditional as well as non-traditional quadruplexes in a given sequence. QUFIND is connected with ENSEMBL and NCBI so that the sequences can be fetched in a real-time manner. The algorithm is designed in such a way that the user is provided with multiple options to customize the base (A, T, G, or C), size of the stem (2-5), loop length (1-30), number of bulges (1-5) as well as the number of mismatches (0-2) enabling the identification of any of the secondary structure as per their interest. QUFIND is designed to predict both CpG islands as well as G4s in a given sequence. Since G4s are very short as compared to the CpG islands, hence, QUFIND can also predict the overlapping G4s within CpG islands. Therefore, the user has the flexibility to identify either overlapping or non-overlapping G4s along with the CpG islands. Additionally, one section of QUFIND is dedicated to comparing the G4s in two viral sequences. The visualization is designed in such a manner that the user is able to see the unique quadruplexes in both the input sequences. The efficiency of QUFIND is calculated on G4s obtained from G4 high throughput sequencing data (n = 1000) or experimentally validated G4s (n = 329). Our results revealed that QUFIND is able to predict G4-quadruplexes obtained from G4-sequencing data with 90.06% prediction accuracy whereas experimentally validated quadruplexes were predicted with 97.26% prediction accuracy.

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

BACKGROUND

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

CONTENT

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

Liver injury in cynomolgus monkeys following intravenous and intrathecal scAAV9 gene therapy delivery

Hepatotoxicity associated with intravenous/intrathecal adeno-associated virus (AAV) gene therapy has been observed in preclinical species and patients. In nonhuman primates, hepatotoxicity following self-complementary AAV9 administration varies from asymptomatic transaminase elevation with minimal to mild microscopic changes to symptomatic elevations of liver function and thromboinflammatory markers with microscopic changes consistent with marked hepatocellular necrosis and deteriorating clinical condition. These transient acute liver injury marker elevations occur from 3-4 days post intravenous administration to ∼2 weeks post intrathecal administration. No transaminase elevation or microscopic changes were observed with intrathecal administration of empty capsids or a "promoterless genome" vector, suggesting that liver injury after cerebrospinal fluid dosing in nonhuman primates is driven by viral transduction and transgene expression. Co-administration of prednisolone after intravenous or intrathecal dosing did not prevent liver enzyme or microscopic changes despite a reduction of T lymphocyte infiltration in liver tissue. Similarly, co-administration of rituximab/everolimus with intrathecal dosing failed to block AAV-driven hepatotoxicity. Self-complementary AAV-induced acute liver injury appears to correlate with high hepatocellular vector load, macrophage activation, and type 1 interferon innate virus-sensing pathway responses. The current work characterizes key aspects pertaining to early AAV-driven hepatotoxicity in cynomolgus macaques, highlighting the usefulness of this nonclinical species in that context.

Hypomethylation of a CpG Site in the CpG Island of the Aquaporin 1 Gene May Be Involved in the Formation of Adult-Onset Non-communicating Hydrocele Testis

Background Water channel aquaporin 1 (AQP1) protein expression is enhanced in the tunica vaginalis of patients with adult-onset non-communicating hydrocele testis and may contribute to the development of non-communicating hydrocele testis. We performed genetic and epigenetic analyses of the AQP1 gene in the tunica vaginalis of patients with adult-onset non-communicating hydrocele testis to elucidate the cause of enhanced AQP1 protein expression. Methodology The genotype was determined for Tag single-nucleotide polymorphisms (SNPs) representing the AQP1 gene and SNPs in the 5'-upstream region of the AQP1 gene. Then, by performing association analysis, the applicability of various genetic models was investigated for each SNP. Moreover, the methylation rate of CpG sites was examined for the CpG island related to the AQP1 gene. Results There was no significant association between each SNP and hydrocele testis for any of the genetic models. The average methylation rate of the 17 CpG sites evaluated was not significantly different between controls and hydrocele testis, but the methylation rate was lower in hydrocele testis than in controls at one CpG site. Conclusions There was a significant decrease in the methylation rate at one of the CpG sites in the CpG island associated with the AQP1 gene in the tunica vaginalis of patients with non-communicating hydrocele testis. This may increase AQP1 protein expression and contribute to the formation of hydrocele testis. SNPs related to the AQP1 gene were not associated with hydrocele testis.

The mitoepigenome responds to stress, suggesting novel mito-nuclear interactions in vertebrates

The mitochondria are central in the cellular response to changing environmental conditions resulting from disease states, environmental exposures or normal physiological processes. Although the influences of environmental stressors upon the nuclear epigenome are well characterized, the existence and role of the mitochondrial epigenome remains contentious. Here, by quantifying the mitochondrial epigenomic response of pineal gland cells to circadian stress, we confirm the presence of extensive cytosine methylation within the mitochondrial genome. Furthermore, we identify distinct epigenetically plastic regions (mtDMRs) which vary in cytosinic methylation, primarily in a non CpG context, in response to stress and in a sex-specific manner. Motifs enriched in mtDMRs contain recognition sites for nuclear-derived DNA-binding factors (ATF4, HNF4A) important in the cellular metabolic stress response, which we found to be conserved across diverse vertebrate taxa. Together, these findings suggest a new layer of mito-nuclear interaction in which the nuclear metabolic stress response could alter mitochondrial transcriptional dynamics through the binding of nuclear-derived transcription factors in a methylation-dependent context.

DNA methylation markers in the detection of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most common primary liver malignancy and has a poor prognosis. Epigenetic modification has been shown to be deregulated during HCC development by dramatically impacting the differentiation, proliferation, and function of cells. One important epigenetic modification is DNA methylation during which methyl groups are added to cytosines without changing the DNA sequence itself. Studies found that methylated DNA markers can be specific for detection of HCC. On the basis of these findings, the utility of methylated DNA markers as novel biomarkers for early-stage HCC has been measured in blood, and indeed superior sensitivity and specificity have been found in several studies when compared to current surveillance methods. However, a variety of factors currently limit the immediate application of these exciting biomarkers. In this review, we provide a detailed rationalisation of the approach and basis for the use of methylation biomarkers for HCC detection and summarise recent studies on methylated DNA markers in HCC focusing on the importance of the aetiological cause of liver disease in the mechanisms leading to cancer.

Diverse genome-wide DNA methylation alterations in canine hepatocellular tumours

BACKGROUND

Canine hepatocellular tumours (HCTs) are common primary liver tumours. However, the exact mechanisms of tumourigenesis remain unclear. Although some genetic mutations have been reported, DNA methylation alterations in canine HCT have not been well studied.

OBJECTIVES

In this study, we aimed to analyse the DNA methylation status of canine HCT.

METHODS

Tissues from 33 hepatocellular carcinomas, 3 hepatocellular adenomas, 1 nodular hyperplasia, 21 non-tumour livers from the patients and normal livers from 5 healthy dogs were used. We analysed the DNA methylation levels of 72,367 cytosine-guanine dinucleotides (CpG sites) in all 63 samples.

RESULTS AND CONCLUSIONS

Although a large fraction of CpG sites that were highly methylated in the normal liver became hypomethylated in tumours from most patients, we also found some patients with less remarkable change or no change in DNA methylation. Hierarchical clustering analysis revealed that 32 of 37 tumour samples differed from normal livers, although the remaining 5 tumour livers fell into the same cluster as normal livers. In addition, the number of hypermethylated genes in tumour livers varied among tumour cases, suggesting various DNA methylation patterns in different tumour groups. However, patient and clinical parameters, such as age, were not associated with DNA methylation status. In conclusion, we found that HCTs undergo aberrant and diverse patterns of genome-wide DNA methylation compared with normal liver tissue, suggesting a complex epigenetic mechanism in canine HCT.

De Novo Mutation Rates in Sticklebacks

Mutation rate is a fundamental parameter in population genetics. Apart from being an important scaling parameter for demographic and phylogenetic inference, it allows one to understand at what rate new genetic diversity is generated and what the expected level of genetic diversity is in a population at equilibrium. However, except for well-established model organisms, accurate estimates of de novo mutation rates are available for a very limited number of organisms from the wild. We estimated mutation rates (µ) in two marine populations of the nine-spined stickleback (Pungitius pungitius) with the aid of several 2- and 3-generational family pedigrees, deep (>50×) whole-genome resequences and a high-quality reference genome. After stringent filtering, we discovered 308 germline mutations in 106 offspring translating to µ = 4.83 × 10-9 and µ = 4.29 × 10-9 per base per generation in the two populations, respectively. Up to 20% of the mutations were shared by full-sibs showing that the level of parental mosaicism was relatively high. Since the estimated µ was 3.1 times smaller than the commonly used substitution rate, recalibration with µ led to substantial increase in estimated divergence times between different stickleback species. Our estimates of the de novo mutation rate should provide a useful resource for research focused on fish population genetics and that of sticklebacks in particular.

A Review of the Epigenetic Clock: Emerging Biomarkers for Asthma and Allergic Disease

DNA methylation (DNAm) is a dynamic, age-dependent epigenetic modification that can be used to study interactions between genetic and environmental factors. Environmental exposures during critical periods of growth and development may alter DNAm patterns, leading to increased susceptibility to diseases such as asthma and allergies. One method to study the role of DNAm is the epigenetic clock-an algorithm that uses DNAm levels at select age-informative Cytosine-phosphate-Guanine (CpG) dinucleotides to predict epigenetic age (EA). The difference between EA and calendar age (CA) is termed epigenetic age acceleration (EAA) and reveals information about the biological capacity of an individual. Associations between EAA and disease susceptibility have been demonstrated for a variety of age-related conditions and, more recently, phenotypes such as asthma and allergic diseases, which often begin in childhood and progress throughout the lifespan. In this review, we explore different epigenetic clocks and how they have been applied, particularly as related to childhood asthma. We delve into how in utero and early life exposures (e.g., smoking, air pollution, maternal BMI) result in methylation changes. Furthermore, we explore the potential for EAA to be used as a biomarker for asthma and allergic diseases and identify areas for further study.

Hypermethylation of RNF125 promotes autophagy-induced oxidative stress in asthma by increasing HMGB1 stability

Asthma is a global chronic airway disease. The expression and role of RNF125, an E3 ubiquitin ligase, in asthma remain uncertain. In this study, we revealed that RNF125 was downregulated in the bronchial epithelium of mice and patients with asthma. Rnf125 hypermethylation was responsible for the low expression of RNF125 in primary airway epithelial cells of mice treated with OVA. Moreover, we demonstrated that RNF125 could attenuate autophagy, oxidative stress, and protect epithelial barrier in vivo and in vitro. Additionally, we identified HMGB1 as a substrate of RNF125, which interacted with the HMG B-box domain of HMGB1 and induced degradation via the ubiquitin proteasome system, reducing autophagy and oxidative stress. Overall, our findings elucidated that hypermethylation of Rnf125 reduced its expression, which promoted autophagy-induced oxidative stress in asthma by increasing HMGB1 stability. These findings offer a theoretical and experimental basis for the pathogenesis of asthma.

Glucocorticoid Receptor Regulates and Interacts with LEDGF/p75 to Promote Docetaxel Resistance in Prostate Cancer Cells

Patients with advanced prostate cancer (PCa) invariably develop resistance to anti-androgen therapy and taxane-based chemotherapy. Glucocorticoid receptor (GR) has been implicated in PCa therapy resistance; however, the mechanisms underlying GR-mediated chemoresistance remain unclear. Lens epithelium-derived growth factor p75 (LEDGF/p75, also known as PSIP1 and DFS70) is a glucocorticoid-induced transcription co-activator implicated in cancer chemoresistance. We investigated the contribution of the GR-LEDGF/p75 axis to docetaxel (DTX)-resistance in PCa cells. GR silencing in DTX-sensitive and -resistant PCa cells decreased LEDGF/p75 expression, and GR upregulation in enzalutamide-resistant cells correlated with increased LEDGF/p75 expression. ChIP-sequencing revealed GR binding sites in the LEDGF/p75 promoter. STRING protein-protein interaction analysis indicated that GR and LEDGF/p75 belong to the same transcriptional network, and immunochemical studies demonstrated their co-immunoprecipitation and co-localization in DTX-resistant cells. The GR modulators exicorilant and relacorilant increased the sensitivity of chemoresistant PCa cells to DTX-induced cell death, and this effect was more pronounced upon LEDGF/p75 silencing. RNA-sequencing of DTX-resistant cells with GR or LEDGF/p75 knockdown revealed a transcriptomic overlap targeting signaling pathways associated with cell survival and proliferation, cancer, and therapy resistance. These studies implicate the GR-LEDGF/p75 axis in PCa therapy resistance and provide a pre-clinical rationale for developing novel therapeutic strategies for advanced PCa.

Toward Identification of Functional Sequences and Variants in Noncoding DNA

Understanding the noncoding part of the genome, which encodes gene regulation, is necessary to identify genetic mechanisms of disease and translate findings from genome-wide association studies into actionable results for treatments and personalized care. Here we provide an overview of the computational analysis of noncoding regions, starting from gene-regulatory mechanisms and their representation in data. Deep learning methods, when applied to these data, highlight important regulatory sequence elements and predict the functional effects of genetic variants. These and other algorithms are used to predict damaging sequence variants. Finally, we introduce rare-variant association tests that incorporate functional annotations and predictions in order to increase interpretability and statistical power.

Methylation is maintained specifically at imprinting control regions but not other DMRs associated with imprinted genes in mice bearing a mutation in the Dnmt1 intrinsically disordered domain

Differential methylation of imprinting control regions in mammals is essential for distinguishing the parental alleles from each other and regulating their expression accordingly. To ensure parent of origin-specific expression of imprinted genes and thereby normal developmental progression, the differentially methylated states that are inherited at fertilization must be stably maintained by DNA methyltransferase 1 throughout subsequent somatic cell division. Further epigenetic modifications, such as the acquisition of secondary regions of differential methylation, are dependent on the methylation status of imprinting control regions and are important for achieving the monoallelic expression of imprinted genes, but little is known about how imprinting control regions direct the acquisition and maintenance of methylation at these secondary sites. Recent analysis has identified mutations that reduce DNA methyltransferase 1 fidelity at some genomic sequences but not at others, suggesting that it may function differently at different loci. We examined the impact of the mutant DNA methyltransferase 1 P allele on methylation at imprinting control regions as well as at secondary differentially methylated regions and non-imprinted sequences. We found that while the P allele results in a major reduction in DNA methylation levels across the mouse genome, methylation is specifically maintained at imprinting control regions but not at their corresponding secondary DMRs. This result suggests that DNA methyltransferase 1 may work differently at imprinting control regions or that there is an alternate mechanism for maintaining methylation at these critical regulatory regions and that maintenance of methylation at secondary DMRs is not solely dependent on the methylation status of the ICR.

Toward a comprehensive catalog of regulatory elements

Regulatory elements are the genomic regions that interact with transcription factors to control cell-type-specific gene expression in different cellular environments. A precise and complete catalog of functional elements encoded by the human genome is key to understanding mammalian gene regulation. Here, we review the current state of regulatory element annotation. We first provide an overview of assays for characterizing functional elements, including genome, epigenome, transcriptome, three-dimensional chromatin interaction, and functional validation assays. We then discuss computational methods for defining regulatory elements, including peak-calling and other statistical modeling methods. Finally, we introduce several high-quality lists of regulatory element annotations and suggest potential future directions.

OGRE: calculate, visualize, and analyze overlap between genomic input regions and public annotations

BACKGROUND

Modern genome sequencing leads to an ever-growing collection of genomic annotations. Combining these elements with a set of input regions (e.g. genes) would yield new insights in genomic associations, such as those involved in gene regulation. The required data are scattered across different databases making a manual approach tiresome, unpractical, and prone to error. Semi-automatic approaches require programming skills in data parsing, processing, overlap calculation, and visualization, which most biomedical researchers lack. Our aim was to develop an automated tool providing all necessary algorithms, benefiting both bioinformaticians and researchers without bioinformatic training.

RESULTS

We developed overlapping annotated genomic regions (OGRE) as a comprehensive tool to associate and visualize input regions with genomic annotations. It does so by parsing regions of interest, mining publicly available annotations, and calculating possible overlaps between them. The user can thus identify location, type, and number of associated regulatory elements. Results are presented as easy to understand visualizations and result tables. We applied OGRE to recent studies and could show high reproducibility and potential new insights. To demonstrate OGRE's performance in terms of running time and output, we have conducted a benchmark and compared its features with similar tools.

CONCLUSIONS

OGRE's functions and built-in annotations can be applied as a downstream overlap association step, which is compatible with most genomic sequencing outputs, and can thus enrich pre-existing analyses pipelines. Compared to similar tools, OGRE shows competitive performance, offers additional features, and has been successfully applied to two recent studies. Overall, OGRE addresses the lack of tools for automatic analysis, local genomic overlap calculation, and visualization by providing an easy to use, end-to-end solution for both biologists and computational scientists.

CpG Methylation Levels in HPA Axis Genes Predict Chronic Pain Outcomes Following Trauma Exposure

Chronic post-traumatic musculoskeletal pain (CPTP) is a common outcome of traumatic stress exposure. Biological factors that influence the development of CPTP are poorly understood, though current evidence indicates that the hypothalamic-pituitary-adrenal (HPA) axis plays a critical role in its development. Little is known about molecular mechanisms underlying this association, including epigenetic mechanisms. Here, we assessed whether peritraumatic DNA methylation levels at 248 5'-C-phosphate-G-3' (CpG) sites in HPA axis genes (FKBP5, NR3C1, CRH, CRHR1, CRHR2, CRHBP, POMC) predict CPTP and whether identified CPTP-associated methylation levels influence expression of those genes. Using participant samples and data collected from trauma survivors enrolled into longitudinal cohort studies (n = 290), we used linear mixed modeling to assess the relationship between peritraumatic blood-based CpG methylation levels and CPTP. A total of 66 (27%) of the 248 CpG sites assessed in these models statistically significantly predicted CPTP, with the three most significantly associated CpG sites originating from the POMC gene region (ie, cg22900229 [β = .124, P < .001], cg16302441 [β = .443, P < .001], cg01926269 [β = .130, P < .001]). Among the genes analyzed, both POMC (z = 2.36, P = .018) and CRHBP (z = 4.89, P < .001) were enriched in CpG sites significantly associated with CPTP. Further, POMC expression was inversely correlated with methylation levels in a CPTP-dependent manner (6-months NRS<4: r = -.59, P < .001; 6-months NRS ≥ 4: r = -.18, P = .2312). Our results suggest that methylation of HPA axis genes including POMC and CRHBP predict risk for and may contribute to vulnerability to CPTP. PERSPECTIVE: Peritraumatic blood levels of CpG methylation sites in HPA axis genes, particularly CpG sites in the POMC gene, predict CPTP development. This data substantially advances our understanding of epigenetic predictors and potential mediators of CPTP, a highly common, morbid, and hard-to-treat form of chronic pain.

Refining the genomic determinants underlying escape from X-chromosome inactivation

X-chromosome inactivation (XCI) epigenetically silences one X chromosome in every cell in female mammals. Although the majority of X-linked genes are silenced, in humans 20% or more are able to escape inactivation and continue to be expressed. Such escape genes are important contributors to sex differences in gene expression, and may impact the phenotypes of X aneuploidies; yet the mechanisms regulating escape from XCI are not understood. We have performed an enrichment analysis of transcription factor binding on the X chromosome, providing new evidence for enriched factors at the transcription start sites of escape genes. The top escape-enriched transcription factors were detected at the RPS4X promoter, a well-described human escape gene previously demonstrated to escape from XCI in a transgenic mouse model. Using a cell line model system that allows for targeted integration and inactivation of transgenes on the mouse X chromosome, we further assessed combinations of RPS4X promoter and genic elements for their ability to drive escape from XCI. We identified a small transgenic construct of only 6 kb capable of robust escape from XCI, establishing that gene-proximal elements are sufficient to permit escape, and highlighting the additive effect of multiple elements that work together in a context-specific fashion.

Mechanisms of human DNA methylation, alteration of methylation patterns in physiological processes and oncology

DNA methylation is one of the epigenetic modifications of the genome, the essence of which is the attachment of a methyl group to nitrogenous bases. In the eukaryote genome, cytosine is methylated in the vast majority of cases. About 98% of cytosines are methylated as part of CpG dinucleotides. They, in turn, form CpG islands, which are clusters of these dinucleotides. Islands located in the regulatory elements of genes are in particular interest. They are assumed to play an important role in the regulation of gene expression in humans. Besides that, cytosine methylation serves the functions of genomic imprinting, transposon suppression, epigenetic memory maintenance, X- chromosome inactivation, and embryonic development. Of particular interest are the enzymatic processes of methylation and demethylation. The methylation process always depends on the work of enzymatic complexes and is very precisely regulated. The methylation process largely depends on the functioning of three groups of enzymes: writers, readers and erasers. Writers include proteins of the DNMT family, readers are proteins containing the MBD, BTB/POZ or SET- and RING-associated domains and erasers are proteins of the TET family. Whereas demethylation can be performed not only by enzymatic complexes, but also passively during DNA replication. Hence, the maintenance of DNA methylation is important. Changes in methylation patterns are observed during embryonic development, aging, and cancers. In both aging and cancer, massive hypomethylation of the genome with local hypermethylation is observed. In this review, we will review the current understanding of the mechanisms of DNA methylation and demethylation in humans, the structure and distribution of CpG islands, the role of methylation in the regulation of gene expression, embryogenesis, aging, and cancer development.

Genomic and Epigenetic Changes Drive Aberrant Skeletal Muscle Differentiation in Rhabdomyosarcoma

Rhabdomyosarcoma (RMS), the most common soft-tissue sarcoma in children and adolescents, represents an aberrant form of skeletal muscle differentiation. Both skeletal muscle development, as well as regeneration of adult skeletal muscle are governed by members of the myogenic family of regulatory transcription factors (MRFs), which are deployed in a highly controlled, multi-step, bidirectional process. Many aspects of this complex process are deregulated in RMS and contribute to tumorigenesis. Interconnected loops of super-enhancers, called core regulatory circuitries (CRCs), define aberrant muscle differentiation in RMS cells. The transcriptional regulation of MRF expression/activity takes a central role in the CRCs active in skeletal muscle and RMS. In PAX3::FOXO1 fusion-positive (PF+) RMS, CRCs maintain expression of the disease-driving fusion oncogene. Recent single-cell studies have revealed hierarchically organized subsets of cells within the RMS cell pool, which recapitulate developmental myogenesis and appear to drive malignancy. There is a large interest in exploiting the causes of aberrant muscle development in RMS to allow for terminal differentiation as a therapeutic strategy, for example, by interrupting MEK/ERK signaling or by interfering with the epigenetic machinery controlling CRCs. In this review, we provide an overview of the genetic and epigenetic framework of abnormal muscle differentiation in RMS, as it provides insights into fundamental mechanisms of RMS malignancy, its remarkable phenotypic diversity and, ultimately, opportunities for therapeutic intervention.

CpG island turnover events predict evolutionary changes in enhancer activity

Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. Here we show that turnover of CpG islands (CGIs), which contribute to enhancer activation, is broadly associated with changes in enhancer activity across mammals, including humans. We integrated maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and found that CGI content in enhancers was strongly associated with increased histone modification levels. CGIs showed widespread turnover across species and species-specific CGIs were strongly enriched for enhancers exhibiting species-specific activity across all tissues and species we examined. Genes associated with enhancers with species-specific CGIs showed concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon. Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

Chicken Erythrocyte: Epigenomic Regulation of Gene Activity

The chicken genome is one-third the size of the human genome and has a similarity of sixty percent when it comes to gene content. Harboring similar genome sequences, chickens' gene arrangement is closer to the human genomic organization than it is to rodents. Chickens have been used as model organisms to study evolution, epigenome, and diseases. The chicken nucleated erythrocyte's physiological function is to carry oxygen to the tissues and remove carbon dioxide. The erythrocyte also supports the innate immune response in protecting the chicken from pathogens. Among the highly studied aspects in the field of epigenetics are modifications of DNA, histones, and their variants. In understanding the organization of transcriptionally active chromatin, studies on the chicken nucleated erythrocyte have been important. Through the application of a variety of epigenomic approaches, we and others have determined the chromatin structure of expressed/poised genes involved in the physiological functions of the erythrocyte. As the chicken erythrocyte has a nucleus and is readily isolated from the animal, the chicken erythrocyte epigenome has been studied as a biomarker of an animal's long-term exposure to stress. In this review, epigenomic features that allow erythroid gene expression in a highly repressive chromatin background are presented.

The Epigenetics of Animal Personality

Animal personality, consistent individual differences in behaviour, is an important concept for understanding how individuals vary in how they cope with environmental challenges. In order to understand the evolutionary significance of animal personality, it is crucial to understand the underlying regulatory mechanisms. Epigenetic marks such as DNA methylation are hypothesised to play a major role in explaining variation in phenotypic changes in response to environmental alterations. Several characteristics of DNA methylation also align well with the concept of animal personality. In this review paper, we summarise the current literature on the role that molecular epigenetic mechanisms may have in explaining personality variation. We elaborate on the potential for epigenetic mechanisms to explain behavioural variation, behavioural development and temporal consistency in behaviour. We then suggest future routes for this emerging field and point to potential pitfalls that may be encountered. We conclude that a more inclusive approach is needed for studying the epigenetics of animal personality and that epigenetic mechanisms cannot be studied without considering the genetic background.