Self-organization of whole-gene expression through coordinated chromatin structural transition

Impact of Gold Nanoparticles and Ionizing Radiation on Whole Chromatin Organization as Detected by Single-Molecule Localization Microscopy

In radiation tumor therapy, irradiation, on one hand, should cause cell death to the tumor. On the other hand, the surrounding non-tumor tissue should be maintained unaffected. Therefore, methods of local dose enhancements are highly interesting. Gold nanoparticles, which are preferentially uptaken by very-fast-proliferating tumor cells, may enhance damaging. However, the results in the literature obtained from cell culture and animal tissue experiments are very contradictory, i.e., only some experiments reveal increased cell killing but others do not. Thus, a better understanding of cellular mechanisms is required. Using the breast cancer cell model SkBr3, the effects of gold nanoparticles in combination with ionizing radiation on chromatin network organization were investigated by Single-Molecule Localization Microscopy (SMLM) and applications of mathematical topology calculations (e.g., Persistent Homology, Principal Component Analysis, etc.). The data reveal a dose and nanoparticle dependent re-organization of chromatin, although colony forming assays do not show a significant reduction of cell survival after the application of gold nanoparticles to the cells. In addition, the spatial organization of γH2AX clusters was elucidated, and characteristic changes were obtained depending on dose and gold nanoparticle application. The results indicate a complex response of ALU-related chromatin and heterochromatin organization correlating to ionizing radiation and gold nanoparticle incorporation. Such complex whole chromatin re-organization is usually associated with changes in genome function and supports the hypothesis that, with the application of gold nanoparticles, not only is DNA damage increasing but also the efficiency of DNA repair may be increased. The understanding of complex chromatin responses might help to improve the gold nanoparticle efficiency in radiation treatment.

Peculiar k-mer Spectra Are Correlated with 3D Contact Frequencies and Breakpoint Regions in the Human Genome

BACKGROUND

It is widely accepted that the 3D chromatin organization in human cell nuclei is not random and recent investigations point towards an interactive relation of epigenetic functioning and chromatin (re-)organization. Although chromatin organization seems to be the result of self-organization of the entirety of all molecules available in the cell nucleus, a general question remains open as to what extent chromatin organization might additionally be predetermined by the DNA sequence and, if so, if there are characteristic differences that distinguish typical regions involved in dysfunction-related aberrations from normal ones, since typical DNA breakpoint regions involved in disease-related chromosome aberrations are not randomly distributed along the DNA sequence.

METHODS

Highly conserved k-mer patterns in intronic and intergenic regions have been reported in eukaryotic genomes. In this article, we search and analyze regions deviating from average spectra (ReDFAS) of k-mer word frequencies in the human genome. This includes all assembled regions, e.g., telomeric, centromeric, genic as well as intergenic regions.

RESULTS

A positive correlation between k-mer spectra and 3D contact frequencies, obtained exemplarily from given Hi-C datasets, has been found indicating a relation of ReDFAS to chromatin organization and interactions. We also searched and found correlations of known functional annotations, e.g., genes correlating with ReDFAS. Selected regions known to contain typical breakpoints on chromosomes 9 and 5 that are involved in cancer-related chromosomal aberrations appear to be enriched in ReDFAS. Since transposable elements like ALUs are often assigned as major players in 3D genome organization, we also studied their impact on our examples but could not find a correlation between ALU regions and breakpoints comparable to ReDFAS.

CONCLUSIONS

Our findings might show that ReDFAS are associated with instable regions of the genome and regions with many chromatin contacts which is in line with current research indicating that chromatin loop anchor points lead to genomic instability.

From Cell States to Cell Fates: Control of Cell State Transitions

We examine the coordinated behavior of thousands of genes in cell fate transitions through genome expression as an integrated dynamical system using the concepts of self-organized criticality and coherent stochastic behavior. To quantify the effects of the collective behavior of genes, we adopted the flux balance approach and developed it in a new tool termed expression flux analysis (EFA). Here we describe this tool and demonstrate how its application to specific experimental genome-wide expression data provides new insights into the dynamics of the cell-fate transitions. Particularly, we show that in cell fate change, specific stochastic perturbations can spread over the entire system to guide distinct cell fate transitions through switching cyclic flux flow in the genome engine. Utilization of EFA enables us to elucidate a unified genomic mechanism for when and how cell-fate change occurs through critical transitions.

Synchronization between Attractors: Genomic Mechanism of Cell-Fate Change

Herein, we provide a brief overview of complex systems theory approaches to investigate the genomic mechanism of cell-fate changes. Cell trajectories across the epigenetic landscape, whether in development, environmental responses, or disease progression, are controlled by extensively coordinated genome-wide gene expression changes. The elucidation of the mechanisms underlying these coherent expression changes is of fundamental importance in cell biology and for paving the road to new therapeutic approaches. In previous studies, we pointed at dynamic criticality as a plausible characteristic of genome-wide transition dynamics guiding cell fate. Whole-genome expression develops an engine-like organization (genome engine) in order to establish an autonomous dynamical system, capable of both homeostasis and transition behaviors. A critical set of genes behaves as a critical point (CP) that serves as the organizing center of cell-fate change. When the system is pushed away from homeostasis, the state change that occurs at the CP makes local perturbation spread over the genome, demonstrating self-organized critical (SOC) control of genome expression. Oscillating-Mode genes (which normally keep genome expression on pace with microenvironment fluctuations), when in the presence of an effective perturbative stimulus, drive the dynamics of synchronization, and thus guide the cell-fate transition.

Self-Regulated Symmetry Breaking Model for Stem Cell Differentiation

In conventional disorder-order phase transitions, a system shifts from a highly symmetric state, where all states are equally accessible (disorder) to a less symmetric state with a limited number of available states (order). This transition may occur by varying a control parameter that represents the intrinsic noise of the system. It has been suggested that stem cell differentiation can be considered as a sequence of such symmetry-breaking events. Pluripotent stem cells, with their capacity to develop into any specialized cell type, are considered highly symmetric systems. In contrast, differentiated cells have lower symmetry, as they can only carry out a limited number of functions. For this hypothesis to be valid, differentiation should emerge collectively in stem cell populations. Additionally, such populations must have the ability to self-regulate intrinsic noise and navigate through a critical point where spontaneous symmetry breaking (differentiation) occurs. This study presents a mean-field model for stem cell populations that considers the interplay of cell-cell cooperativity, cell-to-cell variability, and finite-size effects. By introducing a feedback mechanism to control intrinsic noise, the model can self-tune through different bifurcation points, facilitating spontaneous symmetry breaking. Standard stability analysis showed that the system can potentially differentiate into several cell types mathematically expressed as stable nodes and limit cycles. The existence of a Hopf bifurcation in our model is discussed in light of stem cell differentiation.

Advanced image-free analysis of the nano-organization of chromatin and other biomolecules by Single Molecule Localization Microscopy (SMLM)

The cell as a system of many components, governed by the laws of physics and chemistry drives molecular functions having an impact on the spatial organization of these systems and vice versa. Since the relationship between structure and function is an almost universal rule not only in biology, appropriate methods are required to parameterize the relationship between the structure and function of biomolecules and their networks, the mechanisms of the processes in which they are involved, and the mechanisms of regulation of these processes. Single molecule localization microscopy (SMLM), which we focus on here, offers a significant advantage for the quantitative parametrization of molecular organization: it provides matrices of coordinates of fluorescently labeled biomolecules that can be directly subjected to advanced mathematical analytical procedures without the need for laborious and sometimes misleading image processing. Here, we propose mathematical tools for comprehensive quantitative computer data analysis of SMLM point patterns that include Ripley distance frequency analysis, persistent homology analysis, persistent 'imaging', principal component analysis and co-localization analysis. The application of these methods is explained using artificial datasets simulating different, potentially possible and interpretatively important situations. Illustrative analyses of real complex biological SMLM data are presented to emphasize the applicability of the proposed algorithms. This manuscript demonstrated the extraction of features and parameters quantifying the influence of chromatin (re)organization on genome function, offering a novel approach to study chromatin architecture at the nanoscale. However, the ability to adapt the proposed algorithms to analyze essentially any molecular organizations, e.g., membrane receptors or protein trafficking in the cytosol, offers broad flexibility of use.

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.

Selection and Evaluation of mRNA and miRNA Reference Genes for Expression Studies (qPCR) in Archived Formalin-Fixed and Paraffin-Embedded (FFPE) Colon Samples of DSS-Induced Colitis Mouse Model

The choice of appropriate reference genes is essential for correctly interpreting qPCR data and results. However, the majority of animal studies use a single reference gene without any prior evaluation. Therefore, many qPCR results from rodent studies can be misleading, affecting not only reproducibility but also translatability. In this study, the expression stability of reference genes for mRNA and miRNA in archived FFPE samples of 117 C57BL/6JOlaHsd mice (males and females) from 9 colitis experiments (dextran sulfate sodium; DSS) were evaluated and their expression analysis was performed. In addition, we investigated whether normalization reduced/neutralized the influence of inter/intra-experimental factors which we systematically included in the study. Two statistical algorithms (NormFinder and Bestkeeper) were used to determine the stability of reference genes. Multivariate analysis was made to evaluate the influence of normalization with different reference genes on target gene expression in regard to inter/intra-experimental factors. Results show that archived FFPE samples are a reliable source of RNA and imply that the FFPE procedure does not change the ranking of stability of reference genes obtained in fresh tissues. Multivariate analysis showed that the histological picture is an important factor affecting the expression levels of target genes.

Analysis of Dormancy-Associated Transcriptional Networks Reveals a Shared Quiescence Signature in Lung and Colorectal Cancer

Quiescent cancer cells (QCCs) are a common feature of solid tumors, representing a major obstacle to the long-term success of cancer therapies. We isolated QCCs ex vivo from non-small cell lung cancer (NSCLC) and colorectal cancer (CRC) xenografts with a label-retaining strategy and compared QCCs gene expression profiles to identify a shared "quiescence signature". Principal Component Analysis (PCA) revealed a specific component neatly discriminating quiescent and replicative phenotypes in NSCLC and CRC. The discriminating component showed significant overlapping, with 688 genes in common including ZEB2, a master regulator of stem cell plasticity and epithelial-to-mesenchymal transition (EMT). Gene set enrichment analysis showed that QCCs of both NSCLC and CRC had an increased expression of factors related to stemness/self renewal, EMT, TGF-β, morphogenesis, cell adhesion and chemotaxis, whereas proliferating cells overexpressed Myc targets and factors involved in RNA metabolism. Eventually, we analyzed in depth by means of a complex network approach, both the 'morphogenesis module' and the subset of differentially expressed genes shared by NCSLC and CRC. This allowed us to recognize different gene regulation network wiring for quiescent and proliferating cells and to underpin few genes central for network integration that may represent new therapeutic vulnerabilities. Altogether, our results highlight common regulatory pathways in QCCs of lung and colorectal tumors that may be the target of future therapeutic interventions.

A multivariate statistical test for differential expression analysis

Statistical tests of differential expression usually suffer from two problems. Firstly, their statistical power is often limited when applied to small and skewed data sets. Secondly, gene expression data are usually discretized by applying arbitrary criteria to limit the number of false positives. In this work, a new statistical test obtained from a convolution of multivariate hypergeometric distributions, the Hy-test, is proposed to address these issues. Hy-test has been carried out on transcriptomic data from breast and kidney cancer tissues, and it has been compared with other differential expression analysis methods. Hy-test allows implicit discretization of the expression profiles and is more selective in retrieving both differential expressed genes and terms of Gene Ontology. Hy-test can be adopted together with other tests to retrieve information that would remain hidden otherwise, e.g., terms of (1) cell cycle deregulation for breast cancer and (2) "programmed cell death" for kidney cancer.

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.

Networks and Islands of Genome Nano-architecture and Their Potential Relevance for Radiation Biology : (A Hypothesis and Experimental Verification Hints)

The cell nucleus is a complex biological system in which simultaneous reactions and functions take place to keep the cell as an individualized, specialized system running well. The cell nucleus contains chromatin packed in various degrees of density and separated in volumes of chromosome territories and subchromosomal domains. Between the chromatin, however, there is enough "free" space for floating RNA, proteins, enzymes, ATPs, ions, water molecules, etc. which are trafficking by super- and supra-diffusion to the interaction points where they are required. It seems that this trafficking works somehow automatically and drives the system perfectly. After exposure to ionizing radiation causing DNA damage from single base damage up to chromatin double-strand breaks, the whole system "cell nucleus" responds, and repair processes are starting to recover the fully functional and intact system. In molecular biology, many individual epigenetic pathways of DNA damage response or repair of single and double-strand breaks are described. How these responses are embedded into the response of the system as a whole is often out of the focus of consideration. In this article, we want to follow the hypothesis of chromatin architecture's impact on epigenetic pathways and vice versa. Based on the assumption that chromatin acts like an "aperiodic solid state within a limited volume," functionally determined networks and local topologies ("islands") can be defined that drive the appropriate repair process at a given damage site. Experimental results of investigations of the chromatin nano-architecture and DNA repair clusters obtained by means of single-molecule localization microscopy offer hints and perspectives that may contribute to verifying the hypothesis.

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.

The Eco-Evo Mandala: Simplifying Bacterioplankton Complexity into Ecohealth Signatures

The microbiome emits informative signals of biological organization and environmental pressure that aid ecosystem monitoring and prediction. Are the many signals reducible to a habitat-specific portfolio that characterizes ecosystem health? Does an optimally structured microbiome imply a resilient microbiome? To answer these questions, we applied our novel Eco-Evo Mandala to bacterioplankton data from four habitats within the Great Barrier Reef, to explore how patterns in community structure, function and genetics signal habitat-specific organization and departures from theoretical optimality. The Mandala revealed communities departing from optimality in habitat-specific ways, mostly along structural and functional traits related to bacterioplankton abundance and interaction distributions (reflected by ϵ and λ as power law and exponential distribution parameters), which are not linearly associated with each other. River and reef communities were similar in their relatively low abundance and interaction disorganization (low ϵ and λ) due to their protective structured habitats. On the contrary, lagoon and estuarine inshore reefs appeared the most disorganized due to the ocean temperature and biogeochemical stress. Phylogenetic distances (D) were minimally informative in characterizing bacterioplankton organization. However, dominant populations, such as Proteobacteria, Bacteroidetes, and Cyanobacteria, were largely responsible for community patterns, being generalists with a large functional gene repertoire (high D) that increases resilience. The relative balance of these populations was found to be habitat-specific and likely related to systemic environmental stress. The position on the Mandala along the three fundamental traits, as well as fluctuations in this ecological state, conveys information about the microbiome's health (and likely ecosystem health considering bacteria-based multitrophic dependencies) as divergence from the expected relative optimality. The Eco-Evo Mandala emphasizes how habitat and the microbiome's interaction network topology are first- and second-order factors for ecosystem health evaluation over taxonomic species richness. Unhealthy microbiome communities and unbalanced microbes are identified not by macroecological indicators but by mapping their impact on the collective proportion and distribution of interactions, which regulates the microbiome's ecosystem function.