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

Conserved and divergent gene regulatory programs of the mammalian neocortex

Divergence of cis-regulatory elements drives species-specific traits1, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains unclear. Here we investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset and mouse using single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome and chromosomal conformation profiles from a total of over 200,000 cells. From these data, we show evidence that divergence of transcription factor expression corresponds to species-specific epigenome landscapes. We find that conserved and divergent gene regulatory features are reflected in the evolution of the three-dimensional genome. Transposable elements contribute to nearly 80% of the human-specific candidate cis-regulatory elements in cortical cells. Through machine learning, we develop sequence-based predictors of candidate cis-regulatory elements in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Finally, we show that epigenetic conservation combined with sequence similarity helps to uncover functional cis-regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

Comparative single cell epigenomic analysis of gene regulatory programs in the rodent and primate neocortex

Sequence divergence of cis- regulatory elements drives species-specific traits, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains to be elucidated. We investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset, and mouse with single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome, and chromosomal conformation profiles from a total of over 180,000 cells. For each modality, we determined species-specific, divergent, and conserved gene expression and epigenetic features at multiple levels. We find that cell type-specific gene expression evolves more rapidly than broadly expressed genes and that epigenetic status at distal candidate cis -regulatory elements (cCREs) evolves faster than promoters. Strikingly, transposable elements (TEs) contribute to nearly 80% of the human-specific cCREs in cortical cells. Through machine learning, we develop sequence-based predictors of cCREs in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Lastly, we show that epigenetic conservation combined with sequence similarity helps uncover functional cis -regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

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.

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.

Triazole Modified Tetraiodothyroacetic Acid Conjugated to Polyethylene Glycol, a Thyrointegrin αvβ3 Antagonist as a Radio- and Chemo-Sensitizer in Pancreatic Cancer

Thyroid hormone L thyroxine stimulates pancreatic carcinoma cell proliferation via thyrointegrin αvβ3 receptors, and antagonist tetraiodothyroacetic acid (tetrac) inhibits cancer cell growth. Chemically modified bis-triazole-tetrac conjugated with polyethylene glycol (P-bi-TAT) has higher binding affinity to αvβ3 receptors compared to tetrac. We investigated the antiproliferation effect of P-bi-TAT in pancreatic cancer cells (SUIT2) and its radio- and chemo-sensitizing roles in a mouse model of pancreatic cancer. P-bi-TAT treatment increased tumor-targeted radiation-induced cell death and decreased tumor size. P-bi-TAT acted as a chemo-sensitizer and enhanced the 5-fluorouracil (5FU) effect in decreasing pancreatic tumor weight compared to 5FU monotherapy. Withdrawal of treatment continued the tumor regression; however, the 5FU group showed tumor regrowth. The mechanisms of the anti-cancer activity of P-bi-TAT on SUIT2 cells were assessed by microarray experiments, and genome-wide profiling identified significant alterations of 1348 genes’ expression. Both down-regulated and up-regulated transcripts suggest that a molecular interference at the signaling pathway-associated gene expression is the prevalent mode of P-bi-TAT anti-cancer activity. Our data indicate that non-cytotoxic P-bi-TAT is not only an anti-cancer agent but also a radio-sensitizer and chemo-sensitizer that acts on the extracellular domain of the cell surface αvβ3 receptor.

Targeting Thyrointegrin αvβ3 Using Fluorobenzyl Polyethylene Glycol Conjugated Tetraiodothyroacetic Acid (NP751) in Acute Myeloid Leukemia

Background

Acute myeloid leukemia (AML) is associated with poor long-term survival, even with newer therapeutic agents. Here, we show the results of our preclinical study, in which we evaluated the efficacy of a new thyrointegrin αvβ3 antagonist, named fluorobenzyl polyethylene glycol conjugated tetraiodothyroacetic acid (fb-PMT).

Methods and Results

fb-PMT (NP751) is a potent αvβ3 antagonist of molecular weight of 2,478.9 Da. it represents a conjugate of tetraiodothyroacetic acid (tetrac) and monodisperse polyethylene glycol (PEG36), with a 4-fluorobenzyl group capping the other end of the PEG. fb-PMT effectively suppresses the malignant growth of human acute myeloid leukemia (AML) after successful engraftment in transgenic NSG-S xenograft mouse models of either established human AML cell line or primary AML cells. Daily treatment with fb-PMT (1–10 mg/kg body weight) subcutaneously (s.c.) for 3–4 weeks was associated with marked regression of leukemogenesis and extended survival in both models. The efficiency of the fb-PMT therapy was verified using in vivo imaging system (IVIS) imaging, flow cytometry, and histopathological examination to monitor the engraftment of leukemic cells in the bone marrow and other organs. fb-PMT therapy for 3–4 weeks at 3 and 10 mg/kg daily doses exhibited significant reduction (p < 0.0001) of leukemic cell burden of 74% and >95%, respectively. All fb-PMT-treated mice in the 10 mg/kg treatment arm successfully maintained remission after discontinuing the daily treatment. Comprehensive fb-PMT safety assessments demonstrated excellent safety and tolerability at multiple folds above the anticipated human therapeutic doses. Lastly, our genome-wide microarray screens demonstrated that fb-PMT works through the molecular interference mechanism with multiple signaling pathways contributing to growth and survival of leukemic cells.

Conclusion

Our preclinical findings of the potent anticancer activities of fb-PMT and its favorable safety profiles warrant its clinical investigation for the effective and safe management of AML.

Comprehensive multi-omics integration identifies differentially active enhancers during human brain development with clinical relevance

BACKGROUND

Non-coding regulatory elements (NCREs), such as enhancers, play a crucial role in gene regulation, and genetic aberrations in NCREs can lead to human disease, including brain disorders. The human brain is a complex organ that is susceptible to numerous disorders; many of these are caused by genetic changes, but a multitude remain currently unexplained. Understanding NCREs acting during brain development has the potential to shed light on previously unrecognized genetic causes of human brain disease. Despite immense community-wide efforts to understand the role of the non-coding genome and NCREs, annotating functional NCREs remains challenging.

METHODS

Here we performed an integrative computational analysis of virtually all currently available epigenome data sets related to human fetal brain.

RESULTS

Our in-depth analysis unravels 39,709 differentially active enhancers (DAEs) that show dynamic epigenomic rearrangement during early stages of human brain development, indicating likely biological function. Many of these DAEs are linked to clinically relevant genes, and functional validation of selected DAEs in cell models and zebrafish confirms their role in gene regulation. Compared to enhancers without dynamic epigenomic rearrangement, DAEs are subjected to higher sequence constraints in humans, have distinct sequence characteristics and are bound by a distinct transcription factor landscape. DAEs are enriched for GWAS loci for brain-related traits and for genetic variation found in individuals with neurodevelopmental disorders, including autism.

CONCLUSION

This compendium of high-confidence enhancers will assist in deciphering the mechanism behind developmental genetics of human brain and will be relevant to uncover missing heritability in human genetic brain disorders.

Evolution of mouse circadian enhancers from transposable elements

BACKGROUND

Transposable elements are increasingly recognized as a source of cis-regulatory variation. Previous studies have revealed that transposons are often bound by transcription factors and some have been co-opted into functional enhancers regulating host gene expression. However, the process by which transposons mature into complex regulatory elements, like enhancers, remains poorly understood. To investigate this process, we examined the contribution of transposons to the cis-regulatory network controlling circadian gene expression in the mouse liver, a well-characterized network serving an important physiological function.

RESULTS

ChIP-seq analyses reveal that transposons and other repeats contribute ~ 14% of the binding sites for core circadian regulators (CRs) including BMAL1, CLOCK, PER1/2, and CRY1/2, in the mouse liver. RSINE1, an abundant murine-specific SINE, is the only transposon family enriched for CR binding sites across all datasets. Sequence analyses and reporter assays reveal that the circadian regulatory activity of RSINE1 stems from the presence of imperfect CR binding motifs in the ancestral RSINE1 sequence. These motifs matured into canonical motifs through point mutations after transposition. Furthermore, maturation occurred preferentially within elements inserted in the proximity of ancestral CR binding sites. RSINE1 also acquired motifs that recruit nuclear receptors known to cooperate with CRs to regulate circadian gene expression specifically in the liver.

CONCLUSIONS

Our results suggest that the birth of enhancers from transposons is predicated both by the sequence of the transposon and by the cis-regulatory landscape surrounding their genomic integration site.

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.

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.

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.

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.

Beyond the Exome: The Non-coding Genome and Enhancers in Neurodevelopmental Disorders and Malformations of Cortical Development

The development of the human cerebral cortex is a complex and dynamic process, in which neural stem cell proliferation, neuronal migration, and post-migratory neuronal organization need to occur in a well-organized fashion. Alterations at any of these crucial stages can result in malformations of cortical development (MCDs), a group of genetically heterogeneous neurodevelopmental disorders that present with developmental delay, intellectual disability and epilepsy. Recent progress in genetic technologies, such as next generation sequencing, most often focusing on all protein-coding exons (e.g., whole exome sequencing), allowed the discovery of more than a 100 genes associated with various types of MCDs. Although this has considerably increased the diagnostic yield, most MCD cases remain unexplained. As Whole Exome Sequencing investigates only a minor part of the human genome (1–2%), it is likely that patients, in which no disease-causing mutation has been identified, could harbor mutations in genomic regions beyond the exome. Even though functional annotation of non-coding regions is still lagging behind that of protein-coding genes, tremendous progress has been made in the field of gene regulation. One group of non-coding regulatory regions are enhancers, which can be distantly located upstream or downstream of genes and which can mediate temporal and tissue-specific transcriptional control via long-distance interactions with promoter regions. Although some examples exist in literature that link alterations of enhancers to genetic disorders, a widespread appreciation of the putative roles of these sequences in MCDs is still lacking. Here, we summarize the current state of knowledge on cis-regulatory regions and discuss novel technologies such as massively-parallel reporter assay systems, CRISPR-Cas9-based screens and computational approaches that help to further elucidate the emerging role of the non-coding genome in disease. Moreover, we discuss existing literature on mutations or copy number alterations of regulatory regions involved in brain development. We foresee that the future implementation of the knowledge obtained through ongoing gene regulation studies will benefit patients and will provide an explanation to part of the missing heritability of MCDs and other genetic disorders.