Self-organisation of early stress response in the biology of cancer.

The early stress response by AP-1 (FOS/JUN), supported by upregulation of c-Myc and involved in cell-fate changes and adaptation to hostile environments, is increased in cancer. The review shows the biphasic character of this response with negative feed-back typically lasting a few hours as a feature of the genome regulation by self-organising criticality. It involves  rapid splitting of the pericentromeric heterochromatin clusters, opening of the active chromatin, and a massive transcription acceleration wave. Phylostratigraphic analysis revealed that AP-1 genes evolved in the Cambrian explosion ~500 Mya integrating the protein interaction networks of reproduction including proto-placenta intertwined with cytokine and immunity pathways, sex determination, oocyte maturation, and embryonal stemness. The peak  of this response as part of accelerated cell senescence led by AP-1 and IL-1ß was found in breast cancer cell-line resistant to doxorubicin. Adaptability of aggressive cancer to treatments can be explained by emergent stress response evolutionarily protecting reproduction.

Special Issue "Advances in Genome Regulation in Cancer".

Cancer is globally increasing [...].

The Price of Human Evolution: Cancer-Testis Antigens, the Decline in Male Fertility and the Increase in Cancer.

The increasing frequency of general and particularly male cancer coupled with the reduction in male fertility seen worldwide motivated us to seek a potential evolutionary link between these two phenomena, concerning the reproductive transcriptional modules observed in cancer and the expression of cancer-testis antigens (CTA). The phylostratigraphy analysis of the human genome allowed us to link the early evolutionary origin of cancer via the reproductive life cycles of the unicellulars and early multicellulars, potentially driving soma-germ transition, female meiosis, and the parthenogenesis of polyploid giant cancer cells (PGCCs), with the expansion of the CTA multi-families, very late during their evolution. CTA adaptation was aided by retrovirus domestication in the unstable genomes of mammals, for protecting male fertility in stress conditions, particularly that of humans, as compensation for the energy consumption of a large complex brain which also exploited retrotransposition. We found that the early and late evolutionary branches of human cancer are united by the immunity-proto-placental network, which evolved in the Cambrian and shares stress regulators with the finely-tuned sex determination system. We further propose that social stress and endocrine disruption caused by environmental pollution with organic materials, which alter sex determination in male foetuses and further spermatogenesis in adults, bias the development of PGCC-parthenogenetic cancer by default.

Spatial-Temporal Genome Regulation in Stress-Response and Cell-Fate Change.

Complex functioning of the genome in the cell nucleus is controlled at different levels: (a) the DNA base sequence containing all relevant inherited information; (b) epigenetic pathways consisting of protein interactions and feedback loops; (c) the genome architecture and organization activating or suppressing genetic interactions between different parts of the genome. Most research so far has shed light on the puzzle pieces at these levels. This article, however, attempts an integrative approach to genome expression regulation incorporating these different layers. Under environmental stress or during cell development, differentiation towards specialized cell types, or to dysfunctional tumor, the cell nucleus seems to react as a whole through coordinated changes at all levels of control. This implies the need for a framework in which biological, chemical, and physical manifestations can serve as a basis for a coherent theory of gene self-organization. An international symposium held at the Biomedical Research and Study Center in Riga, Latvia, on 25 July 2022 addressed novel aspects of the abovementioned topic. The present article reviews the most recent results and conclusions of the state-of-the-art research in this multidisciplinary field of science, which were delivered and discussed by scholars at the Riga symposium.

The Role of Mitotic Slippage in Creating a "Female Pregnancy-like System" in a Single Polyploid Giant Cancer Cell.

In our recent work, we observed that triple-negative breast cancer MDA-MB-231 cells respond to doxorubicin (DOX) via "mitotic slippage" (MS), discarding cytosolic damaged DNA during the process that provides their resistance to this genotoxic treatment. We also noted two populations of polyploid giant cells: those budding surviving offspring, versus those reaching huge ploidy by repeated MS and persisting for several weeks. Their separate roles in the recovery from treatment remained unclear. The current study was devoted to characterising the origin and relationship of these two sub-populations in the context of MS. MS was hallmarked by the emergence of nuclear YAP1/OCT4A/MOS/EMI2-positivity featuring a soma-germ transition to the meiotic-metaphase-arrested "maternal germ cell". In silico, the link between modules identified in the inflammatory innate immune response to cytosolic DNA and the reproductive module of female pregnancy (upregulating placenta developmental genes) was observed in polyploid giant cells. Asymmetry of the two subnuclei types, one repairing DNA and releasing buds enriched by CDC42/ACTIN/TUBULIN and the other persisting and degrading DNA in a polyploid giant cell, was revealed. We propose that when arrested in MS, a "maternal cancer germ cell" may be parthenogenetically stimulated by the placental proto-oncogene parathyroid-hormone-like-hormone, increasing calcium, thus creating a "female pregnancy-like" system within a single polyploid giant cancer cell.

The Transcriptome and Proteome Networks of Malignant Tumours Reveal Atavistic Attractors of Polyploidy-Related Asexual Reproduction.

The expression of gametogenesis-related (GG) genes and proteins, as well as whole genome duplications (WGD), are the hallmarks of cancer related to poor prognosis. Currently, it is not clear if these hallmarks are random processes associated only with genome instability or are programmatically linked. Our goal was to elucidate this via a thorough bioinformatics analysis of 1474 GG genes in the context of WGD. We examined their association in protein-protein interaction and coexpression networks, and their phylostratigraphic profiles from publicly available patient tumour data. The results show that GG genes are upregulated in most WGD-enriched somatic cancers at the transcriptome level and reveal robust GG gene expression at the protein level, as well as the ability to associate into correlation networks and enrich the reproductive modules. GG gene phylostratigraphy displayed in WGD+ cancers an attractor of early eukaryotic origin for DNA recombination and meiosis, and one relative to oocyte maturation and embryogenesis from early multicellular organisms. The upregulation of cancer-testis genes emerging with mammalian placentation was also associated with WGD. In general, the results suggest the role of polyploidy for soma-germ transition accessing latent cancer attractors in the human genome network, which appear as pre-formed along the whole Evolution of Life.

A Unified Genomic Mechanism of Cell-Fate Change.

The purpose of our studies is to elucidate the nature of massive control of the whole genome expression with a particular emphasis on cell-fate change. The whole genome expression is coordinated by the emergence of a critical point (CP: a peculiar set of biphasic genes) with the genome acting as an integrated dynamical system. In response to stimuli, the genome expression self-organizes into local sub-, near-, and super-critical states, each exhibiting distinct collective behaviors with its center of mass acting as a local attractor, coexisting with the whole genome attractor (GA). The CP serves as the organizing center of cell-fate change, and its activation makes local perturbation to spread over the genome affecting GA. The activation of CP is in turn elicited by genes with elevated temporal variance (oscillating-mode genes), normally in charge to keep genome expression at pace with microenvironment fluctuations. When oscillation exceeds a given threshold, the CP synchronizes with the GA driving genome expression state transition. The expression synchronization wave invading the entire genome is fostered by the fusion-splitting dynamics of silencing pericentromere-associated heterochromatin domains and the consequent folding-unfolding transitions of transcribing euchromatin domains. The proposed mechanism is a unified step toward a time-evolutional transition theory of biological regulation.

Role of the Circadian Clock "Death-Loop" in the DNA Damage Response Underpinning Cancer Treatment Resistance.

Here, we review the role of the circadian clock (CC) in the resistance of cancer cells to genotoxic treatments in relation to whole-genome duplication (WGD) and telomere-length regulation. The CC drives the normal cell cycle, tissue differentiation, and reciprocally regulates telomere elongation. However, it is deregulated in embryonic stem cells (ESCs), the early embryo, and cancer. Here, we review the DNA damage response of cancer cells and a similar impact on the cell cycle to that found in ESCs­overcoming G1/S, adapting DNA damage checkpoints, tolerating DNA damage, coupling telomere erosion to accelerated cell senescence, and favouring transition by mitotic slippage into the ploidy cycle (reversible polyploidy). Polyploidy decelerates the CC. We report an intriguing positive correlation between cancer WGD and the deregulation of the CC assessed by bioinformatics on 11 primary cancer datasets (rho = 0.83; p < 0.01). As previously shown, the cancer cells undergoing mitotic slippage cast off telomere fragments with TERT, restore the telomeres by ALT-recombination, and return their depolyploidised offspring to telomerase-dependent regulation. By reversing this polyploidy and the CC "death loop", the mitotic cycle and Hayflick limit count are thus again renewed. Our review and proposed mechanism support a life-cycle concept of cancer and highlight the perspective of cancer treatment by differentiation.

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.

Heterochromatin Networks: Topology, Dynamics, and Function (a Working Hypothesis).

Open systems can only exist by self-organization as pulsing structures exchanging matter and energy with the outer world. This review is an attempt to reveal the organizational principles of the heterochromatin supra-intra-chromosomal network in terms of nonlinear thermodynamics. The accessibility of the linear information of the genetic code is regulated by constitutive heterochromatin (CHR) creating the positional information in a system of coordinates. These features include scale-free splitting-fusing of CHR with the boundary constraints of the nucleolus and nuclear envelope. The analysis of both the literature and our own data suggests a radial-concentric network as the main structural organization principle of CHR regulating transcriptional pulsing. The dynamic CHR network is likely created together with nucleolus-associated chromatin domains, while the alveoli of this network, including springy splicing speckles, are the pulsing transcription hubs. CHR contributes to this regulation due to the silencing position variegation effect, stickiness, and flexible rigidity determined by the positioning of nucleosomes. The whole system acts in concert with the elastic nuclear actomyosin network which also emerges by self-organization during the transcriptional pulsing process. We hypothesize that the the transcriptional pulsing, in turn, adjusts its frequency/amplitudes specified by topologically associating domains to the replication timing code that determines epigenetic differentiation memory.

Differentiating cancer cells reveal early large-scale genome regulation by pericentric domains.

Finding out how cells prepare for fate change during differentiation commitment was our task. To address whether the constitutive pericentromere-associated domains (PADs) may be involved, we used a model system with known transcriptome data, MCF-7 breast cancer cells treated with the ErbB3 ligand heregulin (HRG), which induces differentiation and is used in the therapy of cancer. PAD-repressive heterochromatin (H3K9me3), centromere-associated-protein-specific, and active euchromatin (H3K4me3) antibodies, real-time PCR, acridine orange DNA structural test (AOT), and microscopic image analysis were applied. We found a two-step DNA unfolding after 15-20 and 60 min of HRG treatment, respectively. This behavior was consistent with biphasic activation of the early response genes (c-fos - fosL1/myc) and the timing of two transcriptome avalanches reported in the literature. In control, the average number of PADs negatively correlated with their size by scale-free distribution, and centromere clustering in turn correlated with PAD size, both indicating that PADs may create and modulate a suprachromosomal network by fusing and splitting a constant proportion of the constitutive heterochromatin. By 15 min of HRG treatment, the bursting unraveling of PADs from the nucleolus boundary occurred, coinciding with the first step of H3K4me3 chromatin unfolding, confirmed by AOT. The second step after 60 min of HRG treatment was associated with transcription of long noncoding RNA from PADs and peaking of fosL1/c-myc response. We hypothesize that the bursting of PAD clusters under a critical silencing threshold pushes the first transcription avalanche, whereas the destruction of the PAD network enables genome rewiring needed for differentiation repatterning, mediated by early response bivalent genes.

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

"Mitotic Slippage" and Extranuclear DNA in Cancer Chemoresistance: A Focus on Telomeres.

Mitotic slippage (MS), the incomplete mitosis that results in a doubled genome in interphase, is a typical response of TP53-mutant tumors resistant to genotoxic therapy. These polyploidized cells display premature senescence and sort the damaged DNA into the cytoplasm. In this study, we explored MS in the MDA-MB-231 cell line treated with doxorubicin (DOX). We found selective release into the cytoplasm of telomere fragments enriched in telomerase reverse transcriptase (hTERT), telomere capping protein TRF2, and DNA double-strand breaks marked by γH2AX, in association with ubiquitin-binding protein SQSTM1/p62. This occurs along with the alternative lengthening of telomeres (ALT) and DNA repair by homologous recombination (HR) in the nuclear promyelocytic leukemia (PML) bodies. The cells in repeated MS cycles activate meiotic genes and display holocentric chromosomes characteristic for inverted meiosis (IM). These giant cells acquire an amoeboid phenotype and finally bud the depolyploidized progeny, restarting the mitotic cycling. We suggest the reversible conversion of the telomerase-driven telomere maintenance into ALT coupled with IM at the sub-telomere breakage sites introduced by meiotic nuclease SPO11. All three MS mechanisms converging at telomeres recapitulate the amoeba-like agamic life-cycle, decreasing the mutagenic load and enabling the recovery of recombined, reduced progeny for return into the mitotic cycle.

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.

When Three Isn't a Crowd: A Digyny Concept for Treatment-Resistant, Near-Triploid Human Cancers.

Near-triploid human tumors are frequently resistant to radio/chemotherapy through mechanisms that are unclear. We recently reported a tight association of male tumor triploidy with XXY karyotypes based on a meta-analysis of 15 tumor cohorts extracted from the Mitelman database. Here we provide a conceptual framework of the digyny-like origin of this karyotype based on the germline features of malignant tumors and adaptive capacity of digyny, which supports survival in adverse conditions. Studying how the recombinatorial reproduction via diploidy can be executed in primary cancer samples and HeLa cells after DNA damage, we report the first evidence that diploid and triploid cell sub-populations constitutively coexist and inter-change genomes via endoreduplicated polyploid cells generated through genotoxic challenge. We show that irradiated triploid HeLa cells can enter tripolar mitosis producing three diploid sub-subnuclei by segregation and pairwise fusions of whole genomes. Considering the upregulation of meiotic genes in tumors, we propose that the reconstructed diploid sub-cells can initiate pseudo-meiosis producing two "gametes" (diploid "maternal" and haploid "paternal") followed by digynic-like reconstitution of a triploid stemline that returns to mitotic cycling. This process ensures tumor survival and growth by (1) DNA repair and genetic variation, (2) protection against recessive lethal mutations using the third genome.

The Cancer Aneuploidy Paradox: In the Light of Evolution.

Aneuploidy should compromise cellular proliferation but paradoxically favours tumour progression and poor prognosis. Here, we consider this paradox in terms of our most recent observations of chemo/radio-resistant cells undergoing reversible polyploidy. The latter perform the segregation of two parental groups of end-to-end linked dyads by pseudo-mitosis creating tetraploid cells through a dysfunctional spindle. This is followed by autokaryogamy and a homologous pairing preceding a bi-looped endo-prophase. The associated RAD51 and DMC1/γ-H2AX double-strand break repair foci are tandemly situated on the AURKB/REC8/kinetochore doublets along replicated chromosome loops, indicative of recombination events. MOS-associated REC8-positive peri-nucleolar centromere cluster organises a monopolar spindle. The process is completed by reduction divisions (bi-polar or by radial cytotomy including pedogamic exchanges) and by the release of secondary cells and/or the formation of an embryoid. Together this process preserves genomic integrity and chromosome pairing, while tolerating aneuploidy by by-passing the mitotic spindle checkpoint. Concurrently, it reduces the chromosome number and facilitates recombination that decreases the mutation load of aneuploidy and lethality in the chemo-resistant tumour cells. This cancer life-cycle has parallels both within the cycling polyploidy of the asexual life cycles of ancient unicellular protists and cleavage embryos of early multicellulars, supporting the atavistic theory of cancer.

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.

DNA methylation of the Oct4A enhancers in embryonal carcinoma cells after etoposide treatment is associated with alternative splicing and altered pluripotency in reversibly senescent cells.

The epigenetic mechanisms underlying chemoresistance in cancer cells resulting from drug-induced reversible senescence are poorly understood. Chemoresistant ESC-like embryonal carcinoma PA1 cells treated with etoposide (ETO) were previously found to undergo prolonged G2 arrest with transient p53-dependent upregulation of opposing fate regulators, p21CIP1 (senescence) and OCT4A (self-renewal). Here we report on the analysis of the DNA methylation state of the distal enhancer (DE) and proximal enhancer (PE) of the Oct4A gene during this dual response. When compared to non-treated controls the methylation level increased from 1.3% to 12.5% and from 3% to 19.4%, in the DE and PE respectively. It included CpG and non-CpG methylation, which was not chaotic but presented two patterns in each enhancer. Discorrelating with methylation of enhancers, the transcription of Oct4A increased, however, a strong expression of the splicing form Oct4B was also induced, along with down-regulation of the Oct4A partners of in the pluripotency/self-renewal network Sox2 and Lin28. WB demonstrated disjoining of the OCT4A protein from the chromatin-bound fraction. In survival clones, methylation of the DE was considerably erased, while some remnant of methylation of the PE was still observed. The alternative splicing for Oct4B was reduced, Oct4A level insignificantly decreased, while the expression of Sox2 and Lin28 recovered, all three became proportionally above the control. These findings indicate the involvement of the transient patterned methylation of the Oct4A enhancers and alternative splicing in the adaptive regulation of cell fate choice during the p53-dependant dual state of reversible senescence in ESC-like cancer stem cells.

Differential staining of peripheral nuclear chromatin with Acridine orange implies an A-form epichromatin conformation of the DNA.

The chromatin observed by conventional electron microscopy under the nuclear envelope constitutes a single layer of dense 30-35 nm granules, while ∼30 nm fibrils laterally attached to them, form large patches of lamin-associated domains (LADs). This particular surface "epichromatin" can be discerned by specific (H2A+H2B+DNA) conformational antibody at the inner nuclear envelope and around mitotic chromosomes. In order to differentiate the DNA conformation of the peripheral chromatin we applied an Acridine orange (AO) DNA structural test involving RNAse treatment and the addition of AO after acid pre-treatment. MCF-7 cells treated in this way revealed yellow/red patches of LADs attached to a thin green nuclear rim and with mitotic chromosomes outlined in green, topologically corresponding to epichromatin epitope staining by immunofluorescence. Differentially from LADs, the epichromatin was unable to provide metachromatic staining by AO, unless thermally denatured at 94oC. DNA enrichment in GC stretches has been recently reported for immunoprecipitated ∼ 1Kb epichromatin domains. Together these data suggest that certain epichromatin segments assume the relatively hydrophobic DNA A-conformation at the nuclear envelope and surface of mitotic chromosomes, preventing AO side dimerisation.  We hypothesize that epichromatin domains form nucleosome superbeads. Hydrophobic interactions stack these superbeads and align them at the nuclear envelope, while repulsing the hydrophilic LADs. The hydrophobicity of epichromatin explains its location at the surface of mitotic chromosomes and its function in mediating chromosome attachment to the restituting nuclear envelope during telophase.