Genomation: a toolkit to summarize, annotate and visualize genomic intervals

Targeted bisulfite sequencing of Scots pine adaptation-related genes

During environmental changes, epigenetic processes can enable adaptive responses faster than natural selection. In plants, very little is known about the role of DNA methylation during long-term adaptation. Scots pine is a widely distributed coniferous species which must adapt to different environmental conditions throughout its long lifespan. Thus, epigenetic modifications may contribute towards this direction. We provide bisulfite next-generation sequencing data from the putative promoters and exons of eight adaptation-related genes (A3IP2, CCA1, COL1, COL2, FTL2, MFT1, PHYO, and ZTL) in three Scots pine populations located in northern and southern parts of Finland. DNA methylation levels were studied in the two seed tissues: the maternal megagametophyte which contributes to embryo viability, and the biparental embryo which represents the next generation. In most genes, differentially methylated cytosines (DMCs) were in line with our previously demonstrated gene expression differences found in the same Scots pine populations. In addition, we found a strong correlation of total methylation levels between the embryo and megagametophyte tissues of a given individual tree, which indicates that DNA methylation can be inherited from the maternal parent. In conclusion, our results imply that DNA methylation differences may contribute to the adaptation of Scots pine populations in different climatic conditions.

Whole-genome DNA methylation profiling revealed epigenetic regulation of NF-κB signaling pathway involved in response to Vibrio alginolyticus infection in the Pacific oyster, Crassostrea gigas

DNA methylation, an essential epigenetic alteration, is tightly linked to a variety of biological processes, such as immune response. To identify the epigenetic regulatory mechanism in Pacific oyster (Crassostrea gigas), whole-genome bisulfite sequencing (WGBS) was conducted on C. gigas at 0 h, 6 h, and 48 h after infection with Vibrio alginolyticus. At 6 h and 48 h, a total of 11,502 and 14,196 differentially methylated regions (DMRs) were identified (p<0.05, FDR<0.001) compared to 0 h, respectively. Gene ontology (GO) analysis showed that differentially methylated genes (DMGs) were significantly enriched in various biological pathways including immunity, cytoskeleton, epigenetic modification, and metabolic processes. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that transcription machinery (ko03021) is one of the most important pathways. Integrated transcriptome and methylome analyses allowed the identification of 167 and 379 DMG-related DEGs at 6 h and 48 h, respectively. These genes were significantly enriched in immune-related pathways, including nuclear factor kappa B (NF-κB) signaling pathway (ko04064) and tumor necrosis factor (TNF) signaling pathway (ko04668). Interestingly, it's observed that the NF-κB pathway could be activated jointly by TNF Receptor Associated Factor 2 (TRAF2) and Baculoviral IAP Repeat Containing 3 (BIRC3, the homolog of human BIRC2) which were regulated by DNA methylation in response to the challenge posed by V. alginolyticus infection. Through this study, we provided insightful information about the epigenetic regulation of immunity-related genes in the C. gigas, which will be valuable for the understanding of the innate immune system modulation and defense mechanism against bacterial infection in invertebrates.

Epigenetic disruptions in the offspring hypothalamus in response to maternal infection

DNA methylation is an epigenetic modification that can alter gene expression, and the incidence can vary across developmental stages, inflammatory conditions, and sexes. The effects of viral maternal viral infection and sex on the DNA methylation patterns were studied in the hypothalamus of a pig model of immune activation during development. DNA methylation at single-base resolution in regions of high CpG density was measured on 24 individual hypothalamus samples using reduced representation bisulfite sequencing. Differential over- and under-methylated sites were identified and annotated to proximal genes and corresponding biological processes. A total of 120 sites were differentially methylated (FDR-adjusted p-value < 0.05) between maternal infection or sex groups. Among the 66 sites differentially methylated between groups exposed to inflammatory signals and control, most sites were over-methylated in the challenged group and included sites in the promoter regions of genes SIRT3 and NRBP1. Among the 54 differentially methylated sites between females and males, most sites were over-methylated in females and included sites in the promoter region of genes TNC and EIF4G1. The analysis of the genes proximal to the differentially methylated sites suggested that biological processes potentially impacted include immune response, neuron migration and ensheathment, peptide signaling, adaptive thermogenesis, and tissue development. These results suggest that translational studies should consider that the prolonged effect of maternal infection during gestation may be enacted through epigenetic regulatory mechanisms that may differ between sexes.

Integrative analysis of the methylome and transcriptome of tomato fruit (Solanum lycopersicum L.) induced by postharvest handling

Tomato fruit ripening is triggered by the demethylation of key genes, which alters their transcriptional levels thereby initiating and propagating a cascade of physiological events. What is unknown is how these processes are altered when fruit are ripened using postharvest practices to extend shelf-life, as these practices often reduce fruit quality. To address this, postharvest handling-induced changes in the fruit DNA methylome and transcriptome, and how they correlate with ripening speed, and ripening indicators such as ethylene, abscisic acid, and carotenoids, were assessed. This study comprehensively connected changes in physiological events with dynamic molecular changes. Ripening fruit that reached 'Turning' (T) after dark storage at 20°C, 12.5°C, or 5°C chilling (followed by 20°C rewarming) were compared to fresh-harvest fruit 'FHT'. Fruit stored at 12.5°C had the biggest epigenetic marks and alterations in gene expression, exceeding changes induced by postharvest chilling. Fruit physiological and chronological age were uncoupled at 12.5°C, as the time-to-ripening was the longest. Fruit ripening to Turning at 12.5°C was not climacteric; there was no respiratory or ethylene burst, rather, fruit were high in abscisic acid. Clear differentiation between postharvest-ripened and 'FHT' was evident in the methylome and transcriptome. Higher expression of photosynthetic genes and chlorophyll levels in 'FHT' fruit pointed to light as influencing the molecular changes in fruit ripening. Finally, correlative analyses of the -omics data putatively identified genes regulated by DNA methylation. Collectively, these data improve our interpretation of how tomato fruit ripening patterns are altered by postharvest practices, and long-term are expected to help improve fruit quality.

Oncogene EVI1 drives acute myeloid leukemia via a targetable interaction with CTBP2

Acute myeloid leukemia (AML) driven by the activation of EVI1 due to chromosome 3q26/MECOM rearrangements is incurable. Because transcription factors such as EVI1 are notoriously hard to target, insight into the mechanism by which EVI1 drives myeloid transformation could provide alternative avenues for therapy. Applying protein folding predictions combined with proteomics technologies, we demonstrate that interaction of EVI1 with CTBP1 and CTBP2 via a single PLDLS motif is indispensable for leukemic transformation. A 4× PLDLS repeat construct outcompetes binding of EVI1 to CTBP1 and CTBP2 and inhibits proliferation of 3q26/MECOM rearranged AML in vitro and in xenotransplant models. This proof-of-concept study opens the possibility to target one of the most incurable forms of AML with specific EVI1-CTBP inhibitors. This has important implications for other tumor types with aberrant expression of EVI1 and for cancers transformed by different CTBP-dependent oncogenic transcription factors.

Chromatin context-dependent regulation and epigenetic manipulation of prime editing

We set out to exhaustively characterize the impact of the cis-chromatin environment on prime editing, a precise genome engineering tool. Using a highly sensitive method for mapping the genomic locations of randomly integrated reporters, we discover massive position effects, exemplified by editing efficiencies ranging from ∼0% to 94% for an identical target site and edit. Position effects on prime editing efficiency are well predicted by chromatin marks, e.g., positively by H3K79me2 and negatively by H3K9me3. Next, we developed a multiplex perturbational framework to assess the interaction of trans-acting factors with the cis-chromatin environment on editing outcomes. Applying this framework to DNA repair factors, we identify HLTF as a context-dependent repressor of prime editing. Finally, several lines of evidence suggest that active transcriptional elongation enhances prime editing. Consistent with this, we show we can robustly decrease or increase the efficiency of prime editing by preceding it with CRISPR-mediated silencing or activation, respectively.

Comprehensive Transcriptomic Analysis of EWSR1::WT1 Targets Identifies CDK4/6 Inhibitors as an Effective Therapy for Desmoplastic Small Round Cell Tumors

Desmoplastic small round cell tumors (DSRCT) are a type of aggressive, pediatric sarcoma characterized by the EWSR1::WT1 fusion oncogene. Targeted therapies for DSRCT have not been developed, and standard multimodal therapy is insufficient, leading to a 5-year survival rate of only 15% to 25%. Here, we depleted EWSR1::WT1 in DSRCT and established its essentiality in vivo. Transcriptomic analysis revealed that EWSR1::WT1 induces unique transcriptional alterations compared with WT1 and other fusion oncoproteins and that EWSR1::WT1 binding directly mediates gene upregulation. The E-KTS isoform of EWSR1::WT1 played a dominant role in transcription, and it bound to the CCND1 promoter and stimulated DSRCT growth through the cyclin D-CDK4/6-RB axis. Treatment with the CDK4/6 inhibitor palbociclib successfully reduced growth in two DSRCT xenograft models. As palbociclib has been approved by the FDA for the treatment of breast cancer, these findings demonstrate the sensitivity of DSRCT to palbociclib and support immediate clinical investigation of palbociclib for treating this aggressive pediatric cancer.

SIGNIFICANCE

EWSR1::WT1 is essential for desmoplastic small round cell tumors and upregulates the cyclin D-CDK4/6-RB axis that can be targeted with palbociclib, providing a targeted therapeutic strategy for treating this deadly tumor type.

Epigenetic Modifications of White Blood Cell DNA Caused by Transient Fetal Infection with Bovine Viral Diarrhea Virus

Bovine viral diarrhea virus (BVDV) infections cause USD 1.5-2 billion in losses annually. Maternal BVDV after 150 days of gestation causes transient fetal infection (TI) in which the fetal immune response clears the virus. The impact of fetal TI BVDV infections on postnatal growth and white blood cell (WBC) methylome as an index of epigenetic modifications was examined by inoculating pregnant heifers with noncytopathic type 2 BVDV or media (sham-inoculated controls) on Day 175 of gestation to generate TI (n = 11) and control heifer calves (n = 12). Fetal infection in TI calves was confirmed by virus-neutralizing antibody titers at birth and control calves were seronegative. Both control and TI calves were negative for BVDV RNA in WBCs by RT-PCR. The mean weight of the TI calves was less than that of the controls (p < 0.05). DNA methyl seq analysis of WBC DNA demonstrated 2349 differentially methylated cytosines (p ≤ 0.05) including 1277 hypomethylated cytosines, 1072 hypermethylated cytosines, 84 differentially methylated regions based on CpGs in promoters, and 89 DMRs in islands of TI WBC DNA compared to controls. Fetal BVDV infection during late gestation resulted in epigenomic modifications predicted to affect fetal development and immune pathways, suggesting potential consequences for postnatal growth and health of TI cattle.

Life-stage specificity and temporal variations in transcriptomes and DNA methylomes of the reef coral Pocillopora damicornis in response to thermal acclimation

Understanding the acclimation capacity of reef corals across generations to thermal stress and its underlying molecular underpinnings could provide insights into their resilience and adaptive responses to future climate change. Here, we acclimated adult brooding coral Pocillopora damicornis to high temperature (32 °C vs. 29 °C) for three weeks and analyzed the changes in phenotypes, transcriptomes and DNA methylomes of adult corals and their brooded larvae. Results showed that although adult corals did not show noticeable bleaching after thermal exposure, they released fewer but larger larvae. Interestingly, larval cohorts from two consecutive lunar days exhibited contrasting physiological resistance to thermal stress, as evidenced by the divergent responses of area-normalized symbiont densities and photochemical efficiency to thermal stress. RNA-seq and whole-genome bisulfite sequencing revealed that adult and larval corals mounted distinct transcriptional and DNA methylation changes in response to thermal stress. Remarkably, larval transcriptomes and DNA methylomes also varied greatly among lunar days and thermal treatments, aligning well with their physiological metrics. Overall, our study shows that changes in transcriptomes and DNA methylomes in response to thermal acclimation can be highly life stage-specific. More importantly, thermally-acclimated adult corals could produce larval offspring with temporally contrasting photochemical performance and thermal resilience, and such variations in larval phenotypes are associated with differential transcriptomes and DNA methylomes, and are likely to increase the likelihood of reproductive success and plasticity of larval propagules under thermal stress.

Neuronal-specific methylome and hydroxymethylome analysis reveal significant loci associated with alcohol use disorder

Background: Alcohol use disorder (AUD) is a complex condition associated with adverse health consequences that affect millions of individuals worldwide. Epigenetic modifications, including DNA methylation (5 mC), have been associated with AUD and other alcohol-related traits. Epigenome-wide association studies (EWAS) have identified differentially methylated genes associated with AUD in human peripheral and brain tissue. More recently, epigenetic studies of AUD have also evaluated DNA hydroxymethylation (5 hmC) in the human brain. However, most of the epigenetic work in postmortem brain tissue has examined bulk tissue. In this study, we investigated neuronal-specific 5 mC and 5 hmC alterations at CpG sites associated with AUD in the human orbitofrontal cortex (OFC). Methods: Neuronal nuclei from the OFC were evaluated in 34 human postmortem brain samples (10 AUD, 24 non-AUD). Reduced representation oxidative bisulfite sequencing was used to assess 5 mC and 5 hmC at the genome-wide level. Differential 5 mC and 5 hmC were evaluated using the methylKit R package and significance was set at false discovery rate < 0.05 and differential methylation > 2. Functional enrichment analyses were performed, and gene-level convergence was evaluated in an independent dataset that assessed 5 mC and 5 hmC of AUD in bulk cortical tissue. Results: We identified 417 5 mC and 363 5hmC significant differential CpG sites associated with AUD, with 59% in gene promoters. Some of the identified genes have been previously implicated in alcohol consumption, including SYK, DNMT3A for 5 mC, GAD1, DLX1, DLX2, for 5 hmC and GATA4 in both. Convergence with a previous AUD 5 mC and 5 hmC study was observed for 28 genes. We also identified 5 and 35 differential regions for 5 mC and 5 hmC, respectively. Lastly, GWAS enrichment analysis showed an association with AUD for differential 5 mC genes. Discussion: This study reveals neuronal-specific methylome and hydroxymethylome dysregulation associated with AUD, identifying both previously reported and potentially novel gene associations with AUD. Our findings provide new insights into the epigenomic dysregulation of AUD in the human brain.

Identification of epigenetically regulated genes involved in plant-virus interaction and their role in virus-triggered induced resistance

BACKGROUND

Plant responses to a wide range of stresses are known to be regulated by epigenetic mechanisms. Pathogen-related investigations, particularly against RNA viruses, are however scarce. It has been demonstrated that Arabidopsis thaliana plants defective in some members of the RNA-directed DNA methylation (RdDM) or histone modification pathways presented differential susceptibility to the turnip mosaic virus. In order to identify genes directly targeted by the RdDM-related RNA Polymerase V (POLV) complex and the histone demethylase protein JUMONJI14 (JMJ14) during infection, the transcriptomes of infected mutant and control plants were obtained and integrated with available chromatin occupancy data for various epigenetic proteins and marks.

RESULTS

A comprehensive list of virus-responsive gene candidates to be regulated by the two proteins was obtained. Twelve genes were selected for further characterization, confirming their dynamic regulation during the course of infection. Several epigenetic marks on their promoter sequences were found using in silico data, raising confidence that the identified genes are actually regulated by epigenetic mechanisms. The altered expression of six of these genes in mutants of the methyltransferase gene CURLY LEAF and the histone deacetylase gene HISTONE DEACETYLASE 19 suggests that some virus-responsive genes may be regulated by multiple coordinated epigenetic complexes. A temporally separated multiple plant virus infection experiment in which plants were transiently infected with one virus and then infected by a second one was designed to investigate the possible roles of the identified POLV- and JMJ14-regulated genes in wild-type (WT) plants. Plants that had previously been stimulated with viruses were found to be more resistant to subsequent virus challenge than control plants. Several POLV- and JMJ14-regulated genes were found to be regulated in virus induced resistance in WT plants, with some of them poisoned to be expressed in early infection stages.

CONCLUSIONS

A set of confident candidate genes directly regulated by the POLV and JMJ14 proteins during virus infection was identified, with indications that some of them may be regulated by multiple epigenetic modules. A subset of these genes may also play a role in the tolerance of WT plants to repeated, intermittent virus infections.

Genomic and Epigenomic Changes in the Progeny of Cold-Stressed Arabidopsis thaliana Plants

Plants are continuously exposed to various environmental stresses. Because they can not escape stress, they have to develop mechanisms of remembering stress exposures somatically and passing it to the progeny. We studied the Arabidopsis thaliana ecotype Columbia plants exposed to cold stress for 25 continuous generations. Our study revealed that multigenerational exposure to cold stress resulted in the changes in the genome and epigenome (DNA methylation) across generations. Main changes in the progeny were due to the high frequency of genetic mutations rather than epigenetic changes; the difference was primarily in single nucleotide substitutions and deletions. The progeny of cold-stressed plants exhibited the higher rate of missense non-synonymous mutations as compared to the progeny of control plants. At the same time, epigenetic changes were more common in the CHG (C = cytosine, H = cytosine, adenine or thymine, G = guanine) and CHH contexts and favored hypomethylation. There was an increase in the frequency of C to T (thymine) transitions at the CHH positions in the progeny of cold stressed plants; because this type of mutations is often due to the deamination of the methylated cytosines, it can be hypothesized that environment-induced changes in methylation contribute to mutagenesis and may be to microevolution processes and that RNA-dependent DNA methylation plays a crucial role. Our work supports the existence of heritable stress response in plants and demonstrates that genetic changes prevail.

Interplay between coding and non-coding regulation drives the Arabidopsis seed-to-seedling transition

Translation of seed stored mRNAs is essential to trigger germination. However, when RNAPII re-engages RNA synthesis during the seed-to-seedling transition has remained in question. Combining csRNA-seq, ATAC-seq and smFISH in Arabidopsis thaliana we demonstrate that active transcription initiation is detectable during the entire germination process. Features of non-coding regulation such as dynamic changes in chromatin accessible regions, antisense transcription, as well as bidirectional non-coding promoters are widespread throughout the Arabidopsis genome. We show that sensitivity to exogenous ABSCISIC ACID (ABA) during germination depends on proximal promoter accessibility at ABA-responsive genes. Moreover, we provide genetic validation of the existence of divergent transcription in plants. Our results reveal that active enhancer elements are transcribed producing non-coding enhancer RNAs (eRNAs) as widely documented in metazoans. In sum, this study defining the extent and role of coding and non-coding transcription during key stages of germination expands our understanding of transcriptional mechanisms underlying plant developmental transitions.

Cold exposure impacts DNA methylation patterns in cattle sperm

DNA methylation is influenced by various exogenous factors such as nutrition, temperature, toxicants, and stress. Bulls from the Pacific Northwest region of the United States and other northern areas are exposed to extreme cold temperatures during winter. However, the effects of cold exposure on the methylation patterns of bovine sperm remain unclear. To address, DNA methylation profiles of sperm collected during late spring and winter from the same bulls were analyzed using whole genome bisulfite sequencing (WGBS). Bismark (0.22.3) were used for mapping the WGBS reads and R Bioconductor package DSS was used for differential methylation analysis. Cold exposure induced 3,163 differentially methylated cytosines (DMCs) with methylation difference ≥10% and a q-value < 0.05. We identified 438 differentially methylated regions (DMRs) with q-value < 0.05, which overlapped with 186 unique genes. We also identified eight unique differentially methylated genes (DMGs) (Pax6, Macf1, Mest, Ubqln1, Smg9, Ctnnb1, Lsm4, and Peg10) involved in embryonic development, and nine unique DMGs (Prmt6, Nipal1, C21h15orf40, Slc37a3, Fam210a, Raly, Rgs3, Lmbr1, and Gan) involved in osteogenesis. Peg10 and Mest, two paternally expressed imprinted genes, exhibited >50% higher methylation. The differential methylation patterns of six distinct DMRs: Peg10, Smg9 and Mest related to embryonic development and Lmbr1, C21h15orf40 and Prtm6 related to osteogenesis, were assessed by methylation-specific PCR (MS-PCR), which confirmed the existence of variable methylation patterns in those locations across the two seasons. In summary, cold exposure induces differential DNA methylation patterns in genes that appear to affect embryonic development and osteogenesis in the offspring. Our findings suggest the importance of replicating the results of the current study with a larger sample size and exploring the potential of these changes in affecting offspring development.

AAV-mediated delivery of a Sleeping Beauty transposon and an mRNA-encoded transposase for the engineering of therapeutic immune cells

Engineering cells for adoptive therapy requires overcoming limitations in cell viability and, in the efficiency of transgene delivery, the duration of transgene expression and the stability of genomic integration. Here we report a gene-delivery system consisting of a Sleeping Beauty (SB) transposase encoded into a messenger RNA delivered by an adeno-associated virus (AAV) encoding an SB transposon that includes the desired transgene, for mediating the permanent integration of the transgene. Compared with lentiviral vectors and with the electroporation of plasmids of transposon DNA or minicircle DNA, the gene-delivery system, which we named MAJESTIC (for 'mRNA AAV-SB joint engineering of stable therapeutic immune cells'), offers prolonged transgene expression, as well as higher transgene expression, therapeutic-cell yield and cell viability. MAJESTIC can deliver chimeric antigen receptors (CARs) into T cells (which we show lead to strong anti-tumour activity in vivo) and also transduce natural killer cells, myeloid cells and induced pluripotent stem cells with bi-specific CARs, kill-switch CARs and synthetic T-cell receptors.

Bioframe: operations on genomic intervals in Pandas dataframes

MOTIVATION

Genomic intervals are one of the most prevalent data structures in computational genome biology, and used to represent features ranging from genes, to DNA binding sites, to disease variants. Operations on genomic intervals provide a language for asking questions about relationships between features. While there are excellent interval arithmetic tools for the command line, they are not smoothly integrated into Python, one of the most popular general-purpose computational and visualization environments.

RESULTS

Bioframe is a library to enable flexible and performant operations on genomic interval dataframes in Python. Bioframe extends the Python data science stack to use cases for computational genome biology by building directly on top of two of the most commonly-used Python libraries, NumPy and Pandas. The bioframe API enables flexible name and column orders, and decouples operations from data formats to avoid unnecessary conversions, a common scourge for bioinformaticians. Bioframe achieves these goals while maintaining high performance and a rich set of features.

AVAILABILITY AND IMPLEMENTATION

Bioframe is open-source under MIT license, cross-platform, and can be installed from the Python Package Index. The source code is maintained by Open2C on GitHub at https://github.com/open2c/bioframe.

Comparative interactome analysis of α-arrestin families in human and Drosophila

The α-arrestins form a large family of evolutionally conserved modulators that control diverse signaling pathways, including both G-protein-coupled receptor (GPCR)-mediated and non-GPCR-mediated pathways, across eukaryotes. However, unlike β-arrestins, only a few α-arrestin targets and functions have been characterized. Here, using affinity purification and mass spectrometry, we constructed interactomes for 6 human and 12 Drosophila α-arrestins. The resulting high-confidence interactomes comprised 307 and 467 prey proteins in human and Drosophila, respectively. A comparative analysis of these interactomes predicted not only conserved binding partners, such as motor proteins, proteases, ubiquitin ligases, RNA splicing factors, and GTPase-activating proteins, but also those specific to mammals, such as histone modifiers and the subunits of V-type ATPase. Given the manifestation of the interaction between the human α-arrestin, TXNIP, and the histone-modifying enzymes, including HDAC2, we undertook a global analysis of transcription signals and chromatin structures that were affected by TXNIP knockdown. We found that TXNIP activated targets by blocking HDAC2 recruitment to targets, a result that was validated by chromatin immunoprecipitation assays. Additionally, the interactome for an uncharacterized human α-arrestin ARRDC5 uncovered multiple components in the V-type ATPase, which plays a key role in bone resorption by osteoclasts. Our study presents conserved and species-specific protein-protein interaction maps for α-arrestins, which provide a valuable resource for interrogating their cellular functions for both basic and clinical research.

MetageneCluster: a Python package for filtering conflicting signal trends in metagene plots

BACKGROUND

Metagene plots provide a visualization of biological signal trends over subsections of the genome and are used to perform high-level analysis of experimental data by aggregating genome-level data to create an average profile. The generation of metagene plots is useful for summarizing the results of many sequencing-based applications. Despite their prevalence and utility, the standard metagene plot is blind to conflicting signals within data. If multiple distinct trends occur, they can interact destructively, creating a plot that does not accurately represent any of the underlying trends.

RESULTS

We present MetageneCluster, a Python tool to generate a collection of representative metagene plots based on k-means clustering of genomic regions of interest. Clustering the data by similarity allows us to identify patterns within the features of interest. We are then able to summarize each pattern present in the data, rather than averaging across the entire feature space. We show that our method performs well when used to identify conflicting signals in real-world genome-level data.

CONCLUSIONS

Overall, MetageneCluster is a user-friendly tool for the creation of metagene plots that capture distinct patterns in underlying sequence data.

XIST directly regulates X-linked and autosomal genes in naive human pluripotent cells

X chromosome inactivation (XCI) serves as a paradigm for RNA-mediated regulation of gene expression, wherein the long non-coding RNA XIST spreads across the X chromosome in cis to mediate gene silencing chromosome-wide. In female naive human pluripotent stem cells (hPSCs), XIST is in a dispersed configuration, and XCI does not occur, raising questions about XIST's function. We found that XIST spreads across the X chromosome and induces dampening of X-linked gene expression in naive hPSCs. Surprisingly, XIST also targets specific autosomal regions, where it induces repressive chromatin changes and gene expression dampening. Thereby, XIST equalizes X-linked gene dosage between male and female cells while inducing differences in autosomes. The dispersed Xist configuration and autosomal localization also occur transiently during XCI initiation in mouse PSCs. Together, our study identifies XIST as the regulator of X chromosome dampening, uncovers an evolutionarily conserved trans-acting role of XIST/Xist, and reveals a correlation between XIST/Xist dispersal and autosomal targeting.

Genomic and epigenomic insights into the mechanism of cold response in upland cotton (Gossypium hirsutum)

Functional genome research, including gene transcriptional and posttranslational modifications of histones, can benefit greatly from a high-quality genome assembly. Histone modification plays a significant role in modulating the responses to abiotic stress in plants. However, there are limited reports on the involvement of dynamic changes in histone modification in cold stress response in upland cotton. In this study, the genome of an elite accession, YM11, with considerable cold stress tolerance was de novo assembled, which yielded a genome of 2343.06 Mb with a contig N50 of 88.96 Mb, and a total of 73,821 protein-coding gene models were annotated. Comparisons among YM11 and five Gossypium allopolyploid cotton assemblies highlighted a large amount of structural variations and presence/absence variations. We analyzed transcriptome and metabolome changes in YM11 seedlings subjected to cold stress. Using the CUT&Tag method, genome-wide H3K4me3 and H3K9ac modification patterns and effect of histone changes on gene expression were profiled during cold stress. Significant and consistently changing histone modifications and the gene expressions were screened, of which transcription factors (TFs) were highlighted. Our results suggest a positive correlation between the changes in H3K4me3, H3K9ac modifications and cold stress-responsive gene activation. This genome assembly and comprehensive analysis of genome-wide histone modifications and gene expression provide insights into the genomic variation and epigenetic responses to cold stress in upland cotton.

Grandmaternal allergen sensitization reprograms epigenetic and airway responses to allergen in second-generation offspring

Asthma susceptibility is influenced by environmental, genetic, and epigenetic factors. DNA methylation is one form of epigenetic modification that regulates gene expression and is both inherited and modified by environmental exposures throughout life. Prenatal development is a particularly vulnerable time period during which exposure to maternal asthma increases asthma risk in offspring. How maternal asthma affects DNA methylation in offspring and what the consequences of differential methylation are in subsequent generations are not fully known. In this study, we tested the effects of grandmaternal house dust mite (HDM) allergen sensitization during pregnancy on airway physiology and inflammation in HDM-sensitized and challenged second-generation mice. We also tested the effects of grandmaternal HDM sensitization on tissue-specific DNA methylation in allergen-naïve and -sensitized second-generation mice. Descendants of both allergen- and vehicle-exposed grandmaternal founders exhibited airway hyperreactivity after HDM sensitization. However, grandmaternal allergen sensitization significantly potentiated airway hyperreactivity and altered the epigenomic trajectory in second-generation offspring after HDM sensitization compared with HDM-sensitized offspring from vehicle-exposed founders. As a result, biological processes and signaling pathways associated with epigenetic modifications were distinct between lineages. A targeted analysis of pathway-associated gene expression found that Smad3 was significantly dysregulated as a result of grandmaternal allergen sensitization. These data show that grandmaternal allergen exposure during pregnancy establishes a unique epigenetic trajectory that reprograms allergen responses in second-generation offspring and may contribute to asthma risk.NEW & NOTEWORTHY Asthma susceptibility is influenced by environmental, genetic, and epigenetic factors. This study shows that maternal allergen exposure during pregnancy promotes unique epigenetic trajectories in second-generation offspring at baseline and in response to allergen sensitization, which is associated with the potentiation of airway hyperreactivity. These effects are one mechanism by which maternal asthma may influence the inheritance of asthma risk.

Neuronal-specific methylome and hydroxymethylome analysis reveal replicated and novel loci associated with alcohol use disorder

Alcohol use disorder (AUD) is a complex condition associated with adverse health consequences that affect millions of individuals worldwide. Epigenetic modifications, including DNA methylation (5mC), have been associated with AUD and other alcohol-related traits. Epigenome-wide association studies (EWAS) have identified differentially methylated genes associated with AUD in human peripheral and brain tissue. More recently, epigenetic studies of AUD have also evaluated DNA hydroxymethylation (5hmC) in the human brain. However, most of the epigenetic work in postmortem brain tissue has examined bulk tissue. In this study, we investigated neuronal-specific 5mC and 5hmC alterations at CpG sites associated with AUD in the human orbitofrontal cortex (OFC). Neuronal nuclei from the OFC were evaluated in 34 human postmortem brain samples (10 AUD, 24 non-AUD). Reduced representation oxidative bisulfite sequencing was used to assess 5mC and 5hmC at the genome-wide level. Differential 5mC and 5hmC were evaluated using the methylKit R package and significance was set at false discovery rate <0.05 and differential methylation >2. Functional enrichment analyses were performed and replication was evaluated replication in an independent dataset that assessed 5mC and 5hmC of AUD in bulk cortical tissue. We identified 417 5mC and 363 5hmC genome-wide significant differential CpG sites associated with AUD, with 59% in gene promoters. We also identified genes previously implicated in alcohol consumption, such as SYK, CHRM2, DNMT3A, and GATA4, for 5mC and GATA4, and GAD1, GATA4, DLX1 for 5hmC. Replication was observed for 28 CpG sites from a previous AUD 5mC and 5hmC study, including FOXP1. Lastly, GWAS enrichment analysis showed an association with AUD for differential 5mC genes. This study reveals neuronal-specific methylome and hydroxymethylome dysregulation associated with AUD. We replicated previous findings and identified novel associations with AUD for both 5mC and 5hmC marks within the OFC. Our findings provide new insights into the epigenomic dysregulation of AUD in the human brain.

Patient Information Summarization in Clinical Settings: Scoping Review

BACKGROUND

Information overflow, a common problem in the present clinical environment, can be mitigated by summarizing clinical data. Although there are several solutions for clinical summarization, there is a lack of a complete overview of the research relevant to this field.

OBJECTIVE

This study aims to identify state-of-the-art solutions for clinical summarization, to analyze their capabilities, and to identify their properties.

METHODS

A scoping review of articles published between 2005 and 2022 was conducted. With a clinical focus, PubMed and Web of Science were queried to find an initial set of reports, later extended by articles found through a chain of citations. The included reports were analyzed to answer the questions of where, what, and how medical information is summarized; whether summarization conserves temporality, uncertainty, and medical pertinence; and how the propositions are evaluated and deployed. To answer how information is summarized, methods were compared through a new framework "collect-synthesize-communicate" referring to information gathering from data, its synthesis, and communication to the end user.

RESULTS

Overall, 128 articles were included, representing various medical fields. Exclusively structured data were used as input in 46.1% (59/128) of papers, text in 41.4% (53/128) of articles, and both in 10.2% (13/128) of papers. Using the proposed framework, 42.2% (54/128) of the records contributed to information collection, 27.3% (35/128) contributed to information synthesis, and 46.1% (59/128) presented solutions for summary communication. Numerous summarization approaches have been presented, including extractive (n=13) and abstractive summarization (n=19); topic modeling (n=5); summary specification (n=11); concept and relation extraction (n=30); visual design considerations (n=59); and complete pipelines (n=7) using information extraction, synthesis, and communication. Graphical displays (n=53), short texts (n=41), static reports (n=7), and problem-oriented views (n=7) were the most common types in terms of summary communication. Although temporality and uncertainty information were usually not conserved in most studies (74/128, 57.8% and 113/128, 88.3%, respectively), some studies presented solutions to treat this information. Overall, 115 (89.8%) articles showed results of an evaluation, and methods included evaluations with human participants (median 15, IQR 24 participants): measurements in experiments with human participants (n=31), real situations (n=8), and usability studies (n=28). Methods without human involvement included intrinsic evaluation (n=24), performance on a proxy (n=10), or domain-specific tasks (n=11). Overall, 11 (8.6%) reports described a system deployed in clinical settings.

CONCLUSIONS

The scientific literature contains many propositions for summarizing patient information but reports very few comparisons of these proposals. This work proposes to compare these algorithms through how they conserve essential aspects of clinical information and through the "collect-synthesize-communicate" framework. We found that current propositions usually address these 3 steps only partially. Moreover, they conserve and use temporality, uncertainty, and pertinent medical aspects to varying extents, and solutions are often preliminary.

Model-based characterization of the equilibrium dynamics of transcription initiation and promoter-proximal pausing in human cells

In metazoans, both transcription initiation and the escape of RNA polymerase (RNAP) from promoter-proximal pausing are key rate-limiting steps in gene expression. These processes play out at physically proximal sites on the DNA template and appear to influence one another through steric interactions. Here, we examine the dynamics of these processes using a combination of statistical modeling, simulation, and analysis of real nascent RNA sequencing data. We develop a simple probabilistic model that jointly describes the kinetics of transcription initiation, pause-escape, and elongation, and the generation of nascent RNA sequencing read counts under steady-state conditions. We then extend this initial model to allow for variability across cells in promoter-proximal pause site locations and steric hindrance of transcription initiation from paused RNAPs. In an extensive series of simulations, we show that this model enables accurate estimation of initiation and pause-escape rates. Furthermore, we show by simulation and analysis of real data that pause-escape is often strongly rate-limiting and that steric hindrance can dramatically reduce initiation rates. Our modeling framework is applicable to a variety of inference problems, and our software for estimation and simulation is freely available.

Intrinsic epigenetic state of primary osteosarcoma drives metastasis

Osteosarcoma (OS) is the most common primary malignant bone tumor affecting the pediatric population with high potential to metastasize to distal sites, most commonly the lung. Insights into defining molecular features contributing to metastatic potential are lacking. We have mapped the active chromatin landscapes of OS tumors by integrating histone H3 lysine acetylated chromatin (H3K27ac) profiles (n=13), chromatin accessibility profiles (n=11) and gene expression (n=13) to understand the differences in their active chromatin profiles and its impact on molecular mechanisms driving the malignant phenotypes. Primary OS tumors from patients with metastasis (primary met) have a distinct active chromatin landscape compared to primary tumors from patients without metastatic disease (localized). The difference in chromatin activity shapes the transcriptional profile of OS. We identified novel candidate genes involved in OS pathogenesis and metastasis, including PPP1R1B, PREX1 and IGF2BP1, which exhibit increased chromatin activity in primary met along with higher transcript levels. Overall, differential chromatin activity in primary met occurs in proximity of genes regulating actin cytoskeleton organization, cellular adhesion, and extracellular matrix suggestive of their role in facilitating OS metastasis. Furthermore, chromatin profiling of tumors from metastatic lung lesions noted increases in chromatin activity in genes involved in cell migration and key intracellular signaling cascades, including the Wnt pathway. Thus, this data demonstrates that metastatic potential is intrinsically present in primary metastatic tumors and the cellular chromatin profiles further adapt to allow for successful dissemination, migration, and colonization at the distal metastatic site.

Developmental, cytogenetic and epigenetic consequences of removing complex proteins and adding melatonin during in vitro maturation of bovine oocytes

Background

In vitro maturation (IVM) of germinal vesicle intact oocytes prior to in vitro fertilization (IVF) is practiced widely in animals. In human assisted reproduction it is generally reserved for fertility preservation or where ovarian stimulation is contraindicated. Standard practice incorporates complex proteins (CP), in the form of serum and/or albumin, into IVM media to mimic the ovarian follicle environment. However, the undefined nature of CP, together with batch variation and ethical concerns regarding their origin, necessitate the development of more defined formulations. A known component of follicular fluid, melatonin, has multifaceted roles including that of a metabolic regulator and antioxidant. In certain circumstances it can enhance oocyte maturation. At this stage in development, the germinal-vesicle intact oocyte is prone to aneuploidy and epigenetic dysregulation.

Objectives

To determine the developmental, cytogenetic and epigenetic consequences of removing CP and including melatonin during bovine IVM.

Materials and methods

The study comprised a 2 x 2 factorial arrangement comparing (i) the inclusion or exclusion of CP, and (ii) the addition (100 nM) or omission of melatonin, during IVM. Cumulus-oocyte complexes (COCs) were retrieved from stimulated cycles. Following IVM and IVF, putative zygotes were cultured to Day 8 in standard media. RNAseq was performed on isolated cumulus cells, cytogenetic analyses (SNP-based algorithms) on isolated trophectoderm cells, and DNA methylation analysis (reduced representation bisulfite sequencing) on isolated cells of the inner-cell mass.

Results

Removal of CP during IVM led to modest reductions in blastocyst development, whilst added melatonin was beneficial in the presence but detrimental in the absence of CP. The composition of IVM media did not affect the nature or incidence of chromosomal abnormalities but cumulus-cell transcript expression indicated altered metabolism (primarily lipid) in COCs. These effects preceded the establishment of distinct metabolic and epigenetic signatures several days later in expanded and hatching blastocysts.

Conclusions

These findings highlight the importance of lipid, particularly sterol, metabolism by the COC during IVM. They lay the foundation for future studies that seek to develop chemically defined systems of IVM for the generation of transferrable embryos that are both cytogenetically and epigenetically normal.

Sleeping Beauty Transposon Insertions into Nucleolar DNA by an Engineered Transposase Localized in the Nucleolus

Transposons are nature's gene delivery vehicles that can be harnessed for experimental and therapeutic purposes. The Sleeping Beauty (SB) transposon shows efficient transposition and long-term transgene expression in human cells, and is currently under clinical development for gene therapy. SB transposition occurs into the human genome in a random manner, which carries a risk of potential genotoxic effects associated with transposon integration. Here, we evaluated an experimental strategy to manipulate SB's target site distribution by preferentially compartmentalizing the SB transposase to the nucleolus, which contains repetitive ribosomal RNA (rRNA) genes. We generated a fusion protein composed of the nucleolar protein nucleophosmin (B23) and the SB100X transposase, which was found to retain almost full transposition activity as compared to unfused transposase and to be predominantly localized to nucleoli in transfected human cells. Analysis of transposon integration sites generated by B23-SB100X revealed a significant enrichment into the p-arms of chromosomes containing nucleolus organizing regions (NORs), with preferential integration into the p13 and p11.2 cytobands directly neighboring the NORs. This bias in the integration pattern was accompanied by an enrichment of insertions into nucleolus-associated chromatin domains (NADs) at the periphery of nucleolar DNA and into lamina-associated domains (LADs). Finally, sub-nuclear targeting of the transposase resulted in preferential integration into chromosomal domains associated with the Upstream Binding Transcription Factor (UBTF) that plays a critical role in the transcription of 47S rDNA gene repeats of the NORs by RNA Pol I. Future modifications of this technology may allow the development of methods for specific gene insertion for precision genetic engineering.

Origin matters: Using a local reference genome improves measures in population genomics

Genome sequencing enables answering fundamental questions about the genetic basis of adaptation, population structure and epigenetic mechanisms. Yet, we usually need a suitable reference genome for mapping population-level resequencing data. In some model systems, multiple reference genomes are available, giving the challenging task of determining which reference genome best suits the data. Here, we compared the use of two different reference genomes for the three-spined stickleback (Gasterosteus aculeatus), one novel genome derived from a European gynogenetic individual and the published reference genome of a North American individual. Specifically, we investigated the impact of using a local reference versus one generated from a distinct lineage on several common population genomics analyses. Through mapping genome resequencing data of 60 sticklebacks from across Europe and North America, we demonstrate that genetic distance among samples and the reference genomes impacts downstream analyses. Using a local reference genome increased mapping efficiency and genotyping accuracy, effectively retaining more and better data. Despite comparable distributions of the metrics generated across the genome using SNP data (i.e. π, Tajima's D and FST ), window-based statistics using different references resulted in different outlier genes and enriched gene functions. A marker-based analysis of DNA methylation distributions had a comparably high overlap in outlier genes and functions, yet with distinct differences depending on the reference genome. Overall, our results highlight how using a local reference genome decreases reference bias to increase confidence in downstream analyses of the data. Such results have significant implications in all reference-genome-based population genomic analyses.

The NCOR-HDAC3 co-repressive complex modulates the leukemogenic potential of the transcription factor ERG

The ERG (ETS-related gene) transcription factor is linked to various types of cancer, including leukemia. However, the specific ERG domains and co-factors contributing to leukemogenesis are poorly understood. Drug targeting a transcription factor such as ERG is challenging. Our study reveals the critical role of a conserved amino acid, proline, at position 199, located at the 3' end of the PNT (pointed) domain, in ERG's ability to induce leukemia. P199 is necessary for ERG to promote self-renewal, prevent myeloid differentiation in hematopoietic progenitor cells, and initiate leukemia in mouse models. Here we show that P199 facilitates ERG's interaction with the NCoR-HDAC3 co-repressor complex. Inhibiting HDAC3 reduces the growth of ERG-dependent leukemic and prostate cancer cells, indicating that the interaction between ERG and the NCoR-HDAC3 co-repressor complex is crucial for its oncogenic activity. Thus, targeting this interaction may offer a potential therapeutic intervention.

Epigenetic Findings in Twins with Esophageal Atresia

Esophageal atresia (EA) is the most common malformation of the upper gastrointestinal tract. The estimated incidence of EA is 1 in 3500 births. EA is more frequently observed in boys and in twins. The exact cause of isolated EA remains unknown; a multifactorial etiology, including epigenetic gene expression modifications, is considered. The study included six pairs of twins (three pairs of monozygotic twins and three pairs of dizygotic twins) in which one child was born with EA as an isolated defect, while the other twin was healthy. DNA samples were obtained from the blood and esophageal tissue of the child with EA as well as from the blood of the healthy twin. The reduced representation bisulfite sequencing (RRBS) technique was employed for a whole-genome methylation analysis. The analyses focused on comparing the CpG island methylation profiles between patients with EA and their healthy siblings. Hypermethylation in the promoters of 219 genes and hypomethylation in the promoters of 78 genes were observed. A pathway enrichment analysis revealed the statistically significant differences in methylation profile of 10 hypermethylated genes in the Rho GTPase pathway, previously undescribed in the field of EA (ARHGAP36, ARHGAP4, ARHGAP6, ARHGEF6, ARHGEF9, FGD1, GDI1, MCF2, OCRL, and STARD8).

Germinal center output is sustained by HELLS-dependent DNA-methylation-maintenance in B cells

HELLS/LSH (Helicase, Lymphoid Specific) is a SNF2-like chromatin remodelling protein involved in DNA methylation. Its loss-of-function in humans causes humoral immunodeficiency, called ICF4 syndrome (Immunodeficiency, Centromeric Instability, Facial anomalies). Here we show by our newly generated B-cell-specific Hells conditional knockout mouse model that HELLS plays a pivotal role in T-dependent B-cell responses. HELLS deficiency induces accelerated decay of germinal center (GC) B cells and impairs the generation of high affinity memory B cells and circulating antibodies. Mutant GC B cells undergo dramatic DNA hypomethylation and massive de-repression of evolutionary recent retrotransposons, which surprisingly does not directly affect their survival. Instead, they prematurely upregulate either memory B cell markers or the transcription factor ATF4, which is driving an mTORC1-dependent metabolic program typical of plasma cells. Treatment of wild type mice with a DNMT1-specific inhibitor phenocopies the accelerated kinetics, thus pointing towards DNA-methylation maintenance by HELLS being a crucial mechanism to fine-tune the GC transcriptional program and enable long-lasting humoral immunity.

Epigenetic regulation of functional candidate genes for milk production traits in dairy sheep subjected to protein restriction in the prepubertal stage

BACKGROUND

As the prepubertal stage is a crucial point for the proper development of the mammary gland and milk production, this study aims to evaluate how protein restriction at this stage can affect methylation marks in milk somatic cells. Here, 28 Assaf ewes were subjected to 42.3% nutritional protein restriction (14 animals, NPR) or fed standard diets (14 animals, C) during the prepubertal stage. During the second lactation, the milk somatic cells of these ewes were sampled, and the extracted DNA was subjected to whole-genome bisulfite sequencing.

RESULTS

A total of 1154 differentially methylated regions (DMRs) were identified between the NPR and C groups. Indeed, the results of functional enrichment analyses of the genes harboring these DMRs suggested their relevant effects on the development of the mammary gland and lipid metabolism in sheep. The additional analysis of the correlations of the mean methylation levels within these DMRs with fat, protein, and dry extract percentages in the milk and milk somatic cell counts suggested associations between several DMRs and milk production traits. However, there were no phenotypic differences in these traits between the NPR and C groups.

CONCLUSION

In light of the above, the results obtained in the current study might suggest potential candidate genes for the regulation of milk production traits in the sheep mammary gland. Further studies focusing on elucidating the genetic mechanisms affected by the identified DMRs may help to better understand the biological mechanisms modified in the mammary gland of dairy sheep as a response to nutritional challenges and their potential effects on milk production.

Lineage Plasticity and Stemness Phenotypes in Prostate Cancer: Harnessing the Power of Integrated "Omics" Approaches to Explore Measurable Metrics

Prostate cancer (PCa), the most frequent and second most lethal cancer type in men in developed countries, is a highly heterogeneous disease. PCa heterogeneity, therapy resistance, stemness, and lethal progression have been attributed to lineage plasticity, which refers to the ability of neoplastic cells to undergo phenotypic changes under microenvironmental pressures by switching between developmental cell states. What remains to be elucidated is how to identify measurements of lineage plasticity, how to implement them to inform preclinical and clinical research, and, further, how to classify patients and inform therapeutic strategies in the clinic. Recent research has highlighted the crucial role of next-generation sequencing technologies in identifying potential biomarkers associated with lineage plasticity. Here, we review the genomic, transcriptomic, and epigenetic events that have been described in PCa and highlight those with significance for lineage plasticity. We further focus on their relevance in PCa research and their benefits in PCa patient classification. Finally, we explore ways in which bioinformatic analyses can be used to determine lineage plasticity based on large omics analyses and algorithms that can shed light on upstream and downstream events. Most importantly, an integrated multiomics approach may soon allow for the identification of a lineage plasticity signature, which would revolutionize the molecular classification of PCa patients.

Deficiency of the Polycomb Protein RYBP and TET Methylcytosine Oxidases Promotes Extensive CpG Island Hypermethylation and Malignant Transformation

Hypermethylation of CpG islands (CGI) is a common feature of cancer cells and predominantly affects Polycomb-associated genomic regions. Elucidating the underlying mechanisms leading to DNA hypermethylation in human cancer could help identify chemoprevention strategies. Here, we evaluated the role of Polycomb complexes and 5-methylcytosine (5mC) oxidases in protecting CGIs from DNA methylation and observed that four genes coding for components of Polycomb repressive complex 1 (PRC1) are downregulated in tumors. Inactivation of RYBP, a key activator of variant PRC1 complexes, in combination with all three 5mC oxidases (TET proteins) in nontumorigenic bronchial epithelial cells led to widespread hypermethylation of Polycomb-marked CGIs affecting almost 4,000 target genes, which closely resembled the DNA hypermethylation landscape observed in human squamous cell lung tumors. The RYBP- and TET-deficient cells showed methylation-associated aberrant regulation of cancer-relevant pathways, including defects in the Hippo tumor suppressor network. Notably, the quadruple knockout cells acquired a transformed phenotype, including anchorage-independent growth and formation of squamous cell carcinomas in mice. This work provides a mechanism promoting hypermethylation of CGIs and shows that such hypermethylation can lead to cell transformation. The breakdown of a two-pronged protection mechanism can be a route towards genome-wide hypermethylation of CGIs in tumors.

SIGNIFICANCE

Dysfunction of the Polycomb component RYBP in combination with loss of 5-methylcytosine oxidases promotes widespread hypermethylation of CpG islands in bronchial cells and induces tumorigenesis, resembling changes seen in human lung tumors.

Profiling neuronal methylome and hydroxymethylome of opioid use disorder in the human orbitofrontal cortex

Opioid use disorder (OUD) is influenced by genetic and environmental factors. While recent research suggests epigenetic disturbances in OUD, this is mostly limited to DNA methylation (5mC). DNA hydroxymethylation (5hmC) has been widely understudied. We conducted a multi-omics profiling of OUD in a male cohort, integrating neuronal-specific 5mC and 5hmC as well as gene expression profiles from human postmortem orbitofrontal cortex (OUD = 12; non-OUD = 26). Single locus methylomic analysis and co-methylation analysis showed a higher number of OUD-associated genes and gene networks for 5hmC compared to 5mC; these were enriched for GPCR, Wnt, neurogenesis, and opioid signaling. 5hmC marks also showed a higher correlation with gene expression patterns and enriched for GWAS of psychiatric traits. Drug interaction analysis revealed interactions with opioid-related drugs, some used as OUD treatments. Our multi-omics findings suggest an important role of 5hmC and reveal loci epigenetically dysregulated in OFC neurons of individuals with OUD.

Cancer lineage-specific regulation of YAP responsive elements revealed through large-scale functional epigenomic screens

YAP is a key transcriptional co-activator of TEADs, it regulates cell growth and is frequently activated in cancer. In Malignant Pleural Mesothelioma (MPM), YAP is activated by loss-of-function mutations in upstream components of the Hippo pathway, while, in Uveal Melanoma (UM), YAP is activated in a Hippo-independent manner. To date, it is unclear if and how the different oncogenic lesions activating YAP impact its oncogenic program, which is particularly relevant for designing selective anti-cancer therapies. Here we show that, despite YAP being essential in both MPM and UM, its interaction with TEAD is unexpectedly dispensable in UM, limiting the applicability of TEAD inhibitors in this cancer type. Systematic functional interrogation of YAP regulatory elements in both cancer types reveals convergent regulation of broad oncogenic drivers in both MPM and UM, but also strikingly selective programs. Our work reveals unanticipated lineage-specific features of the YAP regulatory network that provide important insights to guide the design of tailored therapeutic strategies to inhibit YAP signaling across different cancer types.

Histone exchange sensors reveal variant specific dynamics in mouse embryonic stem cells

Eviction of histones from nucleosomes and their exchange with newly synthesized or alternative variants is a central epigenetic determinant. Here, we define the genome-wide occupancy and exchange pattern of canonical and non-canonical histone variants in mouse embryonic stem cells by genetically encoded exchange sensors. While exchange of all measured variants scales with transcription, we describe variant-specific associations with transcription elongation and Polycomb binding. We found considerable exchange of H3.1 and H2B variants in heterochromatin and repeat elements, contrasting the occupancy and little exchange of H3.3 in these regions. This unexpected association between H3.3 occupancy and exchange of canonical variants is also evident in active promoters and enhancers, and further validated by reduced H3.1 dynamics following depletion of H3.3-specific chaperone, HIRA. Finally, analyzing transgenic mice harboring H3.1 or H3.3 sensors demonstrates the vast potential of this system for studying histone exchange and its impact on gene expression regulation in vivo.

Comprehensive DNA Methylation Analysis Indicates That Pancreatic Intraepithelial Neoplasia Lesions Are Acinar-Derived and Epigenetically Primed for Carcinogenesis

Pancreatic ductal adenocarcinoma (PDAC) is believed to arise from the accumulation of a series of somatic mutations and is also frequently associated with pancreatic intraepithelial neoplasia (PanIN) lesions. However, there is still debate as to whether the cell type-of-origin of PanINs and PDACs in humans is acinar or ductal. As cell type identity is maintained epigenetically, DNA methylation changes during pancreatic neoplasia can provide a compelling perspective to examine this question. Here, we performed laser-capture microdissection on surgically resected specimens from 18 patients to isolate, with high purity, DNA for whole-genome bisulfite sequencing from four relevant cell types: acini, nonneoplastic ducts, PanIN lesions, and PDAC lesions. Differentially methylated regions (DMR) were identified using two complementary analytical approaches: bsseq, which identifies any DMRs but is particularly useful for large block-like DMRs, and informME, which profiles the potential energy landscape across the genome and is particularly useful for identifying differential methylation entropy. Both global methylation profiles and block DMRs clearly implicated an acinar origin for PanINs. At the gene level, PanIN lesions exhibited an intermediate acinar-ductal phenotype resembling acinar-to-ductal metaplasia. In 97.6% of PanIN-specific DMRs, PanIN lesions had an intermediate methylation level between normal and PDAC, which suggests from an information theory perspective that PanIN lesions are epigenetically primed to progress to PDAC. Thus, epigenomic analysis complements histopathology to define molecular progression toward PDAC. The shared epigenetic lineage between PanIN and PDAC lesions could provide an opportunity for prevention by targeting aberrantly methylated progression-related genes.

SIGNIFICANCE

Analysis of DNA methylation landscapes provides insights into the cell-of-origin of PanIN lesions, clarifies the role of PanIN lesions as metaplastic precursors to human PDAC, and suggests potential targets for chemoprevention.

Molecular and phenotypic effects of early exposure to an environmentally relevant pesticide mixture in the Pacific oyster, Crassostrea gigas

Early life stages are crucial for organism development, especially for those displaying external fertilization, whose gametes and early stages face environmental stressors such as xenobiotics. The pacific oyster, Crassostrea gigas, is considered a model species in ecotoxicology because of its ecological characteristics (benthic, sessile, filter feeding). So far studies have investigated the impact of xenobiotics at embryotoxic, genotoxic and physiological endpoints, sometimes at the multigenerational scale, highlighting the role of epigenetic mechanisms in transmitting alterations induced by exposure to single xenobiotics. However, to date, little is known about the impact of environmentally-mimicking contaminants cocktails. Thus, we examined the impact of an early exposure to environmentally relevant mixture on the Pacific oyster life history. We studied transcriptomic, epigenetic and physiological alterations induced in oysters exposed to 18 pesticides and metals at environmental concentration (nominal sum concentration: 2.85 μg.L-1, measured sum concentration: 3.74 ± 0.013 μg.L-1) during embryo-larval stage (0-48 h post fertilization, hpf). No significant differences in embryo-larval abnormalities at 24 hpf were observed during larval and spat rearing; the swimming behaviour of exposed individuals was disturbed, while they were longer and heavier at specific time points, and exhibited a lower epinephrine-induced metamorphosis rate as well as a higher survival rate in the field. In addition, RNA-seq analyses of gastrula embryos revealed the differential expression of development-related genes (e.g. Hox orthologues and cell cycle regulators) between control and exposed oysters. Whole-genome DNA methylation analyses demonstrated a significant modification of DNA methylation in exposed larvae marked by a demethylation trend. Those findings suggest that early exposure to an environmentally relevant pesticide mixture induces multi-scale latent effects possibly affecting life history traits in the Pacific oyster.

Profiling the immune epigenome across global cattle breeds

BACKGROUND

Understanding the variation between well and poorly adapted cattle breeds to local environments and pathogens is essential for breeding cattle with improved climate and disease-resistant phenotypes. Although considerable progress has been made towards identifying genetic differences between breeds, variation at the epigenetic and chromatin levels remains poorly characterized. Here, we generate, sequence and analyse over 150 libraries at base-pair resolution to explore the dynamics of DNA methylation and chromatin accessibility of the bovine immune system across three distinct cattle lineages.

RESULTS

We find extensive epigenetic divergence between the taurine and indicine cattle breeds across immune cell types, which is linked to the levels of local DNA sequence divergence between the two cattle sub-species. The unique cell type profiles enable the deconvolution of complex cellular mixtures using digital cytometry approaches. Finally, we show distinct sub-categories of CpG islands based on their chromatin and methylation profiles that discriminate between classes of distal and gene proximal islands linked to discrete transcriptional states.

CONCLUSIONS

Our study provides a comprehensive resource of DNA methylation, chromatin accessibility and RNA expression profiles of three diverse cattle populations. The findings have important implications, from understanding how genetic editing across breeds, and consequently regulatory backgrounds, may have distinct impacts to designing effective cattle epigenome-wide association studies in non-European breeds.

Multi-omics data reveals novel impacts of human papillomavirus integration on the epigenomic and transcriptomic signatures of cervical tumorigenesis

Integration of human papilloma virus (HPV) DNA into the human genome may progressively contribute to cervical carcinogenesis. To explore how HPV integration affects gene expression by altering DNA methylation during carcinogenesis, we analyzed a multiomics dataset for cervical cancer. We obtained multiomics data by HPV-capture sequencing, RNA sequencing, and Whole Genome Bisulfite Sequencing from 50 patients with cervical cancer. We detected 985 and 485 HPV-integration sites in matched tumor and adjacent paratumor tissues. Of these, LINC00486 (n = 19), LINC02425 (n = 11), LLPH (n = 11), PROS1 (n = 5), KLF5 (n = 4), LINC00392 (n = 3), MIR205HG (n = 3) and NRG1 (n = 3) were identified as high-frequency HPV-integrated genes, including five novel recurrent genes. Patients at clinical stage II had the highest number of HPV integrations. E6 and E7 genes of HPV16 but not HPV18 showed significantly fewer breakpoints than random distribution. HPV integrations occurring in exons were associated with altered gene expression in tumor tissues but not in paratumor tissues. A list of HPV-integrated genes regulated at transcriptomic or epigenetic level was reported. We also carefully checked the candidate genes with regulation pattern correlated in both levels. HPV fragments integrated at MIR205HG mainly came from the L1 gene of HPV16. RNA expression of PROS1 was downregulated when HPV integrated in its upstream region. RNA expression of MIR205HG was elevated when HPV integrated into its enhancer. The promoter methylation levels of PROS1 and MIR205HG were all negatively correlated with their gene expressions. Further experimental validations proved that upregulation of MIR205HG could promote the proliferative and migrative abilities of cervical cancer cells. Our data provides a new atlas for epigenetic and transcriptomic regulations regarding HPV integrations in cervical cancer genome. We demonstrate that HPV integration may affect gene expression by altering methylation levels of MIR205HG and PROS1. Our study provides novel biological and clinical insights into HPV-induced cervical cancer.

Upregulation of extracellular proteins in a mouse model of Alzheimer's disease

Various risk factors of Alzheimer's disease (AD) are known, such as advanced age, possession of certain genetic variants, accumulation of toxic amyloid-β (Aβ) peptides, and unhealthy lifestyle. An estimate of heritability of AD ranges from 0.13 to 0.25, indicating that its phenotypic variation is accounted for mostly by non-genetic factors. DNA methylation is regarded as an epigenetic mechanism that interfaces the genome with non-genetic factors. The Tg2576 mouse model has been insightful in AD research. These transgenic mice express a mutant form of human amyloid precursor protein linked to familial AD. At 9-13 months of age, these mice show elevated levels of Aβ peptides and cognitive impairment. The current literature lacks integrative multiomics of the animal model. We applied transcriptomics and DNA methylomics to the same brain samples from ~ 11-month-old transgenic mice. We found that genes involved in extracellular matrix structures and functions are transcriptionally upregulated, and genes involved in extracellular protein secretion and localization are differentially methylated in the transgenic mice. Integrative analysis found enrichment of GO terms related to memory and synaptic functionability. Our results indicate a possibility of transcriptional modulation by DNA methylation underlying AD neuropathology.

Chromatin context-dependent regulation and epigenetic manipulation of prime editing

Prime editing is a powerful means of introducing precise changes to specific locations in mammalian genomes. However, the widely varying efficiency of prime editing across target sites of interest has limited its adoption in the context of both basic research and clinical settings. Here, we set out to exhaustively characterize the impact of the cis- chromatin environment on prime editing efficiency. Using a newly developed and highly sensitive method for mapping the genomic locations of a randomly integrated "sensor", we identify specific epigenetic features that strongly correlate with the highly variable efficiency of prime editing across different genomic locations. Next, to assess the interaction of trans -acting factors with the cis -chromatin environment, we develop and apply a pooled genetic screening approach with which the impact of knocking down various DNA repair factors on prime editing efficiency can be stratified by cis -chromatin context. Finally, we demonstrate that we can dramatically modulate the efficiency of prime editing through epigenome editing, i.e. altering chromatin state in a locus-specific manner in order to increase or decrease the efficiency of prime editing at a target site. Looking forward, we envision that the insights and tools described here will broaden the range of both basic research and therapeutic contexts in which prime editing is useful.

Intergenerational impact of dietary protein restriction in dairy ewes on epigenetic marks in the perirenal fat of their suckling lambs

In sheep, nutrition during the prepubertal stage is essential for growth performance and mammary gland development. However, the potential effects of nutrient restriction in a prepuberal stage over the progeny still need to be better understood. Here, the intergenerational effect of maternal protein restriction at prepubertal age (2 months of age) on methylation patterns was evaluated in the perirenal fat of Assaf suckling lambs. In total, 17 lambs from ewes subjected to dietary protein restriction (NPR group, 44% less protein) and 17 lambs from control ewes (C group) were analyzed. These lambs were ranked based on their carcass proportion of perirenal and cavitary fat and classified into HighPCF and LowPCF groups. The perirenal tissue from 4 NPR-LowPCF, 4 NPR-HighPCF, 4 C-LowPCF, and 4 C-HighPCF lambs was subjected to whole-genome bisulfite sequencing and differentially methylated regions (DMRs) were identified. Among other relevant processes, these DMRs were mapped in genes responsible for regulating the transition of brown to white adipose tissue and nonshivering thermoregulation, which might be associated with better adaptation/survival of lambs in the perinatal stage. The current study provides important biological insights about the intergenerational effect on the methylation pattern of an NPR in replacement ewes.

Therapeutic immune cell engineering with an mRNA : AAV- Sleeping Beauty composite system

Adoptive cell therapy has shown clinical success in patients with hematological malignancies. Immune cell engineering is critical for production, research, and development of cell therapy; however, current approaches for generation of therapeutic immune cells face various limitations. Here, we establish a composite gene delivery system for the highly efficient engineering of therapeutic immune cells. This system, termed MAJESTIC ( m RNA A AV-Sleeping-Beauty J oint E ngineering of S table T herapeutic I mmune C ells), combines the merits of mRNA, AAV vector, and transposon into one composite system. In MAJESTIC, the transient mRNA component encodes a transposase that mediates permanent genomic integration of the Sleeping Beauty (SB) transposon, which carries the gene-of-interest and is embedded within the AAV vector. This system can transduce diverse immune cell types with low cellular toxicity and achieve highly efficient and stable therapeutic cargo delivery. Compared with conventional gene delivery systems, such as lentiviral vector, DNA transposon plasmid, or minicircle electroporation, MAJESTIC shows higher cell viability, chimeric antigen receptor (CAR) transgene expression, therapeutic cell yield, as well as prolonged transgene expression. CAR-T cells generated by MAJESTIC are functional and have strong anti-tumor activity in vivo . This system also demonstrates versatility for engineering different cell therapy constructs such as canonical CAR, bi-specific CAR, kill switch CAR, and synthetic TCR; and for CAR delivery into various immune cells, including T cells, natural killer cells, myeloid cells, and induced pluripotent stem cells.

The model diatom Phaeodactylum tricornutum provides insights into the diversity and function of microeukaryotic DNA methyltransferases

Cytosine methylation is an important epigenetic mark involved in the transcriptional control of transposable elements in mammals, plants and fungi. The Stramenopiles-Alveolate-Rhizaria (SAR) lineages are a major group of ecologically important marine microeukaryotes, including the phytoplankton groups diatoms and dinoflagellates. However, little is known about their DNA methyltransferase diversity. Here, we performed an in-silico analysis of DNA methyltransferases found in marine microeukaryotes and showed that they encode divergent DNMT3, DNMT4, DNMT5 and DNMT6 enzymes. Furthermore, we found three classes of enzymes within the DNMT5 family. Using a CRISPR/Cas9 strategy we demonstrated that the loss of the DNMT5a gene correlates with a global depletion of DNA methylation and overexpression of young transposable elements in the model diatom Phaeodactylum tricornutum. The study provides a view of the structure and function of a DNMT family in the SAR supergroup using an attractive model species.

Pig fetal skeletal muscle development is associated with genome-wide DNA hypomethylation and corresponding alterations in transcript and microRNA expression

Fetal myogenesis represents a critical period of porcine skeletal muscle development and requires coordinated expression of thousands of genes. Epigenetic mechanisms, including DNA methylation, drive transcriptional regulation during development; however, these processes are understudied in developing porcine tissues. We performed bisulfite sequencing to assess DNA methylation in pig longissimus dorsi muscle at 41- and 70-days gestation (dg), as well as RNA- and small RNA-sequencing to identify coordinated changes in methylation and expression between myogenic stages. We identified 45 739 differentially methylated regions (DMRs) between stages, and the majority (N = 34 232) were hypomethylated at 70 versus 41 dg. Integration of methylation and transcriptomic data revealed strong associations between differential gene methylation and expression. Differential miRNA methylation was significantly negatively correlated with abundance, and dynamic expression of assayed miRNAs persisted postnatally. Motif analysis revealed significant enrichment of myogenic regulatory factor motifs among hypomethylated regions, suggesting that DNA hypomethylation may function to increase accessibility of muscle-specific transcription factors. We show that developmental DMRs are enriched for GWAS SNPs for muscle- and meat-related traits, demonstrating the potential for epigenetic processes to influence phenotypic diversity. Our results enhance understanding of DNA methylation dynamics of porcine myogenesis and reveal putative cis-regulatory elements governed by epigenetic processes.

Passer, a highly active transposon from a fish genome, as a potential new robust genetic manipulation tool

The discovery of new, active DNA transposons can expand the range of genetic tools and provide more options for genomic manipulation. In this study, a bioinformatics analysis suggested that Passer (PS) transposons, which are members of the pogo superfamily, show signs of recent and current activity in animals and may be active in some species. Cell-based transposition assays revealed that the native PS transposases from Gasterosteus aculeatus and Danio rerio displayed very high activity in human cells relative to the Sleeping Beauty transposon. A typical overproduction inhibition phenomenon was observed for PS, and transposition capacity was decreased by ∼12% with each kilobase increase in the insertion size. Furthermore, PS exhibited a pronounced integration preference for genes and their transcriptional regulatory regions. We further show that two domesticated human proteins derived from PS transposases have lost their transposition activity. Overall, PS may represent an alternative with a potentially efficient genetic manipulation tool for transgenesis and mutagenesis applications.

Epigenetic and transcriptional regulations prime cell fate before division during human pluripotent stem cell differentiation

Stem cells undergo cellular division during their differentiation to produce daughter cells with a new cellular identity. However, the epigenetic events and molecular mechanisms occurring between consecutive cell divisions have been insufficiently studied due to technical limitations. Here, using the FUCCI reporter we developed a cell-cycle synchronised human pluripotent stem cell (hPSC) differentiation system for uncovering epigenome and transcriptome dynamics during the first two divisions leading to definitive endoderm. We observed that transcription of key differentiation markers occurs before cell division, while chromatin accessibility analyses revealed the early inhibition of alternative cell fates. We found that Activator protein-1 members controlled by p38/MAPK signalling are necessary for inducing endoderm while blocking cell fate shifting toward mesoderm, and that enhancers are rapidly established and decommissioned between different cell divisions. Our study has practical biomedical utility for producing hPSC-derived patient-specific cell types since p38/MAPK induction increased the differentiation efficiency of insulin-producing pancreatic beta-cells.

A Cautionary Tale of Sexing by Methylation: Hybrid Bisulfite-Conversion Sequencing of Immunoprecipitated Methylated DNA in Chrysemys picta Turtles with Temperature-Dependent Sex Determination Reveals Contrasting Patterns of Somatic and Gonadal Methylation, but No Unobtrusive Sex Diagnostic

Background: The gonads of Chrysemys picta, a turtle with temperature-dependent sex determination (TSD), exhibit differential DNA methylation between males and females, but whether the same is true in somatic tissues remains unknown. Such differential DNA methylation in the soma would provide a non-lethal sex diagnostic for TSD turtle hatchings who lack visually detectable sexual dimorphism when young. Methods: Here, we tested multiple approaches to study DNA methylation in tail clips of Chrysemys picta hatchlings, to identify differentially methylated candidate regions/sites that could serve as molecular sex markers To detect global differential methylation in the tails we used methylation-sensitive ELISA, and to test for differential local methylation we developed a novel hybrid method by sequencing immunoprecipitated and bisulfite converted DNA (MeDIP-BS-seq) followed by PCR validation of candidate regions/sites after digestion with a methylation-sensitive restriction enzyme. Results: We detected no global differences in methylation between males and females via ELISA. While we detected inter-individual variation in DNA methylation in the tails, this variation was not sexually dimorphic, in contrast with hatchling gonads. Conclusions: Results highlight that differential DNA methylation is tissue-specific and plays a key role in gonadal formation (primary sexual development) and maintenance post-hatching, but not in the somatic tail tissue.