Functional genomic fabrics are remodeled in a mouse model of Chagasic cardiomyopathy and restored following cell therapy

Low-Salt Diet Regulates the Metabolic and Signal Transduction Genomic Fabrics, and Remodels the Cardiac Normal and Chronic Pathological Pathways

Low-salt diet (LSD) is a constant recommendation to hypertensive patients, but the genomic mechanisms through which it improves cardiac pathophysiology are still not fully understood. Our publicly accessible transcriptomic dataset of the left ventricle myocardium of adult male mice subjected to prolonged LSD or normal diet was analyzed from the perspective of the Genomic Fabric Paradigm. We found that LSD shifted the metabolic priorities by increasing the transcription control for fatty acids biosynthesis while decreasing it for steroid hormone biosynthesis. Moreover, LSD remodeled pathways responsible for cardiac muscle contraction (CMC), chronic Chagas (CHA), diabetic (DIA), dilated (DIL), and hypertrophic (HCM) cardiomyopathies, and their interplays with the glycolysis/glucogenesis (GLY), oxidative phosphorylation (OXP), and adrenergic signaling in cardiomyocytes (ASC). For instance, the statistically (p < 0.05) significant coupling between GLY and ASC was reduced by LSD from 13.82% to 2.91% (i.e., -4.75×), and that of ASC with HCM from 10.50% to 2.83% (-3.71×). The substantial up-regulation of the CMC, ASC, and OXP genes, and the significant weakening of the synchronization of the expression of the HCM, CHA, DIA, and DIL genes within their respective fabrics justify the benefits of the LSD recommendation.

Theory and Applications of the (Cardio) Genomic Fabric Approach to Post-Ischemic and Hypoxia-Induced Heart Failure

The genomic fabric paradigm (GFP) characterizes the transcriptome topology by the transcripts’ abundances, the variability of the expression profile, and the inter-coordination of gene expressions in each pathophysiological condition. The expression variability analysis provides an indirect estimate of the cell capability to limit the stochastic fluctuations of the expression levels of key genes, while the expression coordination analysis determines the gene networks in functional pathways. This report illustrates the theoretical bases and the mathematical framework of the GFP with applications to our microarray data from mouse models of post ischemic, and constant and intermittent hypoxia-induced heart failures. GFP analyses revealed the myocardium priorities in keeping the expression of key genes within narrow intervals, determined the statistically significant gene interlinkages, and identified the gene master regulators in the mouse heart left ventricle under normal and ischemic conditions. We quantified the expression regulation, alteration of the expression control, and remodeling of the gene networks caused by the oxygen deprivation and determined the efficacy of the bone marrow mono-nuclear stem cell injections to restore the normal transcriptome. Through the comprehensive assessment of the transcriptome, GFP would pave the way towards the development of personalized gene therapy of cardiac diseases.

Cell-Based Therapies for Heart Failure

Heart failure has reached epidemic proportions with the advances in cardiovascular therapies for ischemic heart diseases and the progressive aging of the world population. Efficient pharmacological therapies are available for treating heart failure, but unfortunately, even with optimized therapy, prognosis is often poor. Their last therapeutic option is, therefore, a heart transplantation with limited organ supply and complications related to immunosuppression. In this setting, cell therapies have emerged as an alternative. Many clinical trials have now been performed using different cell types and injection routes. In this perspective, we will analyze the results of such trials and discuss future perspectives for cell therapies as an efficacious treatment of heart failure.

Retinal Genomic Fabric Remodeling after Optic Nerve Injury

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 three independent characteristics for the expression of each gene: level, variability, and correlation with each other gene. Thus, the 17,657 quantified genes in our study generated a total of 155,911,310 values to analyze. This represents 8830x more data per condition than a traditional transcriptomic analysis. ONC led to a 57% reduction in RGC numbers as detected by retrograde labeling with 1,1'-dioctadecyl-3,3,3,3'-tetramethylindocarbocyanine perchlorate (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 upregulated 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 pairwise 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 provided novel insights beyond the regulation of individual gene expression and disclosed 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.

Transcriptomic uniqueness and commonality of the ion channels and transporters in the four heart chambers

Myocardium transcriptomes of left and right atria and ventricles from four adult male C57Bl/6j mice were profiled with Agilent microarrays to identify the differences responsible for the distinct functional roles of the four heart chambers. Female mice were not investigated owing to their transcriptome dependence on the estrous cycle phase. Out of the quantified 16,886 unigenes, 15.76% on the left side and 16.5% on the right side exhibited differential expression between the atrium and the ventricle, while 5.8% of genes were differently expressed between the two atria and only 1.2% between the two ventricles. The study revealed also chamber differences in gene expression control and coordination. We analyzed ion channels and transporters, and genes within the cardiac muscle contraction, oxidative phosphorylation, glycolysis/gluconeogenesis, calcium and adrenergic signaling pathways. Interestingly, while expression of Ank2 oscillates in phase with all 27 quantified binding partners in the left ventricle, the percentage of in-phase oscillating partners of Ank2 is 15% and 37% in the left and right atria and 74% in the right ventricle. The analysis indicated high interventricular synchrony of the ion channels expressions and the substantially lower synchrony between the two atria and between the atrium and the ventricle from the same side.

Powerful quantifiers for cancer transcriptomics

Every day, investigators find a new link between a form of cancer and a particular alteration in the sequence or/and expression level of a key gene, awarding this gene the title of “biomarker”. The clinician may choose from numerous available panels to assess the type of cancer based on the mutation or expression regulation (“transcriptomic signature”) of “driver” genes. However, cancer is not a “one-gene show” and, together with the alleged biomarker, hundreds other genes are found as mutated or/and regulated in cancer samples. Regardless of the platform, a well-designed transcriptomic study produces three independent features for each gene: Average expression level, expression variability and coordination with expression of each other gene. While the average expression level is used in all studies to identify what genes were up-/down-regulated or turn on/off, the other two features are unfairly ignored. We use all three features to quantify the transcriptomic change during the progression of the disease and recovery in response to a treatment. Data from our published microarray experiments on cancer nodules and surrounding normal tissue from surgically removed tumors prove that the transcriptomic topologies are not only different in histopathologically distinct regions of a tumor but also dynamic and unique for each human being. We show also that the most influential genes in cancer nodules [the Gene Master Regulators (GMRs)] are significantly less influential in the normal tissue. As such, “smart” manipulation of the cancer GMRs expression may selectively kill cancer cells with little consequences on the normal ones. Therefore, we strongly recommend a really personalized approach of cancer medicine and present the experimental procedure and the mathematical algorithm to identify the most legitimate targets (GMRs) for gene therapy.

Biomarkers, Master Regulators and Genomic Fabric Remodeling in a Case of Papillary Thyroid Carcinoma

Publicly available (own) transcriptomic data have been analyzed to quantify the alteration in functional pathways in thyroid cancer, establish the gene hierarchy, identify potential gene targets and predict the effects of their manipulation. The expression data have been generated by profiling one case of papillary thyroid carcinoma (PTC) and genetically manipulated BCPAP (papillary) and 8505C (anaplastic) human thyroid cancer cell lines. The study used the genomic fabric paradigm that considers the transcriptome as a multi-dimensional mathematical object based on the three independent characteristics that can be derived for each gene from the expression data. We found remarkable remodeling of the thyroid hormone synthesis, cell cycle, oxidative phosphorylation and apoptosis pathways. Serine peptidase inhibitor, Kunitz type, 2 (SPINT2) was identified as the Gene Master Regulator of the investigated PTC. The substantial increase in the expression synergism of SPINT2 with apoptosis genes in the cancer nodule with respect to the surrounding normal tissue (NOR) suggests that SPINT2 experimental overexpression may force the PTC cells into apoptosis with a negligible effect on the NOR cells. The predictive value of the expression coordination for the expression regulation was validated with data from 8505C and BCPAP cell lines before and after lentiviral transfection with DDX19B.

Pulmonary Hypertension Remodels the Genomic Fabrics of Major Functional Pathways

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-week-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, as well as neointimal lesions and obliteration of the lumen in small arteries. We found that PH increased the expression of caveolin1, matrix metallopeptidase 2, and numerous inflammatory and cell proliferation genes. The cell cycle, vascular smooth muscle contraction, and oxidative phosphorylation pathways, as well as their interplay, were largely perturbed. Our results also suggest that the upregulated 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 the HM group, indicate metabolic shift toward aerobic glycolysis.

Coordinated Activity of Transcriptional Networks Responding to the Pattern of Action Potential Firing in Neurons

Transcriptional responses to the appropriate temporal pattern of action potential firing are essential for long-term adaption of neuronal properties to the functional activity of neural circuits and environmental experience. However, standard transcriptome analysis methods can be too limited in identifying critical aspects that coordinate temporal coding of action potential firing with transcriptome response. A Pearson correlation analysis was applied to determine how pairs of genes in the mouse dorsal root ganglion (DRG) neurons are coordinately expressed in response to stimulation producing the same number of action potentials by two different temporal patterns. Analysis of 4728 distinct gene-pairs related to calcium signaling, 435,711 pairs of transcription factors, 820 pairs of voltage-gated ion channels, and 86,862 pairs of calcium signaling genes with transcription factors indicated that genes become coordinately activated by distinct action potential firing patterns and this depends on the duration of stimulation. Moreover, a measure of expression variance revealed that the control of transcripts abundances is sensitive to the pattern of stimulation. Thus, action potentials impact intracellular signaling and the transcriptome in dynamic manner that not only alter gene expression levels significantly (as previously reported) but also affects the control of their expression fluctuations and profoundly remodel the transcriptional networks.

The Gene Master Regulators (GMR) Approach Provides Legitimate Targets for Personalized, Time-Sensitive Cancer Gene Therapy

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.

Estrogen Protects Neurotransmission Transcriptome During Status Epilepticus

Women with epilepsy commonly have premature onset of menopause. The decrease in estrogen levels is associated with increased occurrence of neurodegenerative processes and cognitive decline. Previously, we found that estradiol (E2) replacement in ovariectomized (OVX) female rats significantly reduced the seizure-related damage in the sensitive hilar region of hippocampal dentate gyrus (DG). However, the complex mechanisms by which E2 empowers the genomic fabrics of neurotransmission to resist damaging effects of status epilepticus (SE) are still unclear. We determined the protective effects of the estradiol replacement against kainic acid-induced SE-associated transcriptomic alterations in the DG of OVX rats. Without E2 replacement, SE altered expression of 44% of the DG genes. SE affected all major functional pathways, including apoptosis (61%), Alzheimer's disease (47%), cell cycle (59%), long-term potentiation (62%), and depression (55%), as well as synaptic vesicle cycle (62%), glutamatergic (53%), GABAergic (49%), cholinergic (52%), dopaminergic (55%), and serotonergic (49%) neurotransmission. However, in rats with E2 replacement the percentage of significantly affected genes after SE was reduced to the average 11% (from 8% for apoptosis to 32% for GABAergic synapse). Interestingly, while SE down-regulated most of the synaptic receptor genes in oil-injected females it had little effect on these receptors after E2-replacement. Our novel Pathway Protection analysis indicated that the E2-replacement prevented SE-related damage from 50% for GABA to 75% for dopaminergic transmission. The 15% synergistic expression between genes involved in estrogen signaling (ESG) and neurotransmission explains why low E2 levels result in down-regulation of neurotransmission. Interestingly, in animals with E2-replacement, SE switched 131 synergistically expressed ESG-neurotransmission gene pairs into antagonistically expressed gene pairs. Thus, the ESG pathway acts like a buffer against SE-induced alteration of neurotransmission that may contribute to the E2-mediated maintenance of brain function after the SE injury in postmenopausal women. We also show that the long-term potentiation is lost in OVX rats following SE but not in those with E2 replacement. The electrophysiological findings in OVX female rats with SE are corroborated by the high percentage of long-term potentiation regulated genes (62%) in oil-injected while only 13% of genes were regulated following SE in E2-replaced rats.

ACTH and PMX53 recover synaptic transcriptome alterations in a rat model of infantile spasms

We profiled the gene expression in the hypothalamic arcuate nuclei (ARC) of 20 male and 20 female rats to determine the infantile spasms (IS) related transcriptomic alteration of neurotransmission and recovery following two treatments. Rats were prenatally exposed to betamethasone or saline followed by repeated postnatal subjection to NMDA-triggered IS. Rats with spasms were treated with ACTH, PMX53 or saline. Since ACTH, the first line treatment for IS, has inconsistent efficacy and potential harsh side effects, PMX53, a potent complement C5ar1 antagonist, was suggested as a therapeutic alternative given its effects in other epilepsy models. Novel measures that consider all genes and are not affected by arbitrary cut-offs were used, in addition to standard statistical tests, to quantify regulation and recovery of glutamatergic, GABAergic, cholinergic, dopaminergic and serotonergic pathways. Although IS alters expression of ~30% of the ARC genes in both sexes the transcriptomic effects are 3× more severe in males than their female counterparts, as indicated by the Weighted Pathway Regulation measure. Both treatments significantly restored the ARC neurotransmission transcriptome to the non-IS condition with PMX53 performing slightly better, as measured by the Pathway Restoration Efficiency, suggesting these treatments may reduce autistic traits often associated with IS.

Regeneration of neurotransmission transcriptome in a model of epileptic encephalopathy after antiinflammatory treatment

Inflammation is an established etiopathogenesis factor of infantile spasms (IS), a therapy-resistant epileptic syndrome of infancy. We investigated the IS-associated transcriptomic alterations of neurotransmission in rat hypothalamic arcuate nucleus, how they are corrected by antiinflamatory treatments and whether there are sex differences. IS was triggered by repeated intraperitoneal administration of N-methyl-D-aspartic acid following anti-inflammatory treatment (adreno-cortico-tropic-hormone (ACTH) or PMX53) or normal saline vehicle to prenatally exposed to betamethasone young rats. We found that treatments with both ACTH and PMX53 resulted in substantial recovery of the genomic fabrics of all types of synaptic transmission altered by IS. While ACTH represents the first line of treatment for IS, the even higher efficiency of PMX53 (an antagonist of the complement C5a receptor) in restoring the normal transcriptome was not expected. In addition to the childhood epilepsy, the recovery of the neurotransmission genomic fabrics by PMX53 also gives hope for the autism spectrum disorders that share a high comorbidity with IS. Our results revealed significant sex dichotomy in both IS-associated transcriptomic alterations (males more affected) and in the efficiency of PMX53 anti-inflammatory treatment (better for males). Our data further suggest that anti-inflammatory treatments correcting alterations in the inflammatory transcriptome may become successful therapies for refractory epilepsies.