DNA methylation profiles of human active and inactive X chromosomes

Sex Differences of Neutrophil Extracellular Traps on Lipopolysaccharide-Stimulated Human Neutrophils

Objective: Sex differences exist in sepsis, but the commitment of neutrophils to these differences remains unclear. Neutrophil extracellular traps (NETs) function to remove pathogens, yet excessive NETs release can contribute to organ damage. This study explores effects of the gender hormones on endotoxin-induced NETs using neutrophils from both male and female sources. Methods: Blood samples were collected from healthy volunteers. Isolated neutrophils were seeded in collagen-coated cell culture plates, and NETs were induced by lipopolysaccharide (LPS) treatment. After 15 minutes of LPS treatment, 17β-estradiol (0.03-272.4 ng/mL), testosterone enanthate (0.01-10 ng/mL), dimethyl sulfoxide, or ethanol (vehicle control) was added to the plates. These were incubated for three hours at 37°C with 5% CO2. Neutrophil extracellular traps formation was assessed using immunofluorescence staining. Results: Lipopolysaccharide-induced NETs formation was significantly greater in females than in males. In male-derived neutrophils, 17β-estradiol at above the blood concentrations significantly suppressed LPS-induced NETs. No effect was seen while using testosterone enanthate to NETs at any concentration. In female-derived neutrophils, 17β-estradiol, which was near to the highest concentration of non-pregnant women's blood, tended to increase NETs. Testosterone enanthate, which was near to female blood concentration, significantly promoted NETs. Conclusions: Sex differences existed in LPS-induced NETs of human neutrophil. In males, high concentrations of 17β-estradiol administration may have a suppressive effect on excessive NETs during infection. In females, endogenous gender hormones may promote NETs during infection. Sex differences in neutrophils may need to be considered in organ damage owing to NETs excess such as sepsis.

The conneXion between sex and immune responses

There are notable sex-based differences in immune responses to pathogens and self-antigens, with female individuals exhibiting increased susceptibility to various autoimmune diseases, and male individuals displaying preferential susceptibility to some viral, bacterial, parasitic and fungal infections. Although sex hormones clearly contribute to sex differences in immune cell composition and function, the presence of two X chromosomes in female individuals suggests that differential gene expression of numerous X chromosome-linked immune-related genes may also influence sex-biased innate and adaptive immune cell function in health and disease. Here, we review the sex differences in immune system composition and function, examining how hormones and genetics influence the immune system. We focus on the genetic and epigenetic contributions responsible for altered X chromosome-linked gene expression, and how this impacts sex-biased immune responses in the context of pathogen infection and systemic autoimmunity.

Extraembryonic gut endoderm cells undergo programmed cell death during development

Despite a distinct developmental origin, extraembryonic cells in mice contribute to gut endoderm and converge to transcriptionally resemble their embryonic counterparts. Notably, all extraembryonic progenitors share a non-canonical epigenome, raising several pertinent questions, including whether this landscape is reset to match the embryonic regulation and if extraembryonic cells persist into later development. Here we developed a two-colour lineage-tracing strategy to track and isolate extraembryonic cells over time. We find that extraembryonic gut cells display substantial memory of their developmental origin including retention of the original DNA methylation landscape and resulting transcriptional signatures. Furthermore, we show that extraembryonic gut cells undergo programmed cell death and neighbouring embryonic cells clear their remnants via non-professional phagocytosis. By midgestation, we no longer detect extraembryonic cells in the wild-type gut, whereas they persist and differentiate further in p53-mutant embryos. Our study provides key insights into the molecular and developmental fate of extraembryonic cells inside the embryo.

Therapeutic targeting of DNA methylation alterations in cancer

DNA methylation is a critical component of gene regulation and plays an important role in the development of cancer. Hypermethylation of tumor suppressor genes and silencing of DNA repair pathways facilitate uncontrolled cell growth and synergize with oncogenic mutations to perpetuate cancer phenotypes. Additionally, aberrant DNA methylation hinders immune responses crucial for antitumor immunity. Thus, inhibiting dysregulated DNA methylation is a promising cancer therapy. Pharmacologic inhibition of DNA methylation reactivates silenced tumor suppressors and bolster immune responses through induction of viral mimicry. Now, with the advent of immunotherapies and discovery of the immune-modulatory effects of DNA methylation inhibitors, there is great interest in understanding how targeting DNA methylation in combination with other therapies can enhance antitumor immunity. Here, we describe the role of aberrant DNA methylation in cancer and mechanisms by which it promotes tumorigenesis and modulates immune responses. Finally, we review the initial discoveries and ongoing efforts to target DNA methylation as a cancer therapeutic.

Sex-related immunity: could Toll-like receptors be the answer in acute inflammatory response?

An increasing number of studies have highlighted the existence of a sex-specific immune response, wherein men experience a worse prognosis in cases of acute inflammatory diseases. Initially, this sex-dependent inflammatory response was attributed to the influence of sex hormones. However, a growing body of evidence has shifted the focus toward the influence of chromosomes rather than sex hormones in shaping these inflammatory sex disparities. Notably, certain pattern recognition receptors, such as Toll-like receptors (TLRs), and their associated immune pathways have been implicated in driving the sex-specific immune response. These receptors are encoded by genes located on the X chromosome. TLRs are pivotal components of the innate immune system, playing crucial roles in responding to infectious diseases, including bacterial and viral pathogens, as well as trauma-related conditions. Importantly, the TLR-mediated inflammatory responses, as indicated by the production of specific proteins and cytokines, exhibit discernible sex-dependent patterns. In this review, we delve into the subject of sex bias in TLR activation and explore its clinical implications relatively to both the X chromosome and the hormonal environment. The overarching objective is to enhance our understanding of the fundamental mechanisms underlying these sex differences.

Mutation rate heterogeneity at the sub-gene scale due to local DNA hypomethylation

Local mutation rates in human are highly heterogeneous, with known variability at the scale of megabase-sized chromosomal domains, and, on the other extreme, at the scale of oligonucleotides. The intermediate, kilobase-scale heterogeneity in mutation risk is less well characterized. Here, by analyzing thousands of somatic genomes, we studied mutation risk gradients along gene bodies, representing a genomic scale spanning roughly 1-10 kb, hypothesizing that different mutational mechanisms are differently distributed across gene segments. The main heterogeneity concerns several kilobases at the transcription start site and further downstream into 5' ends of gene bodies; these are commonly hypomutated with several mutational signatures, most prominently the ubiquitous C > T changes at CpG dinucleotides. The width and shape of this mutational coldspot at 5' gene ends is variable across genes, and corresponds to variable interval of lowered DNA methylation depending on gene activity level and regulation. Such hypomutated loci, at 5' gene ends or elsewhere, correspond to DNA hypomethylation that can associate with various landmarks, including intragenic enhancers, Polycomb-marked regions, or chromatin loop anchor points. Tissue-specific DNA hypomethylation begets tissue-specific local hypomutation. Of note, direction of mutation risk is inverted for AID/APOBEC3 cytosine deaminase activity, whose signatures are enriched in hypomethylated regions.

DNA methylation differences between the female and male X chromosomes in human brain

The mechanisms of X chromosome inactivation suggest fundamental epigenetic differences between the female and male X chromosomes. However, DNA methylation studies often exclude the X chromosomes. In addition, many previous studies relied on techniques that examine non-randomly selected subsets of positions such as array-based methods, rather than assessing the whole X chromosome. Consequently, our understanding of X chromosome DNA methylation lags behind that of autosomes. Here we addressed this gap of knowledge by studying X chromosome DNA methylation using 89 whole genome bisulfite sequencing (WGBS) maps from neurons and oligodendrocytes. Using this unbiased and comprehensive data, we show that DNA methylation of the female X chromosomes is globally reduced (hypomethylated) across the entire chromosome compared to the male X chromosomes and autosomes. On the other hand, the majority of X-linked promoters were more highly methylated (hypermethylated) in females compared to males, consistent with the role of DNA methylation in X chromosome inactivation and dosage compensation. Remarkably, hypermethylation of female X promoters was limited to a group of previously lowly methylated promoters. The other group of highly methylated promoters were both hyper- and hypo-methylated in females with no obvious association with gene expression. Therefore, X chromosome inactivation by DNA methylation was exclusive to a subset of promoters with distinctive epigenetic feature. Apart from this group of promoters, differentially methylated regions in the female and male X chromosomes were dominated by female hypomethylation. Our study furthers the understanding of X-chromosome dosage regulation by DNA methylation on the chromosomal level as well as on individual gene level.

Mechanisms of sex differences in Alzheimer's disease

Alzheimer's disease (AD) and the mechanisms underlying its etiology and progression are complex and multifactorial. The higher AD risk in women may serve as a clue to better understand these complicated processes. In this review, we examine aspects of AD that demonstrate sex-dependent effects and delve into the potential biological mechanisms responsible, compiling findings from advanced technologies such as single-cell RNA sequencing, metabolomics, and multi-omics analyses. We review evidence that sex hormones and sex chromosomes interact with various disease mechanisms during aging, encompassing inflammation, metabolism, and autophagy, leading to unique characteristics in disease progression between men and women.

Assessing methylation detection for primary human tissue using Nanopore sequencing

DNA methylation most commonly occurs as 5-methylcytosine (5-mC) in the human genome and has been associated with human diseases. Recent developments in single-molecule sequencing technologies (Oxford Nanopore Technologies (ONT) and Pacific Biosciences) have enabled readouts of long, native DNA molecules, including cytosine methylation. ONT recently upgraded their Nanopore sequencing chemistry and kits from R9 to the R10 version, which yielded increased accuracy and sequencing throughput. However the effects on methylation detection have not yet been documented. Here we performed a series of computational analyses to characterize differences in Nanopore-based 5mC detection between the ONT R9 and R10 chemistries. We compared 5mC calls in R9 and R10 for three human genome datasets: a cell line, a frontal cortex brain sample, and a blood sample. We performed an in-depth analysis on CpG islands and homopolymer regions, and documented high concordance for methylation detection among sequencing technologies. The strongest correlation was observed between Nanopore R10 and Illumina bisulfite technologies for cell line-derived datasets. Subtle differences in methylation datasets between technologies can impact analysis tools such as differential methylation calling software. Our findings show that comparisons can be drawn between methylation data from different Nanopore chemistries using guided hypotheses. This work will facilitate comparison among Nanopore data cohorts derived using different chemistries from large scale sequencing efforts, such as the NIH CARD Long Read Initiative.

Role of sex in immune response and epigenetic mechanisms

The functioning of the human immune system is highly dependent on the sex of the individual, which comes by virtue of sex chromosomes and hormonal differences. Epigenetic mechanisms such as X chromosome inactivation, mosaicism, skewing, and dimorphism in X chromosome genes and Y chromosome regulatory genes create a sex-based variance in the immune response between males and females. This leads to differential susceptibility in immune-related disorders like infections, autoimmunity, and malignancies. Various naturally available immunomodulators are also available which target immune pathways containing X chromosome genes.

The Complete Sequence and Comparative Analysis of Ape Sex Chromosomes

Apes possess two sex chromosomes-the male-specific Y and the X shared by males and females. The Y chromosome is crucial for male reproduction, with deletions linked to infertility. The X chromosome carries genes vital for reproduction and cognition. Variation in mating patterns and brain function among great apes suggests corresponding differences in their sex chromosome structure and evolution. However, due to their highly repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the state-of-the-art experimental and computational methods developed for the telomere-to-telomere (T2T) human genome, we produced gapless, complete assemblies of the X and Y chromosomes for five great apes (chimpanzee, bonobo, gorilla, Bornean and Sumatran orangutans) and a lesser ape, the siamang gibbon. These assemblies completely resolved ampliconic, palindromic, and satellite sequences, including the entire centromeres, allowing us to untangle the intricacies of ape sex chromosome evolution. We found that, compared to the X, ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements. This divergence on the Y arises from the accumulation of lineage-specific ampliconic regions and palindromes (which are shared more broadly among species on the X) and from the abundance of transposable elements and satellites (which have a lower representation on the X). Our analysis of Y chromosome genes revealed lineage-specific expansions of multi-copy gene families and signatures of purifying selection. In summary, the Y exhibits dynamic evolution, while the X is more stable. Finally, mapping short-read sequencing data from >100 great ape individuals revealed the patterns of diversity and selection on their sex chromosomes, demonstrating the utility of these reference assemblies for studies of great ape evolution. These complete sex chromosome assemblies are expected to further inform conservation genetics of nonhuman apes, all of which are endangered species.

Differential usage of DNA modifications in neurons, astrocytes, and microglia

BACKGROUND

Cellular identity is determined partly by cell type-specific epigenomic profiles that regulate gene expression. In neuroscience, there is a pressing need to isolate and characterize the epigenomes of specific CNS cell types in health and disease. In this study, we developed an in vivo tagging mouse model (Camk2a-NuTRAP) for paired isolation of neuronal DNA and RNA without cell sorting and then used this model to assess epigenomic regulation, DNA modifications in particular, of gene expression between neurons and glia.

RESULTS

After validating the cell-specificity of the Camk2a-NuTRAP model, we performed TRAP-RNA-Seq and INTACT-whole genome oxidative bisulfite sequencing (WGoxBS) to assess the neuronal translatome and epigenome in the hippocampus of young mice (4 months old). WGoxBS findings were validated with enzymatic methyl-Seq (EM-Seq) and nanopore sequencing. Comparing neuronal data to microglial and astrocytic data from NuTRAP models, microglia had the highest global mCG levels followed by astrocytes and then neurons, with the opposite pattern observed for hmCG and mCH. Differentially modified regions between cell types were predominantly found within gene bodies and distal intergenic regions, rather than proximal promoters. Across cell types there was a negative correlation between DNA modifications (mCG, mCH, hmCG) and gene expression at proximal promoters. In contrast, a negative correlation of gene body mCG and a positive relationship between distal promoter and gene body hmCG with gene expression was observed. Furthermore, we identified a neuron-specific inverse relationship between mCH and gene expression across promoter and gene body regions.

CONCLUSIONS

Neurons, astrocytes, and microglia demonstrate different genome-wide levels of mCG, hmCG, and mCH that are reproducible across analytical methods. However, modification-gene expression relationships are conserved across cell types. Enrichment of differential modifications across cell types in gene bodies and distal regulatory elements, but not proximal promoters, highlights epigenomic patterning in these regions as potentially greater determinants of cell identity. These findings also demonstrate the importance of differentiating between mC and hmC in neuroepigenomic analyses, as up to 30% of what is conventionally interpreted as mCG can be hmCG, which often has a different relationship to gene expression than mCG.

Mosaic loss of the Y chromosome in human neurodegenerative and oncological diseases

The development of new biomarkers for prediction and early detection of human diseases, as well as for monitoring the response to therapy is one of the most relevant areas of modern human genetics and genomics. Until recently, it was believed that the function of human Y chromosome genes was limited to determining sex and controlling spermatogenesis. Thanks to occurance of large databases of the genome-wide association study (GWAS), there has been a transition to the use of large samples for analyzing genetic changes in both normal and pathological conditions. This has made it possible to assess the association of mosaic aneuploidy of the Y chromosome in somatic cells with a shorter lifespan in men compared to women. Based on data from the UK Biobank, an association was found between mosaic loss of the Y chromosome (mLOY) in peripheral blood leukocytes and the age of men over 70, as well as a number of oncological, cardiac, metabolic, neurodegenerative, and psychiatric diseases. As a result, mLOY in peripheral blood cells has been considered a potential marker of biological age in men and as a marker of certain age-related diseases. Currently, numerous associations have been identified between mLOY and genes based on GWAS and transcriptomes in affected tissues. However, the exact cause of mLOY and the impact and consequences of this phenomenon at the whole organism level have not been established. In particular, it is unclear whether aneuploidy of the Y chromosome in blood cells may affect the development of pathologies that manifest in other organs, such as the brain in Alzheimer's disease, or whether it is a neutral biomarker of general genomic instability. This review examines the main pathologies and genetic factors associated with mLOY, as well as the hypotheses regarding their interplay. Special attention is given to recent studies on mLOY in brain cells in Alzheimer's disease.

The inactive X chromosome accumulates widespread epigenetic variability with age

BACKGROUND

Loss of epigenetic control is a hallmark of aging. Among the most prominent roles of epigenetic mechanisms is the inactivation of one of two copies of the X chromosome in females through DNA methylation. Hence, age-related disruption of X-chromosome inactivation (XCI) may contribute to the aging process in women.

METHODS

We analyzed 9,777 CpGs on the X chromosome in whole blood samples from 2343 females and 1688 males (Illumina 450k methylation array) and replicated findings in duplicate using one whole blood and one purified monocyte data set (in total, 991/924 females/males). We used double generalized linear models to detect age-related differentially methylated CpGs (aDMCs), whose mean methylation level differs with age, and age-related variably methylated CpGs (aVMCs), whose methylation level becomes more variable with age.

RESULTS

In females, aDMCs were relatively uncommon (n = 33) and preferentially occurred in regions known to escape XCI. In contrast, many CpGs (n = 987) were found to display an increased variance with age (aVMCs). Of note, the replication rate of aVMCs was also high in purified monocytes (94%), indicating an independence of cell composition. aVMCs accumulated in CpG islands and regions subject to XCI suggesting that they stemmed from the inactive X. In males, carrying an active copy of the X chromosome only, aDMCs (n = 316) were primarily driven by cell composition, while aVMCs replicated well (95%) but were infrequent (n = 37).

CONCLUSIONS

Our results imply that age-related DNA methylation differences at the inactive X chromosome are dominated by the accumulation of variability.

A novel quantitative targeted analysis of X-chromosome inactivation (XCI) using nanopore sequencing

X-chromosome inactivation (XCI) analyses often assist in diagnostics of X-linked traits, however accurate assessment remains challenging with current methods. We developed a novel strategy using amplification-free Cas9 enrichment and Oxford nanopore technologies sequencing called XCI-ONT, to investigate and rigorously quantify XCI in human androgen receptor gene (AR) and human X-linked retinitis pigmentosa 2 gene (RP2). XCI-ONT measures methylation over 116 CpGs in AR and 58 CpGs in RP2, and separate parental X-chromosomes without PCR bias. We show the usefulness of the XCI-ONT strategy over the PCR-based golden standard XCI technique that only investigates one or two CpGs per gene. The results highlight the limitations of using the golden standard technique when the XCI pattern is partially skewed and the advantages of XCI-ONT to rigorously quantify XCI. This study provides a universal XCI-method on DNA, which is highly valuable in clinical and research framework of X-linked traits.

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 .

A systematic review of neuroimaging epigenetic research: calling for an increased focus on development

Epigenetic mechanisms, such as DNA methylation (DNAm), have gained increasing attention as potential biomarkers and mechanisms underlying risk for neurodevelopmental, psychiatric and other brain-based disorders. Yet, surprisingly little is known about the extent to which DNAm is linked to individual differences in the brain itself, and how these associations may unfold across development - a time of life when many of these disorders emerge. Here, we systematically review evidence from the nascent field of Neuroimaging Epigenetics, combining structural or functional neuroimaging measures with DNAm, and the extent to which the developmental period (birth to adolescence) is represented in these studies. We identified 111 articles published between 2011-2021, out of which only a minority (21%) included samples under 18 years of age. Most studies were cross-sectional (85%), employed a candidate-gene approach (67%), and examined DNAm-brain associations in the context of health and behavioral outcomes (75%). Nearly half incorporated genetic data, and a fourth investigated environmental influences. Overall, studies support a link between peripheral DNAm and brain imaging measures, but there is little consistency in specific findings and it remains unclear whether DNAm markers present a cause, correlate or consequence of brain alterations. Overall, there is large heterogeneity in sample characteristics, peripheral tissue and brain outcome examined as well as the methods used. Sample sizes were generally low to moderate (median nall = 98, ndevelopmental = 80), and attempts at replication or meta-analysis were rare. Based on the strengths and weaknesses of existing studies, we propose three recommendations on how advance the field of Neuroimaging Epigenetics. We advocate for: (1) a greater focus on developmentally oriented research (i.e. pre-birth to adolescence); (2) the analysis of large, prospective, pediatric cohorts with repeated measures of DNAm and imaging to assess directionality; and (3) collaborative, interdisciplinary science to identify robust signals, triangulate findings and enhance translational potential.

Differential usage of DNA modifications in neurons, astrocytes, and microglia

Background

Cellular identity is determined partly by cell type-specific epigenomic profiles that regulate gene expression. In neuroscience, there is a pressing need to isolate and characterize the epigenomes of specific CNS cell types in health and disease. This is especially true as for DNA modifications where most data are derived from bisulfite sequencing that cannot differentiate between DNA methylation and hydroxymethylation. In this study, we developed an in vivo tagging mouse model (Camk2a-NuTRAP) for paired isolation of neuronal DNA and RNA without cell sorting and then used this model to assess epigenomic regulation of gene expression between neurons and glia.

Results

After validating the cell-specificity of the Camk2a-NuTRAP model, we performed TRAP-RNA-Seq and INTACT whole genome oxidative bisulfite sequencing to assess the neuronal translatome and epigenome in the hippocampus of young mice (3 months old). These data were then compared to microglial and astrocytic data from NuTRAP models. When comparing the different cell types, microglia had the highest global mCG levels followed by astrocytes and then neurons, with the opposite pattern observed for hmCG and mCH. Differentially modified regions between cell types were predominantly found within gene bodies and distal intergenic regions, with limited differences occurring within proximal promoters. Across cell types there was a negative correlation between DNA modifications (mCG, mCH, hmCG) and gene expression at proximal promoters. In contrast, a negative correlation of mCG with gene expression within the gene body while a positive relationship between distal promoter and gene body hmCG and gene expression was observed. Furthermore, we identified a neuron-specific inverse relationship between mCH and gene expression across promoter and gene body regions.

Conclusions

In this study, we identified differential usage of DNA modifications across CNS cell types, and assessed the relationship between DNA modifications and gene expression in neurons and glia. Despite having different global levels, the general modification-gene expression relationship was conserved across cell types. The enrichment of differential modifications in gene bodies and distal regulatory elements, but not proximal promoters, across cell types highlights epigenomic patterning in these regions as potentially greater determinants of cell identity.

Epigenetic regulation of plastin 3 expression by the macrosatellite DXZ4 and the transcriptional regulator CHD4

Dysregulated Plastin 3 (PLS3) levels associate with a wide range of skeletal and neuromuscular disorders and the most common types of solid and hematopoietic cancer. Most importantly, PLS3 overexpression protects against spinal muscular atrophy. Despite its crucial role in F-actin dynamics in healthy cells and its involvement in many diseases, the mechanisms that regulate PLS3 expression are unknown. Interestingly, PLS3 is an X-linked gene and all asymptomatic SMN1-deleted individuals in SMA-discordant families who exhibit PLS3 upregulation are female, suggesting that PLS3 may escape X chromosome inactivation. To elucidate mechanisms contributing to PLS3 regulation, we performed a multi-omics analysis in two SMA-discordant families using lymphoblastoid cell lines and iPSC-derived spinal motor neurons originated from fibroblasts. We show that PLS3 tissue-specifically escapes X-inactivation. PLS3 is located ∼500 kb proximal to the DXZ4 macrosatellite, which is essential for X chromosome inactivation. By applying molecular combing in a total of 25 lymphoblastoid cell lines (asymptomatic individuals, individuals with SMA, control subjects) with variable PLS3 expression, we found a significant correlation between the copy number of DXZ4 monomers and PLS3 levels. Additionally, we identified chromodomain helicase DNA binding protein 4 (CHD4) as an epigenetic transcriptional regulator of PLS3 and validated co-regulation of the two genes by siRNA-mediated knock-down and overexpression of CHD4. We show that CHD4 binds the PLS3 promoter by performing chromatin immunoprecipitation and that CHD4/NuRD activates the transcription of PLS3 by dual-luciferase promoter assays. Thus, we provide evidence for a multilevel epigenetic regulation of PLS3 that may help to understand the protective or disease-associated PLS3 dysregulation.

Genome-wide DNA methylation profiling after Ayurveda intervention to bronchial asthmatics identifies differential methylation in several transcription factors with immune process related function

BACKGROUND

The Indian traditional medicinal system, Ayurveda, describes several lifestyle practices, processes and medicines as an intervention to treat asthma. Rasayana therapy is one of them and although these treatment modules show improvement in bronchial asthma, their mechanism of action, particularly the effect on DNA methylation, is largely understudied.

OBJECTIVES

Our study aimed at identifying the contribution of DNA methylation changes in modulating bronchial asthma phenotype upon Ayurveda intervention.

MATERIALS AND METHODS

In this study, genome-wide methylation profiling in peripheral blood DNA of healthy controls and bronchial asthmatics before (BT) and after (AT) Ayurveda treatment was performed using array-based profiling of reference-independent methylation status (aPRIMES) coupled to microarray technique.

RESULTS

We identified 4820 treatment-associated DNA methylation signatures (TADS) and 11,643 asthma-associated DNA methylation signatures (AADS), differentially methylated [FDR (≤0.1) adjusted p-values] in AT and HC groups respectively, compared to BT group. Neurotrophin TRK receptor signaling pathway was significantly enriched for differentially methylated genes in bronchial asthmatics, compared to AT and HC subjects. Additionally, we identified over 100 differentially methylated immune-related genes located in the promoter/5'-UTR regions of TADS and AADS. Various immediate-early response and immune regulatory genes with functions such as transcription factor activity (FOXD1, FOXD2, GATA6, HOXA3, HOXA5, MZF1, NFATC1, NKX2-2, NKX2-3, RUNX1, KLF11), G-protein coupled receptor activity (CXCR4, PTGER4), G-protein coupled receptor binding (UCN), DNA binding (JARID2, EBF2, SOX9), SNARE binding (CAPN10), transmembrane signaling receptor activity (GP1BB), integrin binding (ITGA6), calcium ion binding (PCDHGA12), actin binding (TRPM7, PANX1, TPM1), receptor tyrosine kinase binding (PIK3R2), receptor activity (GDNF), histone methyltransferase activity (MLL5), and catalytic activity (TSTA3) were found to show consistent methylation status between AT and HC group in microarray data.

CONCLUSIONS

Our study reports the DNA methylation-regulated genes in bronchial asthmatics showing improvement in symptoms after Ayurveda intervention. DNA methylation regulation in the identified genes and pathways represents the Ayurveda intervention responsive genes and may be further explored as diagnostic, prognostic, and therapeutic biomarkers for bronchial asthma in peripheral blood.

DNA methylation-based profiling of horse archaeological remains for age-at-death and castration

Age profiling of archaeological bone assemblages can inform on past animal management practices, but is limited by the fragmentary nature of the fossil record and the lack of universal skeletal markers for age. DNA methylation clocks offer new, albeit challenging, alternatives for estimating the age-at-death of ancient individuals. Here, we take advantage of the availability of a DNA methylation clock based on 31,836 CpG sites and dental age markers in horses to assess age predictions in 84 ancient remains. We evaluate our approach using whole-genome sequencing data and develop a capture assay providing reliable estimates for only a fraction of the cost. We also leverage DNA methylation patterns to assess castration practice in the past. Our work opens for a deeper characterization of past husbandry and ritual practices and holds the potential to reveal age mortality profiles in ancient societies, once extended to human remains.

Purified regenerating retinal neurons reveal regulatory role of DNA methylation-mediated Na+/K+-ATPase in murine axon regeneration

While embryonic mammalian central nervous system (CNS) axons readily grow and differentiate, only a minority of fully differentiated mature CNS neurons are able to regenerate injured axons, leading to stunted functional recovery after injury and disease. To delineate DNA methylation changes specifically associated with axon regeneration, we used a Fluorescent-Activated Cell Sorting (FACS)-based methodology in a rat optic nerve transection model to segregate the injured retinal ganglion cells (RGCs) into regenerating and non-regenerating cell populations. Whole-genome DNA methylation profiling of these purified neurons revealed genes and pathways linked to mammalian RGC regeneration. Moreover, whole-methylome sequencing of purified uninjured adult and embryonic RGCs identified embryonic molecular profiles reactivated after injury in mature neurons, and others that correlate specifically with embryonic or adult axon growth, but not both. The results highlight the contribution to both embryonic growth and adult axon regeneration of subunits encoding the Na⁺/K⁺-ATPase. In turn, both biochemical and genetic inhibition of the Na⁺/K⁺-ATPase pump significantly reduced RGC axon regeneration. These data provide critical molecular insights into mammalian CNS axon regeneration, pinpoint the Na⁺/K⁺-ATPase as a key regulator of regeneration of injured mature CNS axons, and suggest that successful regeneration requires, in part, reactivation of embryonic signals.

The epigenetic players and the chromatin marks involved in the articular cartilage during osteoarthritis

Epigenetics defines the modifications of the genome that do not involve a change in the nucleotide sequence of DNA. These modifications constitute a mechanism of gene regulation poorly explored in the context of cartilage physiology. They are now intensively studied by the scientific community working on articular cartilage and its related pathology such as osteoarthritis. Indeed, epigenetic regulations can control the expression of crucial gene in the chondrocytes, the only resident cells of cartilage. Some epigenetic changes are considered as a possible cause of the abnormal gene expression and the subsequent alteration of the chondrocyte phenotype (hypertrophy, proliferation, senescence…) as observed in osteoarthritic cartilage. Osteoarthritis is a joint pathology, which results in impaired extracellular matrix homeostasis and leads ultimately to the progressive destruction of cartilage. To date, there is no pharmacological treatment and the exact causes have yet to be defined. Given that the epigenetic modifying enzymes can be controlled by pharmacological inhibitors, it is thus crucial to describe the epigenetic marks that enable the normal expression of extracellular matrix encoding genes, and those associated with the abnormal gene expression such as degradative enzyme or inflammatory cytokines encoding genes. In this review, only the DNA methylation and histone modifications will be detailed with regard to normal and osteoarthritic cartilage. Although frequently referred as epigenetic mechanisms, the regulatory mechanisms involving microRNAs will not be discussed. Altogether, this review will show how this nascent field influences our understanding of the pathogenesis of OA in terms of diagnosis and how controlling the epigenetic marks can help defining epigenetic therapies.

Who's afraid of the X? Incorporating the X and Y chromosomes into the analysis of DNA methylation array data

BACKGROUND

Many human disease phenotypes manifest differently by sex, making the development of methods for incorporating X and Y-chromosome data into analyses vital. Unfortunately, X and Y chromosome data are frequently excluded from large-scale analyses of the human genome and epigenome due to analytical complexity associated with sex chromosome dosage differences between XX and XY individuals, and the impact of X-chromosome inactivation (XCI) on the epigenome. As such, little attention has been given to considering the methods by which sex chromosome data may be included in analyses of DNA methylation (DNAme) array data.

RESULTS

With Illumina Infinium HumanMethylation450 DNAme array data from 634 placental samples, we investigated the effects of probe filtering, normalization, and batch correction on DNAme data from the X and Y chromosomes. Processing steps were evaluated in both mixed-sex and sex-stratified subsets of the analysis cohort to identify whether including both sexes impacted processing results. We found that identification of probes that have a high detection p-value, or that are non-variable, should be performed in sex-stratified data subsets to avoid over- and under-estimation of the quantity of probes eligible for removal, respectively. All normalization techniques investigated returned X and Y DNAme data that were highly correlated with the raw data from the same samples. We found no difference in batch correction results after application to mixed-sex or sex-stratified cohorts. Additionally, we identify two analytical methods suitable for XY chromosome data, the choice between which should be guided by the research question of interest, and we performed a proof-of-concept analysis studying differential DNAme on the X and Y chromosome in the context of placental acute chorioamnionitis. Finally, we provide an annotation of probe types that may be desirable to filter in X and Y chromosome analyses, including probes in repetitive elements, the X-transposed region, and cancer-testis gene promoters.

CONCLUSION

While there may be no single "best" approach for analyzing DNAme array data from the X and Y chromosome, analysts must consider key factors during processing and analysis of sex chromosome data to accommodate the underlying biology of these chromosomes, and the technical limitations of DNA methylation arrays.

XPO1-Mediated EIF1AX Cytoplasmic Relocation Promotes Tumor Migration and Invasion in Endometrial Carcinoma

Dysregulation of eukaryotic translation initiation factor 1A, X-linked (EIF1AX), has been implicated in the pathogenesis of some cancers. However, the role of EIF1AX in endometrial carcinoma (EC) remains unknown. We investigated the EIF1AX expression in EC patients and assessed its tumorigenesis-associated function and nucleocytoplasmic transport mechanism in vitro and in vivo. The results indicated that the cytoplasmic EIF1AX expression showed a gradual increase when going from endometrium normal tissue, simple endometrial hyperplasia, complex endometrial hyperplasia, and endometrial atypical hyperplasia to EC, while vice versa for the nuclear EIF1AX expression. In addition, the cytoplasmic EIF1AX expression was positively correlated with histologic type, high International Federation of Gynecology and Obstetrics (FIGO) grade, advanced FIGO stage, deeper infiltration, high Ki67 index, and shorter recurrence-free survival in EC patients. In vitro, short hairpin RNA-mediated EIF1AX depletion or SV40NLS-mediated EIF1AX import into the nucleus in multiple human EC cells potently suppressed cell migration and invasion, epithelial-mesenchymal transition, and lung metastasis. Moreover, exportin 1 induced the transport of EIF1AX from the nucleus to the cytoplasm that could be inhibited by leptomycin B treatment or the mutation in the EIF1AX location sequence. These results demonstrate that cytoplasmic EIF1AX may play a key role in the incidence and promotion of EC, and thus, targeting EIF1AX or its nucleocytoplasmic transport process may offer an effective new therapeutic approach to EC.

Somatic XIST activation and features of X chromosome inactivation in male human cancers

Expression of the non-coding RNA XIST is essential for initiating X chromosome inactivation (XCI) during early development in female mammals. As the main function of XCI is to enable dosage compensation of chromosome X genes between the sexes, XCI and XIST expression are generally absent in male normal tissues, except in germ cells and in individuals with supernumerary X chromosomes. Via a systematic analysis of public sequencing data of both cancerous and normal tissues, we report that XIST is somatically activated in a subset of male human cancers across diverse lineages. Some of these cancers display hallmarks of XCI, including silencing of gene expression, reduced chromatin accessibility, and increased DNA methylation across chromosome X, suggesting that the developmentally restricted, female-specific program of XCI can be somatically accessed in male cancers.

DMRscaler: a scale-aware method to identify regions of differential DNA methylation spanning basepair to multi-megabase features

BACKGROUND

Pathogenic mutations in genes that control chromatin function have been implicated in rare genetic syndromes. These chromatin modifiers exhibit extraordinary diversity in the scale of the epigenetic changes they affect, from single basepair modifications by DNMT1 to whole genome structural changes by PRM1/2. Patterns of DNA methylation are related to a diverse set of epigenetic features across this full range of epigenetic scale, making DNA methylation valuable for mapping regions of general epigenetic dysregulation. However, existing methods are unable to accurately identify regions of differential methylation across this full range of epigenetic scale directly from DNA methylation data.

RESULTS

To address this, we developed DMRscaler, a novel method that uses an iterative windowing procedure to capture regions of differential DNA methylation (DMRs) ranging in size from single basepairs to whole chromosomes. We benchmarked DMRscaler against several DMR callers in simulated and natural data comparing XX and XY peripheral blood samples. DMRscaler was the only method that accurately called DMRs ranging in size from 100 bp to 1 Mb (pearson's r = 0.94) and up to 152 Mb on the X-chromosome. We then analyzed methylation data from rare-disease cohorts that harbor chromatin modifier gene mutations in NSD1, EZH2, and KAT6A where DMRscaler identified novel DMRs spanning gene clusters involved in development.

CONCLUSION

Taken together, our results show DMRscaler is uniquely able to capture the size of DMR features across the full range of epigenetic scale and identify novel, co-regulated regions that drive epigenetic dysregulation in human disease.

Comprehensive analysis for clarifying transcriptomics landscapes of spread through air spaces in lung adenocarcinoma

Patients with spread through air spaces (STAS) have worse postoperative survival and a higher recurrence rate in lung adenocarcinoma, even in the earliest phases of the disease. At present, the molecular pathogenesis of STAS is not well understood. Therefore, to illustrate the underlying pathogenic mechanism of STAS, we accomplished a comprehensive analysis of a microarray dataset of STAS. Differential expression analysis revealed 841 differentially expressed genes (DEGs) between STAS_positive and STAS_negative groups. Additionally, we acquired two hub genes associated with survival. Gene set variation analysis (GSVA) confirmed that the main differential signaling pathways between the two groups were hypoxia VHL targets, PKC, and pyrimidine metabolism pathways. Analysis of immune activity showed that the increased expression of MHC-class-Ⅰ was observed in the STAS_positive group. These findings provided novel insights for a better knowledge of pathogenic mechanisms and potential therapeutic markers for STAS treatment.

Differential Regulation of Mouse Hippocampal Gene Expression Sex Differences by Chromosomal Content and Gonadal Sex

Common neurological disorders, like Alzheimer's disease (AD), multiple sclerosis (MS), and autism, display profound sex differences in prevalence and clinical presentation. However, sex differences in the brain with health and disease are often overlooked in experimental models. Sex effects originate, directly or indirectly, from hormonal or sex chromosomal mechanisms. To delineate the contributions of genetic sex (XX v. XY) versus gonadal sex (ovaries v. testes) to the epigenomic regulation of hippocampal sex differences, we used the Four Core Genotypes (FCG) mouse model which uncouples chromosomal and gonadal sex. Transcriptomic and epigenomic analyses of ~ 12-month-old FCG mouse hippocampus, revealed genomic context-specific regulatory effects of genotypic and gonadal sex on X- and autosome-encoded gene expression and DNA modification patterns. X-chromosomal epigenomic patterns, classically associated with X-inactivation, were established almost entirely by genotypic sex, independent of gonadal sex. Differences in X-chromosome methylation were primarily localized to gene regulatory regions including promoters, CpG islands, CTCF binding sites, and active/poised chromatin, with an inverse relationship between methylation and gene expression. Autosomal gene expression demonstrated regulation by both genotypic and gonadal sex, particularly in immune processes. These data demonstrate an important regulatory role of sex chromosomes, independent of gonadal sex, on sex-biased hippocampal transcriptomic and epigenomic profiles. Future studies will need to further interrogate specific CNS cell types, identify the mechanisms by which sex chromosomes regulate autosomes, and differentiate organizational from activational hormonal effects.

The emerging role of non-coding RNAs in the epigenetic regulation of pediatric cancers

Epigenetics is a process that involves the regulation of gene expression without altering the sequence of DNA. Numerous studies have documented that epigenetic mechanisms play a critical role in cell growth, differentiation, and cancer over the past decade. The well-known epigenetic modifications are either on DNA or at the histone proteins. Although several studies have focused on regulating gene expression by non-coding RNAs, the current understanding of their biological functions in various human diseases, particularly in cancers, is inadequate. Only about two percent of DNA is involved in coding the protein-coding genes, and leaving the rest 98 percent is non-coding and the scientific community regarded as junk or noise with no known purpose. Most non-coding RNAs are derived from such junk DNA and are known to be involved in various signaling pathways involving cancer initiation, progression, and the development of therapy resistance in many human cancer types. Recent studies have suggested that non-coding RNAs, especially microRNAs, piwi-interactingRNAs, and long non-coding RNAs, play a significant role in controlling epigenetic mechanism(s), indicating the potential effect of epigenetic modulation of non-coding RNAs on cancer progression. In this review article, we briefly presented epigenetic marks' characteristics, crosstalk between epigenetic modifications and microRNAs, piwi-interactingRNAs, and long non-coding RNAs to uncover the effect on the phenotype of pediatric cancers. Further, current knowledge on understanding the RNA epigenetics will help design novel therapeutics that target epigenetic regulatory networks to benefit cancer patients in the clinic.

InterpolatedXY: a two-step strategy to normalize DNA methylation microarray data avoiding sex bias

MOTIVATION

Data normalization is an essential step to reduce technical variation within and between arrays. Due to the different karyotypes and the effects of X chromosome inactivation, females and males exhibit distinct methylation patterns on sex chromosomes; thus, it poses a significant challenge to normalize sex chromosome data without introducing bias. Currently, existing methods do not provide unbiased solutions to normalize sex chromosome data, usually, they just process autosomal and sex chromosomes indiscriminately.

RESULTS

Here, we demonstrate that ignoring this sex difference will lead to introducing artificial sex bias, especially for thousands of autosomal CpGs. We present a novel two-step strategy (interpolatedXY) to address this issue, which is applicable to all quantile-based normalization methods. By this new strategy, the autosomal CpGs are first normalized independently by conventional methods, such as funnorm or dasen; then the corrected methylation values of sex chromosome-linked CpGs are estimated as the weighted average of their nearest neighbors on autosomes. The proposed two-step strategy can also be applied to other non-quantile-based normalization methods, as well as other array-based data types. Moreover, we propose a useful concept: the sex explained fraction of variance, to quantitatively measure the normalization effect.

AVAILABILITY AND IMPLEMENTATION

The proposed methods are available by calling the function 'adjustedDasen' or 'adjustedFunnorm' in the latest wateRmelon package (https://github.com/schalkwyk/wateRmelon), with methods compatible with all the major workflows, including minfi.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

XIST loss impairs mammary stem cell differentiation and increases tumorigenicity through Mediator hyperactivation

X inactivation (XCI) is triggered by upregulation of XIST, which coats the chromosome in cis, promoting formation of a heterochromatic domain (Xi). XIST role beyond initiation of XCI is only beginning to be elucidated. Here, we demonstrate that XIST loss impairs differentiation of human mammary stem cells (MaSCs) and promotes emergence of highly tumorigenic and metastatic carcinomas. On the Xi, XIST deficiency triggers epigenetic changes and reactivation of genes overlapping Polycomb domains, including Mediator subunit MED14. MED14 overdosage results in increased Mediator levels and hyperactivation of the MaSC enhancer landscape and transcriptional program, making differentiation less favorable. We further demonstrate that loss of XIST and Xi transcriptional instability is common among human breast tumors of poor prognosis. We conclude that XIST is a gatekeeper of human mammary epithelium homeostasis, thus unveiling a paradigm in the control of somatic cell identity with potential consequences for our understanding of gender-specific malignancies.

Considerations and challenges for sex-aware drug repurposing

Sex differences are essential factors in disease etiology and manifestation in many diseases such as cardiovascular disease, cancer, and neurodegeneration [33]. The biological influence of sex differences (including genomic, epigenetic, hormonal, immunological, and metabolic differences between males and females) and the lack of biomedical studies considering sex differences in their study design has led to several policies. For example, the National Institute of Health's (NIH) sex as a biological variable (SABV) and Sex and Gender Equity in Research (SAGER) policies to motivate researchers to consider sex differences [204]. However, drug repurposing, a promising alternative to traditional drug discovery by identifying novel uses for FDA-approved drugs, lacks sex-aware methods that can improve the identification of drugs that have sex-specific responses [7, 11, 14, 33]. Sex-aware drug repurposing methods either select drug candidates that are more efficacious in one sex or deprioritize drug candidates based on if they are predicted to cause a sex-bias adverse event (SBAE), unintended therapeutic effects that are more likely to occur in one sex. Computational drug repurposing methods are encouraging approaches to develop for sex-aware drug repurposing because they can prioritize sex-specific drug candidates or SBAEs at lower cost and time than traditional drug discovery. Sex-aware methods currently exist for clinical, genomic, and transcriptomic information [1, 7, 155]. They have not expanded to other data types, such as DNA variation, which has been beneficial in other drug repurposing methods that do not consider sex [114]. Additionally, some sex-aware methods suffer from poorer performance because a disproportionate number of male and female samples are available to train computational methods [7]. However, there is development potential for several different categories (i.e., data mining, ligand binding predictions, molecular associations, and networks). Low-dimensional representations of molecular association and network approaches are also especially promising candidates for future sex-aware drug repurposing methodologies because they reduce the multiple hypothesis testing burden and capture sex-specific variation better than the other methods [151, 159]. Here we review how sex influences drug response, the current state of drug repurposing including with respect to sex-bias drug response, and how model organism study design choices influence drug repurposing validation.

A comparative methylome analysis reveals conservation and divergence of DNA methylation patterns and functions in vertebrates

Background

Cytosine DNA methylation is a heritable epigenetic mark present in most eukaryotic groups. While the patterns and functions of DNA methylation have been extensively studied in mouse and human, their conservation in other vertebrates remains poorly explored. In this study, we interrogated the distribution and function of DNA methylation in primary fibroblasts of seven vertebrate species including bio-medical models and livestock species (human, mouse, rabbit, dog, cow, pig, and chicken).

Results

Our data highlight both divergence and conservation of DNA methylation patterns and functions. We show that the chicken genome is hypomethylated compared to other vertebrates. Furthermore, compared to mouse, other species show a higher frequency of methylation of CpG-rich DNA. We reveal the conservation of large unmethylated valleys and patterns of DNA methylation associated with X-chromosome inactivation through vertebrate evolution and make predictions of conserved sets of imprinted genes across mammals. Finally, using chemical inhibition of DNA methylation, we show that the silencing of germline genes and endogenous retroviruses (ERVs) are conserved functions of DNA methylation in vertebrates.

Conclusions

Our data highlight conserved properties of DNA methylation in vertebrate genomes but at the same time point to differences between mouse and other vertebrate species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01270-x.

Epigenetic Regulation in Uterine Fibroids-The Role of Ten-Eleven Translocation Enzymes and Their Potential Therapeutic Application

Uterine fibroids (UFs) are monoclonal, benign tumors that contain abnormal smooth muscle cells and the accumulation of extracellular matrix (ECM). Although benign, UFs are a major source of gynecologic and reproductive dysfunction, ranging from menorrhagia and pelvic pain to infertility, recurrent miscarriage, and preterm labor. Many risk factors are involved in the pathogenesis of UFs via genetic and epigenetic mechanisms. The latter involving DNA methylation and demethylation reactions provide specific DNA methylation patterns that regulate gene expression. Active DNA demethylation reactions mediated by ten-eleven translocation proteins (TETs) and elevated levels of 5-hydroxymethylcytosine have been suggested to be involved in UF formation. This review paper summarizes the main findings regarding the function of TET enzymes and their activity dysregulation that may trigger the development of UFs. Understanding the role that epigenetics plays in the pathogenesis of UFs may possibly lead to a new type of pharmacological fertility-sparing treatment method.

Use of the FMR1 Gene Methylation Status to Assess the X-Chromosome Inactivation Pattern: A Stepwise Analysis

X-chromosome inactivation (XCI) is a developmental process to compensate the imbalance in the dosage of X-chromosomal genes in females. A skewing of the XCI pattern may suggest a carrier status for an X-linked disease or explain the presence of a severe phenotype. In these cases, it is important to determine the XCI pattern, conventionally using the gold standard Human Androgen-Receptor Assay (HUMARA), based on the analysis of the methylation status at a polymorphic CAG region in the first exon of the human androgen receptor gene (AR). The aim of this study was to evaluate whether the methylation status of the fragile mental retardation protein translational regulator gene (FMR1) can provide an XCI pattern similar to that obtained by HUMARA. A set of 48 female carriers of FMR1 gene normal-sized alleles was examined using two assays: HUMARA and a FMR1 methylation PCR (mPCR). Ranges were defined to establish the XCI pattern using the methylation pattern of the FMR1 gene by mPCR. Overall, a 77% concordance of the XCI patterns was obtained between the two assays, which led us to propose a set of key points and a stepwise analysis towards obtaining an accurate result for the XCI pattern and to minimize the underlying pitfalls.

Mining RNAseq data reveals dynamic metaboloepigenetic profiles in human, mouse and bovine pre-implantation embryos

Metaboloepigenetic regulation has been reported in stem cells, germ cells, and tumor cells. Embryonic metaboloepigenetics, however, have just begun to be described. Here we analyzed RNAseq data to characterize the metaboloepigenetic profiles of human, mouse, and bovine pre-implantation embryos. In embryos, metaboloepigenetic reprogramming was species-specific, varied with the developmental stage and was disrupted with in vitro culture. Metabolic pathways and gene expressions were strongly correlated with early embryo DNA methylation and were changed with in vitro culture. Although the idea that the in vitro environment may influence development is not new, there has been little progress on improving pregnancy rates after decades using in vitro fertilization. Hence, the present data will contribute to understanding how the in vitro manipulation affects the metaboloepigenetic status of early embryos, which can be used to establish culture strategies aimed at improving the in vitro environment and, consequently, pregnancy rates and offspring health.

•

Embryonic metaboloepigenetic reprogramming is stage- and species-specific

•

In vitro culture disrupts the in vivo embryonic metaboloepigenetic reprogramming

•

Metabolic genes and pathways are highly correlated with embryo methylome

Navigating the hydroxymethylome: experimental biases and quality control tools for the tandem bisulfite and oxidative bisulfite Illumina microarrays

Aim: Tandem bisulfite (BS) and oxidative bisulfite (oxBS) conversion on DNA followed by hybridization to Infinium HumanMethylation BeadChips allows nucleotide resolution of 5-hydroxymethylcytosine genome-wide. Here, the authors compared data quality acquired from BS-treated and oxBS-treated samples. Materials & methods: Raw BeadArray data from 417 pairs of samples across 12 independent datasets were included in the study. Probe call rates were compared between paired BS and oxBS treatments controlling for technical variables. Results: oxBS-treated samples had a significantly lower call rate. Among technical variables, DNA-specific extraction kits performed better with higher call rates after oxBS conversion. Conclusion: The authors emphasize the importance of quality control during oxBS conversion to minimize information loss and recommend using a DNA-specific extraction kit for DNA extraction and an oxBSQC package for data preprocessing.

Variation in Genotype and DNA Methylation Patterns Based on Alcohol Use and CVD in the Korean Genome and Epidemiology Study (KoGES)

Alcohol consumption can increase the risk of chronic diseases, such as myocardial infarction, coronary artery disease, hyperlipidemia, and hypertension. We aimed to assess the association between genotype, DNA methylation patterns, alcohol consumption, and chronic diseases in Korean population. We analyzed 8840 subjects for genotypes and 446 for DNA methylation among the 9351 subjects from the Korean Genome and Epidemiology Study (KoGES). We further divided both groups into two sub-groups according to the presence/absence of chronic diseases. We selected genes whose methylation varied significantly with alcohol consumption, and visualized genotype and DNA methylation patterns specific to each group. Genome-wide association study (GWAS) revealed single nucleotide polymorphisms (SNPs) rs2074356 and rs11066280 in HECT domain E3 ubiquitin protein ligase 4 (HECTD4) to be significantly associated with alcohol consumption in both the presence. The rs12229654 genotype also displayed significantly different patterns with alcohol consumption. Furthermore, we retrieved differentially methylated regions (DMRs) from four groups based on sex and chronic diseases and compared them by drinking status. In genotype analysis, cardiovascular diseases (CVDs) showed a higher proportion in drinker than in non-drinker, but not in DMR analysis. Additionally, we analyzed the enriched Gene Ontology terms and Kyoto Gene and Genome Encyclopedia (KEGG) pathways and visualized the network, heatmap, and upset plot. We show that the pattern of DNA methylation associated with CVD is strongly influenced by alcoholism. Overall, this study identified genetic and epigenetic variants influenced by alcohol consumption and chronic diseases.

Epigenetic-Mediated Regulation of Gene Expression for Biological Control and Cancer: Cell and Tissue Structure, Function, and Phenotype

Epigenetic gene regulatory mechanisms play a central role in the biological control of cell and tissue structure, function, and phenotype. Identification of epigenetic dysregulation in cancer provides mechanistic into tumor initiation and progression and may prove valuable for a variety of clinical applications. We present an overview of epigenetically driven mechanisms that are obligatory for physiological regulation and parameters of epigenetic control that are modified in tumor cells. The interrelationship between nuclear structure and function is not mutually exclusive but synergistic. We explore concepts influencing the maintenance of chromatin structures, including phase separation, recognition signals, factors that mediate enhancer-promoter looping, and insulation and how these are altered during the cell cycle and in cancer. Understanding how these processes are altered in cancer provides a potential for advancing capabilities for the diagnosis and identification of novel therapeutic targets.

Methyl-SNP-seq reveals dual readouts of methylome and variome at molecule resolution while enabling target enrichment

Covalent modifications of genomic DNA are crucial for most organisms to survive. Amplicon-based high-throughput sequencing technologies erase all DNA modifications to retain only sequence information for the four canonical nucleobases, necessitating specialized technologies for ascertaining epigenetic information. To also capture base modification information, we developed Methyl-SNP-seq, a technology that takes advantage of the complementarity of the double helix to extract the methylation and original sequence information from a single DNA molecule. More specifically, Methyl-SNP-seq uses bisulfite conversion of one of the strands to identify cytosine methylation while retaining the original four-bases sequence information on the other strand. As both strands are locked together to link the dual readouts on a single paired-end read, Methyl-SNP-seq allows detecting the methylation status of any DNA even without a reference genome. Because one of the strands retains the original four nucleotide composition, Methyl-SNP-seq can also be used in conjunction with standard sequence-specific probes for targeted enrichment and amplification. We show the usefulness of this technology in a broad spectrum of applications ranging from allele-specific methylation analysis in humans to identification of methyltransferase specificity in complex bacterial communities.

Modified cytosines versus cytosine in a DNA polymerase: retrieving thermodynamic and kinetic constants at the single molecule level

DNA methylation plays key roles in various areas, such as gene expression, regulation, epigenetics, and cancers. Since 5-methylcytosine (5mC) is commonly present in methylated DNA, characterizing the binding kinetics and thermodynamics of the nucleotide to the enzymatic pocket can help to understand the DNA replication process. Furthermore, 5-carboxycytosine (5caC) is a form that appears through the iterative oxidation of 5mC, and its effect on the DNA replication process is still not well known. Here, we immobilized a DNA polymerase (DNAP) with an orientation control on a tip of an atomic force microscope (AFM), and observed the interaction between the immobilized deoxyguanosine triphosphate (dGTP) on the surface and the DNAP in the presence of a DNA duplex. The interaction probability increased as the concentration of the DNA strand, and the affinity constant between the DNAP and DNA was obtained by fitting the change. Increasing the concentration of dGTP in solution diminished the interaction probability, and a fitting allowed us to retrieve the affinity constant between dGTP and the DNAP holding the DNA in the reaction pocket. Because the dissociation constant could be obtained through the loading rate dependence of the unbinding force value, both affinity and kinetic constants for cytosine (C), 5mC, and 5caC in the DNAP were compared in the light of the steric and electronic effect of the substituents at 5-position of cytosine.

Associations between infant sex and DNA methylation across umbilical cord blood, artery, and placenta samples

DNA methylation (DNAm) is vulnerable to dysregulation by environmental exposures during epigenetic reprogramming that occurs in embryogenesis. Sexual dimorphism in environmentally induced DNAm dysregulation has been identified and therefore it is important to understand sex-specific DNAm patterns. DNAm at several autosomal sites has been consistently associated with sex in cord blood and placental foetal tissues. However, there is limited research comparing sex-specific DNAm across tissues, particularly differentially methylated regions (DMRs). This study leverages DNAm data measured using the Illumina HumanMethylation450 BeadChip in cord blood (N = 179), placenta (N = 229), and umbilical artery samples (N = 229) in the PRogramming of Intergenerational Stress Mechanisms (PRISM) cohort to identify autosomal DMRs and differentially methylated positions (DMPs). A replication analyses was conducted in an independent cohort (GEO Accession GSE129841). We identified 183, 257, and 419 DMRs and 2119, 2281, and 3405 DMPs (p_Bonferroni < 0.05) in cord blood, placenta, and artery samples, respectively. Thirty-nine DMRs overlapped in all three tissues, overlapping with genes involved in spermatogenesis (NKAPL, PIWIL2 and AURKC) and X-inactivation (LRIF1). In replication analysis, 85% of DMRs overlapped with those identified in PRISM. Overall, DMRs and DMPs had higher methylation levels among females in cord blood and artery samples, but higher methylation levels among males in placenta samples. Further research is necessary to understand biological mechanisms that contribute to differences in sex-specific DNAm signatures across tissues, as well as to determine if sexual dimorphism in the epigenome impacts response to environmental stressors.

Interaction of Male Specific Lethal complex and genomic imbalance on global gene expression in Drosophila

The inverse dosage effect caused by chromosome number variations shows global consequences in genomic imbalance including sexual dimorphism and an X chromosome-specific response. To investigate the relationship of the MSL complex to genomic imbalance, we over-expressed MSL2 in autosomal and sex chromosomal aneuploids, and analyzed the different transcriptomes. Some candidate genes involved in regulatory mechanisms have also been tested during embryogenesis using TSA-FISH. Here we show that the de novo MSL complex assembled on the X chromosomes in females further reduced the global expression level on the basis of 2/3 down-regulation caused by the inverse dosage effect in trisomy through epigenetic modulations rather than induced dosage compensation. Plus, the sexual dimorphism effect in unbalanced genomes was further examined due to the pre-existing of the MSL complex in males. All these results demonstrate the dynamic functions of the MSL complex on global gene expression in different aneuploid genomes.

Skewed allelic expression on X chromosome associated with aberrant expression of XIST on systemic lupus erythematosus lymphocytes

A common feature of autoimmune diseases, including systemic lupus erythematosus (SLE), is an increased prevalence in women. However, the molecular basis for sex disparity in SLE remains poorly understood. To examine the role of X-linked transcription in SLE adaptive immune cells, we performed RNA-seq in T cell and B cell subsets from either healthy donors or patients with SLE. Analyses of allelic expression (AE) profiles identified a pattern of increased allelic imbalance across the entire X chromosome in SLE lymphocytes. X-linked genes exhibiting AE in SLE had an extensive overlap with genes known to escape X chromosome inactivation (XCI). XIST RNA was overexpressed in SLE patients. Differential XIST expression correlated with AE profiles more positively at X-linked genes than the genome-wide background. Analysis of three independent RNA-seq data verified the XIST-associated skewed AE on X chromosome in SLE. Integrative analyses of DNA methylation profiles showed an increased variability of DNA methylation levels at these AE-related X-linked genes. In cultured lymphoblastic cells, knockdown of XIST specifically altered allelic imbalance patterns between X chromosomes. Our study provides genetic evidence that upregulation of XIST accompanied with more skewed allelic expression on X chromosome is associated with the pathogenesis of SLE and may provide mechanistic insights into the increased incidence of SLE in females.

Molecular Mechanisms of Skewed X-Chromosome Inactivation in Female Hemophilia Patients-Lessons from Wide Genome Analyses

Introduction: Hemophilia A (HA) is an X-linked bleeding disorder caused by factor VIII (FVIII) deficiency or dysfunction due to F8 gene mutations. HA carriers are usually asymptomatic because their FVIII levels correspond to approximately half of the concentration found in healthy individuals. However, in rare cases, a carrier may exhibit symptoms of moderate to severe HA primarily due to skewed inactivation of her non-hemophilic X chromosome. Aim: The aim of the study was to investigate X-chromosome inactivation (XCI) patterns in HA carriers, with special emphasis on three karyotypically normal HA carriers presenting with moderate to severe HA phenotype due to skewed XCI, in an attempt to elucidate the molecular mechanism underlying skewed XCI in these symptomatic HA carriers. The study was based on the hypothesis that the presence of a pathogenic mutation on the non-hemophilic X chromosome is the cause of extreme inactivation of that X chromosome. Methods: XCI patterns were studied by PCR analysis of the CAG repeat region in the HUMARA gene. HA carriers that demonstrated skewed XCI were further studied by whole-exome sequencing (WES) followed by X chromosome-targeted bioinformatic analysis. Results: All three HA carriers presenting with the moderate to severe HA phenotype due to skewed XCI were found to carry pathogenic mutations on their non-hemophilic X chromosomes. Patient 1 was diagnosed with a frameshift mutation in the PGK1 gene that was associated with familial XCI skewing in three generations. Patient 2 was diagnosed with a missense mutation in the SYTL4 gene that was associated with familial XCI skewing in two generations. Patient 3 was diagnosed with a nonsense mutation in the NKAP gene that was associated with familial XCI skewing in two generations. Conclusion: Our results indicate that the main reason for skewed XCI in our female HA patients was negative selection against cells with a disadvantage caused by an additional deleterious mutation on the silenced X chromosome, thus complicating the phenotype of a monogenic X-linked disease. Based on our study, we are currently offering the X inactivation test to symptomatic hemophilia carriers and plan to expand this approach to symptomatic carriers of other X-linked diseases, which can be further used in pregnancy planning.

Leukocytes with chromosome Y loss have reduced abundance of the cell surface immunoprotein CD99

Mosaic loss of chromosome Y (LOY) in immune cells is a male-specific mutation associated with increased risk for morbidity and mortality. The CD99 gene, positioned in the pseudoautosomal regions of chromosomes X and Y, encodes a cell surface protein essential for several key properties of leukocytes and immune system functions. Here we used CITE-seq for simultaneous quantification of CD99 derived mRNA and cell surface CD99 protein abundance in relation to LOY in single cells. The abundance of CD99 molecules was lower on the surfaces of LOY cells compared with cells without this aneuploidy in all six types of leukocytes studied, while the abundance of CD proteins encoded by genes located on autosomal chromosomes were independent from LOY. These results connect LOY in single cells with immune related cellular properties at the protein level, providing mechanistic insight regarding disease vulnerability in men affected with mosaic chromosome Y loss in blood leukocytes.

A Novel cis Regulatory Element Regulates Human XIST in a CTCF-Dependent Manner

The long noncoding RNA XIST is the master regulator for the process of X chromosome inactivation (XCI) in mammalian females. Here, we report the existence of a hitherto-uncharacterized cis regulatory element (cRE) within the first exon of human XIST, which determines the transcriptional status of XIST during the initiation and maintenance phases of XCI. In the initiation phase, pluripotency factors bind to this cRE and keep XIST repressed. In the maintenance phase of XCI, the cRE is enriched for CTCF, which activates XIST transcription. By employing a CRISPR-dCas9-KRAB-based interference strategy, we demonstrate that binding of CTCF to the newly identified cRE is critical for regulating XIST in a YY1-dependent manner. Collectively, our study uncovers the combinatorial effect of multiple transcriptional regulators influencing XIST expression during the initiation and maintenance phases of XCI.

Sex Differences Are Here to Stay: Relevance to Prenatal Care

Sex differences exist in the incidence and presentation of many pregnancy complications, including but not limited to pregnancy loss, spontaneous preterm birth, and fetal growth restriction. Sex differences arise very early in development due to differential gene expression from the X and Y chromosomes, and later may also be influenced by the action of gonadal steroid hormones. Though offspring sex is not considered in most prenatal diagnostic or therapeutic strategies currently in use, it may be beneficial to consider sex differences and the associated mechanisms underlying pregnancy complications. This review will cover (i) the prevalence and presentation of sex differences that occur in perinatal complications, particularly with a focus on the placenta; (ii) possible mechanisms underlying the development of sex differences in placental function and pregnancy phenotypes; and (iii) knowledge gaps that should be addressed in the development of diagnostic or risk prediction tools for such complications, with an emphasis on those for which it would be important to consider sex.