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

VAMP7j: A Splice Variant of Human VAMP7 That Modulates Neurite Outgrowth by Regulating L1CAM Transport to the Plasma Membrane

The vesicle-associated membrane protein 7 (VAMP7) is a SNARE protein of the longin family involved in a wide range of subcellular trafficking events, including neurite sprouting and elongation. The expression of the human gene SYBL1, encoding VAMP7, is finely regulated by alternative splicing. Among the minor isoforms identified so far, VAMP7j is the one most expressed and modulated in the human brain. Therefore, we focused on gaining functional evidence on VAMP7j, which lacks a functional SNARE motif but retains both the longin and transmembrane domains. In human SH-SY5Y cells, we found VAMP7j to modulate neuritogenesis by mediating transport of L1CAM toward the plasma membrane, in a fashion regulated by phosphorylation of the longin domain. VAMP7-mediated regulation of L1CAM trafficking seems at least to differentiate humans from rats, with VAMP7j CNS expression being restricted to primates, including humans. Since L1CAM is a central player in neuritogenesis and axon guidance, these findings suggest the species-specific splicing of SYBL1 is among the fine tuners of human neurodevelopmental complexity.

Evolutionarily recent retrotransposons contribute to schizophrenia

Transposable elements (TEs) are mobile genetic elements that constitute half of the human genome. Recent studies suggest that polymorphic non-reference TEs (nrTEs) may contribute to cognitive diseases, such as schizophrenia, through a cis-regulatory effect. The aim of this work is to identify sets of nrTEs putatively linked to an increased risk of developing schizophrenia. To do so, we inspected the nrTE content of genomes from the dorsolateral prefrontal cortex of schizophrenic and control individuals and identified 38 nrTEs that possibly contribute to the emergence of this psychiatric disorder, two of them further confirmed with haplotype-based methods. We then performed in silico functional inferences and found that 9 of the 38 nrTEs act as expression/alternative splicing quantitative trait loci (eQTLs/sQTLs) in the brain, suggesting a possible role in shaping the human cognitive genome structure. To our knowledge, this is the first attempt at identifying polymorphic nrTEs that can contribute to the functionality of the brain. Finally, we suggest that a neurodevelopmental genetic mechanism, which involves evolutionarily young nrTEs, can be key to understanding the ethio-pathogenesis of this complex disorder.

Maternal immune activation alters placental histone-3 lysine-9 tri-methylation, offspring sensorimotor processing, and hypothalamic transposable element expression in a sex-specific manner

Animal models of maternal immune activation (MIA) are central to identifying the biological mechanisms that underly the association between prenatal infection and neuropsychiatric disorder susceptibility. Many studies, however, have limited their scope to protein coding genes and their role in mediating this inherent risk, while much less attention has been directed towards exploring the roles of the epigenome and transposable elements (TEs). In Experiment 1, we demonstrate the ability of MIA to alter the chromatin landscape of the placenta. We induced MIA by injecting 200 μg/kg (i.p.) of lipopolysaccharide (LPS) on gestational day 15 in Sprague-Dawley rats. We found a sex-specific rearrangement of heterochromatin 24-h after exposure to MIA, as evidenced by an increase in histone-3 lysine-9 trimethylation (H3K9me3). In Experiment 2, MIA was associated with long-term sensorimotor processing deficits as indicated by reduced prepulse inhibition (PPI) of the acoustic startle reflex in adult male and female offspring and an increased mechanical allodynia threshold in males. Analyses of gene expression within the hypothalamus-chosen for its involvement in the sex-specific pathogenesis of schizophrenia and the stress response-revealed significantly higher levels of the stress-sensitive genes Gr and Fkbp5. Deleterious TE expression is often a hallmark of neuropsychiatric disease and we found sex-specific increases in the expression of several TEs including IAP, B2 SINE, and LINE-1 ORF1. The data from this study warrant the future consideration of chromatin stability and TEs as part of the mechanism that drives MIA-associated changes in the brain and behavior.

Evolutionarily recent retrotransposons contribute to schizophrenia

Transposable Elements (TEs) are mobile genetic elements that constitute half of the human genome. Recent studies suggest that polymorphic non-reference TEs (nrTEs) may contribute to cognitive diseases, such as schizophrenia, through a cis-regulatory effect. The aim of this work is to identify sets of nrTEs putatively linked to an increased risk of developing schizophrenia. To do so, we inspected the nrTE content of genomes from the Dorsolateral Prefrontal Cortex of schizophrenic and control individuals, and identified 38 nrTEs which possibly contribute to the emergence of this psychiatric disorder. Furthermore, we performed in silico functional inferences and found, for instance, that 9 of the 38 nrTEs act as expression/alternative splicing quantitative trait loci (eQTLs/sQTLs) in the brain, suggesting a possible role in shaping the human cognitive genome structure. Therefore, to our knowledge, this is the first attempt at identifying polymorphic nrTEs that can contribute to the functionality of the brain. Finally, we suggest that a neurodevelopmental genetic mechanism, which involves evolutionarily young nrTEs, can be the key to understanding the ethiopathogenesis of this complex disorder.

The landscape of hervRNAs transcribed from human endogenous retroviruses across human body sites

BACKGROUND

Human endogenous retroviruses (HERVs), the remnants of ancient retroviruses, account for 8% of the human genome, but most have lost their transcriptional abilities under physiological conditions. However, mounting evidence shows that several expressed HERVs do exert biological functions. Here, we systematically characterize physiologically expressed HERVs and examine whether they may give insight into the molecular fundamentals of human development and disease.

RESULTS

We systematically identify 13,889 expressed HERVs across normal body sites and demonstrate that they are expressed in body site-specific patterns and also by sex, ethnicity, and age. Analyzing cis-ERV-related quantitative trait loci, we find that 5435 hervRNAs are regulated by genetic variants. Combining this with a genome-wide association study, we elucidate that the dysregulation of expressed HERVs might be associated with various complex diseases, particularly neurodegenerative and psychiatric diseases. We further find that physiologically activated hervRNAs are associated with histone modifications rather than DNA demethylation.

CONCLUSIONS

Our results present a locus-specific landscape of physiologically expressed hervRNAs, which represent a hidden layer of genetic architecture in development and disease.

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.

Bridging the Gap Between Environmental Adversity and Neuropsychiatric Disorders: The Role of Transposable Elements

Long regarded as “junk DNA,” transposable elements (TEs) have recently garnered much attention for their role in promoting genetic diversity and plasticity. While many processes involved in mammalian development require TE activity, deleterious TE insertions are a hallmark of several psychiatric disorders. Moreover, stressful events including exposure to gestational infection and trauma, are major risk factors for developing psychiatric illnesses. Here, we will provide evidence demonstrating the intersection of stressful events, atypical TE expression, and their epigenetic regulation, which may explain how neuropsychiatric phenotypes manifest. In this way, TEs may be the “bridge” between environmental perturbations and psychopathology.

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

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

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.

Novel open reading frames in human accelerated regions and transposable elements reveal new leads to understand schizophrenia and bipolar disorder

Schizophrenia (SCZ) and bipolar disorder are debilitating neuropsychiatric disorders arising from a combination of environmental and genetic factors. Novel open reading frames (nORFs) are genomic loci that give rise to previously uncharacterized transcripts and protein products. In our previous work, we have shown that nORFs can be biologically regulated and that they may play a role in cancer and rare diseases. More importantly, we have shown that nORFs may emerge in accelerated regions of the genome giving rise to species-specific functions. We hypothesize that nORFs represent a potentially important group of biological factors that may contribute to SCZ and bipolar disorder pathophysiology. Human accelerated regions (HARs) are genomic features showing human-lineage-specific rapid evolution that may be involved in biological regulation and have additionally been found to associate with SCZ genes. Transposable elements (TEs) are another set of genomic features that have been shown to regulate gene expression. As with HARs, their relevance to SCZ has also been suggested. Here, nORFs are investigated in the context of HARs and TEs. This work shows that nORFs whose expression is disrupted in SCZ and bipolar disorder are in close proximity to HARs and TEs and that some of them are significantly associated with SCZ and bipolar disorder genomic hotspots. We also show that nORF encoded proteins can form structures and potentially constitute novel drug targets.

Corticosterone dynamically regulates retrotransposable element expression in the rat hippocampus and C6 cells

The hippocampus is a highly plastic brain region sensitive to environmental stress. It shows dynamic changes in epigenetic marks associated with stress related learning. Previous work has shown that acute stress induces substantial transient changes in histone H3 lysine 9 trimethylation (H3K9me3). Moreover, increased H3K9me3 is enriched in hippocampal gene deserts accumulating within endogenous retroviruses and transposable elements. We have found that in response to acute glucocorticoid treatment, a similar change in global H3K9me3 is observed. However, when localized we found that H3K9me3 is markedly decreased at B2 short interspersed nuclear elements but not within intracisternal-A particle endogenous retroviruses. Further, decreased H3K9me3 valence within B2 elements was associated with increased transcript abundance. These data demonstrate the capacity for acute glucocorticoids to mobilize transposable elements via epigenetic unmasking. Reconciled with previous findings following acute stress, this suggests the capacity for mobile elements to potentially function as novel regulators given their dynamic regulation by stress and glucocorticoids.

Transposable elements and their KZFP controllers are drivers of transcriptional innovation in the developing human brain

Transposable elements (TEs) account for more than 50% of the human genome and many have been co-opted throughout evolution to provide regulatory functions for gene expression networks. Several lines of evidence suggest that these networks are fine-tuned by the largest family of TE controllers, the KRAB-containing zinc finger proteins (KZFPs). One tissue permissive for TE transcriptional activation (termed "transposcription") is the adult human brain, however comprehensive studies on the extent of this process and its potential contribution to human brain development are lacking. To elucidate the spatiotemporal transposcriptome of the developing human brain, we have analyzed two independent RNA-seq data sets encompassing 16 brain regions from eight weeks postconception into adulthood. We reveal a distinct KZFP:TE transcriptional profile defining the late prenatal to early postnatal transition, and the spatiotemporal and cell type-specific activation of TE-derived alternative promoters driving the expression of neurogenesis-associated genes. Long-read sequencing confirmed these TE-driven isoforms as significant contributors to neurogenic transcripts. We also show experimentally that a co-opted antisense L2 element drives temporal protein relocalization away from the endoplasmic reticulum, suggestive of novel TE dependent protein function in primate evolution. This work highlights the widespread dynamic nature of the spatiotemporal KZFP:TE transcriptome and its importance throughout TE mediated genome innovation and neurotypical human brain development. To facilitate interactive exploration of these spatiotemporal gene and TE expression dynamics, we provide the "Brain TExplorer" web application freely accessible for the community.

Retrotransposons as Drivers of Mammalian Brain Evolution

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

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.

LINE retrotransposons characterize mammalian tissue-specific and evolutionarily dynamic regulatory regions

BACKGROUND

To investigate the mechanisms driving regulatory evolution across tissues, we experimentally mapped promoters, enhancers, and gene expression in the liver, brain, muscle, and testis from ten diverse mammals.

RESULTS

The regulatory landscape around genes included both tissue-shared and tissue-specific regulatory regions, where tissue-specific promoters and enhancers evolved most rapidly. Genomic regions switching between promoters and enhancers were more common across species, and less common across tissues within a single species. Long Interspersed Nuclear Elements (LINEs) played recurrent evolutionary roles: LINE L1s were associated with tissue-specific regulatory regions, whereas more ancient LINE L2s were associated with tissue-shared regulatory regions and with those switching between promoter and enhancer signatures across species.

CONCLUSIONS

Our analyses of the tissue-specificity and evolutionary stability among promoters and enhancers reveal how specific LINE families have helped shape the dynamic mammalian regulome.

Genomic and Transcriptomic Survey Provides New Insight into the Organization and Transposition Activity of Highly Expressed LTR Retrotransposons of Sunflower (Helianthus annuus L.)

LTR retrotransposons (RTEs) play a crucial role in plant genome evolution and adaptation. Although RTEs are generally silenced in somatic plant tissues under non-stressed conditions, some expressed RTEs (exRTEs) escape genome defense mechanisms. As our understanding of exRTE organization in plants is rudimentary, we systematically surveyed the genomic and transcriptomic organization and mobilome (transposition) activity of sunflower (Helianthus annuus L.) exRTEs. We identified 44 transcribed RTEs in the sunflower genome and demonstrated their distinct genomic features: more recent insertion time, longer open reading frame (ORF) length, and smaller distance to neighboring genes. We showed that GAG-encoding ORFs are present at significantly higher frequencies in exRTEs, compared with non-expressed RTEs. Most exRTEs exhibit variation in copy number among sunflower cultivars and one exRTE Gagarin produces extrachromosomal circular DNA in seedling, demonstrating recent and ongoing transposition activity. Nanopore direct RNA sequencing of full-length RTE RNA revealed complex patterns of alternative splicing in RTE RNAs, resulting in isoforms that carry ORFs for distinct RTE proteins. Together, our study demonstrates that tens of expressed sunflower RTEs with specific genomic organization shape the hidden layer of the transcriptome, pointing to the evolution of specific strategies that circumvent existing genome defense mechanisms.

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.

Stress, Adaptation, and the Deep Genome: Why Transposons Matter

Stress is a common, if often unpredictable life event. It can be defined from an evolutionary perspective as a force an organism perceives it must adapt to. Thus stress is a useful tool to study adaptation and the adaptive capacity of organisms. The deep genome, long neglected as a pile of “junk” has emerged as a source of regulatory DNA and RNA as well as a potential stockpile of adaptive capacity at the organismal and species levels. Recent work on the regulation of transposable elements (TEs), the principle constituents of the deep genome, by stress has shown that these elements are responsive to host stress and other environmental cues. Further, we have shown that some are likely directly regulated by the glucocorticoid receptor (GR), one of the two major vertebrate stress steroid receptors in a fashion that appears adaptive. On the basis of this and other emerging evidence I argue that the deep genome may represent an adaptive toolkit for organisms to respond to their environments at both individual and evolutionary scales. This argues that genomes may be adapted for what Waddington called “trait adaptability” rather than being purely passive objects of natural selection and single nucleotide level mutation.

Tripartite Combination of Candidate Pandemic Mitigation Agents: Vitamin D, Quercetin, and Estradiol Manifest Properties of Medicinal Agents for Targeted Mitigation of the COVID-19 Pandemic Defined by Genomics-Guided Tracing of SARS-CoV-2 Targets in Human Cells

Genes required for SARS-CoV-2 entry into human cells, ACE2 and FURIN, were employed as baits to build genomic-guided molecular maps of upstream regulatory elements, their expression and functions in the human body, and pathophysiologically relevant cell types. Repressors and activators of the ACE2 and FURIN genes were identified based on the analyses of gene silencing and overexpression experiments as well as relevant transgenic mouse models. Panels of repressors (VDR; GATA5; SFTPC; HIF1a) and activators (HMGA2; INSIG1; RUNX1; HNF4a; JNK1/c-FOS) were then employed to identify existing drugs manifesting in their effects on gene expression signatures of potential coronavirus infection mitigation agents. Using this strategy, vitamin D and quercetin have been identified as putative 2019 coronavirus disease (COVID-19) mitigation agents. Quercetin has been identified as one of top-scoring candidate therapeutics in the supercomputer SUMMIT drug-docking screen and Gene Set Enrichment Analyses (GSEA) of expression profiling experiments (EPEs), indicating that highly structurally similar quercetin, luteolin, and eriodictyol could serve as scaffolds for the development of efficient inhibitors of SARS-CoV-2 infection. In agreement with this notion, quercetin alters the expression of 98 of 332 (30%) of human genes encoding protein targets of SARS-CoV-2, thus potentially interfering with functions of 23 of 27 (85%) of the SARS-CoV-2 viral proteins in human cells. Similarly, Vitamin D may interfere with functions of 19 of 27 (70%) of the SARS-CoV-2 proteins by altering expression of 84 of 332 (25%) of human genes encoding protein targets of SARS-CoV-2. Considering the potential effects of both quercetin and vitamin D, the inference could be made that functions of 25 of 27 (93%) of SARS-CoV-2 proteins in human cells may be altered. GSEA and EPEs identify multiple drugs, smoking, and many disease conditions that appear to act as putative coronavirus infection-promoting agents. Discordant patterns of testosterone versus estradiol impacts on SARS-CoV-2 targets suggest a plausible molecular explanation of the apparently higher male mortality during the coronavirus pandemic. Estradiol, in contrast with testosterone, affects the expression of the majority of human genes (203 of 332; 61%) encoding SARS-CoV-2 targets, thus potentially interfering with functions of 26 of 27 SARS-CoV-2 viral proteins. A hypothetical tripartite combination consisting of quercetin/vitamin D/estradiol may affect expression of 244 of 332 (73%) human genes encoding SARS-CoV-2 targets. Of major concern is the ACE2 and FURIN expression in many human cells and tissues, including immune cells, suggesting that SARS-CoV-2 may infect a broad range of cellular targets in the human body. Infection of immune cells may cause immunosuppression, long-term persistence of the virus, and spread of the virus to secondary targets. Present analyses and numerous observational studies indicate that age-associated vitamin D deficiency may contribute to the high mortality of older adults and the elderly. Immediate availability for targeted experimental and clinical interrogations of potential COVID-19 pandemic mitigation agents, namely vitamin D and quercetin, as well as of the highly selective (Ki, 600 pm) intrinsically specific FURIN inhibitor (a1-antitrypsin Portland (a1-PDX), is considered an encouraging factor. Observations reported in this contribution are intended to facilitate follow-up targeted experimental studies and, if warranted, randomized clinical trials to identify and validate therapeutically viable interventions to combat the COVID-19 pandemic. Specifically, gene expression profiles of vitamin D and quercetin activities and their established safety records as over-the-counter medicinal substances strongly argue that they may represent viable candidates for further considerations of their potential utility as COVID-19 pandemic mitigation agents. In line with the results of present analyses, a randomized interventional clinical trial evaluating effects of estradiol on severity of the coronavirus infection in COVID19+ and presumptive COVID19+ patients and two interventional randomized clinical trials evaluating effects of vitamin D on prevention and treatment of COVID-19 were listed on the ClinicalTrials.gov website.

Mobile genomics: tools and techniques for tackling transposons

Next-generation sequencing approaches have fundamentally changed the types of questions that can be asked about gene function and regulation. With the goal of approaching truly genome-wide quantifications of all the interaction partners and downstream effects of particular genes, these quantitative assays have allowed for an unprecedented level of detail in exploring biological interactions. However, many challenges remain in our ability to accurately describe and quantify the interactions that take place in those hard to reach and extremely repetitive regions of our genome comprised mostly of transposable elements (TEs). Tools dedicated to TE-derived sequences have lagged behind, making the inclusion of these sequences in genome-wide analyses difficult. Recent improvements, both computational and experimental, allow for the better inclusion of TE sequences in genomic assays and a renewed appreciation for the importance of TE biology. This review will discuss the recent improvements that have been made in the computational analysis of TE-derived sequences as well as the areas where such analysis still proves difficult. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

In search of optimal resilience ratios: Differential influences of neurobehavioral factors contributing to stress-resilience spectra

The ability to adapt to stressful circumstances, known as emotional resilience, is a key factor in the maintenance of mental health. Several individual biomarkers of the stress response (e.g., corticosterone) that influence an animal's position along the continuum that ranges from adaptive allostasis to maladaptive allostatic load have been identified. Extending beyond specific biomarkers of stress responses, however, it is also important to consider stress-related responses relative to other relevant responses for a thorough understanding of the underpinnings of adaptive allostasis. In this review, behavioral, neurobiological, developmental and genomic variables are considered in the context of emotional resilience [e.g., explore/exploit behavioral tendencies; DHEA/CORT ratios and relative proportions of protein-coding/nonprotein-coding (transposable) genomic elements]. As complex and multifaceted relationships between pertinent allostasis biomediators are identified, translational applications for optimal resilience are more likely to emerge as effective therapeutic strategies.

Different Families of Retrotransposons and DNA Transposons Are Actively Transcribed and May Have Transposed Recently in Physcomitrium (Physcomitrella) patens

Similarly to other plant genomes of similar size, more than half of the genome of P. patens is covered by Transposable Elements (TEs). However, the composition and distribution of P. patens TEs is quite peculiar, with Long Terminal Repeat (LTR)-retrotransposons, which form patches of TE-rich regions interleaved with gene-rich regions, accounting for the vast majority of the TE space. We have already shown that RLG1, the most abundant TE in P. patens, is expressed in non-stressed protonema tissue. Here we present a non-targeted analysis of the TE expression based on RNA-Seq data and confirmed by qRT-PCR analyses that shows that, at least four LTR-RTs (RLG1, RLG2, RLC4 and tRLC5) and one DNA transposon (PpTc2) are expressed in P. patens. These TEs are expressed during development or under stresses that P. patens frequently faces, such as dehydratation/rehydratation stresses, suggesting that TEs have ample possibilities to transpose during P. patens life cycle. Indeed, an analysis of the TE polymorphisms among four different P. patens accessions shows that different TE families have recently transposed in this species and have generated genetic variability that may have phenotypic consequences, as a fraction of the TE polymorphisms are within or close to genes. Among the transcribed and mobile TEs, tRLC5 is particularly interesting as it concentrates in a single position per chromosome that could coincide with the centromere, and its expression is specifically induced in young sporophyte, where meiosis takes place.

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.

Early life exposures, neurodevelopmental disorders, and transposable elements

Transposable elements make up a much larger portion of the genome than protein-coding genes, yet we know relatively little about their function in the human genome. However, we are beginning to more fully understand their role in brain development, neuroinflammation, and adaptation to environmental insults such as stress. For instance, glucocorticoid receptor activation regulates transposable elements in the brain following acute stress. Early life is a period of substantial brain development during which transposable elements play a role. Environmental exposures and experiences during early life that promote abnormal regulation of transposable elements may lead to a cascade of events that ultimately increase susceptibility to disorders later in life. Recent attention to transposable elements in psychiatric illness has begun to clarify associations indicative of dysregulation of different classes of transposable elements in stress-related and neurodevelopmental illness. Though individual susceptibility or resiliency to psychiatric illness has not been explained by traditional genetic studies, the wide inter-individual variability in transposable element composition in the human genome make TEs attractive candidates to elucidate this differential susceptibility. In this review, we discuss evidence that regulation of transposable elements in the brain are stage-specific, sensitive to environmental factors, and may be impacted by early life perturbations. We further present evidence of associations with stress-related and neurodevelopmental psychiatric illness from a developmental perspective.

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.

Enhanced Molecular Appreciation of Psychiatric Disorders Through High-Dimensionality Data Acquisition and Analytics

The initial diagnosis, molecular investigation, treatment, and posttreatment care of major psychiatric disorders (schizophrenia and bipolar depression) are all still significantly hindered by the current inability to define these disorders in an explicit molecular signaling manner. High-dimensionality data analytics, using large datastreams from transcriptomic, proteomic, or metabolomic investigations, will likely advance both the appreciation of the molecular nature of major psychiatric disorders and simultaneously enhance our ability to more efficiently diagnose and treat these debilitating conditions. High-dimensionality data analysis in psychiatric research has been heterogeneous in aims and methods and limited by insufficient sample sizes, poorly defined case definitions, methodological inhomogeneity, and confounding results. All of these issues combine to constrain the conclusions that can be extracted from them. Here, we discuss possibilities for overcoming methodological challenges through the implementation of transcriptomic, proteomic, or metabolomics signatures in psychiatric diagnosis and offer an outlook for future investigations. To fulfill the promise of intelligent high-dimensionality data-based differential diagnosis in mental disease diagnosis and treatment, future research will need large, well-defined cohorts in combination with state-of-the-art technologies.