Evolutionary framework of the human interactome: Unicellular and multicellular giant clusters

Polyploidy Promotes Hypertranscription, Apoptosis Resistance, and Ciliogenesis in Cancer Cells and Mesenchymal Stem Cells of Various Origins: Comparative Transcriptome In Silico Study

Mesenchymal stem cells (MSC) attract an increasing amount of attention due to their unique therapeutic properties. Yet, MSC can undergo undesirable genetic and epigenetic changes during their propagation in vitro. In this study, we investigated whether polyploidy can compromise MSC oncological safety and therapeutic properties. For this purpose, we compared the impact of polyploidy on the transcriptome of cancer cells and MSC of various origins (bone marrow, placenta, and heart). First, we identified genes that are consistently ploidy-induced or ploidy-repressed through all comparisons. Then, we selected the master regulators using the protein interaction enrichment analysis (PIEA). The obtained ploidy-related gene signatures were verified using the data gained from polyploid and diploid populations of early cardiomyocytes (CARD) originating from iPSC. The multistep bioinformatic analysis applied to the cancer cells, MSC, and CARD indicated that polyploidy plays a pivotal role in driving the cell into hypertranscription. It was evident from the upregulation of gene modules implicated in housekeeping functions, stemness, unicellularity, DNA repair, and chromatin opening by means of histone acetylation operating via DNA damage associated with the NUA4/TIP60 complex. These features were complemented by the activation of the pathways implicated in centrosome maintenance and ciliogenesis and by the impairment of the pathways related to apoptosis, the circadian clock, and immunity. Overall, our findings suggest that, although polyploidy does not induce oncologic transformation of MSC, it might compromise their therapeutic properties because of global epigenetic changes and alterations in fundamental biological processes. The obtained results can contribute to the development and implementation of approaches enhancing the therapeutic properties of MSC by removing polyploid cells from the cell population.

Long-Term Transcriptomic Changes and Cardiomyocyte Hyperpolyploidy after Lactose Intolerance in Neonatal Rats

Many cardiovascular diseases originate from growth retardation, inflammation, and malnutrition during early postnatal development. The nature of this phenomenon is not completely understood. Here we aimed to verify the hypothesis that systemic inflammation triggered by neonatal lactose intolerance (NLI) may exert long-term pathologic effects on cardiac developmental programs and cardiomyocyte transcriptome regulation. Using the rat model of NLI triggered by lactase overloading with lactose and the methods of cytophotometry, image analysis, and mRNA-seq, we evaluated cardiomyocyte ploidy, signs of DNA damage, and NLI-associated long-term transcriptomic changes of genes and gene modules that differed qualitatively (i.e., were switched on or switched off) in the experiment vs. the control. Our data indicated that NLI triggers the long-term animal growth retardation, cardiomyocyte hyperpolyploidy, and extensive transcriptomic rearrangements. Many of these rearrangements are known as manifestations of heart pathologies, including DNA and telomere instability, inflammation, fibrosis, and reactivation of fetal gene program. Moreover, bioinformatic analysis identified possible causes of these pathologic traits, including the impaired signaling via thyroid hormone, calcium, and glutathione. We also found transcriptomic manifestations of increased cardiomyocyte polyploidy, such as the induction of gene modules related to open chromatin, e.g., "negative regulation of chromosome organization", "transcription" and "ribosome biogenesis". These findings suggest that ploidy-related epigenetic alterations acquired in the neonatal period permanently rewire gene regulatory networks and alter cardiomyocyte transcriptome. Here we provided first evidence indicating that NLI can be an important trigger of developmental programming of adult cardiovascular disease. The obtained results can help to develop preventive strategies for reducing the NLI-associated adverse effects of inflammation on the developing cardiovascular system.

Systemic Alterations of Cancer Cells and Their Boost by Polyploidization: Unicellular Attractor (UCA) Model

Using meta-analyses, we introduce a unicellular attractor (UCA) model integrating essential features of the ‘atavistic reversal’, ‘cancer attractor’, ‘somatic mutation’, ‘genome chaos’, and ‘tissue organization field’ theories. The ‘atavistic reversal’ theory is taken as a keystone. We propose a possible mechanism of this reversal, its refinement called ‘gradual atavism’, and evidence for the ‘serial atavism’ model. We showed the gradual core-to-periphery evolutionary growth of the human interactome resulting in the higher protein interaction density and global interactome centrality in the UC center. In addition, we revealed that UC genes are more actively expressed even in normal cells. The modeling of random walk along protein interaction trajectories demonstrated that random alterations in cellular networks, caused by genetic and epigenetic changes, can result in a further gradual activation of the UC center. These changes can be induced and accelerated by cellular stress that additionally activates UC genes (especially during cell proliferation), because the genes involved in cellular stress response and cell cycle are mostly of UC origin. The functional enrichment analysis showed that cancer cells demonstrate the hyperactivation of energetics and the suppression of multicellular genes involved in communication with the extracellular environment (especially immune surveillance). Collectively, these events can unleash selfish cell behavior aimed at survival at all means. All these changes are boosted by polyploidization. The UCA model may facilitate an understanding of oncogenesis and promote the development of therapeutic strategies.

Gradistics: An underappreciated dimension in evolutionary space

The growth of complexity is an unsolved and underappreciated problem. We consider possible causes of this growth, hypotheses testing, molecular mechanisms, complexity measures, cases of simplification, and significance for biomedicine. We focus on a general ability of regulation, which is based on the growing information storage and processing capacities, as the main proxy of complexity. Natural selection is indifferent to complexity. However, complexification can be inferred from the same first principle, on which natural selection is founded. Natural selection depends on potentially unlimited reproduction under limited environmental conditions. Because of the demographic pressure, the simple ecological niches become fulfilled and diversified (due to species splitting and divergence). Diversification increases complexity of biocenoses. After the filling and diversification of simple niches, the more complex niches can arise. This is the 'atomic orbitals' (AO) model. Complexity has many shortcomings but it has an advantage. This advantage is ability to regulatory adaptation, including behavioral, formed in the evolution by means of genetic adaptation. Regulatory adaptation is much faster than genetic one because it is based on the information previously accumulated via genetic adaptation and learning. Regulatory adaptation further increases complexity of biocenoses. This is the 'regulatory advantage' (RA) model. The comparison of both models allows testable predictions. We focus on the animal evolution because of the appearance of higher regulatory level (nervous system), which is absent in other lineages, and relevance to humans (including biomedical aspects).

Polyploidy and Myc Proto-Oncogenes Promote Stress Adaptation via Epigenetic Plasticity and Gene Regulatory Network Rewiring

Polyploid cells demonstrate biological plasticity and stress adaptation in evolution; development; and pathologies, including cardiovascular diseases, neurodegeneration, and cancer. The nature of ploidy-related advantages is still not completely understood. Here, we summarize the literature on molecular mechanisms underlying ploidy-related adaptive features. Polyploidy can regulate gene expression via chromatin opening, reawakening ancient evolutionary programs of embryonality. Chromatin opening switches on genes with bivalent chromatin domains that promote adaptation via rapid induction in response to signals of stress or morphogenesis. Therefore, stress-associated polyploidy can activate Myc proto-oncogenes, which further promote chromatin opening. Moreover, Myc proto-oncogenes can trigger polyploidization de novo and accelerate genome accumulation in already polyploid cells. As a result of these cooperative effects, polyploidy can increase the ability of cells to search for adaptive states of cellular programs through gene regulatory network rewiring. This ability is manifested in epigenetic plasticity associated with traits of stemness, unicellularity, flexible energy metabolism, and a complex system of DNA damage protection, combining primitive error-prone unicellular repair pathways, advanced error-free multicellular repair pathways, and DNA damage-buffering ability. These three features can be considered important components of the increased adaptability of polyploid cells. The evidence presented here contribute to the understanding of the nature of stress resistance associated with ploidy and may be useful in the development of new methods for the prevention and treatment of cardiovascular and oncological diseases.

Life Entrapped in a Network of Atavistic Attractors: How to Find a Rescue

In view of unified cell bioenergetics, cell bioenergetic problems related to cell overenergization can cause excessive disturbances in current cell fate and, as a result, lead to a change of cell-fate. At the onset of the problem, cell overenergization of multicellular organisms (especially overenergization of mitochondria) is solved inter alia by activation and then stimulation of the reversible Crabtree effect by cells. Unfortunately, this apparently good solution can also lead to a much bigger problem when, despite the activation of the Crabtree effect, cell overenergization persists for a long time. In such a case, cancer transformation, along with the Warburg effect, may occur to further reduce or stop the charging of mitochondria by high-energy molecules. Understanding the phenomena of cancer transformation and cancer development has become a real challenge for humanity. To date, many models have been developed to understand cancer-related mechanisms. Nowadays, combining all these models into one coherent universal model of cancer transformation and development can be considered a new challenge. In this light, the aim of this article is to present such a potentially universal model supported by a proposed new model of cellular functionality evolution. The methods of fighting cancer resulting from unified cell bioenergetics and the two presented models are also considered.

Polyploidy as a Fundamental Phenomenon in Evolution, Development, Adaptation and Diseases

DNA replication during cell proliferation is 'vertical' copying, which reproduces an initial amount of genetic information. Polyploidy, which results from whole-genome duplication, is a fundamental complement to vertical copying. Both organismal and cell polyploidy can emerge via premature cell cycle exit or via cell-cell fusion, the latter giving rise to polyploid hybrid organisms and epigenetic hybrids of somatic cells. Polyploidy-related increase in biological plasticity, adaptation, and stress resistance manifests in evolution, development, regeneration, aging, oncogenesis, and cardiovascular diseases. Despite the prevalence in nature and importance for medicine, agri- and aquaculture, biological processes and epigenetic mechanisms underlying these fundamental features largely remain unknown. The evolutionarily conserved features of polyploidy include activation of transcription, response to stress, DNA damage and hypoxia, and induction of programs of morphogenesis, unicellularity, and longevity, suggesting that these common features confer adaptive plasticity, viability, and stress resistance to polyploid cells and organisms. By increasing cell viability, polyploidization can provide survival under stressful conditions where diploid cells cannot survive. However, in somatic cells it occurs at the expense of specific function, thus promoting developmental programming of adult cardiovascular diseases and increasing the risk of cancer. Notably, genes arising via evolutionary polyploidization are heavily involved in cancer and other diseases. Ploidy-related changes of gene expression presumably originate from chromatin modifications and the derepression of bivalent genes. The provided evidence elucidates the role of polyploidy in evolution, development, aging, and carcinogenesis, and may contribute to the development of new strategies for promoting regeneration and preventing cardiovascular diseases and cancer.

Growth of Biological Complexity from Prokaryotes to Hominids Reflected in the Human Genome

The growth of complexity in evolution is a most intriguing phenomenon. Using gene phylostratigraphy, we showed this growth (as reflected in regulatory mechanisms) in the human genome, tracing the path from prokaryotes to hominids. Generally, the different regulatory gene families expanded at different times, yet only up to the Euteleostomi (bony vertebrates). The only exception was the expansion of transcription factors (TF) in placentals; however, we argue that this was not related to increase in general complexity. Surprisingly, although TF originated in the Prokaryota while chromatin appeared only in the Eukaryota, the expansion of epigenetic factors predated the expansion of TF. Signaling receptors, tumor suppressors, oncogenes, and aging- and disease-associated genes (indicating vulnerabilities in terms of complex organization and strongly enrichment in regulatory genes) also expanded only up to the Euteleostomi. The complexity-related gene properties (protein size, number of alternative splicing mRNA, length of untranslated mRNA, number of biological processes per gene, number of disordered regions in a protein, and density of TF–TF interactions) rose in multicellular organisms and declined after the Euteleostomi, and possibly earlier. At the same time, the speed of protein sequence evolution sharply increased in the genes that originated after the Euteleostomi. Thus, several lines of evidence indicate that molecular mechanisms of complexity growth were changing with time, and in the phyletic lineage leading to humans, the most salient shift occurred after the basic vertebrate body plan was fixed with bony skeleton. The obtained results can be useful for evolutionary medicine.

Isolation and Characterization of Human Colon Adenocarcinoma Stem-Like Cells Based on the Endogenous Expression of the Stem Markers

BACKGROUND

Cancer stem cells' (CSCs) self-maintenance is regulated via the pluripotency pathways promoting the most aggressive tumor phenotype. This study aimed to use the activity of these pathways for the CSCs' subpopulation enrichment and separating cells characterized by the OCT4 and SOX2 expression.

METHODS

To select and analyze CSCs, we used the SORE6x lentiviral reporter plasmid for viral transduction of colon adenocarcinoma cells. Additionally, we assessed cell chemoresistance, clonogenic, invasive and migratory activity and the data of mRNA-seq and intrinsic disorder predisposition protein analysis (IDPPA).

RESULTS

We obtained the line of CSC-like cells selected on the basis of the expression of the OCT4 and SOX2 stem cell factors. The enriched CSC-like subpopulation had increased chemoresistance as well as clonogenic and migration activities. The bioinformatic analysis of mRNA seq data identified the up-regulation of pluripotency, development, drug resistance and phototransduction pathways, and the downregulation of pathways related to proliferation, cell cycle, aging, and differentiation. IDPPA indicated that CSC-like cells are predisposed to increased intrinsic protein disorder.

CONCLUSION

The use of the SORE6x reporter construct for CSCs enrichment allows us to obtain CSC-like population that can be used as a model to search for the new prognostic factors and potential therapeutic targets for colon cancer treatment.

Paradoxes of cancer: Survival at the brink

The fundamental understanding of how Cancer initiates, persists and then progresses is evolving. High-resolution technologies, including single-cell mutation and gene expression measurements, are now attainable, providing an ever-increasing insight into the molecular details. However, this higher resolution has shown that somatic mutation theory itself cannot explain the extraordinary resistance of cancer to extinction. There is a need for a more Systems-based framework of understanding cancer complexity, which in particular explains the regulation of gene expression during cell-fate decisions. Cancer displays a series of paradoxes. Here we attempt to approach them from the view-point of adaptive exploration of gene regulatory networks at the edge of order and chaos, where cell-fate is changed by oscillations between alternative regulators of cellular senescence and reprogramming operating through self-organisation. On this background, the role of polyploidy in accessing the phylogenetically pre-programmed "oncofetal attractor" state, related to unicellularity, and the de-selection of unsuitable variants at the brink of cell survival is highlighted. The concepts of the embryological and atavistic theory of cancer, cancer cell "life-cycle", and cancer aneuploidy paradox are dissected under this lense. Finally, we challenge researchers to consider that cancer "defects" are mostly the adaptation tools of survival programs that have arisen during evolution and are intrinsic of cancer. Recognition of these features should help in the development of more successful anti-cancer treatments.

Phylostratic Shift of Whole-Genome Duplications in Normal Mammalian Tissues towards Unicellularity Is Driven by Developmental Bivalent Genes and Reveals a Link to Cancer

Tumours were recently revealed to undergo a phylostratic and phenotypic shift to unicellularity. As well, aggressive tumours are characterized by an increased proportion of polyploid cells. In order to investigate a possible shared causation of these two features, we performed a comparative phylostratigraphic analysis of ploidy-related genes, obtained from transcriptomic data for polyploid and diploid human and mouse tissues using pairwise cross-species transcriptome comparison and principal component analysis. Our results indicate that polyploidy shifts the evolutionary age balance of the expressed genes from the late metazoan phylostrata towards the upregulation of unicellular and early metazoan phylostrata. The up-regulation of unicellular metabolic and drug-resistance pathways and the downregulation of pathways related to circadian clock were identified. This evolutionary shift was associated with the enrichment of ploidy with bivalent genes (p < 10⁻¹⁶). The protein interactome of activated bivalent genes revealed the increase of the connectivity of unicellulars and (early) multicellulars, while circadian regulators were depressed. The mutual polyploidy-c-MYC-bivalent genes-associated protein network was organized by gene-hubs engaged in both embryonic development and metastatic cancer including driver (proto)-oncogenes of viral origin. Our data suggest that, in cancer, the atavistic shift goes hand-in-hand with polyploidy and is driven by epigenetic mechanisms impinging on development-related bivalent genes.

Cell-cycle dependence of transcriptome gene modules: comparison of regression lines

The transcriptome consists of various gene modules that can be mutually dependent, and ignoring these dependencies may lead to misinterpretation. The most important problem is module dependence on cell-cycle activity. Using meta-analysis of over 30 000 single-cell transcriptomes, we show gene module dependencies on cell-cycle signature, which can be consistently observed in various normal and cancer cells. Transcript levels of receptors, plasma membrane, and differentiation-related genes are negatively regressed on cell-cycle signature. Pluripotency, stress response, DNA repair, chromatin remodeling, proteasomal protein degradation, protein network connectivity, and unicellular evolutionary origin are regressed positively. These effects cannot be explained by partial overlap of corresponding gene sets because they remain if the overlapped genes were removed. We propose a visual analysis of gene module-specific regression lines as complement to an uncurated enrichment analysis. The different lines for a same gene module indicate different cell conditions. The approach is tested on several problems (polyploidy, pluripotency, cancer, phylostratigraphy). Intriguingly, we found variation in cell-cycle activity, which is independent of cell progression through the cycle. The upregulation of G2/M checkpoint genes with downregulation of G2/M transition and cytokinesis is revealed in polyploid cells. A temporal increase in cell-cycle activity at transition from pluripotent to more differentiated state is found in human embryonic stem cells. The upregulation of unicellular interactome cluster in human cancers is shown in single cells with control for cell-cycle activity. The greater scatter around regression line in cancer cells suggests greater heterogeneity caused by deviation from a line of normal cells.

Resolution of Complex Issues in Genome Regulation and Cancer Requires Non-Linear and Network-Based Thermodynamics

The apparent lack of success in curing cancer that was evidenced in the last four decades of molecular medicine indicates the need for a global re-thinking both its nature and the biological approaches that we are taking in its solution. The reductionist, one gene/one protein method that has served us well until now, and that still dominates in biomedicine, requires complementation with a more systemic/holistic approach, to address the huge problem of cross-talk between more than 20,000 protein-coding genes, about 100,000 protein types, and the multiple layers of biological organization. In this perspective, the relationship between the chromatin network organization and gene expression regulation plays a fundamental role. The elucidation of such a relationship requires a non-linear thermodynamics approach to these biological systems. This change of perspective is a necessary step for developing successful 'tumour-reversion' therapeutic strategies.

Meta-Analysis of Cancer Triploidy: Rearrangements of Genome Complements in Male Human Tumors Are Characterized by XXY Karyotypes

Triploidy in cancer is associated with poor prognosis, but its origins remain unclear. Here, we attempted to differentiate between random chromosomal and whole-genome origins of cancer triploidy. In silico meta-analysis was performed on 15 male malignant and five benign tumor cohorts (2928 karyotypes) extracted from the Mitelman Database, comparing their ploidy and combinations of sex chromosomes. A distinct near-triploid fraction was observed in all malignant tumor types, and was especially high in seminoma. For all tumor types, X-chromosome doubling, predominantly observed as XXY, correlated strongly with the near-triploid state (r ≈ 0.9, p < 0.001), negatively correlated with near-diploidy, and did not correlate with near-tetraploidy. A smaller near-triploid component with a doubled X-chromosome was also present in three of the five benign tumor types, especially notable in colon adenoma. Principal component analysis revealed a non-random correlation structure shaping the X-chromosome disomy distribution across all tumor types. We suggest that doubling of the maternal genome followed by pedogamic fusion with a paternal genome (a possible mimic of the fertilization aberration, 69, XXY digyny) associated with meiotic reprogramming may be responsible for the observed rearrangements of genome complements leading to cancer triploidy. The relatively frequent loss of the Y-chromosome results as a secondary factor from chromosome instability.