MIR retrotransposon sequences provide insulators to the human genome

Massively parallel characterization of insulator activity across the genome

A key question in regulatory genomics is whether cis-regulatory elements (CREs) are modular elements that can function anywhere in the genome, or whether they are adapted to certain genomic locations. To distinguish between these possibilities we develop MPIRE (Massively Parallel Integrated Regulatory Elements), a technology for recurrently assaying CREs at thousands of defined locations across the genome in parallel. MPIRE allows us to separate the intrinsic activity of CREs from the effects of their genomic environments. We apply MPIRE to assay three insulator sequences at thousands of genomic locations and find that each insulator functions in locations with distinguishable properties. All three insulators can block enhancers, but each insulator blocks specific enhancers at specific locations. However, only ALOXE3 appears to block heterochromatin silencing. We conclude that insulator function is highly context dependent and that MPIRE is a robust method for revealing the context dependencies of CREs.

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.

Mouse retrotransposons: sequence structure, evolutionary age, genomic distribution and function

Retrotransposons are transposable elements that are transposed via transcription and reverse transcription. Their copies have accumulated in the genome of mammals, occupying approximately 40% of mammalian genomic mass. These copies are often involved in numerous phenomena, such as chromatin spatial organization, gene expression, development and disease, and have been recognized as a driving force in evolution. Different organisms have gained specific retrotransposon subfamilies and retrotransposed copies, such as hundreds of Mus-specific subfamilies with diverse sequences and genomic locations. Despite this complexity, basic information is still necessary for present-day genomic and epigenomic studies. Herein, we describe the characteristics of each subfamily of Mus-specific retrotransposons in terms of sequence structure, phylogenetic relationships, evolutionary age, and preference for A or B compartments of chromatin.

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.

TAD evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find 14% of all human TAD boundaries to be shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons compared to species-specific boundaries. CRISPR-Cas9 knockouts of an ultraconserved boundary in a mouse model lead to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in the upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations and showcases the functional importance of TAD evolution.

In silico discovery of repetitive elements as key sequence determinants of 3D genome folding

Natural and experimental genetic variants can modify DNA loops and insulating boundaries to tune transcription, but it is unknown how sequence perturbations affect chromatin organization genome wide. We developed a deep-learning strategy to quantify the effect of any insertion, deletion, or substitution on chromatin contacts and systematically scored millions of synthetic variants. While most genetic manipulations have little impact, regions with CTCF motifs and active transcription are highly sensitive, as expected. Our unbiased screen and subsequent targeted experiments also point to noncoding RNA genes and several families of repetitive elements as CTCF-motif-free DNA sequences with particularly large effects on nearby chromatin interactions, sometimes exceeding the effects of CTCF sites and explaining interactions that lack CTCF. We anticipate that our disruption tracks may be of broad interest and utility as a measure of 3D genome sensitivity, and our computational strategies may serve as a template for biological inquiry with deep learning.

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.

Genetic variation and domestication of horses revealed by 10 chromosome-level genomes and whole-genome resequencing

Understanding the genetic variations of the horse (Equus caballus) genome will improve breeding conservation and welfare. However, genetic variations in long segments, such as structural variants (SVs), remain understudied. We de novo assembled 10 chromosome-level three-dimensional horse genomes, each representing a distinct breed, and analysed horse SVs using a multi-assembly approach. Our findings suggest that SVs with the accumulation of mammalian-wide interspersed repeats related to long interspersed nuclear elements might be a horse-specific mechanism to modulate genome-wide gene regulatory networks. We found that olfactory receptors were commonly loss and accumulated deleterious mutations, but no purge of deleterious mutations occurred during horse domestication. We examined the potential effects of SVs on the spatial structure of chromatin via topologically associating domains (TADs). Breed-specific TADs were significantly enriched by breed-specific SVs. We identified 4199 unique breakpoint-resolved novel insertions across all chromosomes that account for 2.84 Mb sequences missing from the reference genome. Several novel insertions might have potential functional consequences, as 519 appeared to reside within 449 gene bodies. These genes are primarily involved in pathogen recognition, innate immune responses and drug metabolism. Moreover, 37 diverse horses were resequenced. Combining this with public data, we analysed 97 horses through a comparative population genomics approach to identify the genetic basis underlying breed characteristics using Thoroughbreds as a case study. We provide new scientific evidence for horse domestication, an understanding of the genetic mechanism underlying the phenotypic evolution of horses, and a comprehensive genetic variation resource for further genetic studies of horses.

SHIELD: a platform for high-throughput screening of barrier-type DNA elements in human cells

Chromatin boundary elements contribute to the partitioning of mammalian genomes into topological domains to regulate gene expression. Certain boundary elements are adopted as DNA insulators for safe and stable transgene expression in mammalian cells. These elements, however, are ill-defined and less characterized in the non-coding genome, partially due to the lack of a platform to readily evaluate boundary-associated activities of putative DNA sequences. Here we report SHIELD (Site-specific Heterochromatin Insertion of Elements at Lamina-associated Domains), a platform tailored for the high-throughput screening of barrier-type DNA elements in human cells. SHIELD takes advantage of the high specificity of serine integrase at heterochromatin, and exploits the natural heterochromatin spreading inside lamina-associated domains (LADs) for the discovery of potent barrier elements. We adopt SHIELD to evaluate the barrier activity of 1000 DNA elements in a high-throughput manner and identify 8 candidates with barrier activities comparable to the core region of cHS4 element in human HCT116 cells. We anticipate SHIELD could facilitate the discovery of novel barrier DNA elements from the non-coding genome in human cells.

The Pol III transcriptome: Basic features, recurrent patterns, and emerging roles in cancer

The RNA polymerase III (Pol III) transcriptome is universally comprised of short, highly structured noncoding RNA (ncRNA). Through RNA-protein interactions, the Pol III transcriptome actuates functional activities ranging from nuclear gene regulation (7SK), splicing (U6, U6atac), and RNA maturation and stability (RMRP, RPPH1, Y RNA), to cytoplasmic protein targeting (7SL) and translation (tRNA, 5S rRNA). In higher eukaryotes, the Pol III transcriptome has expanded to include additional, recently evolved ncRNA species that effectively broaden the footprint of Pol III transcription to additional cellular activities. Newly evolved ncRNAs function as riboregulators of autophagy (vault), immune signaling cascades (nc886), and translation (Alu, BC200, snaR). Notably, upregulation of Pol III transcription is frequently observed in cancer, and multiple ncRNA species are linked to both cancer progression and poor survival outcomes among cancer patients. In this review, we outline the basic features and functions of the Pol III transcriptome, and the evidence for dysregulation and dysfunction for each ncRNA in cancer. When taken together, recurrent patterns emerge, ranging from shared functional motifs that include molecular scaffolding and protein sequestration, overlapping protein interactions, and immunostimulatory activities, to the biogenesis of analogous small RNA fragments and noncanonical miRNAs, augmenting the function of the Pol III transcriptome and further broadening its role in cancer. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Processing > Processing of Small RNAs RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

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.

Integrative profiling of extrachromosomal circular DNA in placenta and maternal plasma provides insights into the biology of fetal growth restriction and reveals potential biomarkers

Introduction: Fetal growth restriction (FGR) is a placenta-mediated pregnancy complication that predisposes fetuses to perinatal complications. Maternal plasma cell-free DNA harbors DNA originating from placental trophoblasts, which is promising for the prenatal diagnosis and prediction of pregnancy complications. Extrachromosomal circular DNA (eccDNA) is emerging as an ideal biomarker and target for several diseases. Methods: We utilized eccDNA sequencing and bioinformatic pipeline to investigate the characteristics and associations of eccDNA in placenta and maternal plasma, the role of placental eccDNA in the pathogenesis of FGR, and potential plasma eccDNA biomarkers of FGR. Results: Using our bioinformatics pipelines, we identified multi-chromosomal-fragment and single-fragment eccDNA in placenta, but almost exclusively single-fragment eccDNA in maternal plasma. Relative to that in plasma, eccDNA in placenta was larger and substantially more abundant in exons, untranslated regions, promoters, repetitive elements [short interspersed nuclear elements (SINEs)/Alu, SINEs/mammalian-wide interspersed repeats, long terminal repeats/endogenous retrovirus-like elements, and single recognition particle RNA], and transcription factor binding motifs. Placental multi-chromosomal-fragment eccDNA was enriched in confident enhancer regions predicted to pertain to genes in apoptosis, energy, cell growth, and autophagy pathways. Placental eccDNA-associated genes whose abundance differed between the FGR and control groups were associated with immunity-related gene ontology (GO) terms. The combined analysis of plasma and placental eccDNA-associated genes in the FGR and control groups led to the identification of potential biomarkers that were assigned to the GO terms of the epigenetic regulation of gene expression and nutrient-related processes, respectively. Conclusion: Together, our results highlight links between placenta functions and multi-chromosomal-fragment and single-fragment eccDNA. The integrative analysis of placental and plasma eccDNA confirmed the potential of these molecules as disease-specific biomarkers of FGR.

Population genomics reveals mechanisms and dynamics of de novo expressed open reading frame emergence in Drosophila melanogaster

Novel genes are essential for evolutionary innovations and differ substantially even between closely related species. Recently, multiple studies across many taxa showed that some novel genes arise de novo, that is, from previously noncoding DNA. To characterize the underlying mutations that allowed de novo gene emergence and their order of occurrence, homologous regions must be detected within noncoding sequences in closely related sister genomes. So far, most studies do not detect noncoding homologs of de novo genes because of incomplete assemblies and annotations, and long evolutionary distances separating genomes. Here, we overcome these issues by searching for de novo expressed open reading frames (neORFs), the not-yet fixed precursors of de novo genes that emerged within a single species. We sequenced and assembled genomes with long-read technology and the corresponding transcriptomes from inbred lines of Drosophila melanogaster, derived from seven geographically diverse populations. We found line-specific neORFs in abundance but few neORFs shared by lines, suggesting a rapid turnover. Gain and loss of transcription is more frequent than the creation of ORFs, for example, by forming new start and stop codons. Consequently, the gain of ORFs becomes rate limiting and is frequently the initial step in neORFs emergence. Furthermore, transposable elements (TEs) are major drivers for intragenomic duplications of neORFs, yet TE insertions are less important for the emergence of neORFs. However, highly mutable genomic regions around TEs provide new features that enable gene birth. In conclusion, neORFs have a high birth-death rate, are rapidly purged, but surviving neORFs spread neutrally through populations and within genomes.

TAD Evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species, and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find that only 14% of all human TAD boundaries are shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons, compared to species-specific boundaries. CRISPR-Cas9 knockouts of two ultraconserved boundaries in mouse models leads to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations, and showcase the functional importance of TAD evolution.

Human Endogenous Retroviruses in Diseases

Human endogenous retroviruses (HERVs), which are conserved sequences of ancient retroviruses, are widely distributed in the human genome. Although most HERVs have been rendered inactive by evolution, some have continued to exhibit important cytological functions. HERVs in the human genome perform dual functions: on the one hand, they are involved in important physiological processes such as placental development and immune regulation; on the other hand, their aberrant expression is closely associated with the pathological processes of several diseases, such as cancers, autoimmune diseases, and viral infections. HERVs can also regulate a variety of host cellular functions, including the expression of protein-coding genes and regulatory elements that have evolved from HERVs. Here, we present recent research on the roles of HERVs in viral infections and cancers, including the dysregulation of HERVs in various viral infections, HERV-induced epigenetic modifications of histones (such as methylation and acetylation), and the potential mechanisms of HERV-mediated antiviral immunity. We also describe therapies to improve the efficacy of vaccines and medications either by directly or indirectly targeting HERVs, depending on the HERV.

T3E: a tool for characterising the epigenetic profile of transposable elements using ChIP-seq data

BACKGROUND

Despite the advent of Chromatin Immunoprecipitation Sequencing (ChIP-seq) having revolutionised our understanding of the mammalian genome's regulatory landscape, many challenges remain. In particular, because of their repetitive nature, the sequencing reads derived from transposable elements (TEs) pose a real bioinformatics challenge, to the point that standard analysis pipelines typically ignore reads whose genomic origin cannot be unambiguously ascertained.

RESULTS

We show that discarding ambiguously mapping reads may lead to a systematic underestimation of the number of reads associated with young TE families/subfamilies. We also provide evidence suggesting that the strategy of randomly permuting the location of the read mappings (or the TEs) that is often used to compute the background for enrichment calculations at TE families/subfamilies can result in both false positive and negative enrichments. To address these problems, we present the Transposable Element Enrichment Estimator (T3E), a tool that makes use of ChIP-seq data to characterise the epigenetic profile of associated TE families/subfamilies. T3E weights the number of read mappings assigned to the individual TE copies of a family/subfamily by the overall number of genomic loci to which the corresponding reads map, and this is done at the single nucleotide level. In addition, T3E computes ChIP-seq enrichment relative to a background estimated based on the distribution of the read mappings in the input control DNA. We demonstrated the capabilities of T3E on 23 different ChIP-seq libraries. T3E identified enrichments that were consistent with previous studies. Furthermore, T3E detected context-specific enrichments that are likely to pinpoint unexplored TE families/subfamilies with individual TE copies that have been frequently exapted as cis-regulatory elements during the evolution of mammalian regulatory networks.

CONCLUSIONS

T3E is a novel open-source computational tool (available for use at: https://github.com/michelleapaz/T3E ) that overcomes some of the pitfalls associated with the analysis of ChIP-seq data arising from the repetitive mammalian genome and provides a framework to shed light on the epigenetics of entire TE families/subfamilies.

Characterisation of retrotransposon insertion polymorphisms in whole genome sequencing data from individuals with amyotrophic lateral sclerosis

The genetics of an individual is a crucial factor in understanding the risk of developing the neurodegenerative disease amyotrophic lateral sclerosis (ALS). There is still a large proportion of the heritability of ALS, particularly in sporadic cases, to be understood. Among others, active transposable elements drive inter-individual variability, and in humans long interspersed element 1 (LINE1, L1), Alu and SINE-VNTR-Alu (SVA) retrotransposons are a source of polymorphic insertions in the population. We undertook a pilot study to characterise the landscape of non-reference retrotransposon insertion polymorphisms (non-ref RIPs) in 15 control and 15 ALS individuals' whole genomes from Project MinE, an international project to identify potential genetic causes of ALS. The combination of two bioinformatics tools (mobile element locator tool (MELT) and TEBreak) identified on average 1250 Alu, 232 L1 and 77 SVA non-ref RIPs per genome across the 30 analysed. Further PCR validation of individual polymorphic retrotransposon insertions showed a similar level of accuracy for MELT and TEBreak. Our preliminary study did not identify a specific RIP or a significant difference in the total number of non-ref RIPs in ALS compared to control genomes. The use of multiple bioinformatic tools improved the accuracy of non-ref RIP detection and our study highlights the potential importance of studying these elements further in ALS.

MIR retrotransposons link the epigenome and the transcriptome of coding genes in acute myeloid leukemia

DNMT3A and IDH1/2 mutations combinatorically regulate the transcriptome and the epigenome in acute myeloid leukemia; yet the mechanisms of this interplay are unknown. Using a systems approach within topologically associating domains, we find that genes with significant expression-methylation correlations are enriched in signaling and metabolic pathways. The common denominator across these methylation-regulated genes is the density in MIR retrotransposons of their introns. Moreover, a discrete number of CpGs overlapping enhancers are responsible for regulating most of these genes. Established mouse models recapitulate the dependency of MIR-rich genes on the balanced expression of epigenetic modifiers, while projection of leukemic profiles onto normal hematopoiesis ones further consolidates the dependencies of methylation-regulated genes on MIRs. Collectively, MIR elements on genes and enhancers are susceptible to changes in DNA methylation activity and explain the cooperativity of proteins in this pathway in normal and malignant hematopoiesis.

You shall not pass! A Chromatin barrier story in plants

As in other eukaryotes, the plant genome is functionally organized in two mutually exclusive chromatin fractions, a gene-rich and transcriptionally active euchromatin, and a gene-poor, repeat-rich, and transcriptionally silent heterochromatin. In Drosophila and humans, the molecular mechanisms by which euchromatin is preserved from heterochromatin spreading have been extensively studied, leading to the identification of insulator DNA elements and associated chromatin factors (insulator proteins), which form boundaries between chromatin domains with antagonistic features. In contrast, the identity of factors assuring such a barrier function remains largely elusive in plants. Nevertheless, several genomic elements and associated protein factors have recently been shown to regulate the spreading of chromatin marks across their natural boundaries in plants. In this minireview, we focus on recent findings that describe the spreading of chromatin and propose avenues to improve the understanding of how plant chromatin architecture and transitions between different chromatin domains are defined.

Cellular self-organization: An overdrive in Cambrian diversity?

During the early Cambrian period metazoan life forms diverged at an accelerated rate to occupy multiple ecological niches on earth. A variety of explanations have been proposed to address this major evolutionary phenomenon termed the "Cambrian explosion." While most hypotheses address environmental, developmental, and ecological factors that facilitated evolutionary innovations, the biological basis for accelerated emergence of species diversity in the Cambrian period remains largely conjectural. Herein, we posit that morphogenesis by self-organization enables the uncoupling of genomic mutational landscape from phenotypic diversification. Evidence is provided for a two-tiered interpretation of genomic changes in metazoan animals wherein mutations not only impact upon function of individual cells, but also alter the self-organization outcome during developmental morphogenesis. We provide evidence that the morphological impacts of mutations on self-organization could remain repressed if associated with an unmet negative energetic cost. We posit that accelerated morphological diversification in transition to the Cambrian period has occurred by emergence of dormant (i.e., reserved) morphological novelties whose molecular underpinnings were seeded in the Precambrian period.

Impact of Multiple Sclerosis Risk Polymorphism rs7665090 on MANBA Activity, Lysosomal Endocytosis, and Lymphocyte Activation

Deficiencies in Mannosidase β (MANBA) are associated with neurological abnormalities and recurrent infections. The single nucleotide polymorphism located in the 3′UTR of MANBA, rs7665090, was found to be associated with multiple sclerosis (MS) susceptibility. We aimed to study the functional impact of this polymorphism in lymphocytes isolated from MS patients and healthy controls. A total of 152 MS patients and 112 controls were genotyped for rs7665090. MANBA mRNA expression was quantified through qPCR and MANBA enzymatic activity was analyzed. Upon phytohemagglutinin stimulation, immune activation was evaluated by flow cytometry detection of CD69, endocytic function, and metabolic rates with Seahorse XFp Analyzer, and results were stratified by variation in rs7665090. A significantly reduced gene expression (p < 0.0001) and enzymatic activity (p = 0.018) of MANBA were found in lymphocytes of MS patients compared to those of controls. The rs7665090*GG genotype led to a significant β-mannosidase enzymatic deficiency correlated with lysosomal dysfunction, as well as decreased metabolic activation in lymphocytes of MS patients compared to those of rs7665090*GG controls. In contrast, lymphocytes of MS patients and controls carrying the homozygous AA genotype behaved similarly. Our work provides new evidence highlighting the impact of the MS-risk variant, rs7665090, and the role of MANBA in the immunopathology of MS.

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.

TFIIIC-based chromatin insulators through eukaryotic evolution

Eukaryotic chromosomes are divided into domains with distinct structural and functional properties, such as differing levels of chromatin compaction and gene transcription. Domains of relatively compact chromatin and minimal transcription are termed heterochromatic, whereas euchromatin is more open and actively transcribed. Insulators separate these domains and maintain their distinct features. Disruption of insulators can cause diseases such as cancer. Many insulators contain tRNA genes (tDNAs), examples of which have been shown to block the spread of activating or silencing activities. This characteristic of specific tDNAs is conserved through evolution, such that human tDNAs can serve as barriers to the spread of silencing in fission yeast. Here we demonstrate that tDNAs from the methylotrophic fungus Pichia pastoris can function effectively as insulators in distantly-related budding yeast. Key to the function of tDNAs as insulators is TFIIIC, a transcription factor that is also required for their expression. TFIIIC binds additional loci besides tDNAs, some of which have insulator activity. Although the mechanistic basis of TFIIIC-based insulation has been studied extensively in yeast, it is largely uncharacterized in metazoa. Utilising publicly-available genome-wide ChIP-seq data, we consider the extent to which mechanisms conserved from yeast to man may suffice to allow efficient insulation by TFIIIC in the more challenging chromatin environments of metazoa and suggest features that may have been acquired during evolution to cope with new challenges. We demonstrate the widespread presence at human tDNAs of USF1, a transcription factor with well-established barrier activity in vertebrates. We predict that tDNA-based insulators in higher organisms have evolved through incorporation of modules, such as binding sites for factors like USF1 and CTCF that are absent from yeasts, thereby strengthening function and providing opportunities for regulation between cell types.

The Role of Transposable Elements of the Human Genome in Neuronal Function and Pathology

Transposable elements (TEs) have been extensively studied for decades. In recent years, the introduction of whole-genome and whole-transcriptome approaches, as well as single-cell resolution techniques, provided a breakthrough that uncovered TE involvement in host gene expression regulation underlying multiple normal and pathological processes. Of particular interest is increased TE activity in neuronal tissue, and specifically in the hippocampus, that was repeatedly demonstrated in multiple experiments. On the other hand, numerous neuropathologies are associated with TE dysregulation. Here, we provide a comprehensive review of literature about the role of TEs in neurons published over the last three decades. The first chapter of the present review describes known mechanisms of TE interaction with host genomes in general, with the focus on mammalian and human TEs; the second chapter provides examples of TE exaptation in normal neuronal tissue, including TE involvement in neuronal differentiation and plasticity; and the last chapter lists TE-related neuropathologies. We sought to provide specific molecular mechanisms of TE involvement in neuron-specific processes whenever possible; however, in many cases, only phenomenological reports were available. This underscores the importance of further studies in this area.

Emerging roles for endogenous retroviruses in immune epigenetic regulation

In recent years, there has been significant progress toward understanding the transcriptional networks underlying mammalian immune responses, fueled by advances in regulatory genomic technologies. Epigenomic studies profiling immune cells have generated detailed genome-wide maps of regulatory elements that will be key to deciphering the regulatory networks underlying cellular immune responses and autoimmune disorders. Unbiased analyses of these genomic maps have uncovered endogenous retroviruses as an unexpected ally in the regulation of human immune systems. Despite their parasitic origins, studies are finding an increasing number of examples of retroviral sequences having been co-opted for beneficial immune function and regulation by the host cell. Here, we review how endogenous retroviruses have given rise to numerous regulatory elements that shape the epigenetic landscape of host immune responses. We will discuss the implications of these elements on the function, dysfunction, and evolution of innate immunity.

The HPSE Gene Insulator-A Novel Regulatory Element That Affects Heparanase Expression, Stem Cell Mobilization, and the Risk of Acute Graft versus Host Disease

The HPSE gene encodes heparanase (HPSE), a key player in cancer, inflammation, and autoimmunity. We have previously identified a strong HPSE gene enhancer involved in self-regulation of heparanase by negative feedback exerted in a functional rs4693608 single-nucleotide polymorphism (SNP) dependent manner. In the present study, we analyzed the HPSE gene insulator region, located in intron 9 and containing rs4426765, rs28649799, and rs4364254 SNPs. Our results indicate that this region exhibits HPSE regulatory activity. SNP substitutions lead to modulation of a unique DNA-protein complex that affects insulator activity. Analysis of interactions between enhancer and insulator SNPs revealed that rs4693608 has a major effect on HPSE expression and the risk of post-transplantation acute graft versus host disease (GVHD). The C alleles of insulator SNPs rs4364254 and rs4426765 modify the activity of the HPSE enhancer, resulting in altered HPSE expression and increased risk of acute GVHD. Moreover, rs4426765 correlated with HPSE expression in activated mononuclear cells, as well as with CD3 levels and lymphocyte counts following G-CSF mobilization. rs4363084 and rs28649799 were found to be associated with CD34⁺ levels. Our study provides new insight into the mechanism of HPSE gene regulation and its impact on normal and pathological processes in the hematopoietic system.

More than causing (epi)genomic instability: emerging physiological implications of transposable element modulation

Transposable elements (TEs) initially attracted attention because they comprise a major portion of the genomic sequences in plants and animals. TEs may jump around the genome and disrupt both coding genes as well as regulatory sequences to cause disease. Host cells have therefore evolved various epigenetic and functional RNA-mediated mechanisms to mitigate the disruption of genomic integrity by TEs. TE associated sequences therefore acquire the tendencies of attracting various epigenetic modifiers to induce epigenetic alterations that may spread to the neighboring genes. In addition to posting threats for (epi)genome integrity, emerging evidence suggested the physiological importance of endogenous TEs either as cis-acting control elements for controlling gene regulation or as TE-containing functional transcripts that modulate the transcriptome of the host cells. Recent advances in long-reads sequence analysis technologies, bioinformatics and genetic editing tools have enabled the profiling, precise annotation and functional characterization of TEs despite their challenging repetitive nature. The importance of specific TEs in preimplantation embryonic development, germ cell differentiation and meiosis, cell fate determination and in driving species specific differences in mammals will be discussed.

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.

Allelic diversification after transposable element exaptation promoted gsdf as the master sex determining gene of sablefish

Concepts of evolutionary biology suggest that morphological change may occur by rare punctual but rather large changes, or by more steady and gradual transformations. It can therefore be asked whether genetic changes underlying morphological, physiological, and/or behavioral innovations during evolution occur in a punctual manner, whereby a single mutational event has prominent phenotypic consequences, or if many consecutive alterations in the DNA over longer time periods lead to phenotypic divergence. In the marine teleost, sablefish (Anoplopoma fimbria), complementary genomic and genetic studies led to the identification of a sex locus on the Y Chromosome. Further characterization of this locus resulted in identification of the transforming growth factor, beta receptor 1a (tgfbr1a) gene, gonadal somatic cell derived factor (gsdf), as the main candidate for fulfilling the master sex determining (MSD) function. The presence of different X and Y Chromosome copies of this gene indicated that the male heterogametic (XY) system of sex determination in sablefish arose by allelic diversification. The gsdfY gene has a spatio-temporal expression profile characteristic of a male MSD gene. We provide experimental evidence demonstrating a pivotal role of a transposable element (TE) for the divergent function of gsdfY. By insertion within the gsdfY promoter region, this TE generated allelic diversification by bringing cis-regulatory modules that led to transcriptional rewiring and thus creation of a new MSD gene. This points out, for the first time in the scenario of MSD gene evolution by allelic diversification, a single, punctual molecular event in the appearance of a new trigger for male development.

Triangulating molecular evidence to prioritize candidate causal genes at established atopic dermatitis loci

Genome-wide association studies for atopic dermatitis (AD) have identified 25 reproducible loci. We attempt to prioritize candidate causal genes at these loci using extensive molecular resources compiled into a bioinformatics pipeline. We identified a list of 103 molecular resources for AD aetiology, including expression, protein and DNA methylation QTL datasets in skin or immune-relevant tissues which were tested for overlap with GWAS signals. This was combined with functional annotation using regulatory variant prediction, and features such as promoter-enhancer interactions, expression studies and variant fine-mapping. For each gene at each locus, we condensed the evidence into a prioritization score. Across the investigated loci, we detected significant enrichment of genes with adaptive immune regulatory function and epidermal barrier formation among the top prioritized genes. At 8 loci, we were able to prioritize a single candidate gene (IL6R, ADO, PRR5L, IL7R, ETS1, INPP5D, MDM1, TRAF3). In addition, at 6 of the 25 loci, our analysis prioritizes less familiar candidates (SLC22A5, IL2RA, MDM1, DEXI, ADO, STMN3). Our analysis provides support for previously implicated genes at several AD GWAS loci, as well as evidence for plausible additional candidates at others, which may represent potential targets for drug discovery.

Retrotransposons as Drivers of Mammalian Brain Evolution

Retrotransposons, a large and diverse class of transposable elements that are still active in humans, represent a remarkable force of genomic innovation underlying mammalian evolution. Among the features distinguishing mammals from all other vertebrates, the presence of a neocortex with a peculiar neuronal organization, composition and connectivity is perhaps the one that, by affecting the cognitive abilities of mammals, contributed mostly to their evolutionary success. Among mammals, hominids and especially humans display an extraordinarily expanded cortical volume, an enrichment of the repertoire of neural cell types and more elaborate patterns of neuronal connectivity. Retrotransposon-derived sequences have recently been implicated in multiple layers of gene regulation in the brain, from transcriptional and post-transcriptional control to both local and large-scale three-dimensional chromatin organization. Accordingly, an increasing variety of neurodevelopmental and neurodegenerative conditions are being recognized to be associated with retrotransposon dysregulation. We review here a large body of recent studies lending support to the idea that retrotransposon-dependent evolutionary novelties were crucial for the emergence of mammalian, primate and human peculiarities of brain morphology and function.

Investigating the Potential Roles of SINEs in the Human Genome

Short interspersed nuclear elements (SINEs) are nonautonomous retrotransposons that occupy approximately 13% of the human genome. They are transcribed by RNA polymerase III and can be retrotranscribed and inserted back into the genome with the help of other autonomous retroelements. Because they are preferentially located close to or within gene-rich regions, they can regulate gene expression by various mechanisms that act at both the DNA and the RNA levels. In this review, we summarize recent findings on the involvement of SINEs in different types of gene regulation and discuss the potential regulatory functions of SINEs that are in close proximity to genes, Pol III-transcribed SINE RNAs, and embedded SINE sequences within Pol II-transcribed genes in the human genome. These discoveries illustrate how the human genome has exapted some SINEs into functional regulatory elements.

The transcriptional trajectories of pluripotency and differentiation comprise genes with antithetical architecture and repetitive-element content

Background

Extensive molecular differences exist between proliferative and differentiated cells. Here, we conduct a meta-analysis of publicly available transcriptomic datasets from preimplantation and differentiation stages examining the architectural properties and content of genes whose abundance changes significantly across developmental time points.

Results

Analysis of preimplantation embryos from human and mouse showed that short genes whose introns are enriched in Alu (human) and B (mouse) elements, respectively, have higher abundance in the blastocyst compared to the zygote. These highly expressed genes encode ribosomal proteins or metabolic enzymes. On the other hand, long genes whose introns are depleted in repetitive elements have lower abundance in the blastocyst and include genes from signaling pathways. Additionally, the sequences of the genes that are differentially expressed between the blastocyst and the zygote contain distinct collections of pyknon motifs that differ between up- and down-regulated genes. Further examination of the genes that participate in the stem cell-specific protein interaction network shows that their introns are short and enriched in Alu (human) and B (mouse) elements. As organogenesis progresses, in both human and mouse, we find that the primarily short and repeat-rich expressed genes make way for primarily longer, repeat-poor genes. With that in mind, we used a machine learning-based approach to identify gene signatures able to classify human adult tissues: we find that the most discriminatory genes comprising these signatures have long introns that are repeat-poor and include transcription factors and signaling-cascade genes. The introns of widely expressed genes across human tissues, on the other hand, are short and repeat-rich, and coincide with those with the highest expression at the blastocyst stage.

Conclusions

Protein-coding genes that are characteristic of each trajectory, i.e., proliferation/pluripotency or differentiation, exhibit antithetical biases in their intronic and exonic lengths and in their repetitive-element content. While the respective human and mouse gene signatures are functionally and evolutionarily conserved, their introns and exons are enriched or depleted in organism-specific repetitive elements. We posit that these organism-specific repetitive sequences found in exons and introns are used to effect the corresponding genes’ regulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-020-00928-8.

Transposable elements: a jump toward the future of expression vectors

Expression vectors (EVs) are artificial nucleic acid molecules with a modular structure that allows for the transcription of DNA sequences of interest in either cellular or cell-free environments. These vectors have emerged as cross-disciplinary tools with multiple applications in an expanding Life Sciences market. The cis-regulatory sequences (CRSs) that control the transcription in EVs are typically sourced from either viruses or from characterized genes. However, the recent advancement in transposable elements (TEs) technology provides attractive alternatives that may enable a significant improvement in the design of EVs. Commonly known as "jumping genes," due to their ability to move between genetic loci, TEs are constitutive components of both eukaryotic and prokaryotic genomes. TEs harbor native CRSs that allow the regulated transcription of transposition-related genes. However, some TE-related CRSs display striking characteristics, which provides the opportunity to reconsider TEs as lead actors in the design of EVs. In this article, we provide a synopsis of the transcriptional control elements commonly found in EVs together with an extensive discussion of their advantages and limitations. We also highlight the latest findings that may allow for the implementation of TE-derived sequences in the EVs feasible, possibly improving existing vectors. By introducing this new concept of TEs as a source of regulatory sequences, we aim to stimulate a profitable discussion of the potential advantages and benefits of developing a new generation of EVs based on the use of TE-derived control sequences.

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.

Transposable element-derived sequences in vertebrate development

Transposable elements (TEs) are major components of all vertebrate genomes that can cause deleterious insertions and genomic instability. However, depending on the specific genomic context of their insertion site, TE sequences can sometimes get positively selected, leading to what are called "exaptation" events. TE sequence exaptation constitutes an important source of novelties for gene, genome and organism evolution, giving rise to new regulatory sequences, protein-coding exons/genes and non-coding RNAs, which can play various roles beneficial to the host. In this review, we focus on the development of vertebrates, which present many derived traits such as bones, adaptive immunity and a complex brain. We illustrate how TE-derived sequences have given rise to developmental innovations in vertebrates and how they thereby contributed to the evolutionary success of this lineage.

Retroelement-derived RNA and its role in the brain

Comprising ~40% of the human genome, retroelements are mobile genetic elements which are transcribed into RNA, then reverse-transcribed into DNA and inserted into a new site in the genome. Retroelements are referred to as "genetic parasites", residing among host genes and relying on host machinery for transcription and evolutionary propagation. The healthy brain has the highest expression of retroelement-derived sequences compared to other somatic tissue, which leads to the question: how does retroelement-derived RNA influence human traits and cellular states? While the functional importance of upregulating retroelement expression in the brain is an active area of research, RNA species derived from retroelements influence both self- and host gene expression by contributing to chromatin remodeling, alternative splicing, somatic mosaicism and translational repression. Here, we review the emerging evidence that the functional importance of RNA derived from retroelements is multifaceted. Retroelements can influence organismal states through the seeding of epigenetic states in chromatin, the production of structured RNA and even catalytically active ribozymes, the generation of cytoplasmic ssDNA and RNA/DNA hybrids, the production of viral-like proteins, and the generation of somatic mutations. Comparative sequencing suggests that retroelements can contribute to intraspecies variation through these mechanisms to alter transcript identity and abundance. In humans, an increasing number of neurodevelopmental and neurodegenerative conditions are associated with dysregulated retroelements, including Aicardi-Goutieres syndrome (AGS), Rett syndrome (RTT), Amyotrophic Lateral Sclerosis (ALS), Alzheimer's disease (AD), multiple sclerosis (MS), schizophrenia (SZ), and aging. Taken together, these concepts suggest a larger functional role for RNA derived from retroelements. This review aims to define retroelement-derived RNA, discuss how it impacts the mammalian genome, as well as summarize data supporting phenotypic consequences of this unique RNA subset in the brain.

Host Gene Regulation by Transposable Elements: The New, the Old and the Ugly

The human genome has been under selective pressure to evolve in response to emerging pathogens and other environmental challenges. Genome evolution includes the acquisition of new genes or new isoforms of genes and changes to gene expression patterns. One source of genome innovation is from transposable elements (TEs), which carry their own promoters, enhancers and open reading frames and can act as ‘controlling elements’ for our own genes. TEs include LINE-1 elements, which can retrotranspose intracellularly and endogenous retroviruses (ERVs) that represent remnants of past retroviral germline infections. Although once pathogens, ERVs also represent an enticing source of incoming genetic material that the host can then repurpose. ERVs and other TEs have coevolved with host genes for millions of years, which has allowed them to become embedded within essential gene expression programmes. Intriguingly, these host genes are often subject to the same epigenetic control mechanisms that evolved to combat the TEs that now regulate them. Here, we illustrate the breadth of host gene regulation through TEs by focusing on examples of young (The New), ancient (The Old), and disease-causing (The Ugly) TE integrants.

Boundary sequences flanking the mouse tyrosinase locus ensure faithful pattern of gene expression

Control of gene expression is dictated by cell-type specific regulatory sequences that physically organize the structure of chromatin, including promoters, enhancers and insulators. While promoters and enhancers convey cell-type specific activating signals, insulators prevent the cross-talk of regulatory elements within adjacent loci and safeguard the specificity of action of promoters and enhancers towards their targets in a tissue specific manner. Using the mouse tyrosinase (Tyr) locus as an experimental model, a gene whose mutations are associated with albinism, we described the chromatin structure in cells at two distinct transcriptional states. Guided by chromatin structure, through the use of Chromosome Conformation Capture (3C), we identified sequences at the 5' and 3' boundaries of this mammalian gene that function as enhancers and insulators. By CRISPR/Cas9-mediated chromosomal deletion, we dissected the functions of these two regulatory elements in vivo in the mouse, at the endogenous chromosomal context, and proved their mechanistic role as genomic insulators, shielding the Tyr locus from the expression patterns of adjacent genes.

The tuatara genome reveals ancient features of amniote evolution

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2-is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

An atlas of transposable element-derived alternative splicing in cancer

Transposable element (TE)-derived sequences comprise more than half of the human genome, and their presence has been documented to alter gene expression in a number of different ways, including the generation of alternatively spliced transcript isoforms. Alternative splicing has been associated with tumorigenesis for a number of different cancers. The objective of this study was to broadly characterize the role of human TEs in generating alternatively spliced transcript isoforms in cancer. To do so, we screened for the presence of TE-derived sequences co-located with alternative splice sites that are differentially used in normal versus cancer tissues. We analysed a comprehensive set of alternative splice variants characterized for 614 matched normal-tumour tissue pairs across 13 cancer types, resulting in the discovery of 4820 TE-generated alternative splice events distributed among 723 cancer-associated genes. Short interspersed nuclear elements (Alu) and long interspersed nuclear elements (L1) were found to contribute the majority of TE-generated alternative splice sites in cancer genes. A number of cancer-associated genes, including MYH11, WHSC1 and CANT1, were shown to have overexpressed TE-derived isoforms across a range of cancer types. TE-derived isoforms were also linked to cancer-specific fusion transcripts, suggesting a novel mechanism for the generation of transcriptome diversity via trans-splicing mediated by dispersed TE repeats. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Biased visibility in Hi-C datasets marks dynamically regulated condensed and decondensed chromatin states genome-wide

Background

Proximity ligation based techniques, like Hi-C, involve restriction digestion followed by ligation of formaldehyde cross-linked chromatin. Distinct chromatin states can impact the restriction digestion, and hence the visibility in the contact maps, of engaged loci. Yet, the extent and the potential impact of digestion bias remain obscure and under-appreciated in the literature.

Results

Through analysis of 45 Hi-C datasets, lamina-associated domains (LADs), inactive X-chromosome in mammals, and polytene bands in fly, we first established that the DNA in condensed chromatin had lesser accessibility to restriction endonucleases used in Hi-C as compared to that in decondensed chromatin. The observed bias was independent of known systematic biases, was not appropriately corrected by existing computational methods, and needed an additional optimization step. We then repurposed this bias to identify novel condensed domains outside LADs, which were bordered by insulators and were dynamically associated with the polycomb mediated epigenetic and transcriptional states during development.

Conclusions

Our observations suggest that the corrected one-dimensional read counts of existing Hi-C datasets can be reliably repurposed to study the gene-regulatory dynamics associated with chromatin condensation and decondensation, and that the existing Hi-C datasets should be interpreted with cautions.

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.

Transposable elements as genetic accelerators of evolution: contribution to genome size, gene regulatory network rewiring and morphological innovation

In the current era, as a growing number of genome sequence assemblies have been reported in animals, an in-depth analysis of transposable elements (TEs) is one of the most fundamental and essential studies for evolutionary genomics. Although TEs have, in general, been regarded as non-functional junk/selfish DNA, parasitic elements or harmful mutagens, studies have revealed that TEs have had a substantial and sometimes beneficial impact on host genomes in several ways. First, TEs are themselves diverse and thus provide lineage-specific characteristics to the genomes. Second, because TEs constitute a substantial fraction of animal genomes, they are a major contributing factor to evolutionary changes in genome size and composition. Third, host organisms have co-opted many repetitive sequences as genes, cis-regulatory elements and chromatin domain boundaries, which alter gene regulatory networks and in addition are partly involved in morphological evolution, as has been well documented in mammals. Here, I review the impact of TEs on various aspects of the genome, such as genome size and diversity in animals, as well as the evolution of gene networks and genome architecture in mammals. Given that a number of TE families probably remain to be discovered in many non-model organisms, unknown TEs may have contributed to gene networks in a much wider variety of animals than considered previously.

The epigenomic landscape of transposable elements across normal human development and anatomy

Transposable elements (TEs) have deposited functional regulatory elements throughout the human genome. Although most are silenced, certain TEs have been co-opted by the host. However, a comprehensive, multidimensional picture of the contribution of TEs to normal human gene regulation is still lacking. Here, we quantify the epigenomic status of TEs across human anatomy and development using data from the Roadmap Epigenomics Project. We find that TEs encompass a quarter of the human regulatory epigenome, and 47% of elements can be in an active regulatory state. We demonstrate that SINEs are enriched relative to other classes for active and transcribed marks, that TEs encompass a higher proportion of enhancer states in the hematopoietic lineage, and that DNA methylation of Alu elements decreases with age, corresponding with a loss of CpG islands. Finally, we identify TEs that may perform an evolutionarily conserved regulatory function, providing a systematic profile of TE activity in normal human tissue.

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.

Genome-Wide Identification of Microsatellites and Transposable Elements in the Dromedary Camel Genome Using Whole-Genome Sequencing Data

Transposable elements (TEs) along with simple sequence repeats (SSRs) are prevalent in eukaryotic genome, especially in mammals. Repetitive sequences form approximately one-third of the camelid genomes, so study on this part of genome can be helpful in providing deeper information from the genome and its evolutionary path. Here, in order to improve our understanding regarding the camel genome architecture, the whole genome of the two dromedaries (Yazdi and Trodi camels) was sequenced. Totally, 92- and 84.3-Gb sequence data were obtained and assembled to 137,772 and 149,997 contigs with a N50 length of 54,626 and 54,031 bp in Yazdi and Trodi camels, respectively. Results showed that 30.58% of Yazdi camel genome and 30.50% of Trodi camel genome were covered by TEs. Contrary to the observed results in the genomes of cattle, sheep, horse, and pig, no endogenous retrovirus-K (ERVK) elements were found in the camel genome. Distribution pattern of DNA transposons in the genomes of dromedary, Bactrian, and cattle was similar in contrast with LINE, SINE, and long terminal repeat (LTR) families. Elements like RTE-BovB belonging to LINEs family in cattle and sheep genomes are dramatically higher than genome of dromedary. However, LINE1 (L1) and LINE2 (L2) elements cover higher percentage of LINE family in dromedary genome compared to genome of cattle. Also, 540,133 and 539,409 microsatellites were identified from the assembled contigs of Yazdi and Trodi dromedary camels, respectively. In both samples, di-(393,196) and tri-(65,313) nucleotide repeats contributed to about 42.5% of the microsatellites. The findings of the present study revealed that non-repetitive content of mammalian genomes is approximately similar. Results showed that 9.1 Mb (0.47% of whole assembled genome) of Iranian dromedary's genome length is made up of SSRs. Annotation of repetitive content of Iranian dromedary camel genome revealed that 9,068 and 11,544 genes contain different types of TEs and SSRs, respectively. SSR markers identified in the present study can be used as a valuable resource for genetic diversity investigations and marker-assisted selection (MAS) in camel-breeding programs.