The epigenetic landscape of Alu repeats delineates the structural and functional genomic architecture of colon cancer cells.

Cancer cells exhibit multiple epigenetic changes with prominent local DNA hypermethylation and widespread hypomethylation affecting large chromosomal domains. Epigenome studies often disregard the study of repeat elements owing to technical complexity and their undefined role in genome regulation. We have developed NSUMA (Next-generation Sequencing of UnMethylated Alu), a cost-effective approach allowing the unambiguous interrogation of DNA methylation in more than 130,000 individual Alu elements, the most abundant retrotransposon in the human genome. DNA methylation profiles of Alu repeats have been analyzed in colon cancers and normal tissues using NSUMA and whole-genome bisulfite sequencing. Normal cells show a low proportion of unmethylated Alu (1%-4%) that may increase up to 10-fold in cancer cells. In normal cells, unmethylated Alu elements tend to locate in the vicinity of functionally rich regions and display epigenetic features consistent with a direct impact on genome regulation. In cancer cells, Alu repeats are more resistant to hypomethylation than other retroelements. Genome segmentation based on high/low rates of Alu hypomethylation allows the identification of genomic compartments with differential genetic, epigenetic, and transcriptomic features. Alu hypomethylated regions show low transcriptional activity, late DNA replication, and its extent is associated with higher chromosomal instability. Our analysis demonstrates that Alu retroelements contribute to define the epigenetic landscape of normal and cancer cells and provides a unique resource on the epigenetic dynamics of a principal, but largely unexplored, component of the primate genome.

Chainy: an universal tool for standardized relative quantification in real-time PCR.

Summary: Chainy is a cross-platform web tool providing systematic pipelines and steady criteria to process real-time PCR data, including the calculation of efficiencies from raw data by kinetic methods, evaluation of the suitability of multiple references, standardized normalization using one or more references, and group-wise relative quantification statistical testing. We illustrate the utility of Chainy for differential expression and chromatin immunoprecipitation enrichment (ChIP-QPCR) analysis. Availability and Implementation: Chainy is open source and freely available at http://maplab.cat/chainy. Contact: imallona@igtp.cat. Supplementary information: Supplementary data are available at Bioinformatics online.

Truke, a web tool to check for and handle excel misidentified gene symbols.

BACKGROUND: Genomic datasets accompanying scientific publications show a surprisingly high rate of gene name corruption. This error is generated when files and tables are imported into Microsoft Excel and certain gene symbols are automatically converted into dates. RESULTS: We have developed Truke, a fexible Web tool to detect, tag and fix, if possible, such misconversions. Aside, Truke is language and regional locale-aware, providing file format customization (decimal symbol, field sepator, etc.) following user's preferences. CONCLUSIONS: Truke is a data format conversion tool with a unique corrupted gene symbol detection utility. Truke is freely available without registration at http://maplab.cat/truke .

regioneR: an R/Bioconductor package for the association analysis of genomic regions based on permutation tests.

MOTIVATION: Statistically assessing the relation between a set of genomic regions and other genomic features is a common challenging task in genomic and epigenomic analyses. Randomization based approaches implicitly take into account the complexity of the genome without the need of assuming an underlying statistical model. SUMMARY: regioneR is an R package that implements a permutation test framework specifically designed to work with genomic regions. In addition to the predefined randomization and evaluation strategies, regioneR is fully customizable allowing the use of custom strategies to adapt it to specific questions. Finally, it also implements a novel function to evaluate the local specificity of the detected association. AVAILABILITY AND IMPLEMENTATION: regioneR is an R package released under Artistic-2.0 License. The source code and documents are freely available through Bioconductor (http://www.bioconductor.org/packages/regioneR). CONTACT: rmalinverni@carrerasresearch.org.

Muscle cell identity requires Pax7-mediated lineage-specific DNA demethylation.

BACKGROUND: Skeletal muscle stem cells enable the formation, growth, maintenance, and regeneration of skeletal muscle throughout life. The regeneration process is compromised in several pathological conditions, and muscle progenitors derived from pluripotent stem cells have been suggested as a potential therapeutic source for tissue replacement. DNA methylation is an important epigenetic mechanism in the setting and maintenance of cellular identity, but its role in stem cell determination towards the myogenic lineage is unknown. Here we addressed the DNA methylation dynamics of the major genes orchestrating the myogenic determination and differentiation programs in embryonic stem (ES) cells, their Pax7-induced myogenic derivatives, and muscle stem cells in proliferating and differentiating conditions. RESULTS: Our data showed a common muscle-specific DNA demethylation signature required to acquire and maintain the muscle-cell identity. This specific-DNA demethylation is Pax7-mediated, and it is a prime event in muscle stem cells gene activation. Notably, downregulation of the demethylation-related enzyme Apobec2 in ES-derived myogenic precursors reduced myogenin-associated DNA demethylation and dramatically impaired the expression of differentiation markers and, ultimately, muscle differentiation. CONCLUSIONS: Our results underscore DNA demethylation as a key mechanism driving myogenesis and identify specific Pax7-mediated DNA demethylation signatures to acquire and maintain the muscle-cell identity. Additionally, we provide a panel of epigenetic markers for the efficient and safe generation of ES- and induced pluripotent stem cell (iPS)-derived myogenic progenitors for therapeutic applications.

Deconstruction of DNA methylation patterns during myogenesis reveals specific epigenetic events in the establishment of the skeletal muscle lineage.

The progressive restriction of differentiation potential from pluripotent embryonic stem cells (ESCs) to tissue-specific stem cells involves widespread epigenetic reprogramming, including modulation of DNA methylation patterns. Skeletal muscle stem cells are required for the growth, maintenance, and regeneration of skeletal muscle. To investigate the contribution of DNA methylation to the establishment of the myogenic program, we analyzed ESCs, skeletal muscle stem cells in proliferating (myoblasts) and differentiating conditions (myotubes), and mature myofibers. About 1.000 differentially methylated regions were identified during muscle-lineage determination and terminal differentiation, mainly located in gene bodies and intergenic regions. As a whole, myogenic stem cells showed a gain of DNA methylation, while muscle differentiation was accompanied by loss of DNA methylation in CpG-poor regions. Notably, the hypomethylated regions in myogenic stem cells were neighbored by enhancer-type chromatin, suggesting the involvement of DNA methylation in the regulation of cell-type specific enhancers. Interestingly, we demonstrated the hypomethylation of the muscle cell-identity Myf5 super-enhancer only in muscle cells. Furthermore, we observed that upstream stimulatory factor 1 binding to Myf5 super-enhancer occurs upon DNA demethylation in myogenic stem cells. Taken altogether, we characterized the unique DNA methylation signature of skeletal muscle stem cells and highlighted the importance of DNA methylation-mediated regulation of cell identity Myf5 super-enhancer during cellular differentiation.

A knowledgebase of the human Alu repetitive elements.

Alu elements are the most abundant retrotransposons in the human genome with more than one million copies. Alu repeats have been reported to participate in multiple processes related with genome regulation and compartmentalization. Moreover, they have been involved in the facilitation of pathological mutations in many diseases, including cancer. The contribution of Alus and other repeats in genomic regulation is often overlooked because their study poses technical and analytical challenges hardly attainable with conventional strategies. Here we propose the integration of ontology-based semantic methods to query a knowledgebase for the human Alus. The knowledgebase for the human Alus leverages Sequence (SO) and Gene Ontologies (GO) and is devoted to address functional and genetic information in the genomic context of the Alus. For each Alu element, the closest gene and transcript are stored, as well their functional annotation according to GO, the state of the chromatin and the transcription factors binding sites inside the Alu. The model uses Web Ontology Language (OWL) and Semantic Web Rule Language (SWRL). As a case of use and to illustrate the utility of the tool, we have evaluated the epigenetic states of Alu repeats associated with gene promoters according to their transcriptional activity. The ontology is easily extendable, offering a scaffold for the inclusion of new experimental data. The RDF/XML formalization is freely available at http://aluontology.sourceforge.net/.

Epigenetic remodeling in colorectal cancer results in coordinate gene suppression across an entire chromosome band.

We report a new mechanism in carcinogenesis involving coordinate long-range epigenetic gene silencing. Epigenetic silencing in cancer has always been envisaged as a local event silencing discrete genes. However, in this study of silencing in colorectal cancer, we found common repression of the entire 4-Mb band of chromosome 2q.14.2, associated with global methylation of histone H3 Lys9. DNA hypermethylation within the repressed genomic neighborhood was localized to three separate enriched CpG island 'suburbs', with the largest hypermethylated suburb spanning 1 Mb. These data change our understanding of epigenetic gene silencing in cancer cells: namely, epigenetic silencing can span large regions of the chromosome, and both DNA-methylated and neighboring unmethylated genes can be coordinately suppressed by global changes in histone modification. We propose that loss of gene expression can occur through long-range epigenetic silencing, with similar implications as loss of heterozygosity in cancer.

Evaluation of single CpG sites as proxies of CpG island methylation states at the genome scale.

Methylation of a CpG island is a faithful marker of silencing of its associated gene. Different approaches report the methylation status of a CpG island based on the determination of one or a few CpG sites by assuming the homogeneity of methylation along the element. This strategy is frequently applied in both locus-specific and genome-wide studies, but often without a validation of the representativeness of the interrogated CpG site compared with the whole element. We have evaluated the predictive informativeness of the HpaII sites located in CpG islands using data from high-resolution methylome maps, which offer the possibility to assess the methylation homogeneity of each CpG island and to determine the reporter accuracy of single sites as surrogate markers. An excellent correlation was observed between the HpaII and CpG island methylation levels (r > 0.93). At the qualitative level, the predictive sensitivity of HpaII was >95% with >92% specificity for methylated CpG islands and >90% sensitivity with >95% specificity for unmethylated CpG islands. This analysis provides a global validation framework for strategies based on the use of the methylation-sensitive HpaII restriction enzyme.

Changes in the pattern of DNA methylation associate with twin discordance in systemic lupus erythematosus.

Monozygotic (MZ) twins are partially concordant for most complex diseases, including autoimmune disorders. Whereas phenotypic concordance can be used to study heritability, discordance suggests the role of non-genetic factors. In autoimmune diseases, environmentally driven epigenetic changes are thought to contribute to their etiology. Here we report the first high-throughput and candidate sequence analyses of DNA methylation to investigate discordance for autoimmune disease in twins. We used a cohort of MZ twins discordant for three diseases whose clinical signs often overlap: systemic lupus erythematosus (SLE), rheumatoid arthritis, and dermatomyositis. Only MZ twins discordant for SLE featured widespread changes in the DNA methylation status of a significant number of genes. Gene ontology analysis revealed enrichment in categories associated with immune function. Individual analysis confirmed the existence of DNA methylation and expression changes in genes relevant to SLE pathogenesis. These changes occurred in parallel with a global decrease in the 5-methylcytosine content that was concomitantly accompanied with changes in DNA methylation and expression levels of ribosomal RNA genes, although no changes in repetitive sequences were found. Our findings not only identify potentially relevant DNA methylation markers for the clinical characterization of SLE patients but also support the notion that epigenetic changes may be critical in the clinical manifestations of autoimmune disease.

Identification of intergenerational epigenetic inheritance by whole genome DNA methylation analysis in trios.

Genome-wide association studies have identified thousands of loci associated with common diseases and traits. However, a large fraction of heritability remains unexplained. Epigenetic modifications, such as the observed in DNA methylation have been proposed as a mechanism of intergenerational inheritance. To investigate the potential contribution of DNA methylation to the missing heritability, we analysed the methylomes of four healthy trios (two parents and one offspring) using whole genome bisulphite sequencing. Of the 1.5 million CpGs (19%) with over 20% variability between parents in at least one family and compatible with a Mendelian inheritance pattern, only 3488 CpGs (0.2%) lacked correlation with any SNP in the genome, marking them as potential sites for intergenerational epigenetic inheritance. These markers were distributed genome-wide, with some preference to be located in promoters. They displayed a bimodal distribution, being either fully methylated or unmethylated, and were often found at the boundaries of genomic regions with high/low GC content. This analysis provides a starting point for future investigations into the missing heritability of simple and complex traits.

Mycobacterium manresensis induces trained immunity in vitro.

The COVID-19 pandemic posed a global health crisis, with new severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants weakening vaccine-driven protection. Trained immunity could help tackle COVID-19 disease. Our objective was to analyze whether heat-killed Mycobacterium manresensis (hkMm), an environmental mycobacterium, induces trained immunity and confers protection against SARS-CoV-2 infection. To this end, THP-1 cells and primary monocytes were trained with hkMm. The increased secretion of tumor necrosis factor alpha (TNF-α), interleukin (IL)-6, IL-1ß, and IL-10, metabolic activity, and changes in epigenetic marks suggested hkMm-induced trained immunity in vitro. Healthcare workers at risk of SARS-CoV-2 infection were enrolled into the MANRECOVID19 clinical trial (NCT04452773) and were administered Nyaditum resae (NR, containing hkMm) or placebo. No significant differences in monocyte inflammatory responses or the incidence of SARS-CoV-2 infection were found between the groups, although NR modified the profile of circulating immune cell populations. Our results show that M. manresensis induces trained immunity in vitro but not in vivo when orally administered as NR daily for 14 days.

Colorectal Cancer Is Associated with the Presence of Cancer Driver Mutations in Normal Colon.

Although somatic mutations in colorectal cancer are well characterized, little is known about the accumulation of cancer mutations in the normal colon before cancer. Here, we have developed and applied an ultrasensitive, single-molecule mutational test based on CRISPR-DS technology, which enables mutation detection at extremely low frequency (<0.001) in normal colon from patients with and without colorectal cancer. This testing platform revealed that normal colon from patients with and without colorectal cancer carries mutations in common colorectal cancer genes, but these mutations are more abundant in patients with cancer. Oncogenic KRAS mutations were observed in the normal colon of about one third of patients with colorectal cancer but in none of the patients without colorectal cancer. Patients with colorectal cancer also carried more TP53 mutations than patients without cancer and these mutations were more pathogenic and formed larger clones, especially in patients with early-onset colorectal cancer. Most mutations in the normal colon were different from the driver mutations in tumors, suggesting that the occurrence of independent clones with pathogenic KRAS and TP53 mutations is a common event in the colon of individuals who develop colorectal cancer. These results indicate that somatic evolution contributes to clonal expansions in the normal colon and that this process is enhanced in individuals with cancer, particularly in those with early-onset colorectal cancer. SIGNIFICANCE: This work suggests prevalent somatic evolution in the normal colon of patients with colorectal cancer, highlighting the potential of using ultrasensitive gene sequencing to predict disease risk.

Bivalent domains enforce transcriptional memory of DNA methylated genes in cancer cells.

Silencing of multiple cancer-related genes is associated with de novo methylation of linked CpG islands. Additionally, bivalent histone modification profiles characterized by the juxtaposition of active and inactive histone marks have been observed in genes that become hypermethylated in cancer. It is unknown how these ambiguous epigenetic states are maintained and how they interrelate with adjacent genomic regions with different epigenetic landscapes. Here, we present the analysis of a set of neighboring genes, including many frequently silenced in colon cancer cells, in a chromosomal region at 5q35.2 spanning 1.25 Mb. Promoter DNA methylation occurs only at genes maintained at a low transcriptional state and is characterized by the presence of bivalent histone marks, namely trimethylation of lysines 4 and 27 in histone 3. Chemically induced hyperacetylation and DNA demethylation lead to up-regulation of silenced genes in this locus yet do not resolve bivalent domains into a domain-wide active chromatin conformation. In contrast, active genes in the region become down-regulated after drug treatment, accompanied by a partial loss of chromatin domain boundaries and spreading of the inactive histone mark trimethylated lysine 27 in histone 3. Our results demonstrate that bivalent domains mark the promoters of genes that will become DNA methylated in adult tumor cells to enforce transcriptional silence. These bivalent domains not only remain upon drug induced gene reactivation, but also spread over adjacent CpG islands. These results may have important implications in understanding and managing epigenetic therapies of cancer.

HDAC11 is a novel regulator of fatty acid oxidative metabolism in skeletal muscle.

Skeletal muscle is the largest tissue in mammalian organisms and is a key determinant of basal metabolic rate and whole-body energy metabolism. Histone deacetylase 11 (HDAC11) is the only member of the class IV subfamily of HDACs, and it is highly expressed in skeletal muscle, but its role in skeletal muscle physiology has never been investigated. Here, we describe for the first time the consequences of HDAC11 genetic deficiency in skeletal muscle, which results in the improvement of muscle function enhancing fatigue resistance and muscle strength. Loss of HDAC11 had no obvious impact on skeletal muscle structure but increased the number of oxidative myofibers by promoting a glycolytic-to-oxidative muscle fiber switch. Unexpectedly, HDAC11 was localized in muscle mitochondria and its deficiency enhanced mitochondrial content. In particular, we showed that HDAC11 depletion increased mitochondrial fatty acid ß-oxidation through activating the AMP-activated protein kinase-acetyl-CoA carboxylase pathway and reducing acylcarnitine levels in vivo, thus providing a mechanistic explanation for the improved muscle strength and fatigue resistance. Overall, our data reveal a unique role of HDAC11 in the maintenance of muscle fiber-type balance and the mitochondrial lipid oxidation. These findings shed light on the mechanisms governing muscle metabolism and may have implications for chronic muscle metabolic disease management.

Loss of HDAC11 accelerates skeletal muscle regeneration in mice.

Histone deacetylase 11 (HDAC11) is the latest identified member of the histone deacetylase family of enzymes. It is highly expressed in brain, heart, testis, kidney, and skeletal muscle, although its role in these tissues is poorly understood. Here, we investigate for the first time the consequences of HDAC11 genetic impairment on skeletal muscle regeneration, a process principally dependent on its resident stem cells (satellite cells) in coordination with infiltrating immune cells and stromal cells. Our results show that HDAC11 is dispensable for adult muscle growth and establishment of the satellite cell population, while HDAC11 deficiency advances the regeneration process in response to muscle injury. This effect is not caused by differences in satellite cell activation or proliferation upon injury, but rather by an enhanced capacity of satellite cells to differentiate at early regeneration stages in the absence of HDAC11. Infiltrating HDAC11-deficient macrophages could also contribute to this accelerated muscle regenerative process by prematurely producing high levels of IL-10, a cytokine known to promote myoblast differentiation. Altogether, our results show that HDAC11 depletion advances skeletal muscle regeneration and this finding may have potential implications for designing new strategies for muscle pathologies coursing with chronic damage. DATABASE: Data were deposited in NCBI's Gene Expression Omnibus accessible through GEO Series accession number GSE147423.

Methods for DNA methylation analysis and applications in colon cancer.

It is well established that epigenetic events, in an intimate cooperation with genetic events, are involved in every step of tumorigenesis. DNA methylation, which in mammals takes place in the cytosines that precede a guanine (CpG dinucleotide), is the most well-characterized epigenetic mark. The study of aberrant DNA methylation patterns, such as hypermethylation of CpG islands and global genomic hypomethylation, are common issues in the studies on all types of cancer, and as in other areas of molecular oncology, colorectal cancer has become a privileged target. Besides the great variety of technologies available for the analysis of DNA methylation, most methods are based on three principles: methylation-sensitive enzymes, bisulphite conversion of unmethylated cytosines and immunoprecipitation of 5-methylcytosines. By combining each one of these principles with other genomic methodologies, a large range of approaches aimed at the analysis of methylation from one specific CpG site to a large number of sequences on the genome scale and suitable for different research needs have been developed. The goal of this review is to describe the most widely used methylation methods in the study of cancer, as well as the potential clinical applications of DNA methylation biomarkers in colorectal cancer.

Genome-wide tracking of unmethylated DNA Alu repeats in normal and cancer cells.

Methylation of the cytosine is the most frequent epigenetic modification of DNA in mammalian cells. In humans, most of the methylated cytosines are found in CpG-rich sequences within tandem and interspersed repeats that make up to 45% of the human genome, being Alu repeats the most common family. Demethylation of Alu elements occurs in aging and cancer processes and has been associated with gene reactivation and genomic instability. By targeting the unmethylated SmaI site within the Alu sequence as a surrogate marker, we have quantified and identified unmethylated Alu elements on the genomic scale. Normal colon epithelial cells contain in average 25 486 +/- 10 157 unmethylated Alu's per haploid genome, while in tumor cells this figure is 41 995 +/- 17 187 (P = 0.004). There is an inverse relationship in Alu families with respect to their age and methylation status: the youngest elements exhibit the highest prevalence of the SmaI site (AluY: 42%; AluS: 18%, AluJ: 5%) but the lower rates of unmethylation (AluY: 1.65%; AluS: 3.1%, AluJ: 12%). Data are consistent with a stronger silencing pressure on the youngest repetitive elements, which are closer to genes. Further insights into the functional implications of atypical unmethylation states in Alu elements will surely contribute to decipher genomic organization and gene regulation in complex organisms.

Dihydrofolate reductase amplification and sensitization to methotrexate of methotrexate-resistant colon cancer cells.

Gene amplification is one of the most frequent manifestations of genomic instability in human tumors and plays an important role in tumor progression and acquisition of drug resistance. To better understand the factors involved in acquired resistance to cytotoxic drugs via gene amplification, we have analyzed the structure and dynamics of dihydrofolate reductase (DHFR) gene amplification in HT29 cells treated with methotrexate (MTX). Analysis of the DHFR gene amplification process shows that the amplicon exhibits a complex structure that is consistently reproduced in independent treatments. The cytogenetic manifestation of the amplification in advanced stages of the treatment may be in the form of double minutes or as a homogeneously stained region. To get insights into the mechanisms of resistance, we have also investigated the sensitization to MTX of MTX-resistant cells after drug withdrawal and reexposure to MTX. Passive loss of the DHFR amplicon by withdrawal of the drug results in MTX-sensitive cells exhibiting a substantial reduction of their capacity or even an incapacity to generate resistance when submitted to a second cycle of MTX treatment. On a second round of drug administration, the resistant cells generate a different amplicon structure, suggesting that the formation of the amplicon as in the first cycle of treatment is not feasible. These results indicate that DHFR gene amplification is a "wear and tear" process in HT29 cells and that MTX-resistant cells may become responsive to a second round of treatment if left untreated during a sufficient period of time.