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

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.

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.

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.

Inhibiting tumour metastasis by DQA modified paclitaxel plus ligustrazine micelles in treatment of non-small-cell lung cancer

Lung cancer is a kind of malignant tumour characterized as uncontrolled cell growth in lung. These malignant cell growth can spread beyond the lung by process of metastasis into other tissues or parts of the body. In this study, we developed dequalinium (DQA) modified paclitaxel plus ligustrazine micelles to destroy vasculogenic mimicry (VM) channels and inhibit tumour metastasis. In vitro assays showed that the targeting micelles with centralized particle size distribution showed not only vigoroso cytotoxicity on A549 cells but also strong inhibition on VM channels and tumour metastasis. Mechanism studies indicated that the DQA modified paclitaxel plus ligustrazine micelles could down-regulate the expressions of VEGF, MMP2, TGF-β1 and E-cadherin in A549 cells. In vivo assays indicated that the targeting drug-loaded micelles could enhance the accumulation of chemotherapeutic drugs at tumour sites and exhibit strong tumour inhibitory activity with negligible toxicity. Hence, the DQA modified paclitaxel plus ligustrazine micelles developed in this study may provide a potential strategy for treatment of NSCLC.

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 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.

Promoter expression of HERV-K (HML-2) provirus-derived sequences is related to LTR sequence variation and polymorphic transcription factor binding sites

Background

Increased transcription of the human endogenous retrovirus group HERV-K (HML-2) is often seen during disease. Although the mechanism of its tissue-specific activation is unclear, research shows that LTR CpG hypomethylation alone is not sufficient to induce its promoter activity and that the transcriptional milieu of a malignant cell contributes, at least partly, to differential HML-2 expression.

Results

We analyzed the relationship between LTR sequence variation and promoter expression patterns in human breast cancer cell lines, finding them to be positively correlated. In particular, two proviruses (3q12.3 and 11p15.4) displayed increased activity in almost all tumorigenic cell lines sampled. Using a transcription factor binding site prediction algorithm, we identified two unique binding sites in each 5′ LTR that appeared to be associated with inducing promoter activity during neoplasia. Genomic analysis of the homologous proviruses in several non-human primates indicated post-integration genetic drift in two transcription factor binding sites, away from the ancestral sequence and towards the active form. Based on the sequences of 2504 individuals from the 1000 Genomes Project, the active form of the 11p15.4 site was found to be polymorphic within the human population, with an allele frequency of 51%, whereas the activating mutation in the 3q12.3 provirus was fixed in humans but not present in the orthologous provirus in chimpanzees or gorillas.

Conclusions

These data suggest that stage-specific transcription factors at least partly contribute to LTR promoter activity during transformation and that, in some cases, transcription factor binding site polymorphisms may be responsible for the differential HML-2 expression often seen between individuals.

Electronic supplementary material

The online version of this article (10.1186/s12977-018-0441-2) contains supplementary material, which is available to authorized users.

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.

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.

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.