Transposable Elements and DNA Methylation Create in Embryonic Stem Cells Human-Specific Regulatory Sequences Associated with Distal Enhancers and Noncoding RNAs

Identification of novel cattle (Bos taurus) genes and biological insights of their function in pre-implantation embryo development

BACKGROUND

Appropriate regulation of genes expressed in oocytes and embryos is essential for acquisition of developmental competence in mammals. Here, we hypothesized that several genes expressed in oocytes and pre-implantation embryos remain unknown. Our goal was to reconstruct the transcriptome of oocytes (germinal vesicle and metaphase II) and pre-implantation cattle embryos (blastocysts) using short-read and long-read sequences to identify putative new genes.

RESULTS

We identified 274,342 transcript sequences and 3,033 of those loci do not match a gene present in official annotations and thus are potential new genes. Notably, 63.67% (1,931/3,033) of potential novel genes exhibited coding potential. Also noteworthy, 97.92% of the putative novel genes overlapped annotation with transposable elements. Comparative analysis of transcript abundance identified that 1,840 novel genes (recently added to the annotation) or potential new genes were differentially expressed between developmental stages (FDR < 0.01). We also determined that 522 novel or potential new genes (448 and 34, respectively) were upregulated at eight-cell embryos compared to oocytes (FDR < 0.01). In eight-cell embryos, 102 novel or putative new genes were co-expressed (|r|> 0.85, P < 1 × 10-8) with several genes annotated with gene ontology biological processes related to pluripotency maintenance and embryo development. CRISPR-Cas9 genome editing confirmed that the disruption of one of the novel genes highly expressed in eight-cell embryos reduced blastocyst development (ENSBTAG00000068261, P = 1.55 × 10-7).

CONCLUSIONS

Our results revealed several putative new genes that need careful annotation. Many of the putative new genes have dynamic regulation during pre-implantation development and are important components of gene regulatory networks involved in pluripotency and blastocyst formation.

Crosstalk between RNA m6A and DNA methylation regulates transposable element chromatin activation and cell fate in human pluripotent stem cells

Transposable elements (TEs) are parasitic DNA sequences accounting for over half of the human genome. Tight control of the repression and activation states of TEs is critical for genome integrity, development, immunity and diseases, including cancer. However, precisely how this regulation is achieved remains unclear. Here we develop a targeted proteomic proximity labeling approach to capture TE-associated proteins in human embryonic stem cells (hESCs). We find that the RNA N6-methyladenosine (m6A) reader, YTHDC2, occupies genomic loci of the primate-specific TE, LTR7/HERV-H, specifically through its interaction with m6A-modified HERV-H RNAs. Unexpectedly, YTHDC2 recruits the DNA 5-methylcytosine (5mC)-demethylase, TET1, to remove 5mC from LTR7/HERV-H and prevent epigenetic silencing. Functionally, the YTHDC2/LTR7 axis inhibits neural differentiation of hESCs. Our results reveal both an underappreciated crosstalk between RNA m6A and DNA 5mC, the most abundant regulatory modifications of RNA and DNA in eukaryotes, and the fact that in hESCs this interplay controls TE activity and cell fate.

Short tandem repeats are important contributors to silencer elements in T cells

The action of cis-regulatory elements with either activation or repression functions underpins the precise regulation of gene expression during normal development and cell differentiation. Gene activation by the combined activities of promoters and distal enhancers has been extensively studied in normal and pathological contexts. In sharp contrast, gene repression by cis-acting silencers, defined as genetic elements that negatively regulate gene transcription in a position-independent fashion, is less well understood. Here, we repurpose the STARR-seq approach as a novel high-throughput reporter strategy to quantitatively assess silencer activity in mammals. We assessed silencer activity from DNase hypersensitive I sites in a mouse T cell line. Identified silencers were associated with either repressive or active chromatin marks and enriched for binding motifs of known transcriptional repressors. CRISPR-mediated genomic deletions validated the repressive function of distinct silencers involved in the repression of non-T cell genes and genes regulated during T cell differentiation. Finally, we unravel an association of silencer activity with short tandem repeats, highlighting the role of repetitive elements in silencer activity. Our results provide a general strategy for genome-wide identification and characterization of silencer elements.

HERVs and Cancer-A Comprehensive Review of the Relationship of Human Endogenous Retroviruses and Human Cancers

Genomic instability and genetic mutations can lead to exhibition of several cancer hallmarks in affected cells such as sustained proliferative signaling, evasion of growth suppression, activated invasion, deregulation of cellular energetics, and avoidance of immune destruction. Similar biological changes have been observed to be a result of pathogenic viruses and, in some cases, have been linked to virus-induced cancers. Human endogenous retroviruses (HERVs), once external pathogens, now occupy more than 8% of the human genome, representing the merge of genomic and external factors. In this review, we outline all reported effects of HERVs on cancer development and discuss the HERV targets most suitable for cancer treatments as well as ongoing clinical trials for HERV-targeting drugs. We reviewed all currently available reports of the effects of HERVs on human cancers including solid tumors, lymphomas, and leukemias. Our review highlights the central roles of HERV genes, such as gag, env, pol, np9, and rec in immune regulation, checkpoint blockade, cell differentiation, cell fusion, proliferation, metastasis, and cell transformation. In addition, we summarize the involvement of HERV long terminal repeat (LTR) regions in transcriptional regulation, creation of fusion proteins, expression of long non-coding RNAs (lncRNAs), and promotion of genome instability through recombination.

Modulation of human endogenous retroviruses and cytokines expression in peripheral blood mononuclear cells from autistic children and their parents

BACKGROUND

Putative pathogenic effects mediated by human endogenous retroviruses (HERVs) in neurological and psychiatric disorders in humans have been extensively described. HERVs may alter the development of the brain by means of several mechanisms, including modulation of gene expression, alteration of DNA stability, and activation of immune system. We recently demonstrated that autistic children and their mothers share high expression levels of some HERVs and cytokines in peripheral blood mononuclear cells (PBMCs) ex vivo, suggesting a close mother-child association in Autism Spectrum Disorder (ASD).

RESULTS

In the present study, PBMCs from autistic children and their parents were exposed to stimulating factors (Interleukin-2/Phytohaemagglutinin) or drugs, as Valproic acid and Efavirenz. The results show that HERVs and cytokines expression can be modulated in vitro by different stimuli in PBMCs from autistic children and their mothers, while no significant changes were found in PBMCs ASD fathers or in controls individuals. In particular, in vitro exposure to interleukin-2/Phytohaemagglutinin or valproic acid induces the expression of several HERVs and cytokines while Efavirenz inhibits them.

CONCLUSION

Herein we show that autistic children and their mothers share an intrinsic responsiveness to in vitro microenvironmental changes in expressing HERVs and pro-inflammatory cytokines. Remarkably, the antiretroviral drug Efavirenz restores the expression of specific HERV families to values similar to those of the controls, also reducing the expression of proinflammatory cytokines but keeping the regulatory ones high. Our findings open new perspectives to study the role of HERVs in the biological mechanisms underlying Autism.

Sequence divergence and retrotransposon insertion underlie interspecific epigenetic differences in primates

Changes in the epigenome can affect the phenotype without the presence of changes in the genomic sequence. Given the high identity of the human and chimpanzee genome sequences, a substantial portion of their phenotypic divergence likely arises from epigenomic differences between the two species. In this study, the transcriptome and epigenome were determined for induced pluripotent stem cells (iPSCs) generated from human and chimpanzee individuals. The transcriptome and epigenomes for trimethylated histone H3 at lysine-4 (H3K4me3) and lysine-27 (H3K27me3) showed high levels of similarity between the two species. However, there were some differences in histone modifications. Although such regions, in general, did not show significant enrichment of interspecies nucleotide variations, gains in binding motifs for pluripotency-related transcription factors, especially POU5F1 and SOX2, were frequently found in species-specific H3K4me3 regions. We also revealed that species-specific insertions of retrotransposons, including the LTR5_Hs subfamily in human and a newly identified LTR5_Pt subfamily in chimpanzee, created species-specific H3K4me3 regions associated with increased expression of nearby genes. Human iPSCs have more species-specific H3K27me3 regions, resulting in more abundant bivalent domains. Only a limited number of these species-specific H3K4me3 and H3K27me3 regions overlap with species-biased enhancers in cranial neural crest cells, suggesting that differences in the epigenetic state of developmental enhancers appear late in development. Therefore, iPSCs serve as a suitable starting material for studying evolutionary changes in epigenome dynamics during development.

Molecular diversity and phenotypic pleiotropy of ancient genomic regulatory loci derived from human endogenous retrovirus type H (HERVH) promoter LTR7 and HERVK promoter LTR5_Hs and their contemporary impacts on pathophysiology of Modern Humans

Timelines of population-level effects of viruses on humans varied from the evolutionary scale of million years to contemporary spread of viral infections. Correspondingly, these events are exemplified by: (i) emergence of human endogenous retroviruses (HERVs) from ancient germline infections leading to stable integration of viral genomes into human chromosomes; and (ii) wide-spread viral infections reaching a global pandemic state such as the COVID-19 pandemic. Despite significant efforts, understanding of HERV’s roles in governance of genomic regulatory networks, their impacts on primate evolution and development of human-specific physiological and pathological phenotypic traits remains limited. Remarkably, present analyses revealed that expression of a dominant majority of genes (1696 of 1944 genes; 87%) constituting high-confidence down-steam regulatory targets of defined HERV loci was significantly altered in cells infected with the SARS-CoV-2 coronavirus, a pathogen causing the global COVID-19 pandemic. This study focused on defined sub-sets of DNA sequences derived from HERVs that are expressed at specific stages of human preimplantation embryogenesis and exert regulatory actions essential for self-renewal and pluripotency. Evolutionary histories of LTR7/HERVH and LTR5_Hs/HERVK were charted based on evidence of the earliest presence and expansion of highly conserved (HC) LTR sequences. Sequence conservation analyses of most recent releases 17 primate species’ genomes revealed that LTR7/HERVH have entered germlines of primates in Africa after the separation of the New World Monkey lineage, while LTR5_Hs/HERVK successfully colonized primates’ germlines after the segregation of Gibbons’ species. Subsequently, both LTR7 and LTR5_Hs undergo a marked ~ fourfold–fivefold expansion in genomes of Great Apes. Timelines of quantitative expansion of both LTR7 and LTR5_Hs loci during evolution of Great Apes appear to replicate the consensus evolutionary sequence of increasing cognitive and behavioral complexities of non-human primates, which seems particularly striking for LTR7 loci and 11 distinct LTR7 subfamilies. Consistent with previous reports, identified in this study, 351 human-specific (HS) insertions of LTR7 (175 loci) and LTR5_Hs (176 loci) regulatory sequences have been linked to genes implicated in establishment and maintenance of naïve and primed pluripotent states and preimplantation embryogenesis phenotypes. Unexpectedly, HS-LTRs manifest regulatory connectivity to genes encoding markers of 12 distinct cells’ populations of fetal gonads, as well as genes implicated in physiology and pathology of human spermatogenesis, including Y-linked spermatogenic failure, oligo- and azoospermia. Granular interrogations of genes linked with 11 distinct LTR7 subfamilies revealed that mammalian offspring survival (MOS) genes seem to remain one of consistent regulatory targets throughout ~ 30 MYA of the divergent evolution of LTR7 loci. Differential GSEA of MOS versus non-MOS genes identified clearly discernable dominant enrichment patterns of phenotypic traits affected by MOS genes linked with LTR7 (562 MOS genes) and LTR5_Hs (126 MOS genes) regulatory loci across the large panel of genomics and proteomics databases reflecting a broad spectrum of human physiological and pathological traits. GSEA of LTR7-linked MOS genes identified more than 2200 significantly enriched records of human common and rare diseases and gene signatures of 466 significantly enriched records of Human Phenotype Ontology traits, including Autosomal Dominant (92 genes) and Autosomal Recessive (93 genes) Inheritance. LTR7 regulatory elements appear linked with genes implicated in functional and morphological features of central nervous system, including synaptic transmission and protein–protein interactions at synapses, as well as gene signatures differentially regulated in cells of distinct neurodevelopmental stages and morphologically diverse cell types residing and functioning in human brain. These include Neural Stem/Precursor cells, Radial Glia cells, Bergman Glia cells, Pyramidal cells, Tanycytes, Immature neurons, Interneurons, Trigeminal neurons, GABAergic neurons, and Glutamatergic neurons. GSEA of LTR7-linked genes identified significantly enriched gene sets encoding markers of more than 80 specialized types of neurons and markers of 521 human brain regions, most prominently, subiculum and dentate gyrus. Identification and characterization of 1944 genes comprising high-confidence down-steam regulatory targets of LTR7 and/or LTR5_Hs loci validated and extended these observations by documenting marked enrichments for genes implicated in neoplasm metastasis, intellectual disability, autism, multiple cancer types, Alzheimer’s, schizophrenia, and other brain disorders. Overall, genes representing down-stream regulatory targets of ancient retroviral LTRs exert the apparently cooperative and exceedingly broad phenotypic impacts on human physiology and pathology. This is exemplified by altered expression of 93% high-confidence LTR targets in cells infected by contemporary viruses, revealing a convergence of virus-inflicted aberrations on genomic regulatory circuitry governed by ancient retroviral LTR elements and interference with human cells’ differentiation programs.

A retrotransposon storm marks clinical phenoconversion to late-onset Alzheimer's disease

Recent reports have suggested that the reactivation of otherwise transcriptionally silent transposable elements (TEs) might induce brain degeneration, either by dysregulating the expression of genes and pathways implicated in cognitive decline and dementia or through the induction of immune-mediated neuroinflammation resulting in the elimination of neural and glial cells. In the work we present here, we test the hypothesis that differentially expressed TEs in blood could be used as biomarkers of cognitive decline and development of AD. To this aim, we used a sample of aging subjects (age > 70) that developed late-onset Alzheimer’s disease (LOAD) over a relatively short period of time (12–48 months), for which blood was available before and after their phenoconversion, and a group of cognitive stable subjects as controls. We applied our developed and validated customized pipeline that allows the identification, characterization, and quantification of the differentially expressed (DE) TEs before and after the onset of manifest LOAD, through analyses of RNA-Seq data. We compared the level of DE TEs within more than 600,000 TE-mapping RNA transcripts from 25 individuals, whose specimens we obtained before and after their phenotypic conversion (phenoconversion) to LOAD, and discovered that 1790 TE transcripts showed significant expression differences between these two timepoints (logFC ± 1.5, logCMP > 5.3, nominal p value < 0.01). These DE transcripts mapped both over- and under-expressed TE elements. Occurring before the clinical phenoconversion, this TE storm features significant increases in DE transcripts of LINEs, LTRs, and SVAs, while those for SINEs are significantly depleted. These dysregulations end with signs of manifest LOAD. This set of highly DE transcripts generates a TE transcriptional profile that accurately discriminates the before and after phenoconversion states of these subjects. Our findings suggest that a storm of DE TEs occurs before phenoconversion from normal cognition to manifest LOAD in risk individuals compared to controls, and may provide useful blood-based biomarkers for heralding such a clinical transition, also suggesting that TEs can indeed participate in the complex process of neurodegeneration.

Effects of Anticancer Agent P-bi-TAT on Gene Expression Link the Integrin Thyroid Hormone Receptor to Expression of Stemness and Energy Metabolism Genes in Cancer Cells

Chemically modified forms of tetraiodothyroacetic acid (tetrac), an L-thyroxine derivative, have been shown to exert their anticancer activity at plasma membrane integrin αvβ3 of tumor cells. Via a specific hormone receptor on the integrin, tetrac-based therapeutic agents modulate expression of genes relevant to cancer cell proliferation, survival and energy metabolism. P-bi-TAT, a novel bivalent tetrac-containing synthetic compound has anticancer activity in vitro and in vivo against glioblastoma multiforme (GBM) and other types of human cancers. In the current study, microarray analysis was carried out on a primary culture of human GBM cells exposed to P-bi-TAT (10⁻⁶ tetrac equivalent) for 24 h. P-bi-TAT significantly affected expression of a large panel of genes implicated in cancer cell stemness, growth, survival and angiogenesis. Recent interest elsewhere in ATP synthase as a target in GBM cells caused us to focus attention on expression of genes involved in energy metabolism. Significantly downregulated transcripts included multiple energy-metabolism-related genes: electron transport chain genes ATP5A1 (ATP synthase 1), ATP51, ATP5G2, COX6B1 (cytochrome c oxidase subunit 6B1), NDUFA8 (NADH dehydrogenase (ubiquinone) FA8), NDUFV2I and other NDUF genes. The NDUF and ATP genes are also relevant to control of oxidative phosphorylation and transcription. Qualitatively similar actions of P-bi-TAT on expression of subsets of energy-metabolism-linked genes were also detected in established human GBM and pancreatic cancer cell lines. In conclusion, acting at αvβ3 integrin, P-bi-TAT caused downregulation in human cancer cells of expression of a large number of genes involved in electron transport and oxidative phosphorylation. These observations suggest that cell surface thyroid hormone receptors on αvβ3 regulate expression of genes relevant to tumor cell stemness and energy metabolism.

The relationship of lamins with epigenetic factors during aging

The key factor of genome instability during aging is transposon dysregulation. This may be due to senile changes in the expression of lamins, which epigenetically modulate transposons. Lamins directly physically interact with transposons. Epigenetic regulators such as SIRT7, BAF, and microRNA can also serve as intermediaries for their interactions. There is also an inverse regulation, since transposons are sources of miRNAs that affect lamins. We suggest that lamins can be attributed to epigenetic factors, since they are part of the NURD, interact with histone deacetylases and regulate gene expression without changing the nucleotide sequences. The role of lamins in the etiopathogenesis of premature aging syndromes may be associated with interactions with transposons. In various human cells, LINE1 is present in the heterochromatin domains of the genome associated with lamins, while SIRT7 facilitates the interaction of this retroelement with lamins. Both retroelements and the nuclear lamina play an important role in the antiviral response of organisms. This may be due to the role of lamins in protection from both viruses and transposons, since viruses and transposons are evolutionarily related. Transposable elements and lamins are secondary messengers of environmental stressors that can serve as triggers for aging and carcinogenesis. Transposons play a role in the development of cancer, while the microRNAs derived from them, participating in the etiopathogenesis of tumors, are important in human aging. Lamins have similar properties, since lamins are dysregulated in cancer, and microRNAs affecting them are involved in carcinogenesis. Changes in the expression of specif ic microRNAs were also revealed in laminopathies. Identif ication of the epigenetic mechanisms of interaction of lamins with transposons during aging can become the basis for the development of methods of life extension and targeted therapy of age-associated cancer

Mosaic cis-regulatory evolution drives transcriptional partitioning of HERVH endogenous retrovirus in the human embryo

The human endogenous retrovirus type-H (HERVH) family is expressed in the preimplantation embryo. A subset of these elements are specifically transcribed in pluripotent stem cells where they appear to exert regulatory activities promoting self-renewal and pluripotency. How HERVH elements achieve such transcriptional specificity remains poorly understood. To uncover the sequence features underlying HERVH transcriptional activity, we performed a phyloregulatory analysis of the long terminal repeats (LTR7) of the HERVH family, which harbor its promoter, using a wealth of regulatory genomics data. We found that the family includes at least eight previously unrecognized subfamilies that have been active at different timepoints in primate evolution and display distinct expression patterns during human embryonic development. Notably, nearly all HERVH elements transcribed in ESCs belong to one of the youngest subfamilies we dubbed LTR7up. LTR7 sequence evolution was driven by a mixture of mutational processes, including point mutations, duplications, and multiple recombination events between subfamilies, that led to transcription factor binding motif modules characteristic of each subfamily. Using a reporter assay, we show that one such motif, a predicted SOX2/3 binding site unique to LTR7up, is essential for robust promoter activity in induced pluripotent stem cells. Together these findings illuminate the mechanisms by which HERVH diversified its expression pattern during evolution to colonize distinct cellular niches within the human embryo.

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.

Transcriptional dynamics of transposable elements when converting fibroblast cells of Macaca mulatta to neuroepithelial stem cells

BACKGROUND

Transposable elements (TE) account for more than 50% of human genome. It has been reported that some types of TEs are dynamically regulated in the reprogramming of human cell lines. However, it is largely unknown whether some TEs in Macaca mulatta are also regulated during the reprogramming of cell lines of monkey.

RESULTS

Here, we systematically examined the transcriptional activities of TEs during the conversion of Macaca mulatta fibroblast cells to neuroepithelial stem cells (NESCs). Hundreds of TEs were dynamically regulated during the reprogramming of Macaca mulatta fibroblast cells. Furthermore, 48 Long Terminal Repeats (LTRs), as well as some integrase elements, of Macaca endogenous retrovirus 3 (MacERV3) were transiently activated during the early stages of the conversion process, some of which were further confirmed with PCR experiments. These LTRs were potentially bound by critical transcription factors for reprogramming, such as KLF4 and ETV5.

CONCLUSION

These results suggest that the transcription of TEs are delicately regulated during the reprogramming of Macaca mulatta fibroblast cells. Although the family of ERVs activated during the reprogramming of fibroblast cells in Macaca mulatta is different from those in the reprogramming of human fibroblast cells, our results suggest that the activation of some ERVs is a conserved mechanism in primates for converting fibroblast cells to stem cells.

Genome-wide R-loop Landscapes during Cell Differentiation and Reprogramming

DNA:RNA hybrids play key roles in both physiological and disease states by regulating chromatin and genome organization. Their homeostasis during cell differentiation and cell plasticity remains elusive. Using an isogenic human stem cell platform, we systematically characterize R-loops, DNA methylation, histone modifications, and chromatin accessibility in pluripotent cells and their lineage-differentiated derivatives. We confirm that a portion of R-loops formed co-transcriptionally at pluripotency genes in pluripotent stem cells and at lineage-controlling genes in differentiated lineages. Notably, a subset of R-loops maintained after differentiation are associated with repressive chromatin marks on silent pluripotency genes and undesired lineage genes. Moreover, in reprogrammed pluripotent cells, cell-of-origin-specific R-loops are initially present but are resolved with serial passaging. Our analysis suggests a multifaceted role of R-loops in cell fate determination that may serve as an additional layer of modulation on cell fate memory and cell plasticity.

Genomics-Guided Drawing of Molecular and Pathophysiological Components of Malignant Regulatory Signatures Reveals a Pivotal Role in Human Diseases of Stem Cell-Associated Retroviral Sequences and Functionally-Active hESC Enhancers

Repetitive DNA sequences (repeats) colonized two-third of human genome and a majority of repeats comprised of transposable genetic elements (TE). Evolutionary distinct categories of TE represent nucleic acid sequences that are repeatedly copied from and pasted into chromosomes at multiple genomic locations and acquired a multitude of regulatory functions. Here, genomics-guided maps of stemness regulatory signatures were drawn to dissect the contribution of TE to clinical manifestations of malignant phenotypes of human cancers. From patients’ and physicians’ perspectives, the clinical definition of a tumor’s malignant phenotype could be restricted to the early diagnosis of sub-types of malignancies with the increased risk of existing therapy failure and high likelihood of death from cancer. It is the viewpoint from which the understanding of stemness and malignant regulatory signatures is considered in this contribution. Genomics-guided analyses of experimental and clinical observations revealed the pivotal role of human stem cell-associated retroviral sequences (SCARS) in the origin and pathophysiology of clinically-lethal malignancies. SCARS were defined as the evolutionary- and biologically-related family of genomic regulatory sequences, the principal physiological function of which is to create and maintain the stemness phenotype during human preimplantation embryogenesis. For cell differentiation to occur, SCARS expression must be silenced and SCARS activity remains repressed in most terminally-differentiated human cells which are destined to perform specialized functions in the human body. Epigenetic reprogramming, de-repression, and sustained activity of SCARS results in various differentiation-defective phenotypes. One of the most prominent tissue- and organ-specific clinical manifestations of sustained SCARS activities is diagnosed as a pathological condition defined by a consensus of morphological, molecular, and genetic examinations as the malignant growth. Here, contemporary evidence are acquired, analyzed, and reported defining both novel diagnostic tools and druggable molecular targets readily amenable for diagnosis and efficient therapeutic management of clinically-lethal malignancies. These diagnostic and therapeutic approaches are based on monitoring of high-fidelity molecular signals of continuing SCARS activities in conjunction with genomic regulatory networks of thousands’ functionally-active embryonic enhancers affecting down-stream phenotype-altering genetic loci. Collectively, reported herein observations support a model of SCARS-activation triggered singular source code facilitating the intracellular propagation and intercellular (systemic) dissemination of disease states in the human body.

Human-chimpanzee fused cells reveal cis-regulatory divergence underlying skeletal evolution

Gene regulatory divergence is thought to play a central role in determining human-specific traits. However, our ability to link divergent regulation to divergent phenotypes is limited. Here, we utilized human-chimpanzee hybrid induced pluripotent stem cells to study gene expression separating these species. The tetraploid hybrid cells allowed us to separate cis- from trans-regulatory effects, and to control for non-genetic confounding factors. We differentiated these cells into cranial neural crest cells (CNCCs), the primary cell type giving rise to the face. We discovered evidence of lineage-specific selection on the hedgehog signaling pathway, including a human-specific 6-fold down-regulation of EVC2 (LIMBIN), a key hedgehog gene. Inducing a similar down-regulation of EVC2 substantially reduced hedgehog signaling output. Mice and humans lacking functional EVC2 show striking phenotypic parallels to human-chimpanzee craniofacial differences, suggesting that the regulatory divergence of hedgehog signaling may have contributed to the unique craniofacial morphology of humans.

Targeted expression profiling reveals distinct stages of early canine fibroblast reprogramming are regulated by 2-oxoglutarate hydroxylases

Background

Ectopic expression of a defined set of transcription factors allows the reprogramming of mammalian somatic cells to pluripotency. Despite continuous progress in primate and rodent reprogramming, limited attention has been paid to cell reprogramming in domestic and companion species. Previous studies attempting to reprogram canine cells have mostly assessed a small number of presumptive canine induced pluripotent stem cell (iPSC) lines for generic pluripotency attributes. However, why canine cell reprogramming remains extremely inefficient is poorly understood.

Methods

To better characterize the initial steps of pluripotency induction in canine somatic cells, we optimized an experimental system where canine fetal fibroblasts (cFFs) are transduced with the Yamanaka reprogramming factors by Sendai virus vectors. We use quantitative PCR arrays to measure the expression of 80 target genes at various stages of canine cell reprogramming. We ask how cFF reprogramming is influenced by small molecules affecting the epigenomic modification 5-hydroxymethylcytosine, specifically L-ascorbic acid and retinoic acid (AA/RA).

Results

We found that the expression and catalytic output of a class of 2-oxoglutarate-dependent (2-OG) hydroxylases, known as ten-eleven translocation (TET) enzymes, can be modulated in canine cells treated with AA/RA. We further show that AA/RA treatment induces TET1 expression and facilitates early canine reprogramming, evidenced by upregulation of epithelial and pluripotency markers. Using a chemical inhibitor of 2-OG hydroxylases, we demonstrate that 2-OG hydroxylase activity regulates the expression of a subset of genes involved in mesenchymal-to-epithelial transition (MET) and pluripotency in early canine reprogramming. We identify a set of transcription factors depleted in maturing reprogramming intermediates compared to pluripotent canine embryonic stem cells.

Conclusions

Our findings highlight 2-OG hydroxylases have evolutionarily conserved and divergent functions regulating the early reprogramming of canine somatic cells and show reprogramming conditions can be rationally optimized for the generation of maturing canine iPSC.

Impacts of genomic networks governed by human-specific regulatory sequences and genetic loci harboring fixed human-specific neuro-regulatory single nucleotide mutations on phenotypic traits of modern humans

Recent advances in identification and characterization of human-specific regulatory DNA sequences set the stage for the assessment of their global impact on physiology and pathology of modern humans. Gene set enrichment analyses (GSEA) of 8405 genes linked with 35,074 human-specific neuro-regulatory single-nucleotide changes (hsSNCs) revealed numerous significant associations with morphological structures, physiological processes, and pathological conditions of modern humans. Significantly enriched traits include more than 1000 anatomically distinct regions of the adult human brain, many different types of cells and tissues, more than 200 common human disorders, and more than 1000 records of rare diseases. Thousands of genes connected with neuro-regulatory hsSNCs have been identified, which represent essential genetic elements of the autosomal inheritance and offspring survival phenotypes. A total of 1494 hsSNC-linked genes are associated with either autosomal dominant or recessive inheritance, and 2273 hsSNC-linked genes have been associated with premature death, embryonic lethality, as well as pre-, peri-, neo-, and post-natal lethality phenotypes of both complete and incomplete penetrance. Differential GSEA implemented on hsSNC-linked loci and associated genes identify a set of 7990 hsSNC-target genes linked to evolutionary distinct classes of human-specific regulatory sequences (HSRS). Notably, the expression of a majority of these genes (5389 genes; 67%) is regulated by stem cell–associated retroviral sequences (SCARS) and SCARS-regulated genes captured a dominant fraction (91%) of significant phenotypic associations linked with hsSNCs. Interrogations of the MGI database revealed readily available mouse models tailored for precise experimental definitions of functional effects of hsSNCs and SCARS on genes causally affecting thousands of mammalian phenotypes and implicated in hundreds of common and rare human disorders. These observations suggest that a preponderance of human-specific traits evolved under a combinatorial regulatory control of distinct classes of HSRS and neuro-regulatory loci harboring hsSNCs that are fixed in humans, distinct from other primates, and located in differentially accessible chromatin regions during brain development.

Transdifferentiation of Mouse Embryonic Fibroblasts into Dopaminergic Neurons Reactivates LINE-1 Repetitive Elements

In mammals, LINE-1 (L1) retrotransposons constitute between 15% and 20% of the genome. Although only a few copies have retained the ability to retrotranspose, evidence in brain and differentiating pluripotent cells indicates that L1 retrotransposition occurs and creates mosaics in normal somatic tissues. The function of de novo insertions remains to be understood. The transdifferentiation of mouse embryonic fibroblasts to dopaminergic neuronal fate provides a suitable model for studying L1 dynamics in a defined genomic and unaltered epigenomic background. We found that L1 elements are specifically re-expressed and mobilized during the initial stages of reprogramming and that their insertions into specific acceptor loci coincides with higher chromatin accessibility and creation of new transcribed units. Those events accompany the maturation of neuronal committed cells. We conclude that L1 retrotransposition is a non-random process correlating with chromatin opening and lncRNA production that accompanies direct somatic cell reprogramming.

•

L1 activation accompanies induced dopaminergic neuron maturation

•

L1 inhibition impairs the transdifferentiation potential of MEFs

•

L1 retrotransposition creates a lineage-specific genetic mosaicism

•

L1 insertions correlates with open chromatin and lncRNA transcription

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.

A Catalogue of 59,732 Human-Specific Regulatory Sequences Reveals Unique-to-Human Regulatory Patterns Associated with Virus-Interacting Proteins, Pluripotency, and Brain Development

Extensive searches for genomic regions harboring various types of candidate human-specific regulatory sequences (HSRS) identified thousands' HSRS using high-resolution next-generation sequencing technologies and methodologically diverse comparative analyses of human and nonhuman primates' (NHPs) reference genomes. In this study, a comprehensive catalogue of 59,732 genomic loci harboring candidate HSRS has been assembled to facilitate the systematic analyses of genomic sequences that were either inherited from extinct common ancestors (ECAs) or created de novo in human genomes. These analyses identified thousands of candidate HSRS and HSRS-harboring loci that appear inherited from ECAs, yet absent in genomes of our closest evolutionary relatives, chimpanzee and bonobo, presumably due to the incomplete lineage sorting and/or species-specific loss or regulatory DNA. This pattern is particularly prominent for HSRS-harboring loci that have been putatively associated with human-specific gene expression changes in cerebral organoid models. A prominent majority of regions harboring human-specific mutations associated with human-specific expression changes during brain development is highly conserved in chimpanzee, bonobo, and gorilla genomes. Among NHPs, dominant fractions of HSRS-harboring loci associated with human-specific gene expression in both excitatory neurons (347 loci; 67%) and radial glia (683 loci; 72%) are highly conserved in the gorilla genome. Analysis of 4433 genes encoding virus-interacting proteins (VIPs) revealed that 95.9% of human VIPs are components of human-specific regulatory networks that appear to operate in distinct types of human cells from preimplantation embryos to adult dorsolateral prefrontal cortex. These analyses demonstrate that modern humans captured unique genome-wide combinations of regulatory sequences, divergent subsets of which are highly conserved in distinct species of six NHP separated by 30 million years of evolution. Concurrently, this unique-to-human mosaic of genomic regulatory patterns inherited from ECAs was supplemented with 12,486 created de novo HSRS. Genes encoding VIPs appear to represent a principal genomic target during evolution of human-specific regulatory networks, which contribute to fitness of Homo sapiens and affect a functionally diverse spectrum of biological and cellular processes controlled by VIP-containing liquid-liquid phase-separated condensates.

The Concomitant Expression of Human Endogenous Retroviruses and Embryonic Genes in Cancer Cells under Microenvironmental Changes is a Potential Target for Antiretroviral Drugs

In our genomes there are thousands of copies of human endogenous retroviruses (HERVs) originated from the integration of exogenous retroviruses that infected germ line cells millions of years ago, and currently an altered expression of this elements has been associated to the onset, progression and acquisition of aggressiveness features of many cancers. The transcriptional reactivation of HERVs is mainly an effect of their responsiveness to some factors in cell microenvironment, such as nutrients, hormones and cytokines. We have already demonstrated that, under pressure of microenvironmental changes, HERV-K (HML-2) activation is required to maintain human melanoma cell plasticity and CD133+ cancer stem cells survival. In the present study, the transcriptional activity of HERV-K (HML-2), HERV-H, CD133 and the embryonic transcription factors OCT4, NANOG and SOX2 was evaluated during the in vitro treatment with antiretroviral drugs in cells from melanoma, liver and lung cancers exposed to microenvironmental changes. The exposure to stem cell medium induced a phenotype switching with the generation of sphere-like aggregates, characterized by the concomitant increase of HERV-K (HML-2) and HERV-H, CD133 and embryonic genes transcriptional activity. Although with heterogenic response among the different cell lines, the in vitro treatment with antiretroviral drugs affected HERVs transcriptional activity in parallel with the reduction of CD133 and embryonic genes expression, clonogenic activity and cell growth, accompanied by the induction of apoptosis. The responsiveness to antiretroviral drugs treatment of cancer cells with stemness features and expressing HERVs suggests the use of these drugs as innovative approach to treat aggressive tumours in combination with chemotherapeutic/radiotherapy regimens.

Children With Autism Spectrum Disorder and Their Mothers Share Abnormal Expression of Selected Endogenous Retroviruses Families and Cytokines

The Autism Spectrum Disorder (ASD) is a heterogeneous group of neurodevelopmental disorders, only clinically diagnosed since the lack of reliable biomarkers. Autism etiology is probably attributable to the combination of genetic vulnerability and environmental factors, and recently, maternal immune activation has been linked to derailed neurodevelopment, resulting in ASD in the offspring. Human endogenous retroviruses (HERVs) are relics of ancestral infections, stably integrated in the human DNA. Given the HERV persistence in the genome, some of HERVs have been co-opted for physiological functions during evolution, while their reactivation has been associated with several pathological conditions, including cancer, autoimmune, and neurological and psychiatric disorders. Particularly, due to their intrinsic responsiveness to external stimuli, HERVs can modulate the host immune response and in turn HERVs can be activated by the immune effectors. In previous works we demonstrated high expression levels of HERV-H in blood of autistic patients, closely related with the severity of the disease. Moreover, in a preclinical ASD model we proved changes of expression of several ERV families and cytokines from the intrauterine life to the adulthood, and across generations via maternal lineage. Here we analyzed the expression of HEMO and of selected HERVs and cytokines in blood from ASD patients and their parents and corresponding healthy controls, to look for a common molecular trait within family members. ASD patients and their mothers share altered expression of HERV-H and HEMO and of cytokines such as TNF-α, IFN-γ, IL-10. The multivariate regression models showed a mother-child association by HEMO activity and demonstrated in children and mothers an association between HERV-H and HEMO expression and, only in mothers, between HEMO, and TNF-α expression. Furthermore, high diagnostic performance for HERV-H and HEMO was found, suggesting their potential application for the identification of ASD children and their mothers. The present data support the involvement of HERVs in ASD and suggest HERVs and cytokines as ASD-associated traits. Since ASD is a heterogeneous group of neurodevelopmental disorders, a single determinant alone could be not enough to account for the complexity, and HERV/cytokines expression could be considered in a set of biomarkers, easily detectable in blood, and potentially useful for an early diagnosis.

Novel Bioinformatics Approach Identifies Transcriptional Profiles of Lineage-Specific Transposable Elements at Distinct Loci in the Human Dorsolateral Prefrontal Cortex

Expression of transposable elements (TE) is transiently activated during human preimplantation embryogenesis in a developmental stage- and cell type-specific manner and TE-mediated epigenetic regulation is intrinsically wired in developmental genetic networks in human embryos and embryonic stem cells. However, there are no systematic studies devoted to a comprehensive analysis of the TE transcriptome in human adult organs and tissues, including human neural tissues. To investigate TE expression in the human Dorsolateral Prefrontal Cortex (DLPFC), we developed and validated a straightforward analytical approach to chart quantitative genome-wide expression profiles of all annotated TE loci based on unambiguous mapping of discrete TE-encoded transcripts using a de novo assembly strategy. To initially evaluate the potential regulatory impact of DLPFC-expressed TE, we adopted a comparative evolutionary genomics approach across humans, primates, and rodents to document conservation patterns, lineage-specificity, and colocalizations with transcription factor binding sites mapped within primate- and human-specific TE. We identified 654,665 transcripts expressed from 477,507 distinct loci of different TE classes and families, the majority of which appear to have originated from primate-specific sequences. We discovered 4,687 human-specific and transcriptionally active TEs in DLPFC, of which the prominent majority (80.2%) appears spliced. Our analyses revealed significant associations of DLPFC-expressed TE with primate- and human-specific transcription factor binding sites, suggesting potential cross-talks of concordant regulatory functions. We identified 1,689 TEs differentially expressed in the DLPFC of Schizophrenia patients, a majority of which is located within introns of 1,137 protein-coding genes. Our findings imply that identified DLPFC-expressed TEs may affect human brain structures and functions following different evolutionary trajectories. On one side, hundreds of thousands of TEs maintained a remarkably high conservation for ∼8 My of primates’ evolution, suggesting that they are likely conveying evolutionary-constrained primate-specific regulatory functions. In parallel, thousands of transcriptionally active human-specific TE loci emerged more recently, suggesting that they could be relevant for human-specific behavioral or cognitive functions.

Human Endogenous Retrovirus-H and K Expression in Human Mesenchymal Stem Cells as Potential Markers of Stemness

OBJECTIVE

The human endogenous retroviruses (HERVs) are endogenous retroviruses that were inserted into the germ cell DNA of humans over 30 million years ago. Insertion of HERVs into the chromosomal DNA can influence a number of host genes in various modes during human evolution and their proviral long terminal repeats can participate in the transcriptional regulation of various cellular genes. Our aim was to evaluate the pol gene expression of HERV-K and HERV-H in mesenchymal stem cells (MSCs) in relation with the expression of stemness genes such as NANOG, OCT-4, and SOX-2.

METHODS

MSCs were isolated from bone marrow of healthy donors and expanded until the 5th passage in α-MEM with 10% fetal bovine serum. HERV-K, HERV-H pol gene, NANOG, OCT-4, SOX-2, and GAPDH expression was quantified by real-time PCR in MSCs during the expansion.

RESULTS

HERV-K and HERV-H expression was always higher at p1 compared to other passages and this difference reached a high statistical significance when passage p1 was compared with passage 3. In addition, NANOG, OCT-4, and SOX-2 expression at p1 was significantly higher than their expression at p3. Pearson's test demonstrated a strong correlation between the expression of HERV-K and HERV-H and the expression of NANOG, OCT-4, and SOX-2.

CONCLUSIONS

Our findings showed that HERV-K and H were concurrently expressed with pluripotency biomarkers NANOG, OCT-4, and SOX-2.

The evolution of Great Apes has shaped the functional enhancers' landscape in human embryonic stem cells

High-throughput functional assays of enhancer activity have recently enabled the genome-scale definition of molecular, structural, and biochemical features of these genomic regulatory regions. To infer the evolutionary origin of DNA sequences operating as functional enhancers in human embryonic stem cells (hESC), we examined the patterns of evolutionary conservation and divergence in the genome-wide functional enhancers' landscape of hESC. We show that a prominent majority (up to 94%) of DNA sequences identified in hESC as functional enhancers are conserved in humans and our closest evolutionary relatives, Chimpanzee and Bonobo. More than 91% of functional enhancers that are highly conserved in both Chimpanzee and Bonobo, are conserved among other Great Apes and >75% are conserved in the Rhesus genome. In striking contrast, <5% of DNA sequences operating in hESC as functional enhancers are conserved in rodents. Conserved in primates enhancers' sequences are complemented by 1619 sequences of enhancers that are specific to humans. Enhancers that harbor human-specific sequences appear enriched among the invariant enhancer module maintaining activity in different pluripotent states and these regions are associated with pluripotency- and embryonic-lineage-related genes. However, functional enhancers make up only a minority of all conserved in primates or human-specific transcription factor binding sites. Our analyses revealed that sequences that are conserved during ~8 million years of primate evolution dominate the genomic landscape of functional enhancers in both primed and naïve hESC. Collectively, these observations revealed thousands of evolutionarily conserved sequences that function as a core regulatory network in human embryonic stem cells which has recently undergone further extension after divergence of modern humans from our closest relatives, Chimpanzee and Bonobo.

Global long terminal repeat activation participates in establishing the unique gene expression programme of classical Hodgkin lymphoma

Long terminal repeat (LTR) elements are wide-spread in the human genome and have the potential to act as promoters and enhancers. Their expression is therefore under tight epigenetic control. We previously reported in classical Hodgkin Lymphoma (cHL) that a member of the THE1B class of LTR elements acted as a promoter for the proto-oncogene and growth factor receptor gene CSF1R and that expression of this gene is required for cHL tumour survival. However, to which extent and how such elements participate in globally shaping the unique cHL gene expression programme is unknown. To address this question we mapped the genome-wide activation of THE1-LTRs in cHL cells using a targeted next generation sequencing approach (RACE-Seq). Integration of these data with global gene expression data from cHL and control B cell lines showed a unique pattern of LTR activation impacting on gene expression, including genes associated with the cHL phenotype. We also show that global LTR activation is induced by strong inflammatory stimuli. Together these results demonstrate that LTR activation provides an additional layer of gene deregulation in classical Hodgkin lymphoma and highlight the potential impact of genome-wide LTR activation in other inflammatory diseases.

Exaptation at the molecular genetic level

The realization that body parts of animals and plants can be recruited or coopted for novel functions dates back to, or even predates the observations of Darwin. S.J. Gould and E.S. Vrba recognized a mode of evolution of characters that differs from adaptation. The umbrella term aptation was supplemented with the concept of exaptation. Unlike adaptations, which are restricted to features built by selection for their current role, exaptations are features that currently enhance fitness, even though their present role was not a result of natural selection. Exaptations can also arise from nonaptations; these are characters which had previously been evolving neutrally. All nonaptations are potential exaptations. The concept of exaptation was expanded to the molecular genetic level which aided greatly in understanding the enormous potential of neutrally evolving repetitive DNA-including transposed elements, formerly considered junk DNA-for the evolution of genes and genomes. The distinction between adaptations and exaptations is outlined in this review and examples are given. Also elaborated on is the fact that such distinctions are sometimes more difficult to determine; this is a widespread phenomenon in biology, where continua abound and clear borders between states and definitions are rare.

The Decrease in Human Endogenous Retrovirus-H Activity Runs in Parallel with Improvement in ADHD Symptoms in Patients Undergoing Methylphenidate Therapy

Increasing scientific evidence demonstrated the deregulation of human endogenous retroviruses (HERVs) expression in complex diseases, such as cancer, autoimmune, psychiatric, and neurological disorders. The dynamic regulation of HERV activity and their responsiveness to a variety of environmental stimuli designate HERVs as genetic elements that could be modulated by drugs. Methylphenidate (MPH) is widely used in the treatment of attention deficit hyperactivity disorder (ADHD). The aim of this study was to evaluate the time course of human endogenous retrovirus H (HERV-H) expression in peripheral blood mononuclear cells (PBMCs) with respect to clinical response in ADHD patients undergoing MPH therapy. A fast reduction in HERV-H activity in ADHD patients undergoing MPH therapy was observed in parallel with an improvement in clinical symptoms. Moreover, when PBMCs from drug-naïve patients were cultured in vitro, HERV-H expression increased, while no changes in the expression levels were found in ADHD patients undergoing therapy. This suggests that MPH could affect the HERV-H activity and supports the hypothesis that high expression levels of HERV-H could be considered a distinctive trait of ADHD patients.

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.

Functional Dissection of the Enhancer Repertoire in Human Embryonic Stem Cells

Enhancers are genetic elements that regulate spatiotemporal gene expression. Enhancer function requires transcription factor (TF) binding and correlates with histone modifications. However, the extent to which TF binding and histone modifications functionally define active enhancers remains unclear. Here, we combine chromatin immunoprecipitation with a massively parallel reporter assay (ChIP-STARR-seq) to identify functional enhancers in human embryonic stem cells (ESCs) genome-wide in a quantitative unbiased manner. Although active enhancers associate with TFs, only a minority of regions marked by NANOG, OCT4, H3K27ac, and H3K4me1 function as enhancers, with activity markedly changing under naive versus primed culture conditions. We identify an enhancer set associated with functions extending to non-ESC-specific processes. Moreover, although transposable elements associate with putative enhancers, only some exhibit activity. Similarly, within super-enhancers, large tracts are non-functional, with activity restricted to small sub-domains. This catalog of validated enhancers provides a valuable resource for further functional dissection of the regulatory genome.

•

Massively parallel reporter assay assessed over 350,000 genome regions

•

ChIP-STARR-seq catalogs functional enhancers in primed and naive hESCs

•

Identification of transcription factors and transposable elements linked to enhancers

•

ChIP-STARR-seq dissects super-enhancers into small functional units

Single cell expression analysis of primate-specific retroviruses-derived HPAT lincRNAs in viable human blastocysts identifies embryonic cells co-expressing genetic markers of multiple lineages

Chromosome instability and aneuploidies occur very frequently in human embryos, impairing proper embryogenesis and leading to cell cycle arrest, loss of cell viability, and developmental failures in 50–80% of cleavage-stage embryos. This high frequency of cellular extinction events represents a significant experimental obstacle challenging analyses of individual cells isolated from human preimplantation embryos. We carried out single cell expression profiling of 241 individual cells recovered from 32 human embryos during the early and late stages of viable human blastocyst (VHB) differentiation. Classification of embryonic cells was performed solely based on expression patterns of human pluripotency-associated transcripts (HPAT), which represent a family of primate-specific transposable element-derived lincRNAs highly expressed in human embryonic stem cells and regulating nuclear reprogramming and pluripotency induction. We then validated our findings by analyzing transcriptomes of 1,708 individual cells recovered from more than 100 human embryos and 259 mouse cells from more than 40 mouse embryos at different stages of preimplantation embryogenesis. HPAT's expression-guided spatiotemporal reconstruction of human embryonic development inferred from single-cell expression analysis of VHB differentiation enabled identification of telomerase-positive embryonic cells co-expressing key pluripotency regulatory genes and genetic markers of three major lineages. Follow-up validation analyses confirmed the emergence in human embryos prior to lineage segregation of telomerase-positive cells co-expressing genetic markers of multiple lineages. Observations reported in this contribution support the hypothesis of a developmental pathway of creation embryonic lineages and extraembryonic tissues from telomerase-positive pre-lineage cells manifesting multi-lineage precursor phenotype.

Contrasted patterns of evolution of the LINE-1 retrotransposon in perissodactyls: the history of a LINE-1 extinction

Background

LINE-1 (L1) is the dominant autonomously replicating non-LTR retrotransposon in mammals. Although our knowledge of L1 evolution across the tree of life has considerably improved in recent years, what we know of L1 evolution in mammals is biased and comes mostly from studies in primates (mostly human) and rodents (mostly mouse). It is unclear if patterns of evolution that are shared between those two groups apply to other mammalian orders. Here we performed a detailed study on the evolution of L1 in perissodactyls by making use of the complete genome of the domestic horse and of the white rhinoceros. This mammalian order offers an excellent model to study the extinction of L1 since the rhinoceros is one of the few mammalian species to have lost active L1.

Results

We found that multiple L1 lineages, carrying different 5’UTRs, have been simultaneously active during the evolution of perissodactyls. We also found that L1 has continuously amplified and diversified in horse. In rhinoceros, L1 was very prolific early on. Two successful families were simultaneously active until ~20my ago but became extinct suddenly at exactly the same time.

Conclusions

The general pattern of L1 evolution in perissodactyls is very similar to what was previously described in mouse and human, suggesting some commonalities in the way mammalian genomes interact with L1. We confirmed the extinction of L1 in rhinoceros and we discuss several possible mechanisms.

Electronic supplementary material

The online version of this article (10.1186/s13100-018-0117-4) contains supplementary material, which is available to authorized users.

Contribution of transposable elements and distal enhancers to evolution of human-specific features of interphase chromatin architecture in embryonic stem cells

Transposable elements have made major evolutionary impacts on creation of primate-specific and human-specific genomic regulatory loci and species-specific genomic regulatory networks (GRNs). Molecular and genetic definitions of human-specific changes to GRNs contributing to development of unique to human phenotypes remain a highly significant challenge. Genome-wide proximity placement analysis of diverse families of human-specific genomic regulatory loci (HSGRL) identified topologically associating domains (TADs) that are significantly enriched for HSGRL and designated rapidly evolving in human TADs. Here, the analysis of HSGRL, hESC-enriched enhancers, super-enhancers (SEs), and specific sub-TAD structures termed super-enhancer domains (SEDs) has been performed. In the hESC genome, 331 of 504 (66%) of SED-harboring TADs contain HSGRL and 68% of SEDs co-localize with HSGRL, suggesting that emergence of HSGRL may have rewired SED-associated GRNs within specific TADs by inserting novel and/or erasing existing non-coding regulatory sequences. Consequently, markedly distinct features of the principal regulatory structures of interphase chromatin evolved in the hESC genome compared to mouse: the SED quantity is 3-fold higher and the median SED size is significantly larger. Concomitantly, the overall TAD quantity is increased by 42% while the median TAD size is significantly decreased (p = 9.11E-37) in the hESC genome. Present analyses illustrate a putative global role for transposable elements and HSGRL in shaping the human-specific features of the interphase chromatin organization and functions, which are facilitated by accelerated creation of novel transcription factor binding sites and new enhancers driven by targeted placement of HSGRL at defined genomic coordinates. A trend toward the convergence of TAD and SED architectures of interphase chromatin in the hESC genome may reflect changes of 3D-folding patterns of linear chromatin fibers designed to enhance both regulatory complexity and functional precision of GRNs by creating predominantly a single gene (or a set of functionally linked genes) per regulatory domain structures. Collectively, present analyses reveal critical evolutionary contributions of transposable elements and distal enhancers to creation of thousands primate- and human-specific elements of a chromatin folding code, which defines the 3D context of interphase chromatin both restricting and facilitating biological functions of GRNs.

Living Organisms Author Their Read-Write Genomes in Evolution

Evolutionary variations generating phenotypic adaptations and novel taxa resulted from complex cellular activities altering genome content and expression: (i) Symbiogenetic cell mergers producing the mitochondrion-bearing ancestor of eukaryotes and chloroplast-bearing ancestors of photosynthetic eukaryotes; (ii) interspecific hybridizations and genome doublings generating new species and adaptive radiations of higher plants and animals; and, (iii) interspecific horizontal DNA transfer encoding virtually all of the cellular functions between organisms and their viruses in all domains of life. Consequently, assuming that evolutionary processes occur in isolated genomes of individual species has become an unrealistic abstraction. Adaptive variations also involved natural genetic engineering of mobile DNA elements to rewire regulatory networks. In the most highly evolved organisms, biological complexity scales with "non-coding" DNA content more closely than with protein-coding capacity. Coincidentally, we have learned how so-called "non-coding" RNAs that are rich in repetitive mobile DNA sequences are key regulators of complex phenotypes. Both biotic and abiotic ecological challenges serve as triggers for episodes of elevated genome change. The intersections of cell activities, biosphere interactions, horizontal DNA transfers, and non-random Read-Write genome modifications by natural genetic engineering provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Transposable Element Exaptation into Regulatory Regions Is Rare, Influenced by Evolutionary Age, and Subject to Pleiotropic Constraints

Transposable element (TE)-derived sequences make up approximately half of most mammalian genomes, and many TEs have been co-opted into gene regulatory elements. However, we lack a comprehensive tissue- and genome-wide understanding of how and when TEs gain regulatory activity in their hosts. We evaluated the prevalence of TE-derived DNA in enhancers and promoters across hundreds of human and mouse cell lines and primary tissues. Promoters are significantly depleted of TEs in all tissues compared with their overall prevalence in the genome (P < 0.001); enhancers are also depleted of TEs, though not as strongly as promoters. The degree of enhancer depletion also varies across contexts (1.5-3×), with reproductive and immune cells showing the highest levels of TE regulatory activity in humans. Overall, in spite of the regulatory potential of many TE sequences, they are significantly less active in gene regulation than expected from their prevalence. TE age is predictive of the likelihood of enhancer activity; TEs originating before the divergence of amniotes are 9.2 times more likely to have enhancer activity than TEs that integrated in great apes. Context-specific enhancers are more likely to be TE-derived than enhancers active in multiple tissues, and young TEs are more likely to overlap context-specific enhancers than old TEs (86% vs. 47%). Once TEs obtain enhancer activity in the host, they have similar functional dynamics to one another and non-TE-derived enhancers, likely driven by pleiotropic constraints. However, a few TE families, most notably endogenous retroviruses, have greater regulatory potential. Our observations suggest a model of regulatory co-option in which TE-derived sequences are initially repressed, after which a small fraction obtains context-specific enhancer activity, with further gains subject to pleiotropic constraints.

EnHERV: Enrichment analysis of specific human endogenous retrovirus patterns and their neighboring genes

Human endogenous retroviruses (HERVs) are flanked by long terminal repeats (LTRs), which contain the regulation part of the retrovirus. Remaining HERVs constitute 7% to 8% of the present day human genome, and most have been identified as solo LTRs. The HERV sequences have been associated with several molecular functions as well as certain diseases in human, but their roles in human diseases are yet to be established. We designed EnHERV to make accessible the identified endogenous retrovirus repetitive sequences from Repbase Update (a database of eukaryotic repetitive elements) that are present in the human genome. Defragmentation process was done to improve the RepeatMasker annotation output. The defragmented elements were used as core database in EnHERV. EnHERV is available at http://sysbio.chula.ac.th/enherv and can be searched using either gene lists of user interest or HERV characteristics. Besides the search function, EnHERV also provides an enrichment analysis function that allows users to perform enrichment analysis between selected HERV characteristics and user-input gene lists, especially genes with the expression profile of a certain disease. EnHERV will facilitate exploratory studies of specific HERV characteristics that control gene expression patterns related to various disease conditions. Here we analyzed 25 selected HERV groups/names from all four HERV superfamilies, using the sense and anti-sense directions of the HERV and gene expression profiles from 49 specific tissue and disease conditions. We found that intragenic HERVs were associated with down-regulated genes in most cancer conditions and in psoriatic skin tissues and associated with up-regulated genes in immune cells particularly from systemic lupus erythematosus (SLE) patients. EnHERV allowed the analysis of how different types of LTRs were differentially associated with specific gene expression profiles in particular disease conditions for further studies into their mechanisms and functions.

Prenatal Air Pollution Exposures, DNA Methyl Transferase Genotypes, and Associations with Newborn LINE1 and Alu Methylation and Childhood Blood Pressure and Carotid Intima-Media Thickness in the Children's Health Study

BACKGROUND

Although exposure to ambient air pollutants increases cardiovascular disease risk in adults little is known about the effects of prenatal exposure. Genetic variation and epigenetic alterations are two mechanisms that may influence the effects of early-life exposures on cardiovascular phenotypes.

OBJECTIVES

We investigated whether genetic and epigenetic variation modify associations between prenatal air pollution on markers of cardiovascular risk in childhood.

METHODS

We used linear regression analysis to investigate the associations between prenatal pollutants (PM2.5, PM10, NO2, O3), long interspersed nuclear elements (LINE1) and AluYb8 DNA methylation levels measured in newborn blood spot tests, and carotid intima-media thickness (CIMT) and blood pressure (BP) in 459 participants as part of the Children's Health Study. Interaction terms were also included to test for effect modification of these associations by genetic variation in methylation reprogramming genes.

RESULTS

Prenatal exposure to NO2 in the third trimester of pregnancy was associated with higher systolic BP in 11-year-old children. Prenatal exposure to multiple air pollutants in the first trimester was associated with lower DNA methylation in LINE1, whereas later exposure to O3 was associated with higher LINE1 methylation levels in newborn blood spots. The magnitude of associations with prenatal air pollution varied according to genotype for 11 SNPs within DNA methyltransferase 1 (DNMT1), DNA methyltransferase 3 Beta (DNMT3B), Tet methylcytosine dioxygenase 2 (TET2), and Thymine DNA glycosylase (TDG) genes. Although first-trimester O3 exposure was not associated with CIMT and systolic BP overall, associations within strata of DNMT1 or DNMT3B were observed, and the magnitude and the direction of these associations depended on DNMT1 genotypes.

CONCLUSIONS

Genetic and epigenetic variation in DNA methylation reprogramming genes and in LINE1 retrotransposons may play important roles in downstream cardiovascular consequences of prenatal air pollution exposure. Citation: Breton CV, Yao J, Millstein J, Gao L, Siegmund KD, Mack W, Whitfield-Maxwell L, Lurmann F, Hodis H, Avol E, Gilliland FD. 2016. Prenatal air pollution exposures, DNA methyl transferase genotypes, and associations with newborn LINE1 and Alu methylation and childhood blood pressure and carotid intima-media thickness in the Children's Health Study. Environ Health Perspect 124:1905-1912; http://dx.doi.org/10.1289/EHP181.

Mechanistically Distinct Pathways of Divergent Regulatory DNA Creation Contribute to Evolution of Human-Specific Genomic Regulatory Networks Driving Phenotypic Divergence of Homo sapiens

Thousands of candidate human-specific regulatory sequences (HSRS) have been identified, supporting the hypothesis that unique to human phenotypes result from human-specific alterations of genomic regulatory networks. Collectively, a compendium of multiple diverse families of HSRS that are functionally and structurally divergent from Great Apes could be defined as the backbone of human-specific genomic regulatory networks. Here, the conservation patterns analysis of 18,364 candidate HSRS was carried out requiring that 100% of bases must remap during the alignments of human, chimpanzee, and bonobo sequences. A total of 5,535 candidate HSRS were identified that are: (i) highly conserved in Great Apes; (ii) evolved by the exaptation of highly conserved ancestral DNA; (iii) defined by either the acceleration of mutation rates on the human lineage or the functional divergence from non-human primates. The exaptation of highly conserved ancestral DNA pathway seems mechanistically distinct from the evolution of regulatory DNA segments driven by the species-specific expansion of transposable elements. Genome-wide proximity placement analysis of HSRS revealed that a small fraction of topologically associating domains (TADs) contain more than half of HSRS from four distinct families. TADs that are enriched for HSRS and termed rapidly evolving in humans TADs (revTADs) comprise 0.8–10.3% of 3,127 TADs in the hESC genome. RevTADs manifest distinct correlation patterns between placements of human accelerated regions, human-specific transcription factor-binding sites, and recombination rates. There is a significant enrichment within revTAD boundaries of hESC-enhancers, primate-specific CTCF-binding sites, human-specific RNAPII-binding sites, hCONDELs, and H3K4me3 peaks with human-specific enrichment at TSS in prefrontal cortex neurons (P < 0.0001 in all instances). Present analysis supports the idea that phenotypic divergence of Homo sapiens is driven by the evolution of human-specific genomic regulatory networks via at least two mechanistically distinct pathways of creation of divergent sequences of regulatory DNA: (i) recombination-associated exaptation of the highly conserved ancestral regulatory DNA segments; (ii) human-specific insertions of transposable elements.

Evolutionary impact of transposable elements on genomic diversity and lineage-specific innovation in vertebrates

Since their discovery, a growing body of evidence has emerged demonstrating that transposable elements are important drivers of species diversity. These mobile elements exhibit a great variety in structure, size and mechanisms of transposition, making them important putative actors in organism evolution. The vertebrates represent a highly diverse and successful lineage that has adapted to a wide range of different environments. These animals also possess a rich repertoire of transposable elements, with highly diverse content between lineages and even between species. Here, we review how transposable elements are driving genomic diversity and lineage-specific innovation within vertebrates. We discuss the large differences in TE content between different vertebrate groups and then go on to look at how they affect organisms at a variety of levels: from the structure of chromosomes to their involvement in the regulation of gene expression, as well as in the formation and evolution of non-coding RNAs and protein-coding genes. In the process of doing this, we highlight how transposable elements have been involved in the evolution of some of the key innovations observed within the vertebrate lineage, driving the group's diversity and success.

Nothing in Evolution Makes Sense Except in the Light of Genomics: Read-Write Genome Evolution as an Active Biological Process

The 21st century genomics-based analysis of evolutionary variation reveals a number of novel features impossible to predict when Dobzhansky and other evolutionary biologists formulated the neo-Darwinian Modern Synthesis in the middle of the last century. These include three distinct realms of cell evolution; symbiogenetic fusions forming eukaryotic cells with multiple genome compartments; horizontal organelle, virus and DNA transfers; functional organization of proteins as systems of interacting domains subject to rapid evolution by exon shuffling and exonization; distributed genome networks integrated by mobile repetitive regulatory signals; and regulation of multicellular development by non-coding lncRNAs containing repetitive sequence components. Rather than single gene traits, all phenotypes involve coordinated activity by multiple interacting cell molecules. Genomes contain abundant and functional repetitive components in addition to the unique coding sequences envisaged in the early days of molecular biology. Combinatorial coding, plus the biochemical abilities cells possess to rearrange DNA molecules, constitute a powerful toolbox for adaptive genome rewriting. That is, cells possess "Read-Write Genomes" they alter by numerous biochemical processes capable of rapidly restructuring cellular DNA molecules. Rather than viewing genome evolution as a series of accidental modifications, we can now study it as a complex biological process of active self-modification.

Endogenized viral sequences in mammals

Reverse-transcribed RNA molecules compose a significant portion of the human genome. Many of these RNA molecules were retrovirus genomes either infecting germline cells or having done so in a previous generation but retaining transcriptional activity. This mechanism itself accounts for a quarter of the genomic sequence information of mammals for which there is data. We understand relatively little about the causes and consequences of retroviral endogenization. This review highlights functions ascribed to sequences of viral origin endogenized into mammalian genomes and suggests some of the most pressing questions raised by these observations.

Activation of endogenous human stem cell-associated retroviruses (SCARs) and therapy-resistant phenotypes of malignant tumors

Recent reports revealed consistent activation of specific endogenous retroviral elements in human preimplantation embryos and embryonic stem cells. Activity of stem cell associated retroviruses (SCARs) has been implicated in seeding thousands of human-specific regulatory sequences in the hESC genome. Activation of specific SCARs has been demonstrated in patients diagnosed with multiple types of cancer, autoimmune diseases, and neurodegenerative disorders, and appears associated with clinically lethal therapy resistant death-from-cancer phenotypes in a sub-set of cancer patients diagnosed with different types of malignant tumors. A hallmark feature of human-specific SCAR integration sites is deletions of ancestral DNA. Analysis of human-specific genetic loci of SCARs' stemness networks in tumor samples of TCGA cohorts representing 29 cancer types suggests that this approach may facilitate identification of pan-cancer genomic signatures of clinically-lethal disease defined by the presence of somatic non-silent mutations, gene-level copy number changes, and transcripts and proteins' expression of SCAR-regulated host genes. Present analyses indicate that multiple lines of strong circumstantial evidence support the hypothesis that activation of SCARs' networks may play an important role in cancer progression and metastasis, perhaps contributing to the emergence of clinically-lethal therapy-resistant death-from-cancer phenotypes.

Not so bad after all: retroviruses and long terminal repeat retrotransposons as a source of new genes in vertebrates

Viruses and transposable elements, once considered as purely junk and selfish sequences, have repeatedly been used as a source of novel protein-coding genes during the evolution of most eukaryotic lineages, a phenomenon called 'molecular domestication'. This is exemplified perfectly in mammals and other vertebrates, where many genes derived from long terminal repeat (LTR) retroelements (retroviruses and LTR retrotransposons) have been identified through comparative genomics and functional analyses. In particular, genes derived from gag structural protein and envelope (env) genes, as well as from the integrase-coding and protease-coding sequences, have been identified in humans and other vertebrates. Retroelement-derived genes are involved in many important biological processes including placenta formation, cognitive functions in the brain and immunity against retroelements, as well as in cell proliferation, apoptosis and cancer. These observations support an important role of retroelement-derived genes in the evolution and diversification of the vertebrate lineage.

Viruses, stemness, embryogenesis, and cancer: a miracle leap toward molecular definition of novel oncotargets for therapy-resistant malignant tumors?

Recent breakthrough studies documented consistent activation of specific endogenous retroviruses in human embryonic stem cells and preimplantation human embryos and demonstrated the essential role of the sustained retroviral activities for maintenance of pluripotency and embryonic stem cell identity. Present analysis has led to the hypothesis that activation of the human stem cell-associated retroviruses (SCARs), namely LTR7/HERVH and LTR5_Hs/HERVK, is likely associated with the emergence of clinically lethal therapy resistant death-from-cancer phenotypes in a sub-set of cancer patients diagnosed with different types of malignant tumors.