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.

Genomic Fabrics of the Excretory System's Functional Pathways Remodeled in Clear Cell Renal Cell Carcinoma.

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). The GFP characterizes every single gene in each region by these independent variables: the average expression level (AVE), relative expression variability (REV), and expression correlation (COR) with each other gene. While the traditional approach is limited to only AVE analysis, the novel REV analysis identifies the genes whose correct expression level is critical for cell survival and proliferation. The COR analysis determines the real gene networks responsible for functional pathways. The analyses covered the pathways for 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 anti-cancer therapy should go beyond the individual, to his/her major cancer nodules.

Metabolic Deregulation in Pulmonary Hypertension.

The high morbidity and mortality rate of pulmonary arterial hypertension (PAH) is partially explained by metabolic deregulation. The present study complements our previous publication in "Genes" by identifying significant increases of the glucose transporter solute carrier family 2 (Slc2a1), beta nerve growth factor (Ngf), and nuclear factor erythroid-derived 2-like 2 (Nfe2l2) in three standard PAH rat models. PAH was induced by subjecting the animals to hypoxia (HO), or by injecting with monocrotaline in either normal (CM) or hypoxic (HM) atmospheric conditions. The Western blot and double immunofluorescent experiments were complemented with novel analyses of previously published transcriptomic datasets of the animal lungs from the perspective of the Genomic Fabric Paradigm. We found substantial remodeling of the citrate cycle, pyruvate metabolism, glycolysis/gluconeogenesis, and fructose and mannose pathways. According to the transcriptomic distance, glycolysis/gluconeogenesis was the most affected functional pathway in all three PAH models. PAH decoupled the coordinated expression of many metabolic genes, and replaced phosphomannomutase 2 (Pmm2) with phosphomannomutase 1 (Pmm1) in the center of the fructose and mannose metabolism. We also found significant regulation of key genes involved in PAH channelopathies. In conclusion, our data show that metabolic dysregulation is a major PAH pathogenic factor.

Personalized 3-Gene Panel for Prostate Cancer Target Therapy.

Many years and billions spent for research did not yet produce an effective answer to prostate cancer (PCa). Not only each human, but even each cancer nodule in the same tumor, has unique transcriptome topology. The differences go beyond the expression level to the expression control and networking of individual genes. The unrepeatable heterogeneous transcriptomic organization among men makes the quest for universal biomarkers and "fit-for-all" treatments unrealistic. We present a bioinformatics procedure to identify each patient's unique triplet of PCa Gene Master Regulators (GMRs) and predict consequences of their experimental manipulation. The procedure is based on the Genomic Fabric Paradigm (GFP), which characterizes each individual gene by the independent expression level, expression variability and expression coordination with each other gene. GFP can identify the GMRs whose controlled alteration would selectively kill the cancer cells with little consequence on the normal tissue. The method was applied to microarray data on surgically removed prostates from two men with metastatic PCas (each with three distinct cancer nodules), and DU145 and LNCaP PCa cell lines. The applications verified that each PCa case is unique and predicted the consequences of the GMRs' manipulation. The predictions are theoretical and need further experimental validation.

Synaptonuclear messenger PRR7 inhibits c-Jun ubiquitination and regulates NMDA-mediated excitotoxicity.

Elevated c-Jun levels result in apoptosis and are evident in neurodegenerative disorders such as Alzheimer's disease and dementia and after global cerebral insults including stroke and epilepsy. NMDA receptor (NMDAR) antagonists block c-Jun upregulation and prevent neuronal cell death following excitotoxic insults. However, the molecular mechanisms regulating c-Jun abundance in neurons are poorly understood. Here, we show that the synaptic component Proline rich 7 (PRR7) accumulates in the nucleus of hippocampal neurons following NMDAR activity. We find that PRR7 inhibits the ubiquitination of c-Jun by E3 ligase SCF(FBW) (7) (FBW7), increases c-Jun-dependent transcriptional activity, and promotes neuronal death. Microarray assays show that PRR7 abundance is directly correlated with transcripts associated with cellular viability. Moreover, PRR7 knockdown attenuates NMDAR-mediated excitotoxicity in neuronal cultures in a c-Jun-dependent manner. Our results show that PRR7 links NMDAR activity to c-Jun function and provide new insights into the molecular processes that underlie NMDAR-dependent excitotoxicity.

Papillary Thyroid Cancer Remodels the Genetic Information Processing Pathways.

The genetic causes of the differentiated, highly treatable, and mostly non-fatal papillary thyroid cancer (PTC) are not yet fully understood. The mostly accepted PTC etiology blames the altered sequence or/and expression level of certain biomarker genes. However, tumor heterogeneity and the patient's unique set of favoring factors question the fit-for-all gene biomarkers. Publicly accessible gene expression profiles of the cancer nodule and the surrounding normal tissue from a surgically removed PTC tumor were re-analyzed to determine the cancer-induced alterations of the genomic fabrics responsible for major functional pathways. Tumor data were compared with those of standard papillary and anaplastic thyroid cancer cell lines. We found that PTC regulated numerous genes associated with DNA replication, repair, and transcription. Results further indicated that changes of the gene networking in functional pathways and the homeostatic control of transcript abundances also had major contributions to the PTC phenotype occurrence. The purpose to proliferate and invade the entire gland may explain the substantial transcriptomic differences we detected between the cells of the cancer nodule and those spread in homo-cellular cultures (where they need only to survive). In conclusion, the PTC etiology should include the complex molecular mechanisms involved in the remodeling of the genetic information processing pathways.

Plasticity of the myelination genomic fabric.

This study aimed to quantify the influence of the astrocyte proximity on myelination genomic fabric (MYE) of oligodendrocytes, defined as the most interconnected and stably expressed gene web responsible for myelination. Such quantitation is important to evaluate whether astrocyte signaling may contribute to demyelination when impaired and remyelination when properly restored. For this, we compared changes in the gene expression profiles of immortalized precursor oligodendrocytes (Oli-neu), stimulated to differentiate by the proximity of nontouching astrocytes or treatment with db-cAMP. In a previous paper, we reported that the astrocyte proximity upregulated or turned-on a large number of myelination genes and substantially enriched the Ca(2+)-signaling and cytokine receptor regulatory networks of MYE in Oli-neu cells. Here, we introduce the "transcriptomic distance" to evaluate fabric remodeling and "pair-wise relevance" to identify the most influential gene pairs. Together with the prominence gene analysis used to select and rank the fabric genes, these novel analytical tools provide a comprehensively quantitative view of the physio/pathological transformations of the transcriptomic programs of myelinating cells. Applied to our data, the analyses revealed not only that the astrocyte neighborhood is a substantially more powerful regulator of myelination than the differentiating treatment but also the molecular mechanisms of the two differentiating paradigms are different. By inducing a profound remodeling of MYE and regulatory transcriptomic networks, the astrocyte-oligodendrocyte intercommunication may be considered as a major player in both pathophysiology and therapy of neurodegenerative diseases related to myelination.

Gene expression changes associated with myocarditis and fibrosis in hearts of mice with chronic chagasic cardiomyopathy.

Chronic chagasic cardiomyopathy is a leading cause of heart failure in Latin American countries. About 30% of Trypanosoma cruzi-infected individuals develop this severe symptomatic form of the disease, characterized by intense inflammatory response accompanied by fibrosis in the heart. We performed an extensive microarray analysis of hearts from a mouse model of this disease and identified significant alterations in expression of approximately 12% of the sampled genes. Extensive up-regulations were associated with immune-inflammatory responses (chemokines, adhesion molecules, cathepsins, and major histocompatibility complex molecules) and fibrosis (extracellular matrix components, lysyl oxidase, and tissue inhibitor of metalloproteinase 1). Our results indicate potentially relevant factors involved in the pathogenesis of the disease that may provide new therapeutic targets in chronic Chagas disease.

A Personalized Genomics Approach of the Prostate Cancer.

Decades of research identified genomic similarities among prostate cancer patients and proposed general solutions for diagnostic and treatments. However, each human is a dynamic unique with never repeatable transcriptomic topology and no gene therapy is good for everybody. Therefore, we propose the Genomic Fabric Paradigm (GFP) as a personalized alternative to the biomarkers approach. Here, GFP is applied to three (one primary-"A", and two secondary-"B" & "C") cancer nodules and the surrounding normal tissue ("N") from a surgically removed prostate tumor. GFP proved for the first time that, in addition to the expression levels, cancer alters also the cellular control of the gene expression fluctuations and remodels their networking. Substantial differences among the profiled regions were found in the pathways of P53-signaling, apoptosis, prostate cancer, block of differentiation, evading apoptosis, immortality, insensitivity to anti-growth signals, proliferation, resistance to chemotherapy, and sustained angiogenesis. ENTPD2, AP5M1 BAIAP2L1, and TOR1A were identified as the master regulators of the "A", "B", "C", and "N" regions, and potential consequences of ENTPD2 manipulation were analyzed. The study shows that GFP can fully characterize the transcriptomic complexity of a heterogeneous prostate tumor and identify the most influential genes in each cancer nodule.

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.

TWEAKing the Hippocampus: The Effects of TWEAK on the Genomic Fabric of the Hippocampus in a Neuropsychiatric Lupus Mouse Model.

Neuropsychiatric manifestations of systemic lupus erythematosus (SLE), specifically cognitive dysfunction and mood disorders, are widely prevalent in SLE patients, and yet poorly understood. TNF-like weak inducer of apoptosis (TWEAK) has previously been implicated in the pathogenesis of neuropsychiatric lupus (NPSLE), and we have recently shown its effects on the transcriptome of the cortex of the lupus-prone mice model MRL/lpr. As the hippocampus is thought to be an important focus of NPSLE processes, we explored the TWEAK-induced transcriptional changes that occur in the hippocampus, and isolated several genes (Dnajc28, Syne2, transthyretin) and pathways (PI3K-AKT, as well as chemokine-signaling and neurotransmission pathways) that are most differentially affected by TWEAK activation. While the functional roles of these genes and pathways within NPSLE need to be further investigated, an interesting link between neuroinflammation and neurodegeneration appears to emerge, which may prove to be a promising novel direction in NPSLE research.

Remodeling of Neurotransmission, Chemokine, and PI3K-AKT Signaling Genomic Fabrics in Neuropsychiatric Systemic Lupus Erythematosus.

Cognitive dysfunction and mood changes are prevalent and especially taxing issues for patients with systemic lupus erythematosus (SLE). Tumor necrosis factor (TNF)-like weak inducer of apoptosis (TWEAK) and its cognate receptor Fn14 have been shown to play an important role in neurocognitive dysfunction in murine lupus. We profiled and compared gene expression in the cortices of MRL/+, MRL/lpr (that manifest lupus-like phenotype) and MRL/lpr-Fn14 knockout (Fn14ko) adult female mice to determine the transcriptomic impact of TWEAK/Fn14 on cortical gene expression in lupus. We found that the TWEAK/Fn14 pathway strongly affects the expression level, variability and coordination of the genomic fabrics responsible for neurotransmission and chemokine signaling. Dysregulation of the Phosphoinositide 3-kinase (PI3K)-AKT pathway in the MRL/lpr lupus strain compared with the MRL/+ control and Fn14ko mice was particularly prominent and, therefore, promising as a potential therapeutic target, although the complexity of the transcriptomic fabric highlights important considerations in in vivo experimental models.

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.

Heart rhythm genomic fabric in hypoxia.

The molecular mechanisms by which chronic hypoxia, whether constant (CCH) or intermittent (CIH), alters the heart rhythm are still under debate. Expression level, control, maturational profile and intercoordination of 54 genes encoding heart rhythm determinants (HRDs) were analyzed in 36 mice subjected for 1, 2 or 4 weeks of their early life to normal atmospheric conditions or to CCH or CIH. Our analysis revealed a complex network of genes encoding various heart rate, inotropy and development controllers, receptors, ion channels and transporters, ankyrins, epigenetic modulators and intercalated disc components (adherens, cadherins, catenins, desmosomal, gap and tight junction proteins). The network is remodeled during maturation and substantially and differently altered by CIH and CCH. Gene Prominence Analysis that ranks the genes according to their expression stability and networking within functional gene webs, confirmed the HRD status of certain epigenetic modulators and components of the intercalated discs not yet associated with arrhythmia.

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.

Cellular Environment Remodels the Genomic Fabrics of Functional Pathways in Astrocytes.

We profiled the transcriptomes of primary mouse cortical astrocytes cultured alone or co-cultured with immortalized precursor oligodendrocytes (Oli-neu cells). Filters between the cell types prevented formation of hetero-cellular gap junction channels but allowed for free exchange of the two culture media. We previously reported that major functional pathways in the Oli-neu cells are remodeled by the proximity of non-touching astrocytes and that astrocytes and oligodendrocytes form a panglial transcriptomic syncytium in the brain. Here, we present evidence that the astrocyte transcriptome likewise changes significantly in the proximity of non-touching Oli-neu cells. Our results indicate that the cellular environment strongly modulates the transcriptome of each cell type and that integration in a heterocellular tissue changes not only the expression profile but also the expression control and networking of the genes in each cell phenotype. The significant decrease of the overall transcription control suggests that in the co-culture astrocytes are closer to their normal conditions from the brain. The Oli-neu secretome regulates astrocyte genes known to modulate neuronal synaptic transmission and remodels calcium, chemokine, NOD-like receptor, PI3K-Akt, and thyroid hormone signaling, as well as actin-cytoskeleton, autophagy, cell cycle, and circadian rhythm pathways. Moreover, the co-culture significantly changes the gene hierarchy in the astrocytes.

Genomic Fabric Remodeling in Metastatic Clear Cell Renal Cell Carcinoma (ccRCC): A New Paradigm and Proposal for a Personalized Gene Therapy Approach.

Published transcriptomic data from surgically removed metastatic clear cell renal cell carcinoma samples were analyzed from the genomic fabric paradigm (GFP) perspective to identify the best targets for gene therapy. GFP considers the transcriptome as a multi-dimensional mathematical object constrained by a dynamic set of expression controls and correlations among genes. Every gene in the chest wall metastasis, two distinct cancer nodules, and the surrounding normal tissue of the right kidney was characterized by three independent measures: average expression level, relative expression variation, and expression correlation with each other gene. The analyses determined the cancer-induced regulation, control, and remodeling of the chemokine and vascular endothelial growth factor (VEGF) signaling, apoptosis, basal transcription factors, cell cycle, oxidative phosphorylation, renal cell carcinoma, and RNA polymerase pathways. Interestingly, the three cancer regions exhibited different transcriptomic organization, suggesting that the gene therapy should not be personalized only for every patient but also for each major cancer nodule. The gene hierarchy was established on the basis of gene commanding height, and the gene master regulators DAPK3,TASOR, FAM27C and ALG13 were identified in each profiled region. We delineated the molecular mechanisms by which TASOR overexpression and ALG13 silencing would selectively affect the cancer cells with little consequences for the normal cells.

Trypanosoma cruzi Promotes Transcriptomic Remodeling of the JAK/STAT Signaling and Cell Cycle Pathways in Myoblasts.

Chagas disease is responsible for more than 10,000 deaths per year and about 6 to 7 million infected people worldwide. In its chronic stage, patients can develop mega-colon, mega-esophagus, and cardiomyopathy. Differences in clinical outcomes may be determined, in part, by the genetic background of the parasite that causes Chagas disease. Trypanosoma cruzi has a high genetic diversity, and each group of strains may elicit specific pathological responses in the host. Conflicting results have been reported in studies using various combinations of mammalian host-T. cruzi strains. We previously profiled the transcriptomic signatures resulting from infection of L6E9 rat myoblasts with four reference strains of T. cruzi (Brazil, CL, Y, and Tulahuen). The four strains induced similar overall gene expression alterations in the myoblasts, although only 21 genes were equally affected by all strains. Cardiotrophin-like cytokine factor 1 (Clcf1) was one of the genes found to be consistently upregulated by the infection with all four strains of T. cruzi. This cytokine is a member of the interleukin-6 family that binds to glycoprotein 130 receptor and activates the JAK/STAT signaling pathway, which may lead to muscle cell hypertrophy. Another commonly upregulated gene was tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein theta (Ywhaq, 14-3-3 protein Θ), present in the Cell Cycle Pathway. In the present work, we reanalyzed our previous microarray dataset, aiming at understanding in more details the transcriptomic impact that each strain has on JAK/STAT signaling and Cell Cycle pathways. Using Pearson correlation analysis between the expression levels of gene pairs in biological replicas from each pathway, we determined the coordination between such pairs in each experimental condition and the predicted protein interactions between the significantly altered genes by each strain. We found that although these highlighted genes were similarly affected by all four strains, the downstream genes or their interaction partners were not necessarily equally affected, thus reinforcing the idea of the role of parasite background on host cell transcriptome. These new analyses provide further evidence to the mechanistic understanding of how distinct T. cruzi strains lead to diverse remodeling of host cell transcriptome.

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.

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.