Search

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.

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.

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.

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.

Intraventricular hemorrhage (IVH) results in periventricular inflammation, hypomyelination of the white matter, and hydrocephalus in premature infants. No effective therapy exists to prevent these disorders. Peroxisome proliferator activated receptor-γ (PPAR-γ) agonists reduce inflammation, alleviate free radical generation, and enhance microglial phagocytosis, promoting clearance of debris and red blood cells. We hypothesized that activation of PPAR-γ would enhance myelination, reduce hydrocephalus, and promote neurological recovery in newborns with IVH. These hypotheses were tested in a preterm rabbit model of IVH; autopsy brain samples from premature infants with and without IVH were analyzed. We found that IVH augmented PPAR-γ expression in microglia of both preterm human infants and rabbit kits. The treatment with PPAR-γ agonist or PPAR-γ overexpression by adenovirus delivery further elevated PPAR-γ levels in microglia, reduced proinflammatory cytokines, increased microglial phagocytosis, and improved oligodendrocyte progenitor cell (OPC) maturation in kits with IVH. Transcriptomic analyses of OPCs identified previously unrecognized PPAR-γ–induced genes for purinergic signaling, cyclic adenosine monophosphate generation, and antioxidant production, which would reprogram these progenitors toward promoting myelination. RNA-sequencing analyses of microglia revealed PPAR-γ–triggered down-regulation of several proinflammatory genes and transcripts having roles in Parkinson’s disease and amyotrophic lateral sclerosis, contributing to neurological recovery in kits with IVH. Accordingly, PPAR-γ activation enhanced myelination and neurological function in kits with IVH. This also enhanced microglial phagocytosis of red blood cells but did not reduce hydrocephalus. Treatment with PPAR-γ agonist might enhance myelination and neurological recovery in premature infants with IVH.

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.

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.

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.

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.

Simple Summary We applied the genomic fabric principles for personalized gene therapy to a case of clear cell renal cell carcinoma (ccRCC). Despite decades of research, the process of finding the molecular mechanisms responsible for the disease and, more importantly, the therapeutic solution is still a work in progress. We analyzed the transcriptomes of the chest wall metastasis, two distinct cancer nodules, and the cancer-free surrounding tissue in the surgically removed right kidney of a Fuhrman grade 3 metastatic ccRCC patient. The studies revealed that even histopathologically equally classified cancer nodules from the same kidney have different transcriptomic topologies, requiring tailored therapeutic solutions not only for each patient but even for each cancer nodule. We identified death-associated protein kinase 3 (DAPK3); transcription activation suppressor (TASOR); family with sequence similarity 27, member C, long non-coding RNA (FAM27C); and UDP-N-acetylglucosaminyltransferase subunit (ALG13) as the gene master regulators of the four profiled regions and proposed molecular mechanisms by which expression manipulation of TASOR and ALG13 may selectively destroy the cancer cells without affecting many of the normal cells. Abstract 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.

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.

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.

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 applied systems biology approaches to better understand the molecular mechanisms related with the retinal remodeling after induction of RGC degeneration. 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. Enrichment analysis showed that Complement cascade and Notch signaling pathway were the main affected pathways after ONC. To expand our studies of these two pathways, we examined the coordination of gene expressions within each pathway and with the entire transcriptome. ONC increased the number of synergistically coordinated pairs of genes and the number of similar profiles. This study provided novel findings beyond the regulation of individual gene expression and disclosed changes in the control of expression by Complement cascade and Notch signaling functional pathways important for both RGC degeneration and remodeling of the retinal tissue after ONC.

We hypothesize that distinct cell phenotypes are governed by different sets of gene master regulators (GMRs) whose strongly protected (by the homeostatic mechanisms) abundance modulates most cell processes by coordinating the expression of numerous genes from the corresponding functional pathways. Gene Commanding Height (GCH), a composite measure of gene expression control and coordination, is introduced to establish the gene hierarchy in each phenotype. If the hypothesis is true, than one can selectively destroy cancer nodules from a heterogeneous tissue by altering the expression of genes whose GCHs are high in cancer but low in normal cell phenotype. Here, we test the hypothesis and show its utility for the thyroid cancer (TC) gene therapy. First, we prove that malignant and cancer free surrounding areas of a surgically removed papillary TC (PTC) tumor are governed by different GMRs. Second, we show that stable transfection of a gene induces larger transcriptomic alterations in the cells where it has higher GCH than in other cells. For this, we profiled the transcriptomes of the papillary BCPAP and anaplastic 8505C TC cell lines before and after stable transfection with NEMP1, DDX19B, PANK2 or UBALD1. The four genes were selected to have similar expression levels but significantly different GCH scores in the two cell lines before transfection. Indeed, each of the four genes triggered larger alterations in the cells where they had larger GCH. Our results prove the feasibility of a personalized gene therapy approach that selectively targets the cancer cells from a tissue.

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.

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.

Numerous groups race to discover the gene biomarker whose alteration alone is indicative of a particular disease in all humans. Biomarkers are selected from the most frequently altered genes in large population cohorts. However, thousands of other genes are simultaneously affected, and, in each person, the same disease results from a unique, never-repeatable combination of gene alterations. Therefore, our Genomic Fabric Paradigm (GFP) switches the focus from the alteration of one particular gene to the overall change in selected groups of functionally related genes. Biomarkers are of little therapeutic value, their high alterability indicating low protection by the homeostatic mechanisms as for minor players. Instead of these most alterable genes in all patients, GFP identifies in each patient the genes whose highly protected expression governs major functional pathways by controlling the expression of numerous other genes. Smart manipulation of such (commander) genes would have the maximum therapeutic benefit not for everybody but for the treated person. The genomic fabric is defined as the transcriptome associated with the most interconnected and stably expressed network of genes responsible for a particular functional pathway. The fabric exhibits specificity with respect to race/strain, sex, age, tissue/cell type, and lifestyle and environmental factors. It remodels during development, progression of a disease, and in response to external stimuli. GFP is powered by mathematically advanced analytical tools whose application is illustrated by reprocessing data from previously published gene expression experiments.

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.

The Chicken Egg Genotoxicity Assay (CEGA) demonstrated responsiveness to various DNA-reactive chemicals requiring metabolic activation, which implies broad bioactivation capability. To assess potential metabolic competence, expression profiles of metabolic genes in the embryo-chicken fetal liver were determined using microarray technology. Fertilized chicken eggs were injected under the CEGA protocol with vehicle (deionized water [DW]), the activation-dependent carcinogens, diethylnitrosamine (DEN), and N-nitrosodiethanolamine (NDELA) at doses producing no effect on survival. Previously in CEGA, DEN produced DNA damage, whereas NDELA did not. Expressions of 463 genes known to encode for phase I and II of endo- and xenobiotic metabolism were detected on the array. DW did not affect the expression of the selected genes, deregulating less than 1% of them. In contrast, DEN at 2 mg/egg and NDELA at 4 mg/egg produced significant transcriptomic alterations, up-regulating up to 41% and down-regulating over 31% of studied genes. Both nitrosamines modulated the majority of the genes in a similar manner, sharing 64 up-regulated and 93 down-regulated genes with respect to control group, indicating similarity in the regulation of their metabolism by avian liver. Differences in gene expression between DEN and NDELA were documented for several phase I CYP 450 genes that are responsible for nitrosamine biotransformation, as well as for phase II genes that regulate detoxication reactions. These findings could underlie the difference in genotoxicity of DEN and NDELA in CEGA. In conclusion, the analysis of gene expression profiles in embryo-chicken fetal liver dosed with dialkylnitrosamines demonstrated that avian species possess a complex array of inducible genes coding for biotransformation.

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.

Neuropeptide Y (NPY) is an endogenous modulator of neuronal activity affecting both GABAergic and glutamatergic transmission. Previously, we found that oestradiol modifies the number of NPY immunoreactive neurones in the hippocampal dentate gyrus. In the present study, we investigated which oestrogen receptor type is responsible for these changes in the number of NPY-positive neurones. Furthermore, we determined the effects of oestrogen receptor activation on NPY release. Finally, we examined the contribution of oestrogen toward the remodelling of the GABAergic and glutamatergic gene networks in terms of coupling with Npy gene expression in ovariectomised rats. We found that activation of either oestrogen receptor type (ERα or ERβ) increases the number of NPY-immunopositive neurones and enhances NPY release in the dentate gyrus. We also found that, compared to oestrogen-lacking ovariectomised rats, oestrogen replacement increases the probability of synergistic/antagonistic coupling between the Npy and GABAergic synapse genes, whereas the glutamatergic synapse genes are less likely to be coupled with Npy under similar conditions. The data together suggest that oestrogens play a critical role in the regulation of NPY system activity and are also involved in the coupling/uncoupling of the Npy gene with the GABAergic and glutamatergic synapses in the female rat dentate gyrus.

Postnatal glucocorticoids (GCs) are widely used in the prevention of chronic lung disease in premature infants. Their pharmacologic use is associated with neurodevelopmental delay and cerebral palsy. However, the effect of GC dose and preparation (dexamethasone versus betamethasone) on short and long-term neurological outcomes remains undetermined, and the mechanisms of GC-induced brain injury are unclear. We hypothesized that postnatal GC would induce hypomyelination and motor impairment in a preparation- and dose-specific manner, and that GC receptor (GR) inhibition might restore myelination and neurological function in GC-treated animals. Additionally, GC-induced hypomyelination and neurological deficit might be transient. To test our hypotheses, we treated prematurely delivered rabbit pups with high (0.5 mg/kg/day) or low (0.2 mg/kg/day) doses of dexamethasone or betamethasone. Myelin basic protein (MBP), oligodendrocyte proliferation and maturation, astrocytes, transcriptomic profile, and neurobehavioral functions were evaluated. We found that high-dose GC treatment, but not low-dose, reduced MBP expression and impaired motor function at postnatal day 14. High-dose dexamethasone induced astrogliosis, betamethasone did not. Mifepristone, a GR antagonist, reversed dexamethasone-induced myelination, but not astrogliosis. Both GCs inhibited oligodendrocyte proliferation and maturation. Moreover, high-dose dexamethasone altered genes associated with myelination, cell-cycle, GR, and Mitogen-activated protein kinase. Importantly, GC-induced hypomyelination, gliosis, and motor-deficit, observed at day 14, completely recovered by day 21. Hence, high-dose, but not low-dose, postnatal GC causes reversible reductions in myelination and motor functions. GC treatment induces hypomyelination by GR-dependent genomic mechanisms, but astrogliosis by non-genomic mechanisms. GC-induced motor impairment and neurodevelopmental delay might be transient and recover spontaneously in premature infants.

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.

We previously found that, in a mouse model of Chagas cardiomyopathy, 18% of the 9390 quantified unigenes were significantly regulated by Trypanosoma cruzi infection. However, treatment with bone marrow-derived mononuclear cells (MNCs) resulted in 84% transcriptomic recovery. We have applied new algorithms to reanalyze these datasets with respect to specific pathways [Chagas disease (CHAGAS), cardiac muscle contraction (CMC) and chemokine signaling (CCS)]. In addition to the levels of expression of individual genes we also calculated gene expression variability and coordination of expression of each gene with all others. These additional measures revealed changes in the control of transcript abundances and gene networking in CHAGAS and restoration following MNC treatment, not accessible using the conventional approach limited to the average expression levels. Moreover, our weighted pathway regulation analysis incorporated the contributions of all affected genes, eliminating the arbitrary cut-off criteria of fold-change and/or p-value for significantly regulated genes. The new analyses revealed that T. cruzi infection had large transcriptomic consequences for the CMC pathway and triggered a huge cytokine signaling. Remarkably, MNC therapy not only restored normal expression levels of numerous genes, but it also recovered most of the CHAGAS, CMC and CCS fabrics that were altered by the infection.

Summary Borrelia burgdorferi, the agent of Lyme disease, is maintained in nature within an enzootic cycle involving a mammalian reservoir and an Ixodes sp. tick vector. The transmission, survival and pathogenic potential of B. burgdorferi depend on the bacterium's ability to modulate its transcriptome as it transits between vector and reservoir host. Herein, we employed an amplification-microarray approach to define the B. burgdorferi transcriptomes in fed larvae, fed nymphs and in mammalian host-adapted organisms cultivated in dialysis membrane chambers. The results show clearly that spirochetes exhibit unique expression profiles during each tick stage and during cultivation within the mammal; importantly, none of these profiles resembles that exhibited by in vitro grown organisms. Profound shifts in transcript levels were observed for genes encoding known or predicted lipoproteins as well as proteins involved in nutrient uptake, carbon utilization and lipid synthesis. Stage-specific expression patterns of chemotaxis-associated genes also were noted, suggesting that the composition and interactivities of the chemotaxis machinery components vary considerably in the feeding tick and mammal. The results as a whole make clear that environmental sensing by B. burgdorferi directly or indirectly drives an extensive and tightly integrated modulation of cell envelope constituents, chemotaxis/motility machinery, intermediary metabolism and cellular physiology. These findings provide the necessary transcriptional framework for delineating B. burgdorferi regulatory pathways throughout the enzootic cycle as well as defining the contribution(s) of individual genes to spirochete survival in nature and virulence in humans.

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.

Injury of the CA1 subregion induced by a single injection of kainic acid (1 × KA) in juvenile animals (P20) is attenuated in animals with two prior sustained neonatal seizures on P6 and P9. To identify gene candidates involved in the spatially protective effects produced by early-life conditioning seizures we profiled and compared the transcriptomes of CA1 subregions from control, 1 × KA- and 3 × KA-treated animals. More genes were regulated following 3 × KA (9.6%) than after 1 × KA (7.1%). Following 1 × KA, genes supporting oxidative stress, growth, development, inflammation and neurotransmission were upregulated (e.g. Cacng1, Nadsyn1, Kcng1, Aven, S100a4, GFAP, Vim, Hrsp12 and Grik1). After 3 × KA, protective genes were differentially over-expressed [e.g. Cat, Gpx7, Gad1, Hspa12A, Foxn1, adenosine A1 receptor, Ca(2+) adaptor and homeostasis proteins, Cacnb4, Atp2b2, anti-apoptotic Bcl-2 gene members, intracellular trafficking protein, Grasp and suppressor of cytokine signaling (Socs3)]. Distinct anti-inflammatory interleukins (ILs) not observed in adult tissues [e.g. IL-6 transducer, IL-23 and IL-33 or their receptors (IL-F2 )] were also over-expressed. Several transcripts were validated by real-time polymerase chain reaction (QPCR) and immunohistochemistry. QPCR showed that casp 6 was increased after 1 × KA but reduced after 3 × KA; the pro-inflammatory gene Cox1 was either upregulated or unchanged after 1 × KA but reduced by ~70% after 3 × KA. Enhanced GFAP immunostaining following 1 × KA was selectively attenuated in the CA1 subregion after 3 × KA. The observed differential transcriptional responses may contribute to early-life seizure-induced pre-conditioning and neuroprotection by reducing glutamate receptor-mediated Ca(2+) permeability of the hippocampus and redirecting inflammatory and apoptotic pathways. These changes could lead to new genetic therapies for epilepsy.

Steroid hormones affect neuronal excitability. Their final effect on neuronal function depends on many factors, including the hormone type, its active metabolite (or their combination), treatment duration, dose, brain region, sex, hormone receptor sensitivity, and local steroid hormone production in the brain. Understanding of how these distinct mechanisms are orchestrated to regulate neuronal activity may lead to the design of better treatment strategies for patients with epilepsy, as well as other hormone-sensitive brain disorders.

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.

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.

Chronic chagasic cardiomyopathy is a leading cause of heart failure in Latin American countries, being associated with intense inflammatory response and fibrosis. We have previously shown that bone marrow mononuclear cell (BMC) transplantation improves inflammation, fibrosis, and ventricular diameter in hearts of mice with chronic Chagas disease. Here we investigated the transcriptomic recovery induced by BMC therapy by comparing the heart transcriptomes of control, chagasic, and BMC transplanted mice. Out of the 9390 unique genes quantified in all samples, 1702 had their expression altered in chronic chagasic hearts compared to those of normal mice. Major categories of significantly upregulated genes were related to inflammation, fibrosis and immune responses, while genes involved in mitochondrion function were downregulated. When BMC-treated chagasic hearts were compared to infected mice, 96% of the alterations detected in infected hearts were restored to normal levels, although an additional 109 genes were altered by treatment. Transcriptomic recovery, a new measure that considers both resotrative and side effects of treatment, was remarkably high (84%). Immunofluorescence and morphometric analyses confirmed the effects of BMC therapy in the pattern of inflammatory-immune response and expression of adhesion molecules. In conclusion, by using large-scale gene profiling for unbiased assessment of therapeutic efficacy we demonstrate immunomodulatory effects of BMC therapy in chronic chagasic cardiomyopathy and identify potentially relevant factors involved in the pathogenesis of the disease that may provide new therapeutic targets.

We have used cDNA microarrays to compare gene expression profiles in brains from normal mice to those infected with the ANKA strain of Plasmodium berghei ,a model of cerebral malaria. For each of three brains in each group, we computed ratios of all quantifiable genes with a composite reference sample and then computed ratios of gene expression in infected brains compared to untreated controls. Of the almost 12,000 unigenes adequately quantified in all arrays, approximately 3% were significantly downregulated ( P< 0.05, ≥ 50 %f old change) and about 7% were upregulated. Upon inspection of the lists of regulated genes, we identified a high number encoding proteins of importance to normal brain function or associated with neuropathology, including genes that encode for synaptic proteins or genes involved in cerebellar function as well as genes important in certain neurological diseases such as Alzheimer's disease or autism. These results emphasize the important impact of malarial infection on gene expression in the brain and provide potential biomarkers that may provide novel therapeutic targets to ameliorate the neurological sequelae of this infection.

Trypanosoma cruzi infection is a major cause of cardiomyopathy. Previous gene profiling studies of infected mouse hearts have revealed prominent changes in gene expression within many functional pathways. This variety of transcriptomic changes in infected mice raises the question of whether gene expression alterations in whole hearts are due to changes in infected cardiac myocytes or other cells or even to systemic effects of the infection on the heart. We employed microarrays to examine infected cardiac myocyte cultures 48 h post-infection. Statistical comparison of gene expression levels of 7624 well annotated unigenes in four independent cultures of infected and uninfected myocytes detected substantial (or=1.5 absolute fold changes) in 420 (5.5%) of the sampled genes. Major categories of affected genes included those involved in immune response, extracellular matrix and cell adhesion. These findings on infected cardiac myocytes in culture reveal that alterations in cardiac gene expression described in Chagas disease are the consequence of both direct infection of the myocytes themselves as well as resulting from the presence of other cell types in the myocardium and systemic effects of infection.

Changes in gravitational force such as that experienced by astronauts during space flight induce a redistribution of fluids from the caudad to the cephalad portion of the body together with an elimination of normal head-to-foot hydrostatic pressure gradients. To assess brain gene profile changes associated with microgravity and fluid shift, a large-scale analysis of mRNA expression levels was performed in the brains of 2-week control and hindlimb-unloaded (HU) mice using cDNA microarrays. Although to different extents, all functional categories displayed significantly regulated genes indicating that considerable transcriptomic alterations are induced by HU. Interestingly, the TIC class (transport of small molecules and ions into the cells) had the highest percentage of up-regulated genes, while the most down-regulated genes were those of the JAE class (cell junction, adhesion, extracellular matrix). TIC genes comprised 16% of those whose expression was altered, including sodium channel, nonvoltage-gated 1 beta (Scnn1b), glutamate receptor (Grin1), voltage-dependent anion channel 1 (Vdac1), calcium channel beta 3 subunit (Cacnb3) and others. The analysis performed by Gene-MAPP revealed several altered protein classes and functional pathways such as blood coagulation and immune response, learning and memory, ion channels and cell junction. In particular, data indicate that HU causes an alteration in hemostasis which resolves in a shift toward a more hyper-coagulative state with an increased risk of venous thrombosis. Furthermore, HU treatment seems to impact on key steps of synaptic plasticity and learning processes.

Fluid-induced shear stress is widely recognized as an important biophysical signal in cell-cell mechanotransduction. To identify cellular signaling pathways that are regulated by fluid shear stress, we applied the unbiased approach of transcriptional profiling. Our cDNA array analysis detected that 1,165 of the 6,288 sampled unigenes were significantly affected by pulsatile fluid flow. GenMapp 2.1 analysis revealed pathways of genes regulated by shear stress: angiogenesis, blood vessel morphogenesis, regulation of endothelial cell proliferation, and prostaglandin biosynthesis. Individual genes significantly up-/downregulated by shear stress included vascular endothelial growth factor A (Vegf a), cysteine-rich protein 61 (Cyr61), platelet-derived growth factor-alpha (Pdgf a), connective tissue growth factor (Ctgf), Neuropilin 1 (Nrp1), angiotensin II receptor, type 1 a (Agtr1 a) and fibroblast growth factor 1 (Fgf1). Based on these findings, we hypothesize that fluid shear stress-regulated Vegf most likely stimulates MC3T3-E1 cells through autocrine/paracrine release and may provide a powerful recruitment signal for osteoclasts, endothelial cells, and/or stem cells during bone remodeling.

Gap junctions and purinergic P2 receptors (P2Rs) can be regarded as belonging to a common functional unit, given that they are involved in the transmission of calcium signals between cells. We have previously shown that deletion of the Gja1 gene alters expression levels of numerous genes encoding proteins with diverse functions, including purinergic receptors (P2Rs), and have found that genes synergistically or antagonistically expressed in wild-type tissues are more prone to be similarly or oppositely regulated in Cx43-nulls. We have now explored the use of coordination analysis of gene expression as a strategy to identify interlinked genes encoding functionally related proteins and pull-downs to evaluate their interlinkage. Our findings indicate that, in brain and in cultured astrocytes, several of these coexpressed genes encode proteins that are components of P2R signal-transduction pathways and/or directly interact with these receptors, including the gap junction protein connexin43 (Cx43) and Cx45 as well as pannexins. It is proposed that coordination analysis of gene expression may provide a novel unbiased strategy for the identification of proteins belonging to supramolecular complexes.

We have found that deletion of genes encoding the gap junction proteins Cx43, Cx32 and Cx36 alter the expression levels of large numbers of genes in mouse brain located on all chromosomes and encoding proteins from all major functional categories. Gene regulation in Cx32 and Cx43 null brains was more similar than that in the Cx36 null brain, suggesting the possibility of transcriptomic controls exerted by both genes on both astrocytes and oligodendrocytes. In order to explore the nature of expression linkage among the genes, we examined coordinated expression patterns in wild-type and connexin null brains. Coordination with Cx43 in wild-type brain predicted regulation in Cx43 nulls with considerable accuracy. Moreover, interlinkage within gene networks was greatly perturbed in the Cx43 null brain. These findings suggest several principles regarding regulatory transcriptomic networks involving gap junction genes and raise the issue of the underlying cause of connexin null phenotypes as well as mechanisms of regulation.

International audience; The purpose of this paper is to provide a brief overview of current thinking on the role of connexins, in particular Cx43, in growth regulation, and a more detailed discussion as to potential mechanisms involved with an emphasis on gene expression. While the precise molecular mechanism by which connexins can affect the growth of normal or tumor cells remains elusive, a number of exciting reports have expanded our understanding and are presented in some detail. Thus, we will discuss (Section 2): the role of protein-protein interactions in integrating connexins into multiple signal transduction pathways; phosphorylation at specific sites and reversal of growth inhibition; the role of the carboxy-terminal regulatory domain as a signaling molecule. Some of our latest work on the potential functions of endogenously produced carboxy-terminal fragments of Cx43 are also presented (Section 3). Finally, Section 4 will pay tribute to the rapidly emerging realization that connexins such as Cx43 and Cx32 exert important and extensive effects on gene expression, particularly those genes linked to growth regulation.