Enrichment analysis of Alu elements with different spatial chromatin proximity in the human genome

Evolutionary patterns and functional effects of 3D chromatin structures in butterflies with extensive genome rearrangements

Chromosome rearrangements may distort 3D chromatin architectures and thus change gene regulation, yet how 3D chromatin structures evolve in insects is largely unknown. Here, we obtain chromosome-level genomes for four butterfly species, Graphium cloanthus, Graphium sarpedon, Graphium eurypylus with 2n = 30, 40, and 60, respectively, and Papilio bianor with 2n = 60. Together with large-scale Hi-C data, we find that inter-chromosome rearrangements very rarely disrupted the pre-existing 3D chromatin structure of ancestral chromosomes. However, some intra-chromosome rearrangements changed 3D chromatin structures compared to the ancestral configuration. We find that new TADs and subTADs have emerged across the rearrangement sites where their adjacent compartments exhibit uniform types. Two intra-chromosome rearrangements altered Rel and lft regulation, potentially contributing to wing patterning differentiation and host plant choice. Notably, butterflies exhibited chromatin loops between Hox gene cluster ANT-C and BX-C, unlike Drosophila. Our CRISPR-Cas9 experiments in butterflies confirm that knocking out the CTCF binding site of the loops in BX-C affected the phenotypes regulated by Antp in ANT-C, resulting in legless larva. Our results reveal evolutionary patterns of insect 3D chromatin structures and provide evidence that 3D chromatin structure changes can play important roles in the evolution of traits.

Spatial-Temporal Genome Regulation in Stress-Response and Cell-Fate Change

Complex functioning of the genome in the cell nucleus is controlled at different levels: (a) the DNA base sequence containing all relevant inherited information; (b) epigenetic pathways consisting of protein interactions and feedback loops; (c) the genome architecture and organization activating or suppressing genetic interactions between different parts of the genome. Most research so far has shed light on the puzzle pieces at these levels. This article, however, attempts an integrative approach to genome expression regulation incorporating these different layers. Under environmental stress or during cell development, differentiation towards specialized cell types, or to dysfunctional tumor, the cell nucleus seems to react as a whole through coordinated changes at all levels of control. This implies the need for a framework in which biological, chemical, and physical manifestations can serve as a basis for a coherent theory of gene self-organization. An international symposium held at the Biomedical Research and Study Center in Riga, Latvia, on 25 July 2022 addressed novel aspects of the abovementioned topic. The present article reviews the most recent results and conclusions of the state-of-the-art research in this multidisciplinary field of science, which were delivered and discussed by scholars at the Riga symposium.

The single-cell expression profile of transposable elements and transcription factors in human early biparental and uniparental embryonic development

Transposable elements (TEs) and transcription factors (TFs) are involved in the precise regulation of gene expression during the preimplantation stage. Activation of TEs is a key event for mammalian embryonic genome activation and preimplantation early embryonic development. TFs are involved in the regulation of drastic changes in gene expression patterns, but an inventory of the interplay between TEs and TFs during normal/abnormal human embryonic development is still lacking. Here we used single-cell RNA sequencing data generated from biparental and uniparental embryos to perform an integrative analysis of TE and TF expression. Our results showed that endogenous retroviruses (ERVs) are mainly expressed during the minor embryonic genome activation (EGA) process of early embryos, while Alu is gradually expressed in the middle and later stages. Some important ERVs (e.g., LTR5_Hs, MLT2A1) and Alu TEs are expressed at significantly lower levels in androgenic embryos. Integrative analysis revealed that the expression of the transcription factors CTCF and POU5F1 is correlated with the differential expression of ERV TEs. Comparative coexpression network analysis further showed distinct expression levels of important TFs (e.g., LEUTX and ZSCAN5A) in dizygotic embryos vs. parthenogenetic and androgenic embryos. This systematic investigation of TE and TF expression in human early embryonic development by single-cell RNA sequencing provides valuable insights into mammalian embryonic development.

Comparative 3D genome architecture in vertebrates

BACKGROUND

The three-dimensional (3D) architecture of the genome has a highly ordered and hierarchical nature, which influences the regulation of essential nuclear processes at the basis of gene expression, such as gene transcription. While the hierarchical organization of heterochromatin and euchromatin can underlie differences in gene expression that determine evolutionary differences among species, the way 3D genome architecture is affected by evolutionary forces within major lineages remains unclear. Here, we report a comprehensive comparison of 3D genomes, using high resolution Hi-C data in fibroblast cells of fish, chickens, and 10 mammalian species.

RESULTS

This analysis shows a correlation between genome size and chromosome length that affects chromosome territory (CT) organization in the upper hierarchy of genome architecture, whereas lower hierarchical features, including local transcriptional availability of DNA, are selected through the evolution of vertebrates. Furthermore, conservation of topologically associating domains (TADs) appears strongly associated with the modularity of expression profiles across species. Additionally, LINE and SINE transposable elements likely contribute to heterochromatin and euchromatin organization, respectively, during the evolution of genome architecture.

CONCLUSIONS

Our analysis uncovers organizational features that appear to determine the conservation and transcriptional regulation of functional genes across species. These findings can guide ongoing investigations of genome evolution by extending our understanding of the mechanisms shaping genome architecture.

A hypothetical model of trans-acting R-loops-mediated promoter-enhancer interactions by Alu elements

Enhancers modulate gene expression by interacting with promoters. Models of enhancer-promoter interactions (EPIs) in the literature involve the activity of many components, including transcription factors and nucleic acid. However, the role that sequence similarity plays in EPIs remains largely unexplored. Herein, we report that Alu-derived sequences dominate sequence similarity between enhancers and promoters. After rejecting alternative DNA:DNA and DNA:RNA triplex models, we propose that enhancer-associated RNAs (eRNAs) may directly contact their targeted promoters by forming trans-acting R-loops at those Alu sequences. We show how the characteristic distribution of functional genomic data, such as RNA-DNA proximate ligation reads, binding of transcription factors, and RNA-binding proteins, all align with the Alu sequences of EPIs. We also show that these aligned Alu sequences may be subject to the constraint of coevolution, further implying the functional significance of these R-loop hybrids. Finally, our results imply that eRNA and Alu elements associate in a manner previously unrecognized in EPIs and the evolution of gene regulation networks in mammals.

K-mer Content Changes with Node Degree in Promoter-Enhancer Network of Mouse ES Cells

Maps of Hi-C contacts between promoters and enhancers can be analyzed as networks, with cis-regulatory regions as nodes and their interactions as edges. We checked if in the published promoter-enhancer network of mouse embryonic stem (ES) cells the differences in the node type (promoter or enhancer) and the node degree (number of regions interacting with a given promoter or enhancer) are reflected by sequence composition or sequence similarity of the interacting nodes. We used counts of all k-mers (k = 4) to analyze the sequence composition and the Euclidean distance between the k-mer count vectors (k-mer distance) as the measure of sequence (dis)similarity. The results we obtained with 4-mers are interpretable in terms of dinucleotides. Promoters are GC-rich as compared to enhancers, which is known. Enhancers are enriched in scaffold/matrix attachment regions (S/MARs) patterns and depleted of CpGs. Furthermore, we show that promoters are more similar to their interacting enhancers than vice-versa. Most notably, in both promoters and enhancers, the GC content and the CpG count increase with the node degree. As a consequence, enhancers of higher node degree become more similar to promoters, whereas higher degree promoters become less similar to enhancers. We confirmed the key results also for human keratinocytes.

Analyzing a putative enhancer of optic disc morphology

Background

Genome-wide association studies have identified the CDC7-TGFBR3 intergenic region on chromosome 1 to be strongly associated with optic disc area size. The mechanism of its function remained unclear until new data on eQTL markers emerged from the Genotype-Tissue Expression project. The target region was found to contain a strong silencer of the distal (800 kb) Transcription Factor (TF) gene GFI1 (Growth Factor Independent Transcription Repressor 1) specifically in neuroendocrine cells (pituitary gland). GFI1 has also been reported to be involved in the development of sensory neurons and hematopoiesis. Therefore, GFI1, being a developmental gene, is likely to affect optic disc area size by altering the expression of the associated genes via long-range interactions.

Results

Distribution of haplotypes in the putative enhancer region has been assessed using the data on four continental supergroups generated by the 1000 Genomes Project. The East Asian (EAS) populations were shown to manifest a highly homogenous unimodal haplotype distribution pattern within the region with the major haplotype occurring with the frequency of 0.9. Another European specific haplotype was observed with the frequency of 0.21. The major haplotype appears to be involved in silencing GFI1repressor gene expression, which might be the cause of increased optic disc area characteristic of the EAS populations.

The enhancer/eQTL region overlaps AluJo element, which implies that this particular regulatory element is primate-specific and confined to few tissues.

Conclusion

Population specific distribution of GFI1 enhancer alleles may predispose certain ethnic groups to glaucoma.

Architectural proteins for the formation and maintenance of the 3D genome

Eukaryotic genomes are densely packaged into hierarchical three-dimensional (3D) structures that contain information about gene regulation and many other biological processes. With the development of imaging and sequencing-based technologies, 3D genome studies have revealed that the high-order chromatin structure is composed of hierarchical levels, including chromosome territories, A/B compartments, topologically associated domains, and chromatin loops. However, how this chromatin architecture is formed and maintained is not completely clear. In this review, we introduce experimental methods to investigate the 3D genome, review major architectural proteins that regulate 3D chromatin organization in mammalian cells, such as CTCF (CCCTC-binding factor), cohesin, lamins, and transcription factors, and discuss relevant mechanisms such as phase separation.

TFIIIC Binding to Alu Elements Controls Gene Expression via Chromatin Looping and Histone Acetylation

How repetitive elements, epigenetic modifications, and architectural proteins interact ensuring proper genome expression remains poorly understood. Here, we report regulatory mechanisms unveiling a central role of Alu elements (AEs) and RNA polymerase III transcription factor C (TFIIIC) in structurally and functionally modulating the genome via chromatin looping and histone acetylation. Upon serum deprivation, a subset of AEs pre-marked by the activity-dependent neuroprotector homeobox Protein (ADNP) and located near cell-cycle genes recruits TFIIIC, which alters their chromatin accessibility by direct acetylation of histone H3 lysine-18 (H3K18). This facilitates the contacts of AEs with distant CTCF sites near promoter of other cell-cycle genes, which also become hyperacetylated at H3K18. These changes ensure basal transcription of cell-cycle genes and are critical for their re-activation upon serum re-exposure. Our study reveals how direct manipulation of the epigenetic state of AEs by a general transcription factor regulates 3D genome folding and expression.

Physicochemical Foundations of Life that Direct Evolution: Chance and Natural Selection are not Evolutionary Driving Forces

The current framework of evolutionary theory postulates that evolution relies on random mutations generating a diversity of phenotypes on which natural selection acts. This framework was established using a top-down approach as it originated from Darwinism, which is based on observations made of complex multicellular organisms and, then, modified to fit a DNA-centric view. In this article, it is argued that based on a bottom-up approach starting from the physicochemical properties of nucleic and amino acid polymers, we should reject the facts that (i) natural selection plays a dominant role in evolution and (ii) the probability of mutations is independent of the generated phenotype. It is shown that the adaptation of a phenotype to an environment does not correspond to organism fitness, but rather corresponds to maintaining the genome stability and integrity. In a stable environment, the phenotype maintains the stability of its originating genome and both (genome and phenotype) are reproduced identically. In an unstable environment (i.e., corresponding to variations in physicochemical parameters above a physiological range), the phenotype no longer maintains the stability of its originating genome, but instead influences its variations. Indeed, environment- and cellular-dependent physicochemical parameters define the probability of mutations in terms of frequency, nature, and location in a genome. Evolution is non-deterministic because it relies on probabilistic physicochemical rules, and evolution is driven by a bidirectional interplay between genome and phenotype in which the phenotype ensures the stability of its originating genome in a cellular and environmental physicochemical parameter-depending manner.

The role of environmental exposures and the epigenome in health and disease

The genetic material of every organism exists within the context of regulatory networks that govern gene expression, collectively called the epigenome. Epigenetics has taken center stage in the study of diseases such as cancer and diabetes, but its integration into the field of environmental health is still emerging. As the Environmental Mutagenesis and Genomics Society (EMGS) celebrates its 50th Anniversary this year, we have come together to review and summarize the seminal advances in the field of environmental epigenomics. Specifically, we focus on the role epigenetics may play in multigenerational and transgenerational transmission of environmentally induced health effects. We also summarize state of the art techniques for evaluating the epigenome, environmental epigenetic analysis, and the emerging field of epigenome editing. Finally, we evaluate transposon epigenetics as they relate to environmental exposures and explore the role of noncoding RNA as biomarkers of environmental exposures. Although the field has advanced over the past several decades, including being recognized by EMGS with its own Special Interest Group, recently renamed Epigenomics, we are excited about the opportunities for environmental epigenetic science in the next 50 years. Environ. Mol. Mutagen. 61:176-192, 2020. © 2019 Wiley Periodicals, Inc.

New insights on human essential genes based on integrated analysis and the construction of the HEGIAP web-based platform

Essential genes are those whose loss of function compromises organism viability or results in profound loss of fitness. Recent gene-editing technologies have provided new opportunities to characterize essential genes. Here, we present an integrated analysis that comprehensively and systematically elucidates the genetic and regulatory characteristics of human essential genes. First, we found that essential genes act as ‘hubs’ in protein–protein interaction networks, chromatin structure and epigenetic modification. Second, essential genes represent conserved biological processes across species, although gene essentiality changes differently among species. Third, essential genes are important for cell development due to their discriminate transcription activity in embryo development and oncogenesis. In addition, we developed an interactive web server, the Human Essential Genes Interactive Analysis Platform (http://sysomics.com/HEGIAP/), which integrates abundant analytical tools to enable global, multidimensional interpretation of gene essentiality. Our study provides new insights that improve the understanding of human essential genes.

The modular organization of the cerebral cortex: Evolutionary significance and possible links to neurodevelopmental conditions

The recognition of discernible anatomical regularities that appear to self-organize during development makes apparent the modular organization of the cerebral cortex. The metabolic cost engendered in sustaining interneuronal communications has emphasized the viability of short connections among neighboring neurons. This pattern of connectivity establishes a microcircuit which is repeated in parallel throughout the cerebral cortex. This canonical circuit is contained within the smallest module of information processing of the cerebral cortex; one which Vernon Mountcastle called the minicolumn. Plasticity within the brain is accounted, in part, by the presence of weak linkages that allow minicolumns to process information from a variety of sources and to quickly adapt to environmental exigencies without a need for genetic change. Recent research suggests that interlaminar correlated firing between minicolumns during the decision phase of target selection provides for the emergence of some executive functions. Bottlenecks of information processing within this modular minicolumnar organization may account for a variety of mental disorders observed in neurodevelopmental conditions.

The Transcription Factor ThPOK Orchestrates Stochastic Interchromosomal Interactions Required for IFNB1 Virus-Inducible Gene Expression

Virus infection induces stochastic activation of the interferon-β gene. Three previously identified Alu-like DNA elements called NRCs (NF-κB reception centers) function by capturing and delivering NF-κB to the IFNB1 enhancer via stochastic interchromosomal interactions. We show that the transcription factor ThPOK binds cooperatively with NF-κB to NRCs and mediates their physical proximity with the IFNB1 gene via its ability to oligomerize when bound to DNA. ThPOK knockdown significantly decreased the frequency of interchromosomal interactions, NF-κB DNA binding to the IFNB1 enhancer, and virus-induced IFNB1 gene activation. We also demonstrate that cooperative DNA binding between ThPOK and NF-κB on the same face of the double DNA helix is required for interchromosomal interactions and distinguishes NRCs from various other Alu elements bearing κB sites. These studies show how DNA binding cooperativity of stereospecifically aligned transcription factors provides the necessary ultrasensitivity for the all-or-none stochastic cell responses to virus infection.

ALUminating the Path of Atherosclerosis Progression: Chaos Theory Suggests a Role for Alu Repeats in the Development of Atherosclerotic Vascular Disease

Atherosclerosis (ATH) and coronary artery disease (CAD) are chronic inflammatory diseases with an important genetic background; they derive from the cumulative effect of multiple common risk alleles, most of which are located in genomic noncoding regions. These complex diseases behave as nonlinear dynamical systems that show a high dependence on their initial conditions; thus, long-term predictions of disease progression are unreliable. One likely possibility is that the nonlinear nature of ATH could be dependent on nonlinear correlations in the structure of the human genome. In this review, we show how chaos theory analysis has highlighted genomic regions that have shared specific structural constraints, which could have a role in ATH progression. These regions were shown to be enriched with repetitive sequences of the Alu family, genomic parasites that have colonized the human genome, which show a particular secondary structure and are involved in the regulation of gene expression. Here, we show the impact of Alu elements on the mechanisms that regulate gene expression, especially highlighting the molecular mechanisms via which the Alu elements alter the inflammatory response. We devote special attention to their relationship with the long noncoding RNA (lncRNA); antisense noncoding RNA in the INK4 locus (ANRIL), a risk factor for ATH; their role as microRNA (miRNA) sponges; and their ability to interfere with the regulatory circuitry of the (nuclear factor kappa B) NF-κB response. We aim to characterize ATH as a nonlinear dynamic system, in which small initial alterations in the expression of a number of repetitive elements are somehow amplified to reach phenotypic significance.

Differential hypomethylation of the repetitive Tol2/Alu-rich sequences in the genome of Bodianus species (Labriformes, Labridae)

Representatives of the order Labriformes show karyotypes of extreme conservatism together with others with high chromosomal diversification. However, the cytological characterization of epigenetic modifications remains unknown for the majority of the species. In the family Labridae, the most abundant fishes on tropical reefs, the genomes of the genus Bodianus Bloch, 1790 have been characterized by the occurrence of a peculiar chromosomal region, here denominated BOD. This region is exceptionally decondensed, heterochromatic, argentophilic, GC-neutral and, in contrast to classical secondary constrictions, shows no signals of hybridization with 18S rDNA probes. In order to characterize the BOD region, the methylation pattern, the distribution of Alu and Tol2 retrotransposons and of 18S and 5S rDNA sites, respectively, were analyzed by Fluorescence In Situ Hybridization (FISH) on metaphase chromosomes of two Bodianus species, B. insularis Gomon & Lubbock, 1980 and B. pulchellus (Poey, 1860). Immunolocalization of the 5-methylcytosine revealed hypermethylated chromosomal regions, dispersed along the entire length of the chromosomes of both species, while the BOD regions exhibited a hypomethylated pattern. Hypomethylation of the BOD region is associated with the precise co-location of Tol2 and Alu elements, suggesting their active participation in the regulatory epigenetic process. This evidence underscores a probable differential methylation action during the cell cycle, as well as the role of Tol2/Alu elements in functional processes of fish genomes.

Challenges for Super-Resolution Localization Microscopy and Biomolecular Fluorescent Nano-Probing in Cancer Research

Understanding molecular interactions and regulatory mechanisms in tumor initiation, progression, and treatment response are key requirements towards advanced cancer diagnosis and novel treatment procedures in personalized medicine. Beyond decoding the gene expression, malfunctioning and cancer-related epigenetic pathways, investigations of the spatial receptor arrangements in membranes and genome organization in cell nuclei, on the nano-scale, contribute to elucidating complex molecular mechanisms in cells and tissues. By these means, the correlation between cell function and spatial organization of molecules or molecular complexes can be studied, with respect to carcinogenesis, tumor sensitivity or tumor resistance to anticancer therapies, like radiation or antibody treatment. Here, we present several new applications for bio-molecular nano-probes and super-resolution, laser fluorescence localization microscopy and their potential in life sciences, especially in biomedical and cancer research. By means of a tool-box of fluorescent antibodies, green fluorescent protein (GFP) tagging, or specific oligonucleotides, we present tumor relevant re-arrangements of Erb-receptors in membranes, spatial organization of Smad specific ubiquitin protein ligase 2 (Smurf2) in the cytosol, tumor cell characteristic heterochromatin organization, and molecular re-arrangements induced by radiation or antibody treatment. The main purpose of this article is to demonstrate how nano-scaled distance measurements between bio-molecules, tagged by appropriate nano-probes, can be applied to elucidate structures and conformations of molecular complexes which are characteristic of tumorigenesis and treatment responses. These applications open new avenues towards a better interpretation of the spatial organization and treatment responses of functionally relevant molecules, at the single cell level, in normal and cancer cells, offering new potentials for individualized medicine.

Discovery of rare, diagnostic AluYb8/9 elements in diverse human populations

BACKGROUND

Polymorphic human Alu elements are excellent tools for assessing population structure, and new retrotransposition events can contribute to disease. Next-generation sequencing has greatly increased the potential to discover Alu elements in human populations, and various sequencing and bioinformatics methods have been designed to tackle the problem of detecting these highly repetitive elements. However, current techniques for Alu discovery may miss rare, polymorphic Alu elements. Combining multiple discovery approaches may provide a better profile of the polymorphic Alu mobilome. AluYb8/9 elements have been a focus of our recent studies as they are young subfamilies (~2.3 million years old) that contribute ~30% of recent polymorphic Alu retrotransposition events. Here, we update our ME-Scan methods for detecting Alu elements and apply these methods to discover new insertions in a large set of individuals with diverse ancestral backgrounds.

RESULTS

We identified 5,288 putative Alu insertion events, including several hundred novel AluYb8/9 elements from 213 individuals from 18 diverse human populations. Hundreds of these loci were specific to continental populations, and 23 non-reference population-specific loci were validated by PCR. We provide high-quality sequence information for 68 rare AluYb8/9 elements, of which 11 have hallmarks of an active source element. Our subfamily distribution of rare AluYb8/9 elements is consistent with previous datasets, and may be representative of rare loci. We also find that while ME-Scan and low-coverage, whole-genome sequencing (WGS) detect different Alu elements in 41 1000 Genomes individuals, the two methods yield similar population structure results.

CONCLUSION

Current in-silico methods for Alu discovery may miss rare, polymorphic Alu elements. Therefore, using multiple techniques can provide a more accurate profile of Alu elements in individuals and populations. We improved our false-negative rate as an indicator of sample quality for future ME-Scan experiments. In conclusion, we demonstrate that ME-Scan is a good supplement for next-generation sequencing methods and is well-suited for population-level analyses.

Combining Low Temperature Fluorescence DNA-Hybridization, Immunostaining, and Super-Resolution Localization Microscopy for Nano-Structure Analysis of ALU Elements and Their Influence on Chromatin Structure

Immunostaining and fluorescence in situ hybridization (FISH) are well established methods for specific labelling of chromatin in the cell nucleus. COMBO-FISH (combinatorial oligonucleotide fluorescence in situ hybridization) is a FISH method using computer designed oligonucleotide probes specifically co-localizing at given target sites. In combination with super resolution microscopy which achieves spatial resolution far beyond the Abbe Limit, it allows new insights into the nano-scaled structure and organization of the chromatin of the nucleus. To avoid nano-structural changes of the chromatin, the COMBO-FISH labelling protocol was optimized omitting heat treatment for denaturation of the target. As an example, this protocol was applied to ALU elements—dispersed short stretches of DNA which appear in different kinds in large numbers in primate genomes. These ALU elements seem to be involved in gene regulation, genomic diversity, disease induction, DNA repair, etc. By computer search, we developed a unique COMBO-FISH probe which specifically binds to ALU consensus elements and combined this DNA–DNA labelling procedure with heterochromatin immunostainings in formaldehyde-fixed cell specimens. By localization microscopy, the chromatin network-like arrangements of ALU oligonucleotide repeats and heterochromatin antibody labelling sites were simultaneously visualized and quantified. This novel approach which simultaneously combines COMBO-FISH and immunostaining was applied to chromatin analysis on the nanoscale after low-linear-energy-transfer (LET) radiation exposure at different doses. Dose-correlated curves were obtained from the amount of ALU representing signals, and the chromatin re-arrangements during DNA repair after irradiation were quantitatively studied on the nano-scale. Beyond applications in radiation research, the labelling strategy of immunostaining and COMBO-FISH with localization microscopy will also offer new potentials for analyses of subcellular elements in combination with other specific chromatin targets.

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.