Genomic landscape of human, bat, and ex vivo DNA transposon integrations

Comprehensive Identification and Characterization of HML-9 Group in Chimpanzee Genome

Endogenous retroviruses (ERVs) are related to long terminal repeat (LTR) retrotransposons, comprising gene sequences of exogenous retroviruses integrated into the host genome and inherited according to Mendelian law. They are considered to have contributed greatly to the evolution of host genome structure and function. We previously characterized HERV-K HML-9 in the human genome. However, the biological function of this type of element in the genome of the chimpanzee, which is the closest living relative of humans, largely remains elusive. Therefore, the current study aims to characterize HML-9 in the chimpanzee genome and to compare the results with those in the human genome. Firstly, we report the distribution and genetic structural characterization of the 26 proviral elements and 38 solo LTR elements of HML-9 in the chimpanzee genome. The results showed that the distribution of these elements displayed a non-random integration pattern, and only six elements maintained a relatively complete structure. Then, we analyze their phylogeny and reveal that the identified elements all cluster together with HML-9 references and with those identified in the human genome. The HML-9 integration time was estimated based on the 2-LTR approach, and the results showed that HML-9 elements were integrated into the chimpanzee genome between 14 and 36 million years ago and into the human genome between 18 and 49 mya. In addition, conserved motifs, cis-regulatory regions, and enriched PBS sequence features in the chimpanzee genome were predicted based on bioinformatics. The results show that pathways significantly enriched for ERV LTR-regulated genes found in the chimpanzee genome are closely associated with disease development, including neurological and neurodevelopmental psychiatric disorders. In summary, the identification, characterization, and genomics of HML-9 presented here not only contribute to our understanding of the role of ERVs in primate evolution but also to our understanding of their biofunctional significance.

Transposable elements and their role in aging

Transposable elements (TEs) are an important part of eukaryotic genomes. The role of somatic transposition in aging, carcinogenesis, and other age-related diseases has been determined. This review discusses the fundamental properties of TEs and their complex interactions with cellular processes, which are crucial for understanding the diverse effects of their activity on the genetics and epigenetics of the organism. The interactions of TEs with recombination, replication, repair, and chromosomal regulation; the ability of TEs to maintain a balance between their own activity and repression, the involvement of TEs in the creation of new or alternative genes, the expression of coding/non-coding RNA, and the role in DNA damage and modification of regulatory networks are reviewed. The contribution of the derepressed TEs to age-dependent effects in individual cells/tissues in different organisms was assessed. Conflicting information about TE activity under stress as well as theories of aging mechanisms related to TEs is discussed. On the one hand, transposition activity in response to stressors can lead to organisms acquiring adaptive innovations of great importance for evolution at the population level. On the other hand, the TE expression can cause decreased longevity and stress tolerance at the individual level. The specific features of TE effects on aging processes in germline and soma and the ways of their regulation in cells are highlighted. Recent results considering somatic mutations in normal human and animal tissues are indicated, with the emphasis on their possible functional consequences. In the context of aging, the correlation between somatic TE activation and age-related changes in the number of proteins required for heterochromatin maintenance and longevity regulation was analyzed. One of the original features of this review is a discussion of not only effects based on the TEs insertions and the associated consequences for the germline cell dynamics and somatic genome, but also the differences between transposon- and retrotransposon-mediated structural genome changes and possible phenotypic characteristics associated with aging and various age-related pathologies. Based on the analysis of published data, a hypothesis about the influence of the species-specific features of number, composition, and distribution of TEs on aging dynamics of different animal genomes was formulated.

Noncanonical DNA structures are drivers of genome evolution

In addition to the canonical right-handed double helix, other DNA structures, termed 'non-B DNA', can form in the genomes across the tree of life. Non-B DNA regulates multiple cellular processes, including replication and transcription, yet its presence is associated with elevated mutagenicity and genome instability. These discordant cellular roles fuel the enormous potential of non-B DNA to drive genomic and phenotypic evolution. Here we discuss recent studies establishing non-B DNA structures as novel functional elements subject to natural selection, affecting evolution of transposable elements (TEs), and specifying centromeres. By highlighting the contributions of non-B DNA to repeated evolution and adaptation to changing environments, we conclude that evolutionary analyses should include a perspective of not only DNA sequence, but also its structure.

Non-B DNA: a major contributor to small- and large-scale variation in nucleotide substitution frequencies across the genome

Approximately 13% of the human genome can fold into non-canonical (non-B) DNA structures (e.g. G-quadruplexes, Z-DNA, etc.), which have been implicated in vital cellular processes. Non-B DNA also hinders replication, increasing errors and facilitating mutagenesis, yet its contribution to genome-wide variation in mutation rates remains unexplored. Here, we conducted a comprehensive analysis of nucleotide substitution frequencies at non-B DNA loci within noncoding, non-repetitive genome regions, their ±2 kb flanking regions, and 1-Megabase windows, using human-orangutan divergence and human single-nucleotide polymorphisms. Functional data analysis at single-base resolution demonstrated that substitution frequencies are usually elevated at non-B DNA, with patterns specific to each non-B DNA type. Mirror, direct and inverted repeats have higher substitution frequencies in spacers than in repeat arms, whereas G-quadruplexes, particularly stable ones, have higher substitution frequencies in loops than in stems. Several non-B DNA types also affect substitution frequencies in their flanking regions. Finally, non-B DNA explains more variation than any other predictor in multiple regression models for diversity or divergence at 1-Megabase scale. Thus, non-B DNA substantially contributes to variation in substitution frequencies at small and large scales. Our results highlight the role of non-B DNA in germline mutagenesis with implications to evolution and genetic diseases.

Endogenous Retroviruses Activity as a Molecular Signature of Neurodevelopmental Disorders

Human endogenous retroviruses (HERVs) are genetic elements resulting from relics of ancestral infection of germline cells, now recognized as cofactors in the etiology of several complex diseases. Here we present a review of findings supporting the role of the abnormal HERVs activity in neurodevelopmental disorders. The derailment of brain development underlies numerous neuropsychiatric conditions, likely starting during prenatal life and carrying on during subsequent maturation of the brain. Autism spectrum disorders, attention deficit hyperactivity disorders, and schizophrenia are neurodevelopmental disorders that arise clinically during early childhood or adolescence, currently attributed to the interplay among genetic vulnerability, environmental risk factors, and maternal immune activation. The role of HERVs in human embryogenesis, their intrinsic responsiveness to external stimuli, and the interaction with the immune system support the involvement of HERVs in the derailed neurodevelopmental process. Although definitive proofs that HERVs are involved in neurobehavioral alterations are still lacking, both preclinical models and human studies indicate that the abnormal expression of ERVs could represent a neurodevelopmental disorders-associated biological trait in affected individuals and their parents.

Human endogenous retroviruses role in cancer cell stemness

Cancer incidence and mortality, metastasis, drug resistance and recurrence are still the critical issues of oncological diseases. In this scenario, increasing scientific evidences demonstrate that the activation of human endogenous retroviruses (HERVs) is involved in the aggressiveness of tumors such as melanoma, breast, germ cell, renal, ovarian, liver and haematological cancers. In their dynamic regulation, HERVs have also proved to be important determinants of pluripotency in human embryonic stem cells (ESC) and of the reprogramming process of induced pluripotent stem cells (iPSCs). In many types of tumors, essential characteristics of aggressiveness have been associated with the achievement of stemness features, often accompanied with the identification of defined subpopulations, termed cancer stem cells (CSCs), which possess stem cell-like properties and sustain tumorigenesis. Indeed, CSCs show high self-renewal capacity with a peculiar potential in tumor initiation, progression, metastasis, heterogeneity, recurrence, radiotherapy and drug resistance. However, HERVs role in CSCs biology is still not fully elucidated. In this regard, CD133 is a widely recognized marker of CSCs, and our group demonstrated, for the first time, the requirement of HERV-K activation to expand and maintain a CD133+ melanoma cell subpopulation with stemness features in response to microenvironmental modifications. The review will discuss HERVs expression as cancer hallmark, with particular focus on their role in the regulation of cancer stemness features and the potential involvement as targets for therapy.

Coevolution between transposable elements and recombination

One of the most striking patterns of genome structure is the tight, typically negative, association between transposable elements (TEs) and meiotic recombination rates. While this is a highly recurring feature of eukaryotic genomes, the mechanisms driving correlations between TEs and recombination remain poorly understood, and distinguishing cause versus effect is challenging. Here, we review the evidence for a relation between TEs and recombination, and discuss the underlying evolutionary forces. Evidence to date suggests that overall TE densities correlate negatively with recombination, but the strength of this correlation varies across element types, and the pattern can be reversed. Results suggest that heterogeneity in the strength of selection against ectopic recombination and gene disruption can drive TE accumulation in regions of low recombination, but there is also strong evidence that the regulation of TEs can influence local recombination rates. We hypothesize that TE insertion polymorphism may be important in driving within-species variation in recombination rates in surrounding genomic regions. Furthermore, the interaction between TEs and recombination may create positive feedback, whereby TE accumulation in non-recombining regions contributes to the spread of recombination suppression. Further investigation of the coevolution between recombination and TEs has important implications for our understanding of the evolution of recombination rates and genome structure.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Human Endogenous Retroviruses in Neurological Diseases

The causes of multiple sclerosis and amyotrophic lateral sclerosis have long remained elusive. A new category of pathogenic components, normally dormant within human genomes, has been identified: human endogenous retroviruses (HERVs). These represent ∼8% of the human genome, and environmental factors have reproducibly been shown to trigger their expression. The resulting production of envelope (Env) proteins from HERV-W and HERV-K appears to engage pathophysiological pathways leading to the pathognomonic features of MS and ALS, respectively. Pathogenic HERV elements may thus provide a missing link in understanding these complex diseases. Moreover, their neutralization may represent a promising strategy to establish novel and more powerful therapeutic approaches.

Profiling of Human Molecular Pathways Affected by Retrotransposons at the Level of Regulation by Transcription Factor Proteins

Endogenous retroviruses and retrotransposons also termed retroelements (REs) are mobile genetic elements that were active until recently in human genome evolution. REs regulate gene expression by actively reshaping chromatin structure or by directly providing transcription factor binding sites (TFBSs). We aimed to identify molecular processes most deeply impacted by the REs in human cells at the level of TFBS regulation. By using ENCODE data, we identified ~2 million TFBS overlapping with putatively regulation-competent human REs located in 5-kb gene promoter neighborhood (~17% of all TFBS in promoter neighborhoods; ~9% of all RE-linked TFBS). Most of REs hosting TFBS were highly diverged repeats, and for the evolutionary young (0–8% diverged) elements we identified only ~7% of all RE-linked TFBS. The gene-specific distributions of RE-linked TFBS generally correlated with the distributions for all TFBS. However, several groups of molecular processes were highly enriched in the RE-linked TFBS regulation. They were strongly connected with the immunity and response to pathogens, with the negative regulation of gene transcription, ubiquitination, and protein degradation, extracellular matrix organization, regulation of STAT signaling, fatty acids metabolism, regulation of GTPase activity, protein targeting to Golgi, regulation of cell division and differentiation, development and functioning of perception organs and reproductive system. By contrast, the processes most weakly affected by the REs were linked with the conservative aspects of embryo development. We also identified differences in the regulation features by the younger and older fractions of the REs. The regulation by the older fraction of the REs was linked mainly with the immunity, cell adhesion, cAMP, IGF1R, Notch, Wnt, and integrin signaling, neuronal development, chondroitin sulfate and heparin metabolism, and endocytosis. The younger REs regulate other aspects of immunity, cell cycle progression and apoptosis, PDGF, TGF beta, EGFR, and p38 signaling, transcriptional repression, structure of nuclear lumen, catabolism of phospholipids, and heterocyclic molecules, insulin and AMPK signaling, retrograde Golgi-ER transport, and estrogen signaling. The immunity-linked pathways were highly represented in both categories, but their functional roles were different and did not overlap. Our results point to the most quickly evolving molecular pathways in the recent and ancient evolution of human genome.

Gonadal morphogenesis and establishment of the germline in the phyllostomid bat Sturnira lilium

In vertebrates such as the mouse and the human, primordial germ cells (PGCs) arise at the base of the allantois and are carried to the epithelium of the posterior intestine, to later migrate to the primordial gonads. In the case of bats, almost nothing is known about this process. To clarify the dynamics of PGCs during gonadal morphogenesis in the phyllostomid bat Sturnira lilium, the proteins for the Ddx4, Sox9 and Mis genes were detected in the gonads of embryos at different stages of development. We identified 15 stages (St) of embryonic development in Sturnira lilium. We found that the formation of the genital ridge and the establishment of the undifferentiated gonad take place between stages 11 and 14. The onset of morphological differentiation in the gonad is first detected in the male gonads at St17. The first PGCs in meiosis are detected in the ovary at St19, whereas in the testicles, the PGCs were in mitotic arrest. Structural changes leading to testicular and ovarian development in Sturnira lilium are observed to be similar to those described for the mouse; however, differences will be established concerning the time taken for these processes to occur.

Integration site selection by retroviruses and transposable elements in eukaryotes

Transposable elements and retroviruses are found in most genomes, can be pathogenic and are widely used as gene-delivery and functional genomics tools. Exploring whether these genetic elements target specific genomic sites for integration and how this preference is achieved is crucial to our understanding of genome evolution, somatic genome plasticity in cancer and ageing, host-parasite interactions and genome engineering applications. High-throughput profiling of integration sites by next-generation sequencing, combined with large-scale genomic data mining and cellular or biochemical approaches, has revealed that the insertions are usually non-random. The DNA sequence, chromatin and nuclear context, and cellular proteins cooperate in guiding integration in eukaryotic genomes, leading to a remarkable diversity of insertion site distribution and evolutionary strategies.

Chromatin states shape insertion profiles of the piggyBac, Tol2 and Sleeping Beauty transposons and murine leukemia virus

DNA transposons and retroviruses are versatile tools in functional genomics and gene therapy. To facilitate their application, we conducted a genome-wide insertion site profiling of the piggyBac (PB), Tol2 and Sleeping Beauty (SB) transposons and the murine leukemia virus (MLV) in mouse embryonic stem cells (ESCs). PB and MLV preferred highly expressed genes, whereas Tol2 and SB preferred weakly expressed genes. However, correlations with DNase I hypersensitive sites were different for all vectors, indicating that chromatin accessibility is not the sole determinant. Therefore, we analysed various chromatin states. PB and MLV highly correlated with Cohesin, Mediator and ESC-specific transcription factors. Notably, CTCF sites were correlated with PB but not with MLV, suggesting MLV prefers smaller promoter-enhancer loops, whereas PB insertion encompasses larger chromatin loops termed topologically associating domains. Tol2 also correlated with Cohesin and CTCF. However, correlations with ESC-specific transcription factors were weaker, suggesting that Tol2 prefers transcriptionally weak chromatin loops. Consistently, Tol2 insertions were associated with bivalent histone modifications characteristic of silent and inducible loci. SB showed minimum preference to all chromatin states, suggesting the least adverse effect on adjacent genes. These results will be useful for vector selection for various applications.

Hsp90 and Physiological Stress Are Linked to Autonomous Transposon Mobility and Heritable Genetic Change in Nematodes

Transposable elements (TEs) have been recognized as potentially powerful drivers of genomic evolutionary change, but factors affecting their mobility and regulation remain poorly understood. Chaperones such as Hsp90 buffer environmental perturbations by regulating protein conformation, but are also part of the PIWI-interacting RNA pathway, which regulates genomic instability arising from mobile TEs in the germline. Stress-induced mutagenesis from TE movement could thus arise from functional trade-offs in the dual roles of Hsp90. We examined the functional constraints of Hsp90 and its role as a regulator of TE mobility by exposing nematodes (Caenorhabditis elegans and Caenorhabditis briggsae) to environmental stress, with and without RNAi-induced silencing of Hsp90. TE excision frequency increased with environmental stress intensity at multiple loci in several strains of each species. These effects were compounded by RNAi-induced knockdown of Hsp90. Mutation frequencies at the unc-22 marker gene in the offspring of animals exposed to environmental stress and Hsp90 RNAi mirrored excision frequency in response to these treatments. Our results support a role for Hsp90 in the suppression of TE mobility, and demonstrate that that the Hsp90 regulatory pathway can be overwhelmed with moderate environmental stress. By compromising genomic stability in germline cells, environmentally induced mutations arising from TE mobility and insertion can have permanent and heritable effects on both the phenotype and the genotype of subsequent generations.

TE-domestication and horizontal transfer in a putative Nef-AP1mu mimic of HLA-A cytoplasmic domain re-trafficking

Genes of the major histocompatibility complex (MHC; also called HLA in human) are polymorphic elements in the genomes of sharks to humans. Class-I and class-II MHC loci appear responsible for much of the genetic linkage to myriad disease states via the capacity to bind short (~8-15 a.a.) peptides of a given pathogen's proteome, or in some cases, the altered proteomes of cancerous cells, and even (in autoimmunity) certain nominal 'self' peptides (Janeway, 2004).(1) Unfortunately, little is known about how the canonical structure of the MHC-I/-II peptide-presenting gene evolved, particularly since beyond ~500 Mya (sharks) no paralogs exist.(2,3) We previously reported that HLA-A isotype alleles with the α1-helix, R65 motif, are wide-spread in phylogeny, but that the α 2-helix, H151R motif, has apparently segregated out of most species. Surprisingly, an uncharacterized orf in T. syrichta (Loc-103275158) encoded R151, but within a truncated A-23 like gene containing 5'- and 3'- footprints of the transposon (TE), tigger-1; the extant tarsier A-23 allele is totally missing exon-3 and part-of exon-4; together, suggesting TE-mediated inactivation of an intact/ancestral A-23 allele (Murray, 2015a).(4) The unique Loc-103275158 orf encodes a putative 15-exon transcript with no apparent paralogs throughout phylogeny. However, an HLA-A11 like gene in M. leucophaeus with a shortened C-terminal domain, and an HLA-A like orf in C. atys with two linked α1/α2/α3 domains, both contain a second transmembrane segment, which is conserved in Loc-103275158. Thus, we could model the putative protein with its Nef-like tail domain docked to its MHC-I like α3 domain (i.e., on the same side of a membrane). This modeled tertiary structure is strikingly similar to the solved structure of the Nef:MHC-I CD:AP1mu transporter (Jia, 2012).(5) Nef:AP1mu binds the CD of MHC-I in trafficking MHC-I away from the trans-golgi and into the endocytic pathway in HIV-1 infected cells. The CD loop of the Loc-103275158 provisional protein conserved the nominal MHC-I CD tyrosine phosphorylation site, and it has an N-terminal SH3 domain that we docked in one conformation to its internal Nef-like domain. Here, we suggest that phosphorylation of the protein's CD-loop signals an exchange between the internal Nef-like domain and a lentiviral-Nef for binding the N-terminal SH3 domain - freeing the Nef-like domain to bind MHC-I CD. Since the 5'-tigger sequence encodes part of the pseudo α1/α2 MHC-I domain, and the 3'-tigger part of the Nef-like domain, we speculate that transposition proceeded phylogenetically disparate horizontal transfers, involving adjacent 5'- and 3'- parasitic footprints, which we also found in the Loc-103275158 orf.

Integration and Fixation Preferences of Human and Mouse Endogenous Retroviruses Uncovered with Functional Data Analysis

Endogenous retroviruses (ERVs), the remnants of retroviral infections in the germ line, occupy ~8% and ~10% of the human and mouse genomes, respectively, and affect their structure, evolution, and function. Yet we still have a limited understanding of how the genomic landscape influences integration and fixation of ERVs. Here we conducted a genome-wide study of the most recently active ERVs in the human and mouse genome. We investigated 826 fixed and 1,065 in vitro HERV-Ks in human, and 1,624 fixed and 242 polymorphic ETns, as well as 3,964 fixed and 1,986 polymorphic IAPs, in mouse. We quantitated >40 human and mouse genomic features (e.g., non-B DNA structure, recombination rates, and histone modifications) in ±32 kb of these ERVs' integration sites and in control regions, and analyzed them using Functional Data Analysis (FDA) methodology. In one of the first applications of FDA in genomics, we identified genomic scales and locations at which these features display their influence, and how they work in concert, to provide signals essential for integration and fixation of ERVs. The investigation of ERVs of different evolutionary ages (young in vitro and polymorphic ERVs, older fixed ERVs) allowed us to disentangle integration vs. fixation preferences. As a result of these analyses, we built a comprehensive model explaining the uneven distribution of ERVs along the genome. We found that ERVs integrate in late-replicating AT-rich regions with abundant microsatellites, mirror repeats, and repressive histone marks. Regions favoring fixation are depleted of genes and evolutionarily conserved elements, and have low recombination rates, reflecting the effects of purifying selection and ectopic recombination removing ERVs from the genome. In addition to providing these biological insights, our study demonstrates the power of exploiting multiple scales and localization with FDA. These powerful techniques are expected to be applicable to many other genomic investigations.

The effects of chromatin organization on variation in mutation rates in the genome

The variation in local rates of mutations can affect both the evolution of genes and their function in normal and cancer cells. Deciphering the molecular determinants of this variation will be aided by the elucidation of distinct types of mutations, as they differ in regional preferences and in associations with genomic features. Chromatin organization contributes to regional variation in mutation rates, but its contribution differs among mutation types. In both germline and somatic mutations, base substitutions are more abundant in regions of closed chromatin, perhaps reflecting error accumulation late in replication. By contrast, a distinctive mutational state with very high levels of insertions and deletions (indels) and substitutions is enriched in regions of open chromatin. These associations indicate an intricate interplay between the nucleotide sequence of DNA and its dynamic packaging into chromatin, and have important implications for current biomedical research. This Review focuses on recent studies showing associations between chromatin state and mutation rates, including pairwise and multivariate investigations of germline and somatic (particularly cancer) mutations.