The 3D folding of metazoan genomes correlates with the association of similar repetitive elements

Development and evolution of Drosophila chromatin landscape in a 3D genome context

Little is known about how the epigenomic states change during development and evolution in a 3D genome context. Here we use Drosophila pseudoobscura with complex turnover of sex chromosomes as a model to address this, by collecting massive epigenomic and Hi-C data from five developmental stages and three adult tissues. We reveal that over 60% of the genes and transposable elements (TE) exhibit at least one developmental transition of chromatin state. Transitions on specific but not housekeeping enhancers are associated with specific chromatin loops and topologically associated domain borders (TABs). While evolutionarily young TEs are generally silenced, old TEs more often have been domesticated as interacting TABs or specific enhancers. But on the recently evolved X chromosome, young TEs are instead often active and recruited as TABs, due to acquisition of dosage compensation. Overall we characterize how Drosophila epigenomic landscapes change during development and in response to chromosome evolution, and highlight the important roles of TEs in genome organization and regulation.

Loops are geometric catalysts for DNA integration

The insertion of DNA elements within genomes underpins both genetic diversity and disease when unregulated. Most of DNA insertions are not random and the physical mechanisms underlying the integration site selection are poorly understood. Here, we perform Molecular Dynamics simulations to study the insertion of DNA elements, such as viral DNA or transposons, into naked DNA or chromatin substrates. More specifically, we explore the role of loops within the polymeric substrate and discover that they act as 'geometric catalysts' for DNA integration by reducing the energy barrier for substrate deformation. Additionally, we discover that the 1D pattern and 3D conformation of loops have a marked effect on the distribution of integration sites. Finally, we show that loops may compete with nucleosomes to attract DNA integrations. These results may be tested in vitro and they may help to understand patterns of DNA insertions with implications in genome evolution and engineering.

The role of transposable elements in human evolution and methods for their functional analysis: current status and future perspectives

Transposable elements (TEs) are mobile DNA sequences that can insert themselves into various locations within the genome, causing mutations that may provide advantages or disadvantages to individuals and species. The insertion of TEs can result in genetic variation that may affect a wide range of human traits including genetic disorders. Understanding the role of TEs in human biology is crucial for both evolutionary and medical research. This review discusses the involvement of TEs in human traits and disease susceptibility, as well as methods for functional analysis of TEs.

ZNF689 deficiency promotes intratumor heterogeneity and immunotherapy resistance in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is an aggressive disease characterized by remarkable intratumor heterogeneity (ITH), which poses therapeutic challenges. However, the clinical relevance and key determinant of ITH in TNBC are poorly understood. Here, we comprehensively characterized ITH levels using multi-omics data across our center's cohort (n = 260), The Cancer Genome Atlas cohort (n = 134), and four immunotherapy-treated cohorts (n = 109). Our results revealed that high ITH was associated with poor patient survival and immunotherapy resistance. Importantly, we identified zinc finger protein 689 (ZNF689) deficiency as a crucial determinant of ITH formation. Mechanistically, the ZNF689-TRIM28 complex was found to directly bind to the promoter of long interspersed element-1 (LINE-1), inducing H3K9me3-mediated transcriptional silencing. ZNF689 deficiency reactivated LINE-1 retrotransposition to exacerbate genomic instability, which fostered ITH. Single-cell RNA sequencing, spatially resolved transcriptomics and flow cytometry analysis confirmed that ZNF689 deficiency-induced ITH inhibited antigen presentation and T-cell activation, conferring immunotherapy resistance. Pharmacological inhibition of LINE-1 significantly reduced ITH, enhanced antitumor immunity, and eventually sensitized ZNF689-deficient tumors to immunotherapy in vivo. Consistently, ZNF689 expression positively correlated with favorable prognosis and immunotherapy response in clinical samples. Altogether, our study uncovers a previously unrecognized mechanism underlying ZNF689 deficiency-induced ITH and suggests LINE-1 inhibition combined with immunotherapy as a novel treatment strategy for TNBC.

Transposable elements in mammalian chromatin organization

Transposable elements (TEs) are mobile DNA elements that comprise almost 50% of mammalian genomic sequence. TEs are capable of making additional copies of themselves that integrate into new positions in host genomes. This unique property has had an important impact on mammalian genome evolution and on the regulation of gene expression because TE-derived sequences can function as cis-regulatory elements such as enhancers, promoters and silencers. Now, advances in our ability to identify and characterize TEs have revealed that TE-derived sequences also regulate gene expression by both maintaining and shaping 3D genome architecture. Studies are revealing how TEs contribute raw sequence that can give rise to the structures that shape chromatin organization, and thus gene expression, allowing for species-specific genome innovation and evolutionary novelty.

The Conserved Chromatin Remodeler SMARCAD1 Interacts with TFIIIC and Architectural Proteins in Human and Mouse

In vertebrates, SMARCAD1 participates in transcriptional regulation, heterochromatin maintenance, DNA repair, and replication. The molecular basis underlying its involvement in these processes is not well understood. We identified the RNA polymerase III general transcription factor TFIIIC as an interaction partner of native SMARCAD1 in mouse and human models using endogenous co-immunoprecipitations. TFIIIC has dual functionality, acting as a general transcription factor and as a genome organizer separating chromatin domains. We found that its partnership with SMARCAD1 is conserved across different mammalian cell types, from somatic to pluripotent cells. Using purified proteins, we confirmed that their interaction is direct. A gene expression analysis suggested that SMARCAD1 is dispensable for TFIIIC function as an RNA polymerase III transcription factor in mouse ESCs. The distribution of TFIIIC and SMARCAD1 in the ESC genome is distinct, and unlike in yeast, SMARCAD1 is not enriched at active tRNA genes. Further analysis of SMARCAD1-binding partners in pluripotent and differentiated mammalian cells reveals that SMARCAD1 associates with several factors that have key regulatory roles in chromatin organization, such as cohesin, laminB, and DDX5. Together, our work suggests for the first time that the SMARCAD1 enzyme participates in genome organization in mammalian nuclei through interactions with architectural proteins.

Young SINEs in pig genomes impact gene regulation, genetic diversity, and complex traits

Transposable elements (TEs) are a major source of genetic polymorphisms and play a role in chromatin architecture, gene regulatory networks, and genomic evolution. However, their functional role in pigs and contributions to complex traits are largely unknown. We created a catalog of TEs (n = 3,087,929) in pigs and found that young SINEs were predominantly silenced by histone modifications, DNA methylation, and decreased accessibility. However, some transcripts from active young SINEs showed high tissue-specificity, as confirmed by analyzing 3570 RNA-seq samples. We also detected 211,067 dimorphic SINEs in 374 individuals, including 340 population-specific ones associated with local adaptation. Mapping these dimorphic SINEs to genome-wide associations of 97 complex traits in pigs, we found 54 candidate genes (e.g., ANK2 and VRTN) that might be mediated by TEs. Our findings highlight the important roles of young SINEs and provide a supplement for genotype-to-phenotype associations and modern breeding in pigs.

Native functions of short tandem repeats

Over a third of the human genome is comprised of repetitive sequences, including more than a million short tandem repeats (STRs). While studies of the pathologic consequences of repeat expansions that cause syndromic human diseases are extensive, the potential native functions of STRs are often ignored. Here, we summarize a growing body of research into the normal biological functions for repetitive elements across the genome, with a particular focus on the roles of STRs in regulating gene expression. We propose reconceptualizing the pathogenic consequences of repeat expansions as aberrancies in normal gene regulation. From this altered viewpoint, we predict that future work will reveal broader roles for STRs in neuronal function and as risk alleles for more common human neurological diseases.

Analysis of endometrial carcinoma TCGA reveals differences in DNA methylation in tumors from Black and White women

OBJECTIVE

Racial disparities exist in cancer patients both in incidence and death rates. In endometrial cancer, Black patients are reported to have higher incidence of aggressive endometrial cancer subtypes and higher death rates than White women. To date, diagnostic and prognostic biomarkers associated with race-specific methylation driven genes have yet to be identified. The objective of this study was to explore DNA methylation patterns in endometrial tumor samples from White and Black women.

METHODS

Differentially methylated CpGs (DMCs) and differentially methylated regions (DMRs) were identified in White tumor samples compared to Black tumor samples using Endometrial Carcinoma (EC) methylation and clinical data from The Cancer Genome Atlas (TCGA). Survival analysis was performed using survival R package and results were visualized using Kaplan-Meier plots. To access the correlation between changes in methylation and gene expression, we downloaded raw RNA-sequencing by Expectation-Maximization (RSEM) counts data from The Cancer Genome Atlas (TCGA) using TCGABiolinks package (v2.18.0).

RESULTS

Our analysis revealed 704 differentially methylated CpGs in tumors from Black and White women. These differentially methylated genes showed strong negative correlation with gene expression and statistically significant enrichment in regulatory regions such as DNase I hypersensitivity sites (DHSs) and transcription factor binding sites (TFBSs). Increased variability in methylation occurred in genes such as the insulin signaling pathway in Black tumor samples.

CONCLUSION

By using two independent statistical method based on means (DMR, DMCs) and variances (DVCs) on the endometrial carcinoma TCGA data, we showed that endometrial tumors from Black women are hypomethylated and more hypervariable than tumors from White women. In-depth investigation of these methylation driven markers can aid in successful management of endometrial cancer disparities and improved overall survival in Black and White populations.

Homologous pairing in short double-stranded DNA-grafted colloidal microspheres

Homologous pairing (HP), i.e., the pairing of similar or identical double-stranded DNA, is an insufficiently understood fundamental biological process. HP is now understood to also occur without protein mediation, but crucial mechanistic details remain poorly established. Unfortunately, systematic studies of sequence dependence are not practical due to the enormous number of nucleotide permutations and multiple possible conformations involved in existing biophysical strategies even when using as few as 150 basepairs. Here, we show that HP can occur in DNA as short as 18 basepairs in a colloidal microparticle-based system. Exemplary systematic studies include resolving opposing reports of the impact of % AT composition, validating the impact of nucleotide order and triplet framework and revealing isotropic bendability to be crucial for HP. These studies are enabled by statistical analysis of crystal size and fraction within coexisting fluid-crystal phases of double-stranded DNA-grafted colloidal microspheres, where crystallization is predicated by HP.

Roles of transposable elements in the regulation of mammalian transcription

Transposable elements (TEs) comprise about half of the mammalian genome. TEs often contain sequences capable of recruiting the host transcription machinery, which they use to express their own products and promote transposition. However, the regulatory sequences carried by TEs may affect host transcription long after the TEs have lost the ability to transpose. Recent advances in genome analysis and engineering have facilitated systematic interrogation of the regulatory activities of TEs. In this Review, we discuss diverse mechanisms by which TEs contribute to transcription regulation. Notably, TEs can donate enhancer and promoter sequences that influence the expression of host genes, modify 3D chromatin architecture and give rise to novel regulatory genes, including non-coding RNAs and transcription factors. We discuss how TEs spur regulatory evolution and facilitate the emergence of genetic novelties in mammalian physiology and development. By virtue of their repetitive and interspersed nature, TEs offer unique opportunities to dissect the effects of mutation and genomic context on the function and evolution of cis-regulatory elements. We argue that TE-centric studies hold the key to unlocking general principles of transcription regulation and evolution.

3D or Not 3D: Shaping the Genome during Development

One of the most fundamental questions in developmental biology is how one fertilized cell can give rise to a fully mature organism and how gene regulation governs this process. Precise spatiotemporal gene expression is required for development and is believed to be achieved through a complex interplay of sequence-specific information, epigenetic modifications, trans-acting factors, and chromatin folding. Here we review the role of chromatin folding during development, the mechanisms governing 3D genome organization, and how it is established in the embryo. Furthermore, we discuss recent advances and debated questions regarding the contribution of the 3D genome to gene regulation during organogenesis. Finally, we describe the mechanisms that can reshape the 3D genome, including disease-causing structural variations and the emerging view that transposable elements contribute to chromatin organization.

Recent Bioinformatic Progress to Identify Epigenetic Changes Associated to Transposable Elements

Transposable elements (TEs) are recognized for their great impact on the functioning and evolution of their host genomes. They are associated to various deleterious effects, which has led to the evolution of regulatory epigenetic mechanisms to control their activity. Despite these negative effects, TEs are also important actors in the evolution of genomes by promoting genetic diversity and new regulatory elements. Consequently, it is important to study the epigenetic modifications associated to TEs especially at a locus-specific level to determine their individual influence on gene functioning. To this aim, this short review presents the current bioinformatic tools to achieve this task.

CENP-B-mediated DNA loops regulate activity and stability of human centromeres

Chromosome inheritance depends on centromeres, epigenetically specified regions of chromosomes. While conventional human centromeres are known to be built of long tandem DNA repeats, much of their architecture remains unknown. Using single-molecule techniques such as AFM, nanopores, and optical tweezers, we find that human centromeric DNA exhibits complex DNA folds such as local hairpins. Upon binding to a specific sequence within centromeric regions, the DNA-binding protein CENP-B compacts centromeres by forming pronounced DNA loops between the repeats, which favor inter-chromosomal centromere compaction and clustering. This DNA-loop-mediated organization of centromeric chromatin participates in maintaining centromere position and integrity upon microtubule pulling during mitosis. Our findings emphasize the importance of DNA topology in centromeric regulation and stability.

An RNA Polymerase III General Transcription Factor Engages in Cell Type-Specific Chromatin Looping

Transcription factors (TFs) bind DNA in a sequence-specific manner and are generally cell type-specific factors and/or developmental master regulators. In contrast, general TFs (GTFs) are part of very large protein complexes and serve for RNA polymerases’ recruitment to promoter sequences, generally in a cell type-independent manner. Whereas, several TFs have been proven to serve as anchors for the 3D genome organization, the role of GTFs in genome architecture have not been carefully explored. Here, we used ChIP-seq and Hi-C data to depict the role of TFIIIC, one of the RNA polymerase III GTFs, in 3D genome organization. We find that TFIIIC genome occupancy mainly occurs at specific regions, which largely correspond to Alu elements; other characteristic classes of repetitive elements (REs) such as MIR, FLAM-C and ALR/alpha are also found depending on the cell’s developmental origin. The analysis also shows that TFIIIC-enriched regions are involved in cell type-specific DNA looping, which does not depend on colocalization with the master architectural protein CTCF. This work extends previous knowledge on the role of TFIIIC as a bona fide genome organizer whose action participates in cell type-dependent 3D genome looping via binding to REs.

Principles and functions of pericentromeric satellite DNA clustering into chromocenters

Simple non-coding tandem repeats known as satellite DNA are observed widely across eukaryotes. These repeats occupy vast regions at the centromere and pericentromere of chromosomes but their contribution to cellular function has remained incompletely understood. Here, we review the literature on pericentromeric satellite DNA and discuss its organization and functions across eukaryotic species. We specifically focus on chromocenters, DNA-dense nuclear foci that contain clustered pericentromeric satellite DNA repeats from multiple chromosomes. We first discuss chromocenter formation and the roles that epigenetic modifications, satellite DNA transcripts and sequence-specific satellite DNA-binding play in this process. We then review the newly emerging functions of chromocenters in genome encapsulation, the maintenance of cell fate and speciation. We specifically highlight how the rapid divergence of satellite DNA repeats impacts reproductive isolation between closely related species. Together, we underline the importance of this so-called 'junk DNA' in fundamental biological processes.

Nucleic acids movement and its relation to genome dynamics of repetitive DNA: Is cellular and intercellular movement of DNA and RNA molecules related to the evolutionary dynamic genome components?: Is cellular and intercellular movement of DNA and RNA molecules related to the evolutionary dynamic genome components?

There is growing evidence of evolutionary genome plasticity. The evolution of repetitive DNA elements, the major components of most eukaryotic genomes, involves the amplification of various classes of mobile genetic elements, the expansion of satellite DNA, the transfer of fragments or entire organellar genomes and may have connections with viruses. In addition to various repetitive DNA elements, a plethora of large and small RNAs migrate within and between cells during individual development as well as during evolution and contribute to changes of genome structure and function. Such migration of DNA and RNA molecules often results in horizontal gene transfer, thus shaping the whole genomic network of interconnected species. Here, we propose that a high evolutionary dynamism of repetitive genome components is often related to the migration/movement of DNA or RNA molecules. We speculate that the cytoplasm is probably an ideal compartment for such evolutionary experiments.

Neuron-specific chromosomal megadomain organization is adaptive to recent retrotransposon expansions

Regulatory mechanisms associated with repeat-rich sequences and chromosomal conformations in mature neurons remain unexplored. Here, we map cell-type specific chromatin domain organization in adult mouse cerebral cortex and report strong enrichment of Endogenous Retrovirus 2 (ERV2) repeat sequences in the neuron-specific heterochromatic B₂^NeuN+ megabase-scaling subcompartment. Single molecule long-read sequencing and comparative Hi-C chromosomal contact mapping in wild-derived SPRET/EiJ (Mus spretus) and laboratory inbred C57BL/6J (Mus musculus) reveal neuronal reconfigurations tracking recent ERV2 expansions in the murine germline, with significantly higher B₂^NeuN+ contact frequencies at sites with ongoing insertions in Mus musculus. Neuronal ablation of the retrotransposon silencer Kmt1e/Setdb1 triggers B₂^NeuN+ disintegration and rewiring with open chromatin domains enriched for cellular stress response genes, along with severe neuroinflammation and proviral assembly with infiltration of dendrites . We conclude that neuronal megabase-scale chromosomal architectures include an evolutionarily adaptive heterochromatic organization which, upon perturbation, results in transcriptional dysregulation and unleashes ERV2 proviruses with strong neuronal tropism.

Three-dimensional genome organization via triplex-forming RNAs

An increasing number of long noncoding RNAs (lncRNAs) have been proposed to act as nuclear organization factors during interphase. Direct RNA-DNA interactions can be achieved by the formation of triplex helix structures where a single-stranded RNA molecule hybridizes by complementarity into the major groove of double-stranded DNA. However, whether and how these direct RNA-DNA associations influence genome structure in interphase chromosomes remain poorly understood. Here we theorize that RNA organizes the genome in space via a triplex-forming mechanism. To test this theory, we apply a computational modeling approach of chromosomes that combines restraint-based modeling with polymer physics. Our models suggest that colocalization of triplex hotspots targeted by lncRNAs could contribute to large-scale chromosome compartmentalization cooperating, rather than competing, with architectural transcription factors such as CTCF.

Eukaryotic Genomes Show Strong Evolutionary Conservation of k-mer Composition and Correlation Contributions between Introns and Intergenic Regions

Several strongly conserved DNA sequence patterns in and between introns and intergenic regions (IIRs) consisting of short tandem repeats (STRs) with repeat lengths <3 bp have already been described in the kingdom of Animalia. In this work, we expanded the search and analysis of conserved DNA sequence patterns to a wider range of eukaryotic genomes. Our aims were to confirm the conservation of these patterns, to support the hypothesis on their functional constraints and/or the identification of unknown patterns. We pairwise compared genomic DNA sequences of genes, exons, CDS, introns and intergenic regions of 34 Embryophyta (land plants), 30 Protista and 29 Fungi using established k-mer-based (alignment-free) comparison methods. Additionally, the results were compared with values derived for Animalia in former studies. We confirmed strong correlations between the sequence structures of IIRs spanning over the entire domain of Eukaryotes. We found that the high correlations within introns, intergenic regions and between the two are a result of conserved abundancies of STRs with repeat units ≤2 bp (e.g., (AT)n). For some sequence patterns and their inverse complementary sequences, we found a violation of equal distribution on complementary DNA strands in a subset of genomes. Looking at mismatches within the identified STR patterns, we found specific preferences for certain nucleotides stable over all four phylogenetic kingdoms. We conclude that all of these conserved patterns between IIRs indicate a shared function of these sequence structures related to STRs.

Nuclear Dynamics and Chromatin Structure: Implications for Pancreatic Cancer

Changes in nuclear shape have been extensively associated with the dynamics and functionality of cancer cells. In most normal cells, nuclei have a regular ellipsoid shape and minimal variation in nuclear size; however, an irregular nuclear contour and abnormal nuclear size is often observed in cancer, including pancreatic cancer. Furthermore, alterations in nuclear morphology have become the 'gold standard' for tumor staging and grading. Beyond the utility of altered nuclear morphology as a diagnostic tool in cancer, the implications of altered nuclear structure for the biology and behavior of cancer cells are profound as changes in nuclear morphology could impact cellular responses to physical strain, adaptation during migration, chromatin organization, and gene expression. Here, we aim to highlight and discuss the factors that regulate nuclear dynamics and their implications for pancreatic cancer biology.

Toward the massive genome of Proteus anguinus-illuminating longevity, regeneration, convergent evolution, and metabolic disorders

Deciphering the genetic code of organisms with unusual phenotypes can help answer fundamental biological questions and provide insight into mechanisms relevant to human biomedical research. The cave salamander Proteus anguinus (Urodela: Proteidae), also known as the olm, is an example of a species with unique morphological and physiological adaptations to its subterranean environment, including regenerative abilities, resistance to prolonged starvation, and a life span of more than 100 years. However, the structure and sequence of the olm genome is still largely unknown owing to its enormous size, estimated at nearly 50 gigabases. An international Proteus Genome Research Consortium has been formed to decipher the olm genome. This perspective provides the scientific and biomedical rationale for exploring the olm genome and outlines potential outcomes, challenges, and methodological approaches required to analyze and annotate the genome of this unique amphibian.

Chromatin Conformation in Development and Disease

Chromatin domains and loops are important elements of chromatin structure and dynamics, but much remains to be learned about their exact biological role and nature. Topological associated domains and functional loops are key to gene expression and hold the answer to many questions regarding developmental decisions and diseases. Here, we discuss new findings, which have linked chromatin conformation with development, differentiation and diseases and hypothesized on various models while integrating all recent findings on how chromatin architecture affects gene expression during development, evolution and disease.

Complete Chloroplast Genome Sequence of Fortunella venosa (Champ. ex Benth.) C.C.Huang (Rutaceae): Comparative Analysis, Phylogenetic Relationships, and Robust Support for Its Status as an Independent Species

Fortunella venosa (Rutaceae) is an endangered species endemic to China and its taxonomic status has been controversial. The genus Fortunella contains a variety of important economic plants with high value in food, medicine, and ornamental. However, the placement of Genus Fortunella into Genus Citrus has led to controversy on its taxonomy and Systematics. In this present research, the Chloroplast genome of F. venosa was sequenced using the second-generation sequencing, and its structure and phylogenetic relationship analyzed. The results showed that the Chloroplast genome size of F. venosa was 160,265 bp, with a typical angiosperm four-part ring structure containing a large single copy region (LSC) (87,597 bp), a small single copy region (SSC) (18,732 bp), and a pair of inverted repeat regions (IRa\IRb) (26,968 bp each). There are 134 predicted genes in Chloroplast genome, including 89 protein-coding genes, 8 rRNAs, and 37 tRNAs. The GC-content of the whole Chloroplast genome was 43%, with the IR regions having a higher GC content than the LSC and the SSC regions. There were no rearrangements present in the Chloroplast genome; however, the IR regions showed obvious contraction and expansion. A total of 108 simple sequence repeats (SSRs) were present in the entire chloroplast genome and the nucleotide polymorphism was high in LSC and SSC. In addition, there is a preference for codon usage with the non-coding regions being more conserved than the coding regions. Phylogenetic analysis showed that species of Fortunella are nested in the genus of Citrus and the independent species status of F. venosa is supported robustly, which is significantly different from F. japonica. These findings will help in the development of DNA barcodes that can be useful in the study of the systematics and evolution of the genus Fortunella and the family Rutaceae.

Order and disorder: abnormal 3D chromatin organization in human disease

A precise three-dimensional (3D) organization of chromatin is central to achieve the intricate transcriptional patterns that are required to form complex organisms. Growing evidence supports an important role of 3D chromatin architecture in development and delineates its alterations as prominent causes of disease. In this review, we discuss emerging concepts on the fundamental forces shaping genomes in space and on how their disruption can lead to pathogenic phenotypes. We describe the molecular mechanisms underlying a wide range of diseases, from the systemic effects of coding mutations on 3D architectural factors, to the more tissue-specific phenotypes resulting from genetic and epigenetic modifications at specific loci. Understanding the connection between the 3D organization of the genome and its underlying biological function will allow a better interpretation of human pathogenesis.

Losing DNA methylation at repetitive elements and breaking bad

Background

DNA methylation is an epigenetic chromatin mark that allows heterochromatin formation and gene silencing. It has a fundamental role in preserving genome stability (including chromosome stability) by controlling both gene expression and chromatin structure. Therefore, the onset of an incorrect pattern of DNA methylation is potentially dangerous for the cells. This is particularly important with respect to repetitive elements, which constitute the third of the human genome.

Main body

Repetitive sequences are involved in several cell processes, however, due to their intrinsic nature, they can be a source of genome instability. Thus, most repetitive elements are usually methylated to maintain a heterochromatic, repressed state. Notably, there is increasing evidence showing that repetitive elements (satellites, long interspersed nuclear elements (LINEs), Alus) are frequently hypomethylated in various of human pathologies, from cancer to psychiatric disorders. Repetitive sequences’ hypomethylation correlates with chromatin relaxation and unscheduled transcription. If these alterations are directly involved in human diseases aetiology and how, is still under investigation.

Conclusions

Hypomethylation of different families of repetitive sequences is recurrent in many different human diseases, suggesting that the methylation status of these elements can be involved in preservation of human health. This provides a promising point of view towards the research of therapeutic strategies focused on specifically tuning DNA methylation of DNA repeats.

Nucleosome-induced homology recognition in chromatin

One of the least understood properties of chromatin is the ability of its similar regions to recognize each other through weak interactions. Theories based on electrostatic interactions between helical macromolecules suggest that the ability to recognize sequence homology is an innate property of the non-ideal helical structure of DNA. However, this theory does not account for the nucleosomal packing of DNA. Can homologous DNA sequences recognize each other while wrapped up in the nucleosomes? Can structural homology arise at the level of nucleosome arrays? Here, we present a theoretical model for the recognition potential well between chromatin fibres sliding against each other. This well is different from the one predicted for bare DNA; the minima in energy do not correspond to literal juxtaposition, but are shifted by approximately half the nucleosome repeat length. The presence of this potential well suggests that nucleosome positioning may induce mutual sequence recognition between chromatin fibres and facilitate the formation of chromatin nanodomains. This has implications for nucleosome arrays enclosed between CTCF-cohesin boundaries, which may form stiffer stem-like structures instead of flexible entropically favourable loops. We also consider switches between chromatin states, e.g. through acetylation/deacetylation of histones, and discuss nucleosome-induced recognition as a precursory stage of genetic recombination.

The emerging role of chromatin remodelers in neurodevelopmental disorders: a developmental perspective

Neurodevelopmental disorders (NDDs), including intellectual disability (ID) and autism spectrum disorders (ASD), are a large group of disorders in which early insults during brain development result in a wide and heterogeneous spectrum of clinical diagnoses. Mutations in genes coding for chromatin remodelers are overrepresented in NDD cohorts, pointing towards epigenetics as a convergent pathogenic pathway between these disorders. In this review we detail the role of NDD-associated chromatin remodelers during the developmental continuum of progenitor expansion, differentiation, cell-type specification, migration and maturation. We discuss how defects in chromatin remodelling during these early developmental time points compound over time and result in impaired brain circuit establishment. In particular, we focus on their role in the three largest cell populations: glutamatergic neurons, GABAergic neurons, and glia cells. An in-depth understanding of the spatiotemporal role of chromatin remodelers during neurodevelopment can contribute to the identification of molecular targets for treatment strategies.

Toward incorporating epigenetics into regulation of gene expression in the parasite Trichomonas vaginalis

Trichomonas vaginalis is an extracellular parasite that colonizes the human urogenital tract, causing a highly prevalent sexually transmitted infection. The parasite must change its transcriptional profile in order to establish and maintain infection. However, few core regulatory elements and transcription factors have been identified to date and little is known about other mechanisms that may control these rapid changes in gene expression during parasite infection. In the last years, epigenetic mechanisms involved in the regulation of gene expression have been gaining major attention in this parasite. In this review, we summarize and discuss the major advances of the last few years with regard to epigenetics (DNA methylation, post-translational histone modifications, and histone variants) in the parasite T. vaginalis. These studies can shed light into our current understanding of this parasite's biology with far-reaching implications for the prognosis and treatment of trichomoniasis.

In situ genome sequencing resolves DNA sequence and structure in intact biological samples

Understanding genome organization requires integration of DNA sequence and three-dimensional spatial context; however, existing genome-wide methods lack either base pair sequence resolution or direct spatial localization. Here, we describe in situ genome sequencing (IGS), a method for simultaneously sequencing and imaging genomes within intact biological samples. We applied IGS to human fibroblasts and early mouse embryos, spatially localizing thousands of genomic loci in individual nuclei. Using these data, we characterized parent-specific changes in genome structure across embryonic stages, revealed single-cell chromatin domains in zygotes, and uncovered epigenetic memory of global chromosome positioning within individual embryos. These results demonstrate how IGS can directly connect sequence and structure across length scales from single base pairs to whole organisms.

3D Genome of macaque fetal brain reveals evolutionary innovations during primate corticogenesis

Elucidating the regulatory mechanisms of human brain evolution is essential to understanding human cognition and mental disorders. We generated multi-omics profiles and constructed a high-resolution map of 3D genome architecture of rhesus macaque during corticogenesis. By comparing the 3D genomes of human, macaque, and mouse brains, we identified many human-specific chromatin structure changes, including 499 topologically associating domains (TADs) and 1,266 chromatin loops. The human-specific loops are significantly enriched in enhancer-enhancer interactions, and the regulated genes show human-specific expression changes in the subplate, a transient zone of the developing brain critical for neural circuit formation and plasticity. Notably, many human-specific sequence changes are located in the human-specific TAD boundaries and loop anchors, which may generate new transcription factor binding sites and chromatin structures in human. Collectively, the presented data highlight the value of comparative 3D genome analyses in dissecting the regulatory mechanisms of brain development and evolution.

Homology-mediated inter-chromosomal interactions in hexaploid wheat lead to specific subgenome territories following polyploidization and introgression

BACKGROUND

Polyploidization and introgression are major events driving plant genome evolution and influencing crop breeding. However, the mechanisms underlying the higher-order chromatin organization of subgenomes and alien chromosomes are largely unknown.

RESULTS

We probe the three-dimensional chromatin architecture of Aikang 58 (AK58), a widely cultivated allohexaploid wheat variety in China carrying the 1RS/1BL translocation chromosome. The regions involved in inter-chromosomal interactions, both within and between subgenomes, have highly similar sequences. Subgenome-specific territories tend to be connected by subgenome-dominant homologous transposable elements (TEs). The alien 1RS chromosomal arm, which was introgressed from rye and differs from its wheat counterpart, has relatively few inter-chromosome interactions with wheat chromosomes. An analysis of local chromatin structures reveals topologically associating domain (TAD)-like regions covering 52% of the AK58 genome, the boundaries of which are enriched with active genes, zinc-finger factor-binding motifs, CHH methylation, and 24-nt small RNAs. The chromatin loops are mostly localized around TAD boundaries, and the number of gene loops is positively associated with gene activity.

CONCLUSIONS

The present study reveals the impact of the genetic sequence context on the higher-order chromatin structure and subgenome stability in hexaploid wheat. Specifically, we characterized the sequence homology-mediated inter-chromosome interactions and the non-canonical role of subgenome-biased TEs. Our findings may have profound implications for future investigations of the interplay between genetic sequences and higher-order structures and their consequences on polyploid genome evolution and introgression-based breeding of crop plants.

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.

Both p53 codon 72 Arg/Arg and pro/Arg genotypes in glioblastoma multiforme are associated with a better prognosis in bevacizumab treatment

BACKGROUND

It has previously been shown that bevacizumab, when added to chemotherapy, improved overall survival in several cancers. In glioblastoma multiforme (GBM), bevacizumab increased progression-free survival and it is widely used for tumor recurrence, though it has failed to improve overall survival (OS) in controlled trials. However, an effective biomarker for predicting the prognosis of bevacizumab treatment has yet to be identified. This study, therefore, aimed to retrospectively analyze the polymorphisms of p53 codon 72 and the clinical characteristics of GBM specimens from Taiwanese patients.

METHODS

The polymorphisms of p53 codon 72 in 99 patients with GBM treated at Taichung Veterans General Hospital in Taiwan from 2007 to 2017 were analyzed using direct DNA sequencing and PCR-RFLP analysis.

RESULTS

We found that among these GBM patients, the distribution of codon 72 polymorphisms was 28.3% for proline homozygotes (Pro/Pro), 38.4% for arginine homozygotes (Arg/Arg), and 33.3% for proline/arginine heterozygotes (Pro/Arg). Although the polymorphisms of p53 codon 72 were not directly associated with the overall survival of GBM, both the Arg/Arg and Arg/Pro genotypes were associated with significant benefits in terms of overall survival in patients treated with CCRT plus bevacizumab compared to patients treated with CCRT alone.

CONCLUSIONS

This pilot study suggests that both the Arg/Arg and Arg/Pro genotypes of p53 codon 72 polymorphism may have value as independent prognostic or predictive parameters for bevacizumab treatment response and failure. Relatedly, the results of the study further demonstrate the utility of stratifying GBM patients according to bevacizumab sensitivity.

90 YEARS OF PROGESTERONE: Molecular mechanisms of progesterone receptor action on the breast cancer genome

Gene regulation by steroid hormones has been at the forefront in elucidating the intricacies of transcriptional regulation in eukaryotes ever since the discovery by Karlson and Clever that the insect steroid hormone ecdysone induces chromatin puffs in giant chromosomes. After the successful cloning of the hormone receptors toward the end of the past century, detailed mechanistic insight emerged in some model systems, in particular the MMTV provirus. With the arrival of next generation DNA sequencing and the omics techniques, we have gained even further insight into the global cellular response to steroid hormones that in the past decades also extended to the function of the 3D genome topology. More recently, advances in high resolution microcopy, single cell genomics and the new vision of liquid-liquid phase transitions in the context of nuclear space bring us closer than ever to unravelling the logic of gene regulation and its complex integration of global cellular signaling networks. Using the function of progesterone and its cellular receptor in breast cancer cells, we will briefly summarize the history and describe the present extent of our knowledge on how regulatory proteins deal with the chromatin structure to gain access to DNA sequences and interpret the genomic instructions that enable cells to respond selectively to external signals by reshaping their gene regulatory networks.

Genome survey sequencing and identification of genomic SSR markers for Rhododendron micranthum

Rhododendron micranthum is an evergreen shrub species widely distributed in China that has high ornamental and medicinal value. However, there is a lack of molecular and genomic data for this plant, which severely restricts the development of its relevant research. The objective of the present study was to conduct a first genomic survey of R. micranthum and determine its whole-genome sequencing scheme. Next-generation sequencing (Illumina Hi-Seq Xten) was used to measure the genome size of R. micranthum, K-mer analysis were employed to investigate its genomic profile. Finally, we conducted bioinformatics methods to performed SSR (simple sequence repeat) prediction based on the genomic data. The genome size of R. micranthum was estimated to be 554.22 Mb. The heterozygosity ratio was 0.93%, and the sequence repeat ratio was calculated to be 49.17%. The clean reads of R. micranthum were assembled into 2281551 scaffolds with a N50 value of 916 bp. A total of 479724 SSR molecular markers were identified in the R. micranthum genome, and 871656 pairs of primers designed for application. Among of them, 100 primer pairs were validated, and 71 primer pairs were successfully amplified. In summary, the R. micranthum genome is complex with high heterozygosity and low repeated sequences. In future whole-genome research in R. micranthum, higher-depth '2+3' (Illumina+PacBio) sequencing may yield better assembly results.

Structure encoding in DNA

It is proposed that transposons and related long non-coding RNA define the fine structure of body parts. Although morphogens have long been known to direct the formation of many gross structures in early embryonic development, they do not have the necessary precision to define a structure down to the individual cellular level. Using the distinction between procedural and declarative knowledge in information processing as an analogy, it is hypothesized that DNA encodes fine structure in a manner that is different from the genetic code for proteins. The hypothesis states that repeated or near-repeated sequences that are in transposons and non-coding RNA define body part structures. As the cells in a body part go through the epigenetic process of differentiation, the action of methylation serves to inactivate all but the relevant structure definitions and some associated cell type genes. The transposons left active will then physically modify the DNA sequence in the heterochromatin to establish the local context in the three-dimensional body part structure. This brings the encoded definition of the cell type to the histone. The histone code for that cell type starts the regulatory cascade that turns on the genes associated with that particular type of cell, transforming it from a multipotent cell to a fully differentiated cell. This mechanism creates structures in the musculoskeletal system, the organs of the body, the major parts of the brain, and other systems.

Enhancer-gene maps in the human and zebrafish genomes using evolutionary linkage conservation

The spatiotemporal expression of genes is controlled by enhancer sequences that bind transcription factors. Identifying the target genes of enhancers remains difficult because enhancers regulate gene expression over long genomic distances. To address this, we used an evolutionary approach to build two genome-wide maps of predicted enhancer-gene associations in the human and zebrafish genomes. Evolutionary conserved sequences were linked to their predicted target genes using PEGASUS, a bioinformatics method that relies on evolutionary conservation of synteny. The analysis of these maps revealed that the number of predicted enhancers linked to a gene correlate with its expression breadth. Comparison of both maps identified hundreds of putative vertebrate ancestral regulatory relationships from which we could determine that predicted enhancer-gene distances scale with genome size despite strong positional conservation. The two maps represent a resource for further studies, including the prioritization of sequence variants in whole genome sequence of patients affected by genetic diseases.

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.

Integrating transposable elements in the 3D genome

Chromosome organisation is increasingly recognised as an essential component of genome regulation, cell fate and cell health. Within the realm of transposable elements (TEs) however, the spatial information of how genomes are folded is still only rarely integrated in experimental studies or accounted for in modelling. Whilst polymer physics is recognised as an important tool to understand the mechanisms of genome folding, in this commentary we discuss its potential applicability to aspects of TE biology. Based on recent works on the relationship between genome organisation and TE integration, we argue that existing polymer models may be extended to create a predictive framework for the study of TE integration patterns. We suggest that these models may offer orthogonal and generic insights into the integration profiles (or "topography") of TEs across organisms. In addition, we provide simple polymer physics arguments and preliminary molecular dynamics simulations of TEs inserting into heterogeneously flexible polymers. By considering this simple model, we show how polymer folding and local flexibility may generically affect TE integration patterns. The preliminary discussion reported in this commentary is aimed to lay the foundations for a large-scale analysis of TE integration dynamics and topography as a function of the three-dimensional host genome.

New Aspects of Magnesium Function: A Key Regulator in Nucleosome Self-Assembly, Chromatin Folding and Phase Separation

Metal cations are associated with many biological processes. The effects of these cations on nucleic acids and chromatin were extensively studied in the early stages of nucleic acid and chromatin research. The results revealed that some monovalent and divalent metal cations, including Mg²⁺, profoundly affect the conformations and stabilities of nucleic acids, the folding of chromatin fibers, and the extent of chromosome condensation. Apart from these effects, there have only been a few reports on the functions of these cations. In 2007 and 2013, however, Mg²⁺-implicated novel phenomena were found: Mg²⁺ facilitates or enables both self-assembly of identical double-stranded (ds) DNA molecules and self-assembly of identical nucleosomes in vitro. These phenomena may be deeply implicated in the heterochromatin domain formation and chromatin-based phase separation. Furthermore, a recent study showed that elevation of the intranuclear Mg²⁺ concentration causes unusual differentiation of mouse ES (embryonic stem) cells. All of these phenomena seem to be closely related to one another. Mg²⁺ seems to be a key regulator of chromatin dynamics and chromatin-based biological processes.

MIR sequences recruit zinc finger protein ZNF768 to expressed genes

Mammalian-wide interspersed repeats (MIRs) are retrotransposed elements of mammalian genomes. Here, we report the specific binding of zinc finger protein ZNF768 to the sequence motif GCTGTGTG (N20) CCTCTCTG in the core region of MIRs. ZNF768 binding is preferentially associated with euchromatin and promoter regions of genes. Binding was observed for genes expressed in a cell type-specific manner in human B cell line Raji and osteosarcoma U2OS cells. Mass spectrometric analysis revealed binding of ZNF768 to Elongator components Elp1, Elp2 and Elp3 and other nuclear factors. The N-terminus of ZNF768 contains a heptad repeat array structurally related to the C-terminal domain (CTD) of RNA polymerase II. This array evolved in placental animals but not marsupials and monotreme species, displays species-specific length variations, and possibly fulfills CTD related functions in gene regulation. We propose that the evolution of MIRs and ZNF768 has extended the repertoire of gene regulatory mechanisms in mammals and that ZNF768 binding is associated with cell type-specific gene expression.

The Impact of Centromeres on Spatial Genome Architecture

The development of new technologies and experimental techniques is enabling researchers to see what was once unable to be seen. For example, the centromere was first seen as the mediator between spindle fiber and chromosome during mitosis and meiosis. Although this continues to be its most prominent role, we now know that the centromere functions beyond cellular division with important roles in genome organization and chromatin regulation. Here we aim to share the structures and functions of centromeres in various organisms beginning with the diversity of their DNA sequence anatomies. We zoom out to describe their position in the nucleus and ultimately detail the different ways they contribute to genome organization and regulation at the spatial level.

Generative modeling of multi-mapping reads with mHi-C advances analysis of Hi-C studies

Current Hi-C analysis approaches are unable to account for reads that align to multiple locations, and hence underestimate biological signal from repetitive regions of genomes. We developed and validated mHi-C, a multi-read mapping strategy to probabilistically allocate Hi-C multi-reads. mHi-C exhibited superior performance over utilizing only uni-reads and heuristic approaches aimed at rescuing multi-reads on benchmarks. Specifically, mHi-C increased the sequencing depth by an average of 20% resulting in higher reproducibility of contact matrices and detected interactions across biological replicates. The impact of the multi-reads on the detection of significant interactions is influenced marginally by the relative contribution of multi-reads to the sequencing depth compared to uni-reads, cis-to-trans ratio of contacts, and the broad data quality as reflected by the proportion of mappable reads of datasets. Computational experiments highlighted that in Hi-C studies with short read lengths, mHi-C rescued multi-reads can emulate the effect of longer reads. mHi-C also revealed biologically supported bona fide promoter-enhancer interactions and topologically associating domains involving repetitive genomic regions, thereby unlocking a previously masked portion of the genome for conformation capture studies.

Complete chloroplast genome sequence of Betula platyphylla: gene organization, RNA editing, and comparative and phylogenetic analyses

BACKGROUND

Betula platyphylla is a common tree species in northern China that has high economic and medicinal value. Our laboratory has been devoted to genome research on B. platyphylla for approximately 10 years. As primary organelle genomes, the complete genome sequences of chloroplasts are important to study the divergence of species, RNA editing and phylogeny. In this study, we sequenced and analyzed the complete chloroplast (cp) genome sequence of B. platyphylla.

RESULTS

The complete cp genome of B. platyphylla was 160,518 bp in length, which included a pair of inverted repeats (IRs) of 26,056 bp that separated a large single copy (LSC) region of 89,397 bp and a small single copy (SSC) region of 19,009 bp. The annotation contained a total of 129 genes, including 84 protein-coding genes, 37 tRNA genes and 8 rRNA genes. There were 3 genes using alternative initiation codons. Comparative genomics showed that the sequence of the Fagales species cp genome was relatively conserved, but there were still some high variation regions that could be used as molecular markers. The IR expansion event of B. platyphylla resulted in larger cp genomes and rps19 pseudogene formation. The simple sequence repeat (SSR) analysis showed that there were 105 SSRs in the cp genome of B. platyphylla. RNA editing sites recognition indicated that at least 80 RNA editing events occurred in the cp genome. Most of the substitutions were C to U, while a small proportion of them were not. In particular, three editing loci on the rRNA were converted to more than two other bases that had never been reported. For synonymous conversion, most of them increased the relative synonymous codon usage (RSCU) value of the codons. The phylogenetic analysis suggested that B. platyphylla had a closer evolutionary relationship with B. pendula than B. nana.

CONCLUSIONS

In this study, we not only obtained and annotated the complete cp genome sequence of B. platyphylla, but we also identified new RNA editing sites and predicted the phylogenetic relationships among Fagales species. These findings will facilitate genomic, genetic engineering and phylogenetic studies of this important species.

Analysis of 3D genomic interactions identifies candidate host genes that transposable elements potentially regulate

BACKGROUND

The organization of chromatin in the nucleus plays an essential role in gene regulation. About half of the mammalian genome comprises transposable elements. Given their repetitive nature, reads associated with these elements are generally discarded or randomly distributed among elements of the same type in genome-wide analyses. Thus, it is challenging to identify the activities and properties of individual transposons. As a result, we only have a partial understanding of how transposons contribute to chromatin folding and how they impact gene regulation.

RESULTS

Using PCR and Capture-based chromosome conformation capture (3C) approaches, collectively called 4Tran, we take advantage of the repetitive nature of transposons to capture interactions from multiple copies of endogenous retrovirus (ERVs) in the human and mouse genomes. With 4Tran-PCR, reads are selectively mapped to unique regions in the genome. This enables the identification of transposable element interaction profiles for individual ERV families and integration events specific to particular genomes. With this approach, we demonstrate that transposons engage in long-range intra-chromosomal interactions guided by the separation of chromosomes into A and B compartments as well as topologically associated domains (TADs). In contrast to 4Tran-PCR, Capture-4Tran can uniquely identify both ends of an interaction that involve retroviral repeat sequences, providing a powerful tool for uncovering the individual transposable element insertions that interact with and potentially regulate target genes.

CONCLUSIONS

4Tran provides new insight into the manner in which transposons contribute to chromosome architecture and identifies target genes that transposable elements can potentially control.

Structure, clustering and functional insights of repeats configurations in the upstream promoter region of the human coding genes

Background

Repetitive DNA sequences (Repeats) are significant regions in the human genome that have a specific genomic distribution, structure, and several binding sites for genome architecture and function. In consequence, the possible configurations of Repeats in specific and dynamic regions like the gene promoters could define footprints for molecular mechanisms, pathways, and cell function beyond their density in the genome. Here we explored the distribution of Repeats in the upstream promoter region of the human coding genes with the aim to identify specific configurations, clusters and functional meaning of those elements. Our method includes structural descriptions, hierarchical clustering, pathway association, and functional enrichment analysis.

Results

We report here several configurations of Repeats in the upstream promoter region (UPR), which define 2729 patterns for the 80% of the human coding genes. There are 47 types of Repeats in these configurations, where the most frequent were Alu, Low_complexity, MIR, Simple_repeat, LINE/L2, LINE/L1, hAT-Charlie, and ERV1. The distribution, length, and the high frequency of Repeats in the UPR defines several patterns and clusters, where the minimum frequency of configuration among Repeats was higher than 0.7. We found those clusters associated with cellular pathways and ontologies; thus, it was plausible to determine groups of Repeats to specific functional insights, for example, pathways for Genetic Information Processing or Metabolism shows particular groups of Repeats with specific configurations.

Conclusion

Based on these findings, we propose that specific configurations of repetitive elements describe frequent patterns in the upstream promoter for sets of human coding genes, which those correlated to specific and essential cell pathways and functions.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5196-6) contains supplementary material, which is available to authorized users.