Many years and $$$ spent for research did not yet produced a universally effective gene therapy of prostate cancer (PCa). Our studies indicated that not only each human, but even each cancer nodule in the same tumor has unique and dynamic gene expression profile, control and coordination. The tumor heterogeneity of the transcriptome topology is a strong argument in favor of personalized gene therapies, tailored on patients’ primary tumor unique characteristics. Here, we propose a bioinformatics procedure by which to identify the Gene Master Regulators (GMR) of cancer cells from transcriptomic data and predict consequences of their experimental manipulation. The procedure, based on our Genomic Fabric Paradigm (GFP), can determine the most important gene in each cancer nodule whose controlled alteration would selectively kill the cancer cells. In this report, the method is applied to our microarray data on two men PCs (each with three distinct cancer nodules) and two standard human PCa cell lines (DU145 and LNCaP). We expect the industry to produce ready-to-use CRISPR constructs for all genes allowing the clinical oncologist to prescribe the adequate personalized CRISPR cocktail for his/her patient. Thus, the GMR approach may provide the most effective, yet affordable solution in fighting cancer.

Prostate cancer is a leading cause of death among men but its genomic characterization and best therapeutic strategy are still under debate. The Genomic Fabric Paradigm (GFP) considers the transcriptome as a multi-dimensional mathematical object subjected to a dynamic set of expression correlations among the genes. Here, GFP is applied to gene expression profiles of three (one primary, and two secondary) cancer nodules and the surrounding normal tissue from a surgically removed prostate tumor. GFP was used to determine the regulation and rewiring of the P53 signaling, apoptosis, prostate cancer and several other pathways involved in survival and proliferation of the cancer cells. Genes responsible for the block of differentiation, evading apoptosis, immortality, insensitivity to anti-growth signals, proliferation, resistance to chemotherapy and sustained angiogenesis were found as differently regulated in the three cancer nodules with respect to the normal tissue. The analysis indicates that even histo-pathologically equally graded cancer nodules from the same tumor have substantially different transcriptomic organizations, raising legitimate questions about the validity of meta-analyses comparing large populations of healthy and cancer humans. The study suggests that GFP may provide a personalized alternative to the biomarkers’ approach of cancer genomics.

Myocardium transcriptomes of mouse left and right atria and ventricles were profiled separately to identify the differences that might be responsible for the distinct functional roles of the four heart chambers. In total, 16,886 distinct unigenes have been quantified in all 16 samples collected from four adult male mice from the same litter. 15.76% of the quantified genes on the left and 16.5% on the right exhibited differential expression between the corresponding atrium and ventricle of the same side, while 5.8% in atria and 1.2% in ventricles were differently expressed between the left and the right. Beyond the differentially expressed genes, the study revealed distinct expression control and coordination of ion channels and genes within the cardiac muscle contraction, oxidative phosphorylation, glycolysis/glucogenesis, calcium and adrenergic signaling pathways. Interestingly, while expression of Ank2 (encoding ankyrin-B) oscillates in phase with all its binding partners in the left ventricle, the percentage of synergistically expressed partners of Ank2 is 15% and 37% in the left and right atria and 74% in the right ventricle. The analysis revealed also the high interventricular synchrony of the expression of ion channels.

The dynamic and never exactly repeatable tumor transcriptomic profile of people affected by the same form of cancer requires a personalized and time-sensitive approach of the gene therapy. The Gene Master Regulators (GMRs) were defined as genes whose highly controlled expression by the homeostatic mechanisms commands the cell phenotype by modulating major functional pathways through expression correlation with their genes. The Gene Commanding Height (GCH), a measure that combines the expression control and expression correlation with all other genes, is used to establish the gene hierarchy in each cell phenotype. We developed the experimental protocol, the mathematical algorithm and the computer software to identify the GMRs from transcriptomic data in surgically removed tumors, biopsies or blood from cancer patients. The GMR approach is illustrated with applications to our microarray data on human kidney, thyroid and prostate cancer samples, and on thyroid, prostate and blood cancer cell lines. We proved experimentally that each patient has his/her own GMRs, that cancer nuclei and surrounding normal tissue are governed by different GMRs, and that manipulating the expression has larger consequences for genes with higher GCH. Therefore, we launch the hypothesis that silencing the GMR may selectively kill the cancer cells from a tissue.

Keywords: Adrenergic signaling in cardiomyocytes; cardiac muscle contraction; Chagas disease; diabetic cardiomyopathy; dilated cardiomyopathy; glycerolipid metabolism; glycolisis/glucogenesis; hypertrophic cardiomyopathy; purine metabolism; steroid hormone biosynthesis.

Clear cell renal cell carcinoma (ccRCC) is the most frequent form of kidney cancer. Metastatic stages of ccRCC reduce the five-year survival rate to 15%. In this report we analyze the ccRCC-induced remodeling of the five KEGG-constructed excretory functional pathways in a surgically removed right kidney and its metastasis in the chest wall from the perspective of the Genomic Fabric Paradigm (GFP). GFP characterizes every single gene in each region by the independent variables: average expression level (AVE), relative expression variability (REV) and expression correlation (COR) with each other gene. While traditional approach is limited to only AVE analysis, the novel REV analysis indicates genes whose correct expression level is critical for cell survival and proliferation. COR analysis determines the real gene networks responsible for functional pathways. The analyses covered the pathways: aldosterone-regulated sodium reabsorption, collecting duct acid secretion, endocrine and other factor-regulated sodium reabsorption, proximal tubule bicarbonate reclamation, and vasopressin-regulated water reabsorption. The present study confirms the conclusion of our previously published articles on prostate and kidney cancers that even equally graded cancer nodules from the same tumor have different transcriptomic topologies. Therefore, the personalization of the anti-cancer therapy should go beyond the individual to his/her major cancer nodules.

We profiled the transcriptomes of primary mouse cortical astrocytes cultured alone or co-cultured with immortalized precursor oligodendrocytes. The experimental set-up (insert systems) prevented formation of gap junction channels but allowed free exchange of the two culture media. The study complements our previously published reports that the genomic fabrics of major functional pathways in oligodendrocytes are substantially remodeled by the proximity of non-touching astrocytes. Here, we present new analysis indicating that the transcriptomic landscape of astrocytes likewise changes significantly in the proximity of non-touching oligodendrocytes. The research was stimulated by the reported transcriptomic similarity between the brains of Cx43KO and Cx32KO mice, both substantially different from that of the Cx36KO mice. Since the three connexins are expressed in different cell types (Cx43 in astrocytes, Cx32 in oligodendrocytes and Cx36 in neurons), altogether these findings support the idea of a “panglial transcriptomic syncytium” in the mouse brain. Going further, our results suggest that integration in a heterocellular tissue modulates not only the expression profile but also the expression control and networking of the genes in each cell phenotype.

Pulmonary hypertension (PH) is a serious disorder with high morbidity and mortality rate. We analyzed the right ventricular systolic pressure (RVSP), right ventricular hypertrophy (RVH), lung histology and transcriptomes of six weeks old male rats with PH induced by: 1) hypoxia (HO), 2) administration of monocrotaline (CM) or 3) administration of monocrotaline and exposure to hypoxia (HM). The results in PH rats were compared to those in control rats (CO). After four weeks exposure, increased RVSP and RVH, pulmonary arterial wall thickening, and alteration of the lung transcriptome were observed in all PH groups. The HM group exhibited the largest alterations and also neointimal lesions and obliteration of lumen in small arteries. We found that the PH increased the expression of caveolin1, matrix metallopeptidase 2 and numerous inflammatory and cell proliferation genes. The cell-cycle, vascular smooth muscle contraction and the oxidative phosphorylation pathways, as well as their interplay were largely perturbed. Our results also suggest that the up-regulated Rhoa (ras homolog family member A) mediates its action through expression coordination with several ATPases. The upregulation of antioxidant genes and the extensive mitochondrial damage observed especially in HM group, indicate metabolic shift towards aerobic glycolysis.

The high morbidity and mortality rate of the pulmonary arterial hypertension (PAH) is partially explained by the metabolic deregulation associated with the disease. The present study identified statistically significant increase of the glucose transporter solute carrier family 2 (Slc2a1), beta nerve growth factor (Ngf), and nuclear factor erythroid derived 2, like 2 (Nfe2l2) on three standard PAH rat models and their healthy counterpart. PAH was induced by subjecting the animals to hypoxia (HO group) or by injecting them with monocrotaline in either normal (CM) or hypoxic (HM) atmospheric conditions. The Western blot and double immunofluorescent experiments were complemented with the gene expression profiling of the animal lungs analyzed from the perspective of the Genomic Fabric Paradigm. We found substantial remodeling of the genomic fabrics of the citrate cycle, pyruvate metabolism, glycolysis/gluconeogenesis and fructose and manose pathways. According to the novel transcriptomic distance criterion, the most comprehensive measure of the transcriptomic alteration, glycolysis/gluconeogenesis was the most affected functional pathway in all three PAH models. PAH decoupled the coordinated expression of many metabolic genes and replaced Pmm2 with Pmm1 in the center of the fructose and mannose metabolism. Our data show that metabolic dysregulation is a major pathogenic factor of the PAH.

Published transcriptomic data from surgically removed metastatic clear cell renal cell carcinoma (ccRCC) samples were re-analyzed from the Genomic Fabric Perspective that considers the transcriptome a multi-dimensional mathematical object, constrained by a dynamic set of expression correlations among genes. Every gene in the chest wall metastasis (MET), two primary tumors (PTA, PTB) and the surrounding normal tissue (NOR) of the right kidney was characterized by three independent measures: average expression level (AVE), relative expression variation (REV) and expression correlation (COR) with each other gene. AVE was used to determine the regulation of the genomic fabrics of ccRCC, apoptosis, chemokine and VEGF signaling pathways. REV quantified the alteration of the transcripts’ abundances control, while COR determined the remodeling of the transcriptomic networks of chemokine signaling and oxidative phosphorylation genes. The gene hierarchy was established in based on Gene Commanding Height and the Gene Master Regulators (GMR) TASOR (PTA), FAM27C (PTB) and ALG13 (MET) and DAPK3 (NOR) were identified in each profiled region. We predict that TASOR overexpression would block transcription in PTA but not in PTB, while slightly stimulating it in NOR. Silencing of ALG3 would slow-down the cell-cycle in all three cancer regions with practically no effect in NOR.

The genetic causes of the differentiated, highly treatable and mostly non-fatal papillary thyroid cancer (PTC) are not yet fully known. The widely accepted PTC etiology blames altered sequence or/and expression level of certain biomarker genes. However, our previous studies found that the PTC biomarkers played significantly lower roles than the personalized PTC Gene Master Regulators. Our publicly accessible gene expression data from the cancer nodule and the surrounding normal tissue of surgically removed PTC tumor were re-analyzed to determine the transcriptomic effects on the genomic fabrics responsible for the DNA replication, repair and transcription. Tumor results were compared to the gene expression profiles of the papillary (BCPAP) and the anaplastic (8505C) human thyroid cancer cell lines. The analyses were carried out from the Genomic Fabric Paradigm (GFP) perspective that provides the most theoretically possible comprehensive characterization of the transcriptome and its alterations in disease. The study indicated that, in addition to regulating numerous genes, changes in the homeostatic control of transcript abundances and remodeling of the gene networks had major contributions to the PTC occurrence, resilience and proliferation. Therefore, the molecular mechanisms responsible for gene expression control and inter-coordination should also be considered when designing personalized anti-cancer gene therapies.

Glaucoma is a multifactorial neurodegenerative disease, characterized by degeneration of the retinal ganglion cells (RGCs). There has been little progress in developing efficient strategies for neuroprotection in glaucoma. We profiled the retina transcriptome of Lister Hooded rats at 2 weeks after optic nerve crush (ONC) and analyzed the data from the Genomic Fabric Paradigm (GFP) to bring additional insights into the molecular mechanisms of the retinal remodeling after induction of RGC degeneration. GFP considers for the expression of each gene 3 independent characteristics: level, variability and correlation with each other gene. Thus, the 17,657 quantified genes our study generated a total of 155,911,310 values to analyze. This represents 8,830x more data per condition than a traditional transcriptomic analysis. ONC led to a 57% reduction in RGC numbers as detected by retrograde labeling with DiI. We observed a higher Relative Expression Variability after ONC. Gene expression stability was used as a measure of transcription control and disclosed a robust reduction in the number of very stably expressed genes. Predicted Protein-Protein interaction (PPI) analysis with STRING revealed axon and neuron projection as mostly decreased processes, consistent with RGC degeneration. Conversely, immune response PPIs were found among up-regulated genes. Enrichment analysis showed that Complement Cascade and Notch Signaling Pathway, as well as Oxidative Stress and Kit Receptor Pathway were affected after ONC. To expand our studies of altered molecular pathways, we examined the pair-wise coordination of gene expressions within each pathway and within the entire transcriptome using Pearson correlations. ONC increased the number of synergistically coordinated pairs of genes and the number of similar profiles mainly in Complement Cascade and Notch Signaling Pathway. This deep bioinformatic study provides novel insights beyond the regulation of individual gene expression and discloses changes in the control of expression of Complement Cascade and Notch Signaling functional pathways that may be relevant for both RGC degeneration and remodeling of the retinal tissue after ONC.