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

Genomic Characterization of Cisplatin Response Uncovers Priming of Cisplatin-Induced Genes in a Resistant Cell Line

Cisplatin is a chemotherapy drug that kills cancer cells by damaging their DNA. In human cells, this damage is repaired primarily by nucleotide excision repair. While cisplatin is generally effective, many cancers exhibit initial or acquired resistance to it. Here, we studied cisplatin resistance in a defined cell line system. We conducted a comprehensive genomic characterization of the cisplatin-sensitive A2780 ovarian cancer cell line compared to A2780cis, its resistant derivative. The resistant cells acquired less damage, but had similar repair kinetics. Genome-wide mapping of nucleotide excision repair showed a shift in the resistant cells from global genome towards transcription-coupled repair. By mapping gene expression changes following cisplatin treatment, we identified 56 upregulated genes that have higher basal expression in the resistant cell line, suggesting they are primed for a cisplatin response. More than half of these genes are novel to cisplatin- or damage-response. Six out of seven primed genes tested were upregulated in response to cisplatin in additional cell lines, making them attractive candidates for future investigation. These novel candidates for cisplatin resistance could prove to be important prognostic markers or targets for tailored combined therapy in the future.

Impact of Epigenomic Hypermethylation at TP53 on Allogeneic Hematopoietic Cell Transplantation Outcomes for Myelodysplastic Syndromes

Myelodysplastic syndromes (MDS) are a heterogeneous group of hematopoietic stem cell disorders for which allogeneic hematopoietic cell transplantation (HCT) is currently the sole curative treatment. Epigenetic lesions are considered a major pathogenetic determinant in many cancers, and in MDS, epigenetic-based interventions have emerged as life-prolonging therapies. However, the impact of epigenomic aberrations on HCT outcomes among patients with MDS are not well understood. We hypothesized that epigenomic signatures in MDS patients before undergoing HCT serve as a novel prognostic indicator of the risk of post-HCT MDS relapse. To evaluate these epigenomic signatures in MDS patients, we analyzed reduced representation bisulfite sequencing profiles in a matched case-control population of 94 patients. Among these HCT recipients, 47 patients with MDS who relapsed post-HCT (cases) were matched 1:1 to patients with MDS who did not relapse (controls). Only patients with wild-type TP53, RAS pathway, and JAK2 mutations were included in this study to promote the discovery of novel factors. Cases were matched with controls based on conditioning regimen intensity, age, sex, Revised International Prognostic Scoring System, Karnofsky Performance Status, graft type, and donor type. Pre-HCT whole-blood samples from cases and matched controls were obtained from the Center for International Blood and Marrow Transplant Research repository. We comprehensively identified differentially methylated regions (DMRs) by comparing the methylation patterns among matched cases and controls. Our findings show that cases displayed more hyper-DMRs pretransplantation compared with controls, even after adjusting for pre-HCT use of hypomethylating agents. Hyper-DMRs specific to cases were mapped to the transcription start site of 218 unique genes enriched in 5 different signaling pathways that may serve as regions of interest and factors to consider as prognostic determinants of post-HCT relapse in MDS patients. Interestingly, although patients selected for this cohort were wild-type for the TP53 gene, cases showed significantly greater levels of methylation at TP53 compared with controls. These findings indicate that previously identified prognostic genes for MDS, such as TP53, may affect disease relapse not only through genetic mutation, but also through epigenetic methylation mechanisms.

Intergenic RNA mainly derives from nascent transcripts of known genes

BACKGROUND

Eukaryotic genomes undergo pervasive transcription, leading to the production of many types of stable and unstable RNAs. Transcription is not restricted to regions with annotated gene features but includes almost any genomic context. Currently, the source and function of most RNAs originating from intergenic regions in the human genome remain unclear.

RESULTS

We hypothesize that many intergenic RNAs can be ascribed to the presence of as-yet unannotated genes or the "fuzzy" transcription of known genes that extends beyond the annotated boundaries. To elucidate the contributions of these two sources, we assemble a dataset of more than 2.5 billion publicly available RNA-seq reads across 5 human cell lines and multiple cellular compartments to annotate transcriptional units in the human genome. About 80% of transcripts from unannotated intergenic regions can be attributed to the fuzzy transcription of existing genes; the remaining transcripts originate mainly from putative long non-coding RNA loci that are rarely spliced. We validate the transcriptional activity of these intergenic RNAs using independent measurements, including transcriptional start sites, chromatin signatures, and genomic occupancies of RNA polymerase II in various phosphorylation states. We also analyze the nuclear localization and sensitivities of intergenic transcripts to nucleases to illustrate that they tend to be rapidly degraded either on-chromatin by XRN2 or off-chromatin by the exosome.

CONCLUSIONS

We provide a curated atlas of intergenic RNAs that distinguishes between alternative processing of well-annotated genes from independent transcriptional units based on the combined analysis of chromatin signatures, nuclear RNA localization, and degradation pathways.

Glutaminolysis-induced mTORC1 activation drives non-alcoholic steatohepatitis progression

BACKGROUND & AIMS

A holistic insight on the relationship between obesity and metabolic dysfunction-associated fatty liver disease is an unmet clinical need. Omics investigations can be used to investigate the multifaceted role of altered mitochondrial pathways to promote nonalcoholic steatohepatitis, a major risk factor for liver disease-associated death. There are no specific treatments but remission via surgery might offer an opportunity to examine the signaling processes that govern the complex spectrum of chronic liver diseases observed in extreme obesity. We aim to assess the emerging relationship between metabolism, methylation and liver disease.

METHODS

We tailed the flow of information, before and after steatohepatitis remission, from biochemical, histological, and multi-omics analyses in liver biopsies from patients with extreme obesity and successful bariatric surgery. Functional studies were performed in HepG2 cells and primary hepatocytes.

RESULTS

The reversal of hepatic mitochondrial dysfunction and the control of oxidative stress and inflammatory responses revealed the regulatory role of mitogen-activated protein kinases. The reversible metabolic rearrangements leading to steatohepatitis increased the glutaminolysis-induced production of α-ketoglutarate and the hyperactivation of mammalian target of rapamycin complex 1. These changes were crucial for the adenosine monophosphate-activated protein kinase/mammalian target of rapamycin-driven pathways that modulated hepatocyte survival by coordinating apoptosis and autophagy. The signaling activity of α-ketoglutarate and the associated metabolites also affected methylation-related epigenomic remodeling enzymes. Integrative analysis of hepatic transcriptome signatures and differentially methylated genomic regions distinguished patients with and without steatohepatitis.

CONCLUSION

We provide evidence supporting the multifaceted potential of the increased glutaminolysis-induced α-ketoglutarate production and the mammalian target of rapamycin complex 1 dysregulation as a conceivable source of the inefficient adaptive responses leading to steatohepatitis.

LAY SUMMARY

Steatohepatitis is a frequent and threatening complication of extreme obesity without specific treatment. Omics technologies can be used to identify therapeutic targets. We highlight increased glutaminolysis-induced α-ketoglutarate production as a potential source of signals promoting and exacerbating steatohepatitis.

Chromatin-associated MRN complex protects highly transcribing genes from genomic instability

MRN-MDC1 plays a central role in the DNA damage response (DDR) and repair. Using proteomics of isolated chromatin fragments, we identified DDR factors, such as MDC1, among those highly associating with a genomic locus upon transcriptional activation. Purification of MDC1 in the absence of exogenous DNA damage revealed interactions with factors involved in gene expression and RNA processing, in addition to DDR factors. ChIP-seq showed that MRN subunits, MRE11 and NBS1, colocalized throughout the genome, notably at TSSs and bodies of actively transcribing genes, which was dependent on the RNAPII transcriptional complex rather than transcription per se. Depletion of MRN increased RNAPII abundance at MRE11/NBS1-bound genes. Prolonged MRE11 or NBS1 depletion induced single-nucleotide polymorphisms across actively transcribing MRN target genes. These data suggest that association of MRN with the transcriptional machinery constitutively scans active genes for transcription-induced DNA damage to preserve the integrity of the coding genome.

Transcription factor chromatin profiling genome-wide using uliCUT&RUN in single cells and individual blastocysts

Determining chromatin-associated protein localization across the genome has provided insight into the functions of DNA-binding proteins and their connections to disease. However, established protocols requiring large quantities of cell or tissue samples currently limit applications for clinical and biomedical research in this field. Furthermore, most technologies have been optimized to assess abundant histone protein localization, prohibiting the investigation of nonhistone protein localization in low cell numbers. We recently described a protocol to profile chromatin-associated protein localization in as low as one cell: ultra-low-input cleavage under targets and release using nuclease (uliCUT&RUN). Optimized from chromatin immunocleavage and CUT&RUN, uliCUT&RUN is a tethered enzyme-based protocol that utilizes a combination of recombinant protein, antibody recognition and stringent purification to selectively target proteins of interest and isolate the associated DNA. Performed in native conditions, uliCUT&RUN profiles protein localization to chromatin with low input and high precision. Compared with other profiling technologies, uliCUT&RUN can determine nonhistone protein chromatin occupancies in low cell numbers, permitting the investigation into the molecular functions of a range of DNA-binding proteins within rare samples. From sample preparation to sequencing library submission, the uliCUT&RUN protocol takes <2 d to perform, with the accompanying data analysis timeline dependent on experience level.

Altered chromatin conformation and transcriptional regulation in watermelon following genome doubling

Polyploidy has played a crucial role in plant evolution, development and function. Synthetic autopolyploid represents an ideal system to investigate the effects of polyploidization on transcriptional regulation. In this study, we deciphered the impact of genome duplication at phenotypic and molecular levels in watermelon. Overall, 88% of the genes in tetraploid watermelon followed a >1:1 dosage effect, and accordingly, differentially expressed genes were largely upregulated. In addition, a great number of hypomethylated regions (1688) were identified in an isogenic tetraploid watermelon. These differentially methylated regions were localized in promoters and intergenic regions and near transcriptional start sites of the identified upregulated genes, which enhances the importance of methylation in gene regulation. These changes were reflected in sophisticated higher-order chromatin structures. The genome doubling caused switching of 108 A and 626 B compartments that harbored genes associated with growth, development and stress responses.

PAX8 and MECOM are interaction partners driving ovarian cancer

The transcription factor PAX8 is critical for the development of the thyroid and urogenital system. Comprehensive genomic screens furthermore indicate an additional oncogenic role for PAX8 in renal and ovarian cancers. While a plethora of PAX8-regulated genes in different contexts have been proposed, we still lack a mechanistic understanding of how PAX8 engages molecular complexes to drive disease-relevant oncogenic transcriptional programs. Here we show that protein isoforms originating from the MECOM locus form a complex with PAX8. These include MDS1-EVI1 (also called PRDM3) for which we map its interaction with PAX8 in vitro and in vivo. We show that PAX8 binds a large number of genomic sites and forms transcriptional hubs. At a subset of these, PAX8 together with PRDM3 regulates a specific gene expression module involved in adhesion and extracellular matrix. This gene module correlates with PAX8 and MECOM expression in large scale profiling of cell lines, patient-derived xenografts (PDXs) and clinical cases and stratifies gynecological cancer cases with worse prognosis. PRDM3 is amplified in ovarian cancers and we show that the MECOM locus and PAX8 sustain in vivo tumor growth, further supporting that the identified function of the MECOM locus underlies PAX8-driven oncogenic functions in ovarian cancer.

Multi-omics phenotyping of the gut-liver axis reveals metabolic perturbations from a low-dose pesticide mixture in rats

Health effects of pesticides are not always accurately detected using the current battery of regulatory toxicity tests. We compared standard histopathology and serum biochemistry measures and multi-omics analyses in a subchronic toxicity test of a mixture of six pesticides frequently detected in foodstuffs (azoxystrobin, boscalid, chlorpyrifos, glyphosate, imidacloprid and thiabendazole) in Sprague-Dawley rats. Analysis of water and feed consumption, body weight, histopathology and serum biochemistry showed little effect. Contrastingly, serum and caecum metabolomics revealed that nicotinamide and tryptophan metabolism were affected, which suggested activation of an oxidative stress response. This was not reflected by gut microbial community composition changes evaluated by shotgun metagenomics. Transcriptomics of the liver showed that 257 genes had their expression changed. Gene functions affected included the regulation of response to steroid hormones and the activation of stress response pathways. Genome-wide DNA methylation analysis of the same liver samples showed that 4,255 CpG sites were differentially methylated. Overall, we demonstrated that in-depth molecular profiling in laboratory animals exposed to low concentrations of pesticides allows the detection of metabolic perturbations that would remain undetected by standard regulatory biochemical measures and which could thus improve the predictability of health risks from exposure to chemical pollutants.

Heritability of DNA methylation in threespine stickleback (Gasterosteus aculeatus)

Epigenetic mechanisms underlying phenotypic change are hypothesized to contribute to population persistence and adaptation in the face of environmental change. To date, few studies have explored the heritability of intergenerationally stable methylation levels in natural populations, and little is known about the relative contribution of cis- and trans-regulatory changes to methylation variation. Here, we explore the heritability of DNA methylation, and conduct methylation quantitative trait loci (meQTLs) analysis to investigate the genetic architecture underlying methylation variation between marine and freshwater ecotypes of threespine stickleback (Gasterosteus aculeatus). We quantitatively measured genome-wide DNA methylation in fin tissue using reduced representation bisulfite sequencing of F1 and F2 crosses, and their marine and freshwater source populations. We identified cytosines (CpG sites) that exhibited stable methylation levels across generations. We found that additive genetic variance explained an average of 24-35% of the methylation variance, with a number of CpG sites possibly autonomous from genetic control. We also detected both cis- and trans-meQTLs, with only trans-meQTLs overlapping with previously identified genomic regions of high differentiation between marine and freshwater ecotypes. Finally, we identified the genetic architecture underlying two key CpG sites that were differentially methylated between ecotypes. These findings demonstrate a potential role for DNA methylation in facilitating adaptation to divergent environments and improve our understanding of the heritable basis of population epigenomic variation.

Exons and introns exhibit transcriptional strand asymmetry of dinucleotide distribution, damage formation and DNA repair

Recent cancer sequencing efforts have uncovered asymmetry in DNA damage induced mutagenesis between the transcribed and non-transcribed strands of genes. Here, we investigate the major type of damage induced by ultraviolet (UV) radiation, the cyclobutane pyrimidine dimers (CPDs), which are formed primarily in TT dinucleotides. We reveal that a transcriptional asymmetry already exists at the level of TT dinucleotide frequency and therefore also in CPD damage formation. This asymmetry is conserved in vertebrates and invertebrates and is completely reversed between introns and exons. We show the asymmetry in introns is linked to the transcription process itself, and is also found in enhancer elements. In contrast, the asymmetry in exons is not correlated to transcription, and is associated with codon usage preferences. Reanalysis of nucleotide excision repair, normalizing repair to the underlying TT frequencies, we show repair of CPDs is more efficient in exons compared to introns, contributing to the maintenance and integrity of coding regions. Our results highlight the importance of considering the primary sequence of the DNA in determining DNA damage sensitivity and mutagenic potential.

Integrated analysis of behavioral, epigenetic, and gut microbiome analyses in AppNL-G-F, AppNL-F, and wild type mice

Epigenetic mechanisms occurring in the brain as well as alterations in the gut microbiome composition might contribute to Alzheimer’s disease (AD). Human amyloid precursor protein knock-in (KI) mice contain the Swedish and Iberian mutations (App^NL-F) or those two and also the Arctic mutation (App^NL-G-F). In this study, we assessed whether behavioral and cognitive performance in 6-month-old App^NL-F, App^NL-G-F, and C57BL/6J wild-type (WT) mice was associated with the gut microbiome, and whether the genotype modulates this association. The genotype effects observed in behavioral tests were test-dependent. The biodiversity and composition of the gut microbiome linked to various aspects of mouse behavioral and cognitive performance but differences in genotype modulated these relationships. These genotype-dependent associations include members of the Lachnospiraceae and Ruminococcaceae families. In a subset of female mice, we assessed DNA methylation in the hippocampus and investigated whether alterations in hippocampal DNA methylation were associated with the gut microbiome. Among other differentially methylated regions, we identified a 1 Kb region that overlapped ing 3′UTR of the Tomm40 gene and the promoter region of the Apoe gene that and was significantly more methylated in the hippocampus of App^NL-G-F than WT mice. The integrated gut microbiome hippocampal DNA methylation analysis revealed a positive relationship between amplicon sequence variants (ASVs) within the Lachnospiraceae family and methylation at the Apoe gene. Hence, these microbes may elicit an impact on AD-relevant behavioral and cognitive performance via epigenetic changes in AD-susceptibility genes in neural tissue or that such changes in the epigenome can elicit alterations in intestinal physiology that affect the growth of these taxa in the gut microbiome.

Characterizing RNA stability genome-wide through combined analysis of PRO-seq and RNA-seq data

BACKGROUND

The concentrations of distinct types of RNA in cells result from a dynamic equilibrium between RNA synthesis and decay. Despite the critical importance of RNA decay rates, current approaches for measuring them are generally labor-intensive, limited in sensitivity, and/or disruptive to normal cellular processes. Here, we introduce a simple method for estimating relative RNA half-lives that is based on two standard and widely available high-throughput assays: Precision Run-On sequencing (PRO-seq) and RNA sequencing (RNA-seq).

RESULTS

Our method treats PRO-seq as a measure of transcription rate and RNA-seq as a measure of RNA concentration, and estimates the rate of RNA decay required for a steady-state equilibrium. We show that this approach can be used to assay relative RNA half-lives genome-wide, with good accuracy and sensitivity for both coding and noncoding transcription units. Using a structural equation model (SEM), we test several features of transcription units, nearby DNA sequences, and nearby epigenomic marks for associations with RNA stability after controlling for their effects on transcription. We find that RNA splicing-related features are positively correlated with RNA stability, whereas features related to miRNA binding and DNA methylation are negatively correlated with RNA stability. Furthermore, we find that a measure based on U1 binding and polyadenylation sites distinguishes between unstable noncoding and stable coding transcripts but is not predictive of relative stability within the mRNA or lincRNA classes. We also identify several histone modifications that are associated with RNA stability.

CONCLUSION

We introduce an approach for estimating the relative half-lives of individual RNAs. Together, our estimation method and systematic analysis shed light on the pervasive impacts of RNA stability on cellular RNA concentrations.

Interspecific Variation in One-Carbon Metabolism within the Ovarian Follicle, Oocyte, and Preimplantation Embryo: Consequences for Epigenetic Programming of DNA Methylation

One-carbon (1C) metabolism provides methyl groups for the synthesis and/or methylation of purines and pyrimidines, biogenic amines, proteins, and phospholipids. Our understanding of how 1C pathways operate, however, pertains mostly to the (rat) liver. Here we report that transcripts for all bar two genes (i.e., BHMT, MAT1A) encoding enzymes in the linked methionine-folate cycles are expressed in all cell types within the ovarian follicle, oocyte, and blastocyst in the cow, sheep, and pig; as well as in rat granulosa cells (GCs) and human KGN cells (a granulosa-like tumor cell line). Betaine-homocysteine methyltransferase (BHMT) protein was absent in bovine theca and GCs, as was activity of this enzyme in GCs. Mathematical modeling predicted that absence of this enzyme would lead to more volatile S-adenosylmethionine-mediated transmethylation in response to 1C substrate (e.g., methionine) or cofactor provision. We tested the sensitivity of bovine GCs to reduced methionine (from 50 to 10 µM) and observed a diminished flux of 1C units through the methionine cycle. We then used reduced-representation bisulfite sequencing to demonstrate that this reduction in methionine during bovine embryo culture leads to genome-wide alterations to DNA methylation in >1600 genes, including a cohort of imprinted genes linked to an abnormal fetal-overgrowth phenotype. Bovine ovarian and embryonic cells are acutely sensitive to methionine, but further experimentation is required to determine the significance of interspecific variation in BHMT expression.

Genome-wide analysis of allele-specific expression of genes in the model diatom Phaeodactylum tricornutum

Recent advances in next generation sequencing technologies have allowed the discovery of widespread autosomal allele-specific expression (aASE) in mammals and plants with potential phenotypic effects. Extensive numbers of genes with allele-specific expression have been described in the diatom Fragilariopsis cylindrus in association with adaptation to external cues, as well as in Fistulifera solaris in the context of natural hybridization. However, the role of aASE and its extent in diatoms remain elusive. In this study, we investigate allele-specific expression in the model diatom Phaeodactylum tricornutum by the re-analysis of previously published whole genome RNA sequencing data and polymorphism calling. We found that 22% of P. tricornutum genes show moderate bias in allelic expression while 1% show nearly complete monoallelic expression. Biallelic expression associates with genes encoding components of protein metabolism while moderately biased genes associate with functions in catabolism and protein transport. We validated candidate genes by pyrosequencing and found that moderate biases in allelic expression were less stable than monoallelically expressed genes that showed consistent bias upon experimental validations at the population level and in subcloning experiments. Our approach provides the basis for the analysis of aASE in P. tricornutum and could be routinely implemented to test for variations in allele expression under different environmental conditions.

Evolution of DNA methylome from precancerous lesions to invasive lung adenocarcinomas

The evolution of DNA methylome and methylation intra-tumor heterogeneity (ITH) during early carcinogenesis of lung adenocarcinoma has not been systematically studied. We perform reduced representation bisulfite sequencing of invasive lung adenocarcinoma and its precursors, atypical adenomatous hyperplasia, adenocarcinoma in situ and minimally invasive adenocarcinoma. We observe gradual increase of methylation aberrations and significantly higher level of methylation ITH in later-stage lesions. The phylogenetic patterns inferred from methylation aberrations resemble those based on somatic mutations suggesting parallel methylation and genetic evolution. De-convolution reveal higher ratio of T regulatory cells (Tregs) versus CD8 + T cells in later-stage diseases, implying progressive immunosuppression with neoplastic progression. Furthermore, increased global hypomethylation is associated with higher mutation burden, copy number variation burden and AI burden as well as higher Treg/CD8 ratio, highlighting the potential impact of methylation on chromosomal instability, mutagenesis and tumor immune microenvironment during early carcinogenesis of lung adenocarcinomas.

Liquid Biopsy of Methylation Biomarkers in Cell-Free DNA

Liquid biopsies, in particular, analysis of cell-free DNA (cfDNA), have emerged as a promising noninvasive diagnostic approach in oncology. Abnormal distribution of DNA methylation is one of the hallmarks of many cancers and methylation changes occur early during carcinogenesis. Systemic analysis of cfDNA methylation profiles is being developed for cancer early detection, monitoring for minimal residual disease (MRD), predicting treatment response and prognosis, and tracing the tissue origin. This review highlights the advantages and disadvantages of ctDNA profiling for noninvasive diagnosis of early-stage cancers and explores recent advances in the clinical application of ctDNA methylation assays. We also summarize the technologies for ctDNA methylation analysis and provide a brief overview of the bioinformatic approaches for analyzing DNA methylation sequencing data.

Long-read transcriptome sequencing reveals abundant promoter diversity in distinct molecular subtypes of gastric cancer

BACKGROUND

Deregulated gene expression is a hallmark of cancer; however, most studies to date have analyzed short-read RNA sequencing data with inherent limitations. Here, we combine PacBio long-read isoform sequencing (Iso-Seq) and Illumina paired-end short-read RNA sequencing to comprehensively survey the transcriptome of gastric cancer (GC), a leading cause of global cancer mortality.

RESULTS

We performed full-length transcriptome analysis across 10 GC cell lines covering four major GC molecular subtypes (chromosomal unstable, Epstein-Barr positive, genome stable and microsatellite unstable). We identify 60,239 non-redundant full-length transcripts, of which > 66% are novel compared to current transcriptome databases. Novel isoforms are more likely to be cell line and subtype specific, expressed at lower levels with larger number of exons, with longer isoform/coding sequence lengths. Most novel isoforms utilize an alternate first exon, and compared to other alternative splicing categories, are expressed at higher levels and exhibit higher variability. Collectively, we observe alternate promoter usage in 25% of detected genes, with the majority (84.2%) of known/novel promoter pairs exhibiting potential changes in their coding sequences. Mapping these alternate promoters to TCGA GC samples, we identify several cancer-associated isoforms, including novel variants of oncogenes. Tumor-specific transcript isoforms tend to alter protein coding sequences to a larger extent than other isoforms. Analysis of outcome data suggests that novel isoforms may impart additional prognostic information.

CONCLUSIONS

Our results provide a rich resource of full-length transcriptome data for deeper studies of GC and other gastrointestinal malignancies.

YAP contributes to DNA methylation remodeling upon mouse embryonic stem cell differentiation

The Yes-associated protein (YAP), one of the major effectors of the Hippo pathway together with its related protein WW-domain-containing transcription regulator 1 (WWTR1; also known as TAZ), mediates a range of cellular processes from proliferation and death to morphogenesis. YAP and WW-domain-containing transcription regulator 1 (WWTR1; also known as TAZ) regulate a large number of target genes, acting as coactivators of DNA-binding transcription factors or as negative regulators of transcription by interacting with the nucleosome remodeling and histone deacetylase complexes. YAP is expressed in self-renewing embryonic stem cells (ESCs), although it is still debated whether it plays any crucial roles in the control of either stemness or differentiation. Here we show that the transient downregulation of YAP in mouse ESCs perturbs cellular homeostasis, leading to the inability to differentiate properly. Bisulfite genomic sequencing revealed that this transient knockdown caused a genome-wide alteration of the DNA methylation remodeling that takes place during the early steps of differentiation, suggesting that the phenotype we observed might be due to the dysregulation of some of the mechanisms involved in regulation of ESC exit from pluripotency. By gene expression analysis, we identified two molecules that could have a role in the altered genome-wide methylation profile: the long noncoding RNA ephemeron, whose rapid upregulation is crucial for the transition of ESCs into epiblast, and the methyltransferase-like protein Dnmt3l, which, during the embryo development, cooperates with Dnmt3a and Dnmt3b to contribute to the de novo DNA methylation that governs early steps of ESC differentiation. These data suggest a new role for YAP in the governance of the epigenetic dynamics of exit from pluripotency.

Temporal genome-wide DNA methylation signature of post-smolt Pacific salmon challenged with Piscirickettsia salmonis

Piscirickettsiosis is the most important bacterial disease in the Chilean salmon industry, which has sorted several efforts to its control, generating enormous economic losses. Epigenetic alterations, such as DNA methylation, can play a relevant role in the modulation of the metazoans response to pathogens. Bacterial disease may activate global and local immune responses generating intricate responses with significant biological impact in the host. However, it is scarcely understood how bacterial infections influence fish epigenetic alterations. In the present study, we utilized Pacific salmon and Piscirickettsiosis as model, to gain understanding into the dynamics of DNA methylation among fish-bacterial infection interactions. A genome-wide analysis of DNA methylation patterns in female spleen tissue of Pacific salmon was achieved by reduced representation bisulphite sequencing from a time course design. We determined 2,251, 1,918, and 2,516 differentially methylated regions DMRs among infected and control Pacific salmon in 1 dpi, 5 dpi, and 15 dpi, respectively. The mean methylation difference per DMR among control and infected groups was of ~35%, with an oscillatory pattern of hypo, hyper, and hypomethylation across the disease. DMCs, among the control and infected group, showed that they were statistically enriched in intergenic regions and depleted in exons. Functional annotation of the DMR genes demonstrated three KEGG principal categories, associated directly with the host response to pathogens infections. Our results provide the first evidence of epigenetic variation in fish provoked by bacterial infection and demonstrate that this variation can be modulated across the disease.

DNA methylation during human adipogenesis and the impact of fructose

BACKGROUND

Increased adipogenesis and altered adipocyte function contribute to the development of obesity and associated comorbidities. Fructose modified adipocyte metabolism compared to glucose, but the regulatory mechanisms and consequences for obesity are unknown. Genome-wide methylation and global transcriptomics in SGBS pre-adipocytes exposed to 0, 2.5, 5, and 10 mM fructose, added to a 5-mM glucose-containing medium, were analyzed at 0, 24, 48, 96, 192, and 384 h following the induction of adipogenesis.

RESULTS

Time-dependent changes in DNA methylation compared to baseline (0 h) occurred during the final maturation of adipocytes, between 192 and 384 h. Larger percentages (0.1% at 192 h, 3.2% at 384 h) of differentially methylated regions (DMRs) were found in adipocytes differentiated in the glucose-containing control media compared to adipocytes differentiated in fructose-supplemented media (0.0006% for 10 mM, 0.001% for 5 mM, and 0.005% for 2.5 mM at 384 h). A total of 1437 DMRs were identified in 5237 differentially expressed genes at 384 h post-induction in glucose-containing (5 mM) control media. The majority of them inversely correlated with the gene expression, but 666 regions were positively correlated to the gene expression.

CONCLUSIONS

Our studies demonstrate that DNA methylation regulates or marks the transformation of morphologically differentiating adipocytes (seen at 192 h), to the more mature and metabolically robust adipocytes (as seen at 384 h) in a genome-wide manner. Lower (2.5 mM) concentrations of fructose have the most robust effects on methylation compared to higher concentrations (5 and 10 mM), suggesting that fructose may be playing a signaling/regulatory role at lower concentrations of fructose and as a substrate at higher concentrations.

GeneDMRs: An R Package for Gene-Based Differentially Methylated Regions Analysis

DNA methylation in gene or gene body could influence gene transcription. Moreover, methylation in gene regions along with CpG island regions could modulate the transcription to undetectable gene expression levels. Therefore, it is necessary to investigate the methylation levels within the gene, gene body, CpG island regions, and their overlapped regions and then identify the gene-based differentially methylated regions (GeneDMRs). In this study, R package GeneDMRs aims to facilitate computing gene-based methylation rate using next-generation sequencing-based methylome data. The user-friendly GeneDMRs package is presented to analyze the methylation levels in each gene/promoter/exon/intron/CpG island/CpG island shore or each overlapped region (e.g., gene-CpG island/promoter-CpG island/exon-CpG island/intron-CpG island/gene-CpG island shore/promoter-CpG island shore/exon-CpG island shore/intron-CpG island shore). GeneDMRs can also interpret complex interplays between methylation levels and gene expression differences or similarities across physiological conditions or disease states. We used the public reduced representation bisulfite sequencing data of mouse (GSE62392) for evaluating software and revealing novel biologically significant results to supplement the previous research. In addition, the whole-genome bisulfite sequencing data of cattle (GSE106538) given the much larger size was used for further evaluation.

Movement of transposable elements contributes to cichlid diversity

African cichlid fishes are a prime model for studying speciation mechanisms. Despite the development of extensive genomic resources, it has been difficult to determine which sources of genetic variation are responsible for cichlid phenotypic variation. One of their most variable phenotypes is visual sensitivity, with some of the largest spectral shifts among vertebrates. These shifts arise primarily from differential expression of seven cone opsin genes. By mapping expression quantitative trait loci (eQTL) in intergeneric crosses of Lake Malawi cichlids, we previously identified four causative genetic variants that correspond to indels in the promoters of either key transcription factors or an opsin gene. In this comprehensive study, we show that these indels are the result of the movement of transposable elements (TEs) that correlate with opsin expression variation across the Malawi flock. In tracking the evolutionary history of these particular indels, we found they are endemic to Lake Malawi, suggesting that these TEs are recently active and are segregating within the Malawi cichlid lineage. However, an independent indel has arisen at a similar genomic location in one locus outside of the Malawi flock. The convergence in TE movement suggests these loci are primed for TE insertion and subsequent deletions. Increased TE mobility may be associated with interspecific hybridization, which disrupts mechanisms of TE suppression. This might provide a link between cichlid hybridization and accelerated regulatory variation. Overall, our study suggests that TEs may be an important driver of key regulatory changes, facilitating rapid phenotypic change and possibly speciation in African cichlids.

Genome-Wide Assessment of DNA Methylation in Chicken Cardiac Tissue Exposed to Different Incubation Temperatures and CO2 Levels

Temperature and CO₂ concentration during incubation have profound effects on broiler chick development, and numerous studies have identified significant effects on hatch heart weight (HW) as a result of differences in these parameters. Early life environment has also been shown to affect broiler performance later in life; it has thus been suggested that epigenetic mechanisms may mediate long-term physiological changes induced by environmental stimuli. DNA methylation is an epigenetic modification that can confer heritable changes in gene expression. Using reduced-representation bisulfite sequencing (RRBS), we assessed DNA methylation patterns in cardiac tissue of 84 broiler hatchlings incubated at two egg shell temperatures (EST; 37.8°C and 38.9°C) and three CO₂ concentrations (0.1%, 0.4%, and 0.8%) from day 8 of incubation onward. We assessed differential methylation between EST treatments and identified 2,175 differentially methylated (DM) CpGs (1,121 hypermethylated, 1,054 hypomethylated at 38.9° vs. 37.8°) in 269 gene promoters and 949 intragenic regions. DM genes (DMGs) were associated with heart developmental processes, including cardiomyocyte proliferation and differentiation. We identified enriched binding motifs among DM loci, including those for transcription factors associated with cell proliferation and heart development among hypomethylated CpGs that suggest increased binding ability at higher EST. We identified 9,823 DM CpGs between at least two CO₂ treatments, with the greatest difference observed between 0.8 and 0.1% CO₂ that disproportionately impacted genes involved in cardiac muscle development and response to low oxygen levels. Using HW measurements from the same chicks, we performed an epigenome-wide association study (EWAS) for HW, and identified 23 significantly associated CpGs, nine of which were also DM between ESTs. We found corresponding differences in transcript abundance between ESTs in three DMGs (ABLIM2, PITX2, and THRSP). Hypomethylation of an exonic CpG in PITX2 at 38.9°C was associated with increased expression, and suggests increased cell proliferation in broiler hatchlings incubated at higher temperatures. Overall, these results identified numerous epigenetic associations between chick incubation factors and heart development that may manifest in long-term differences in animal performance.

The Gene-Regulatory Footprint of Aging Highlights Conserved Central Regulators

Many genes and pathways have been linked to aging, yet our understanding of underlying molecular mechanisms is still lacking. Here, we measure changes in the transcriptome, histone modifications, and DNA methylome in three metabolic tissues of adult and aged mice. Transcriptome and methylome changes dominate the liver aging footprint, whereas heart and muscle globally increase chromatin accessibility, especially in aging pathways. In mouse and human data from multiple tissues and regulatory layers, age-related transcription factor expression changes and binding site enrichment converge on putative aging modulators, including ZIC1, CXXC1, HMGA1, MECP2, SREBF1, SREBF2, ETS2, ZBTB7A, and ZNF518B. Using Mendelian randomization, we establish possible epidemiological links between expression of some of these transcription factors or their targets, including CXXC1, ZNF518B, and BBC3, and longevity. We conclude that conserved modulators are at the core of the molecular footprint of aging, and variation in tissue-specific expression of some may affect human longevity.

Multi-Omics Investigation of Innate Navitoclax Resistance in Triple-Negative Breast Cancer Cells

Cancer cells employ various defense mechanisms against drug-induced cell death. Investigating multi-omics landscapes of cancer cells before and after treatment can reveal resistance mechanisms and inform new therapeutic strategies. We assessed the effects of navitoclax, a BCL2 family inhibitor, on the transcriptome, methylome, chromatin structure, and copy number variations of MDA-MB-231 triple-negative breast cancer (TNBC) cells. Cells were sampled before treatment, at 72 h of exposure, and after 10-day drug-free recovery from treatment. We observed transient alterations in the expression of stress response genes that were accompanied by corresponding changes in chromatin accessibility. Most of these changes returned to baseline after the recovery period. We also detected lasting alterations in methylation states and genome structure that suggest permanent changes in cell population composition. Using single-cell analyses, we identified 2350 genes significantly upregulated in navitoclax-resistant cells and derived an 18-gene navitoclax resistance signature. We assessed the navitoclax-response-predictive function of this signature in four additional TNBC cell lines in vitro and in silico in 619 cell lines treated with 251 different drugs. We observed a drug-specific predictive value in both experiments, suggesting that this signature could help guiding clinical biomarker studies involving navitoclax.

Experimental Parasite Infection Causes Genome-Wide Changes in DNA Methylation

Parasites are arguably among the strongest drivers of natural selection, constraining hosts to evolve resistance and tolerance mechanisms. Although, the genetic basis of adaptation to parasite infection has been widely studied, little is known about how epigenetic changes contribute to parasite resistance and eventually, adaptation. Here, we investigated the role of host DNA methylation modifications to respond to parasite infections. In a controlled infection experiment, we used the three-spined stickleback fish, a model species for host-parasite studies, and their nematode parasite Camallanus lacustris. We showed that the levels of DNA methylation are higher in infected fish. Results furthermore suggest correlations between DNA methylation and shifts in key fitness and immune traits between infected and control fish, including respiratory burst and functional trans-generational traits such as the concentration of motile sperm. We revealed that genes associated with metabolic, developmental, and regulatory processes (cell death and apoptosis) were differentially methylated between infected and control fish. Interestingly, genes such as the neuropeptide FF receptor 2 and the integrin alpha 1 as well as molecular pathways including the Th1 and Th2 cell differentiation were hypermethylated in infected fish, suggesting parasite-mediated repression mechanisms of immune responses. Altogether, we demonstrate that parasite infection contributes to genome-wide DNA methylation modifications. Our study brings novel insights into the evolution of vertebrate immunity and suggests that epigenetic mechanisms are complementary to genetic responses against parasite-mediated selection.

Contrasting Roles of Transcription Factors Spineless and EcR in the Highly Dynamic Chromatin Landscape of Butterfly Wing Metamorphosis

Development requires highly coordinated changes in chromatin accessibility in order for proper gene regulation to occur. Here, we identify factors associated with major, discrete changes in chromatin accessibility during butterfly wing metamorphosis. By combining mRNA sequencing (mRNA-seq), assay for transposase-accessible chromatin using sequencing (ATAC-seq), and machine learning analysis of motifs, we show that distinct sets of transcription factors are predictive of chromatin opening at different developmental stages. Our data suggest an important role for nuclear hormone receptors early in metamorphosis, whereas PAS-domain transcription factors are strongly associated with later chromatin opening. Chromatin immunoprecipitation sequencing (ChIP-seq) validation of select candidate factors showed spineless binding to be a major predictor of opening chromatin. Surprisingly, binding of ecdysone receptor (EcR), a candidate accessibility factor in Drosophila, was not predictive of opening but instead marked persistent sites. This work characterizes the chromatin dynamics of insect wing metamorphosis, identifies candidate chromatin remodeling factors in insects, and presents a genome assembly of the model butterfly Junonia coenia.

Single-base cytosine methylation analysis in fruits of three Capsicum species

Single-base cytosine methylation analysis across fruits of Capsicum annuum, C. chinense and C. frutescens showed global average methylation ranging from 82.8-89.1%, 77.6-83.9%, and 22.4-25% at CG, CHG and CHH contexts, respectively. High gene-body methylation at CG and CHG was observed across Capsicum species. The C. annuum showed the highest proportion (>80%) of mCs at different genomic regions compared to C. chinense and C. frutescens. Cytosine methylation for transposable-elements were lower in C. frutescens compared to C. annuum and C. chinense. A total of 510,165 CG, 583112 CHG and 277,897 CHH DMRs were identified across three Capsicum species. The differentially methylated regions (DMRs) distribution analysis revealed C. frutescens as more hypo-methylated compared to C. annuum and C. chinense, and also the presence of more intergenic DMRs in Capsicum genome. At CG and CHG context, gene expression and promoter methylation showed inverse correlations. Furthermore, the observed correlation between methylation and expression of genes suggested the potential role of methylation in Capsicum fruit development/ripening.

Genome‑wide analysis of DNA methylation and gene expression changes in an ovalbumin‑induced asthma mouse model

The aim of the present study was to establish an integrated network of DNA methylation and RNA expression in an ovalbumin (OVA)-induced asthma model, and to investigate the epigenetically-regulated genes involved in asthma development. Genome-wide CpG-DNA methylation profiling was conducted through the use of a methylated DNA immunoprecipitation microarray and RNA sequencing was performed using three lung samples from mice with OVA-induced asthma. A total of 35,401 differentially methylated regions (DMRs) were identified between mice with OVA-induced asthma and control mice. Of these, 3,060 were located in promoter regions and 370 of the genes containing these DMRs demonstrated an inverse correlation between methylation and gene expression. Kyoto Encyclopedia of Genes and Genomes pathway analysis identified that 368 genes were upregulated or downregulated in OVA-induced asthma samples, including genes involved in ‘chemokine signalling pathway’, ‘focal adhesion’, ‘leukocyte transendothelial migration’ and ‘vascular smooth muscle contraction signaling’ pathways. Integrated network analysis identified four hub genes, consisting of three upregulated genes [forkhead box O1 (FOXO1), SP1 transcription factor (SP1) and amyloid β precursor protein (APP)], and one downregulated gene [RUNX family transcription factor 1 (RUNX1)], all of which demonstrated an association between DNA methylation and gene expression. These genes were highly interconnected nodes in the Ingenuity Pathway Analysis module and were functionally significant. A total of four interconnected hub genes, FOXO1, RUNX1, SP1 and APP, were identified from the integrated DNA methylation and gene expression networks involved in asthma development. These results suggested that modulating these four genes could effectively control the development of asthma.

Characterization of Global DNA Methylation in Different Gene Regions Reveals Candidate Biomarkers in Pigs with High and Low Levels of Boar Taint

DNA methylation of different gene components, including different exons and introns, or different lengths of exons and introns is associated with differences in gene expression. To investigate the methylation of porcine gene components associated with the boar taint (BT) trait, this study used reduced representation bisulfite sequencing (RRBS) data from nine porcine testis samples in three BT groups (low, medium and high BT). The results showed that the methylation levels of the first exons and first introns were lower than those of the other exons and introns. The first exons/introns of CpG island regions had even lower levels of methylation. A total of 123 differentially methylated promoters (DMPs), 194 differentially methylated exons (DMEs) and 402 differentially methylated introns (DMIs) were identified, of which 80 DMPs (DMP-CpGis), 112 DMEs (DME-CpGis) and 166 DMIs (DMI-CpGis) were discovered in CpG islands. Importantly, GPX1 contained one each of DMP, DME, DMI, DMP-CpGi, DME-CpGi and DMI-CpGi. Gene-GO term relationships and pathways analysis showed DMP-CpGi-related genes are mainly involved in methylation-related biological functions. In addition, gene-gene interaction networks consisted of nodes that were hypo-methylated GPX1, hypo-methylated APP, hypo-methylated ATOX1, hyper-methylated ADRB2, hyper-methylated RPS6KA1 and hyper-methylated PNMT. They could be used as candidate biomarkers for reducing boar taint in pigs, after further validation in large cohorts.

High-Resolution In Vivo Identification of miRNA Targets by Halo-Enhanced Ago2 Pull-Down

The identification of microRNA (miRNA) targets by Ago2 crosslinking-immunoprecipitation (CLIP) methods has provided major insights into the biology of this important class of non-coding RNAs. However, these methods are technically challenging and not easily applicable to an in vivo setting. To overcome these limitations and facilitate the investigation of miRNA functions in vivo, we have developed a method based on a genetically engineered mouse harboring a conditional Halo-Ago2 allele expressed from the endogenous Ago2 locus. By using a resin conjugated to the HaloTag ligand, Ago2-miRNA-mRNA complexes can be purified from cells and tissues expressing the endogenous Halo-Ago2 allele. We demonstrate the reproducibility and sensitivity of this method in mouse embryonic stem cells, developing embryos, adult tissues, and autochthonous mouse models of human brain and lung cancers. This method and the datasets we have generated will facilitate the characterization of miRNA-mRNA networks in vivo under physiological and pathological conditions.

Global changes in nitration levels and DNA binding profile of Trypanosoma cruzi histones induced by incubation with host extracellular matrix

Adhesion of T. cruzi trypomastigotes to components of the extracellular matrix (ECM) is an important step in mammalian host cell invasion. We have recently described a significant increase in the tyrosine nitration levels of histones H2A and H4 when trypomastigotes are incubated with components of the ECM. In this work, we used chromatin immunoprecipitation (ChIP) with an anti-nitrotyrosine antibody followed by mass spectrometry to identify nitrated DNA binding proteins in T. cruzi and to detect alterations in nitration levels induced upon parasite incubation with the ECM. Histone H1, H2B, H2A and H3 were detected among the 9 most abundant nitrated DNA binding proteins using this proteomic approach. One nitrated tyrosine residue (Y29) was identified in Histone H2B in the MS/MS spectrum. In addition, we observed a significant increase in the nitration levels of histones H1, H2B, H2A and H4 upon parasite incubation with ECM. Finally, we used ChIP-Seq to map global changes in the DNA binding profile of nitrated proteins. We observed a significant change in the binding pattern of nitrated proteins to DNA after parasite incubation with ECM. This work provides the first global profile of nitrated DNA binding proteins in T. cruzi and additional evidence for modification in the nitration profile of histones upon parasite incubation with ECM. Our data also indicate that the parasite interaction with the ECM induces alterations in chromatin structure, possibly affecting nuclear functions.

The cooperative action of CSB, CSA, and UVSSA target TFIIH to DNA damage-stalled RNA polymerase II

The response to DNA damage-stalled RNA polymerase II (RNAPIIo) involves the assembly of the transcription-coupled repair (TCR) complex on actively transcribed strands. The function of the TCR proteins CSB, CSA and UVSSA and the manner in which the core DNA repair complex, including transcription factor IIH (TFIIH), is recruited are largely unknown. Here, we define the assembly mechanism of the TCR complex in human isogenic knockout cells. We show that TCR is initiated by RNAPIIo-bound CSB, which recruits CSA through a newly identified CSA-interaction motif (CIM). Once recruited, CSA facilitates the association of UVSSA with stalled RNAPIIo. Importantly, we find that UVSSA is the key factor that recruits the TFIIH complex in a manner that is stimulated by CSB and CSA. Together these findings identify a sequential and highly cooperative assembly mechanism of TCR proteins and reveal the mechanism for TFIIH recruitment to DNA damage-stalled RNAPIIo to initiate repair.

The response to DNA damage-stalled RNA polymerase II leads to the assembly of the transcription-coupled repair (TCR) complex on actively transcribed strands. Here, the authors reveal the complex assembly mechanism of the TCR complex in human cells.

Reactivation of Myc transcription in the mouse heart unlocks its proliferative capacity

It is unclear why some tissues are refractory to the mitogenic effects of the oncogene Myc. Here we show that Myc activation induces rapid transcriptional responses followed by proliferation in some, but not all, organs. Despite such disparities in proliferative response, Myc is bound to DNA at open elements in responsive (liver) and non-responsive (heart) tissues, but fails to induce a robust transcriptional and proliferative response in the heart. Using heart as an exemplar of a non-responsive tissue, we show that Myc-driven transcription is re-engaged in mature cardiomyocytes by elevating levels of the positive transcription elongation factor (P-TEFb), instating a large proliferative response. Hence, P-TEFb activity is a key limiting determinant of whether the heart is permissive for Myc transcriptional activation. These data provide a greater understanding of how Myc transcriptional activity is determined and indicate modification of P-TEFb levels could be utilised to drive regeneration of adult cardiomyocytes for the treatment of heart myopathies.

Robust mapping of polyadenylated and non-polyadenylated RNA 3' ends at nucleotide resolution by 3'-end sequencing

3'-end poly(A)+ sequencing is an efficient and economical method for global measurement of mRNA levels and alternative poly(A) site usage. A common method involves oligo(dT)19V reverse-transcription (RT)-based library preparation and high-throughput sequencing with a custom primer ending in (dT)19. While the majority of library products have the first sequenced nucleotide reflect the bona fide poly(A) site (pA), a substantial fraction of sequencing reads arise from various mis-priming events. These can result in incorrect pA site calls anywhere from several nucleotides downstream to several kilobases upstream from the bona fide pA site. While these mis-priming events can be mitigated by increasing annealing stringency (e.g. increasing temperature from 37 °C to 42 °C), they still persist at an appreciable level (∼10%) and computational methods must be used to prevent artifactual calls. Here we present a bioinformatics workflow for precise mapping of poly(A)+ 3' ends and handling of artifacts due to oligo(dT) mis-priming and sample polymorphisms. We test pA site calling with three different read mapping programs (STAR, BWA, and BBMap), and show that the way in which each handles terminal mismatches and soft clipping has a substantial impact on identifying correct pA sites, with BWA requiring the least post-processing to correct artifacts. We demonstrate the use of this pipeline for mapping pA sites in the model eukaryote S. cerevisiae, and further apply this technology to non-polyadenylated transcripts by employing in vitro polyadenylation prior to library prep (IVP-seq). As proof of principle, we show that a fraction of tRNAs harbor CCU 3' tails instead of the canonical CCA tail, and globally identify 3' ends of splicing intermediates arising from inefficiently spliced transcripts.

RNA-guided retargeting of Sleeping Beauty transposition in human cells

An ideal tool for gene therapy would enable efficient gene integration at predetermined sites in the human genome. Here we demonstrate biased genome-wide integration of the Sleeping Beauty (SB) transposon by combining it with components of the CRISPR/Cas9 system. We provide proof-of-concept that it is possible to influence the target site selection of SB by fusing it to a catalytically inactive Cas9 (dCas9) and by providing a single guide RNA (sgRNA) against the human Alu retrotransposon. Enrichment of transposon integrations was dependent on the sgRNA, and occurred in an asymmetric pattern with a bias towards sites in a relatively narrow, 300 bp window downstream of the sgRNA targets. Our data indicate that the targeting mechanism specified by CRISPR/Cas9 forces integration into genomic regions that are otherwise poor targets for SB transposition. Future modifications of this technology may allow the development of methods for specific gene insertion for precision genetic engineering.

ChromID identifies the protein interactome at chromatin marks

Chromatin modifications regulate genome function by recruiting protein factors to the genome. However, the protein composition at distinct chromatin modifications remains to be fully characterized. Here, we use natural protein domains as modular building blocks to develop engineered chromatin readers (eCRs) selective for DNA methylation and histone tri-methylation at H3K4, H3K9 a H3K27 residues. We first demonstrate their utility as selective chromatin binders in living cells by stably expressing eCRs in mouse embryonic stem cells and measuring their subnuclear localisation, genomic distribution and histone modification–binding preference. By fusing eCRs to the biotin ligase BASU, we establish ChromID, a method for identifying the chromatin-dependent protein interactome based on proximity biotinylation, and apply it to distinct chromatin modifications in mouse stem cells. Using a synthetic dual-modification reader, we also uncover the protein composition at bivalent promoters marked by H3K4me3 and H3K27me3. These results highlight the ability of ChromID to obtain a detailed view of protein interaction networks on chromatin.

Two different epigenetic information channels in wild three-spined sticklebacks are involved in salinity adaptation

Epigenetic inheritance has been proposed to contribute to adaptation and acclimation via two information channels: (i) inducible epigenetic marks that enable transgenerational plasticity and (ii) noninducible epigenetic marks resulting from random epimutations shaped by selection. We studied both postulated channels by sequencing methylomes and genomes of Baltic three-spined sticklebacks (Gasterosteus aculeatus) along a salinity cline. Wild populations differing in salinity tolerance revealed differential methylation (pop-DMS) at genes enriched for osmoregulatory processes. A two-generation experiment demonstrated that 62% of these pop-DMS were noninducible by salinity manipulation, suggesting that they are the result of either direct selection or associated genomic divergence at cis- or trans-regulatory sites. Two-thirds of the remaining inducible pop-DMS increased in similarity to patterns detected in wild populations from corresponding salinities. The level of similarity accentuated over consecutive generations, indicating a mechanism of transgenerational plasticity. While we can attribute natural DNA methylation patterns to the two information channels, their interplay with genomic variation in salinity adaptation is still unresolved.

RADICL-seq identifies general and cell type-specific principles of genome-wide RNA-chromatin interactions

Mammalian genomes encode tens of thousands of noncoding RNAs. Most noncoding transcripts exhibit nuclear localization and several have been shown to play a role in the regulation of gene expression and chromatin remodeling. To investigate the function of such RNAs, methods to massively map the genomic interacting sites of multiple transcripts have been developed; however, these methods have some limitations. Here, we introduce RNA And DNA Interacting Complexes Ligated and sequenced (RADICL-seq), a technology that maps genome-wide RNA–chromatin interactions in intact nuclei. RADICL-seq is a proximity ligation-based methodology that reduces the bias for nascent transcription, while increasing genomic coverage and unique mapping rate efficiency compared with existing methods. RADICL-seq identifies distinct patterns of genome occupancy for different classes of transcripts as well as cell type–specific RNA-chromatin interactions, and highlights the role of transcription in the establishment of chromatin structure.

Mammalian genomes encode tens of thousands of ncRNAs that have important roles in regulation of gene expression and chromatin organization. Here, the authors present RADICLseq to map RNA-chromatin interactions in intact nuclei to shed light on these fine-tuned processes.

Epigenetic reprogramming at estrogen-receptor binding sites alters 3D chromatin landscape in endocrine-resistant breast cancer

Endocrine therapy resistance frequently develops in estrogen receptor positive (ER+) breast cancer, but the underlying molecular mechanisms are largely unknown. Here, we show that 3-dimensional (3D) chromatin interactions both within and between topologically associating domains (TADs) frequently change in ER+ endocrine-resistant breast cancer cells and that the differential interactions are enriched for resistance-associated genetic variants at CTCF-bound anchors. Ectopic chromatin interactions are preferentially enriched at active enhancers and promoters and ER binding sites, and are associated with altered expression of ER-regulated genes, consistent with dynamic remodelling of ER pathways accompanying the development of endocrine resistance. We observe that loss of 3D chromatin interactions often occurs coincidently with hypermethylation and loss of ER binding. Alterations in active A and inactive B chromosomal compartments are also associated with decreased ER binding and atypical interactions and gene expression. Together, our results suggest that 3D epigenome remodelling is a key mechanism underlying endocrine resistance in ER+ breast cancer.

Acute PFOA exposure promotes epigenomic alterations in mouse kidney tissues

Perfluorooctanoic acid (PFOA), a manufactured perfluorochemical is a common surfactant and environmental pollutant found in various consumer products and water sources. Epidemiological studies have demonstrated its association with kidney dysfunction. However, the mechanisms that trigger kidney dysfunction following PFOA exposure is a gap in the field. The work presented explores the potential epigenetic indicators of kidney disease due to exposure to PFOA. In this study, 30 days old CD-1 mice were exposed to 1, 5, 10, or 20 mg/kg/day of PFOA for 10 days. Following acute oral exposure, epigenetic alterations and expression levels of various markers of fibroblast activation were evaluated in kidney tissues. We noted that PFOA-exposed mice exhibited differential methylation yielding 879 differentially methylated regions compared to vehicle. The mRNA expression revealed significant increase in Dnmt1 with decreased Rasal1 expression at higher levels of PFOA exposure suggestive of Rasal1 hypermethylation (an early indicator of fibroblast activation in kidney). Like Dnmt1, we also observed significant increase in Hdac1, 3 and 4. These are class I & II HDACs which are known to be critically altered in some renal diseases. Further, the mRNA expression levels of TGF-β and α-SMA significantly increased compared to vehicle. The KEGG and Go enrichment pathway analysis of reduced representation bisulfite data also revealed pathways implicated in renal fibrosis. Our study shows clear evidence of epigenetic alterations (DNA methylation and HDAC expression changes) in tissues from mouse kidney following PFOA exposure. Our results also suggest that epigenetic alterations in kidney promote the expression of early markers of fibroblast activation.

In-depth characterization of the placental imprintome reveals novel differentially methylated regions across birth weight categories

Imprinted genes play a pivotal role in placental processes underlying fetal development, and much interest centers on discerning whether these loci, via changes in DNA methylation and/or gene expression, inform disruptions in appropriate fetal growth. In this study, we comprehensively profiled DNA methylation across the placental imprintome and assessed the relationship with gene expression levels and aberrant fetal growth.Placental DNA methylation across 153 imprinted loci, including imprint control regions (ICR) and surrounding non-ICR regions, was surveyed using the Nimblegen TruSeq bisulfite sequencing platform among participants enrolled in the Rhode Island Child Health Study (RICHS, n = 163). Methylation and gene expression associations were assessed using eQTM analysis. Differential methylation analysis contrasting small (SGA) and large for gestational age (LGA) infants against appropriate for gestational age (AGA) infants was assessed using the DMRcate R package.We identified 34 SGA-related differentially methylated regions (DMRs) and 9 LGA-related DMRs (FDR<0.05), and these BW-DMRs predominated in promoter and intronic regions. We observed overall hypomethylation among SGA-DMRs overlapping maternally expressed (paternally imprinted) genes while no parent-of-origin effect was observed among LGA DMRs. Three BW-DMRs, mapping to GABRG3, IGF1R and MEST, were common to SGA and LGA placenta. We did not observe significant correlations between BW-DMR-associated CpG methylation and gene expression levels.We report the first in-depth characterization of the placental imprintome in a population-wide setting. Our findings reveal growth-related differences in methylation without concomitant expression differences in regions that extend beyond typically interrogated imprinted loci, highlighting potentially novel placental biomarkers of growth and development.

Preclinical Evaluation of a Cell-Based Gene Therapy Using the Sleeping Beauty Transposon System in Choroidal Neovascularization

Age-related macular degeneration (AMD) is a progressive retinal disorder characterized by imbalanced pro- and antiangiogenic signals. The aim of this study was to evaluate the effect of ex vivo cell-based gene therapy with stable expression of human pigment epithelium-derived factor (PEDF) release using the non-viral Sleeping Beauty (SB100X) transposon system delivered by miniplasmids free of antibiotic resistance markers (pFAR4). Retinal pigment epithelial (RPE) cells and iris pigment epithelial (IPE) cells were co-transfected with pFAR4-inverted terminal repeats (ITRs) CMV-PEDF-BGH and pFAR4-CMV-SB100X-SV40 plasmids. Laser-induced choroidal neovascularization (CNV) was performed in rats, and transfected primary cells (transfected RPE [tRPE] and transfected IPE [tIPE] cells) were injected into the subretinal space. The leakage and CNV areas, vascular endothelial growth factor (VEGF), PEDF protein expression, metalloproteinases 2 and 9 (MMP-2/9), and microglial/macrophage markers were measured. Injection with tRPE/IPE cells significantly reduced the leakage area at 7 and 14 days and the CNV area at 7 days. There was a significant increase in PEDF and the PEDF/VEGF ratio with tRPE cells and a reduction in the MMP-2 activity. Our data demonstrated that ex vivo non-viral gene therapy reduces CNV and could be an effective and safe therapeutic option for angiogenic retinal diseases.

Genome-wide epigenetic isolation by environment in a widespread Anolis lizard

Epigenetic changes can provide a pathway for organisms to respond to local environmental conditions by influencing gene expression. However, we still know little about the spatial distribution of epigenetic variation in natural systems, how it relates to the distribution of genetic variation and the environmental structure of the landscape, and the processes that generate and maintain it. Studies examining spatial patterns of genetic and epigenetic variation can provide valuable insights into how ecological and population processes contribute to epigenetic divergence across heterogeneous landscapes. Here, we perform a comparative analysis of spatial genetic and epigenetic variation based on 8,459 single nucleotide polymorphisms (SNPs) and 8,580 single methylation variants (SMVs) from eight populations of the Puerto Rican crested anole, Anolis cristatellus, an abundant lizard in the adaptive radiations of anoles on the Greater Antilles that occupies a diverse range of habitats. Using generalized dissimilarity modelling and multiple matrix regression, we found that genome-wide epigenetic differentiation is strongly correlated with environmental divergence, even after controlling for the underlying genetic structure. We also detected significant associations between key environmental variables and 96 SMVs, including 42 located in promoter regions or gene bodies. Our results suggest an environmental basis for population-level epigenetic differentiation in this system and contribute to better understanding how environmental gradients structure epigenetic variation in nature.

HOT or not: examining the basis of high-occupancy target regions

High-occupancy target (HOT) regions are segments of the genome with unusually high number of transcription factor binding sites. These regions are observed in multiple species and thought to have biological importance due to high transcription factor occupancy. Furthermore, they coincide with house-keeping gene promoters and consequently associated genes are stably expressed across multiple cell types. Despite these features, HOT regions are solely defined using ChIP-seq experiments and shown to lack canonical motifs for transcription factors that are thought to be bound there. Although, ChIP-seq experiments are the golden standard for finding genome-wide binding sites of a protein, they are not noise free. Here, we show that HOT regions are likely to be ChIP-seq artifacts and they are similar to previously proposed ‘hyper-ChIPable’ regions. Using ChIP-seq data sets for knocked-out transcription factors, we demonstrate presence of false positive signals on HOT regions. We observe sequence characteristics and genomic features that are discriminatory of HOT regions, such as GC/CpG-rich k-mers, enrichment of RNA–DNA hybrids (R-loops) and DNA tertiary structures (G-quadruplex DNA). The artificial ChIP-seq enrichment on HOT regions could be associated to these discriminatory features. Furthermore, we propose strategies to deal with such artifacts for the future ChIP-seq studies.

Circulating Tumor Cell Clustering Shapes DNA Methylation to Enable Metastasis Seeding

The ability of circulating tumor cells (CTCs) to form clusters has been linked to increased metastatic potential. Yet biological features and vulnerabilities of CTC clusters remain largely unknown. Here, we profile the DNA methylation landscape of single CTCs and CTC clusters from breast cancer patients and mouse models on a genome-wide scale. We find that binding sites for stemness- and proliferation-associated transcription factors are specifically hypomethylated in CTC clusters, including binding sites for OCT4, NANOG, SOX2, and SIN3A, paralleling embryonic stem cell biology. Among 2,486 FDA-approved compounds, we identify Na⁺/K⁺ ATPase inhibitors that enable the dissociation of CTC clusters into single cells, leading to DNA methylation remodeling at critical sites and metastasis suppression. Thus, our results link CTC clustering to specific changes in DNA methylation that promote stemness and metastasis and point to cluster-targeting compounds to suppress the spread of cancer.

Cell-Free DNA Methylation Profiling Analysis-Technologies and Bioinformatics

Analysis of circulating nucleic acids in bodily fluids, referred to as “liquid biopsies”, is rapidly gaining prominence. Studies have shown that cell-free DNA (cfDNA) has great potential in characterizing tumor status and heterogeneity, as well as the response to therapy and tumor recurrence. DNA methylation is an epigenetic modification that plays an important role in a broad range of biological processes and diseases. It is well known that aberrant DNA methylation is generalizable across various samples and occurs early during the pathogenesis of cancer. Methylation patterns of cfDNA are also consistent with their originated cells or tissues. Systemic analysis of cfDNA methylation profiles has emerged as a promising approach for cancer detection and origin determination. In this review, we will summarize the technologies for DNA methylation analysis and discuss their feasibility for liquid biopsy applications. We will also provide a brief overview of the bioinformatic approaches for analysis of DNA methylation sequencing data. Overall, this review provides informative guidance for the selection of experimental and computational methods in cfDNA methylation-based studies.

Using R and Bioconductor in Clinical Genomics and Transcriptomics

Bioinformatics pipelines are essential in the analysis of genomic and transcriptomic data generated by next-generation sequencing (NGS). Recent guidelines emphasize the need for rigorous validation and assessment of robustness, reproducibility, and quality of NGS analytic pipelines intended for clinical use. Software tools written in the R statistical language and, in particular, the set of tools available in the Bioconductor repository are widely used in research bioinformatics; and these frameworks offer several advantages for use in clinical bioinformatics, including the breath of available tools, modular nature of software packages, ease of installation, enforcement of interoperability, version control, and short learning curve. This review provides an introduction to R and Bioconductor software, its advantages and limitations for clinical bioinformatics, and illustrative examples of tools that can be used in various steps of NGS analysis.

Genome-wide DNA methylation analysis using next-generation sequencing to reveal candidate genes responsible for boar taint in pigs

Boar taint (BT) is an offensive flavor observed in non-castrated male pigs that reduces the carcass price. Surgical castration effectively avoids the taint but is associated with animal welfare concerns. The functional annotation of farm animal genomes for understanding the biology of complex traits can be used in the selection of breeding animals to achieve favorable phenotypic outcomes. The characterization of pig epigenomes/methylation changes between animals with high and low BT and genome-wide epigenetic markers that can predict BT are lacking. Reduced representation bisulfite sequencing of DNA methylation patterns based on next-generation sequencing is an efficient technology to identify candidate epigenetic biomarkers associated with BT. Three different BT levels were analyzed using reduced representation bisulfite sequencing data to calculate the methylation levels of cytosine and guanine dinucleotide (CpG) sites. The co-analysis of differentially methylated CpG sites identified by this study and differentially expressed genes identified by a previous study found 32 significant co-located genes. The joint analysis of GO terms and pathways revealed that methylation and gene expression of seven candidate genes were associated with BT; in particular, FASN plays a key role in fatty acid biosynthesis, and PEMT might be involved in estrogen regulation and the development of BT. This study is the first to report the genome-wide DNA methylation profiles of BT in pigs using next-generation sequencing and summarize candidate genes associated with epigenetic markers of BT, which could contribute to the understanding of the functional biology of BT traits and selective breeding of pigs against BT based on epigenetic biomarkers.