Advanced Molecular Solutions for Cancer Therapy-The Good, the Bad, and the Ugly of the Biomarker Paradigm.

Identifying the most effective actionable molecules whose "smart" manipulation might selectively kill/slow down/stop the proliferation of cancer cells, with few side effects on the normal cells of the tissue, was for decades the single major objective of countless investigators [...].

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.

Oestrogen treatment restores dentate gyrus development in premature newborns by IGF1 regulation.

Prematurely-born infants cared for in the neonatal units suffer from memory and learning deficits. Prematurity diminishes neurogenesis and synaptogenesis in the hippocampal dentate gyrus (DG). This dysmaturation of neurons is attributed to elevated PSD95, NMDR2A, and IGF1 levels. Since oestrogen treatment plays key roles in the development and plasticity of DG, we hypothesized that 17ß-estradiol (E2) treatment would ameliorate neurogenesis and synaptogenesis in the DG, reversing cognitive deficits in premature newborns. Additionally, E2-induced recovery would be mediated by IGF1 signalling. These hypotheses were tested in a rabbit model of prematurity and nonmaternal care, in which premature kits were gavage-fed and reared by laboratory personnel. We compared E2- and vehicle-treated preterm kits for morphological, molecular, and behavioural parameters. We also treated kits with oestrogen degrader, RAD1901, and assessed IGF1 signalling. We found that E2 treatment increased the number of Tbr2+ and DCX+ neuronal progenitors and increased the density of glutamatergic synapses in the DG. E2 treatment restored PSD95 and NMDAR2A levels and cognitive function in preterm kits. Transcriptomic analyses showed that E2 treatment contributed to recovery by influencing interactions between IGF1R and neurodegenerative, as well as glutamatergic genes. ERα expression was reduced on completion of E2 treatment at D7, followed by D30 elevation. E2-induced fluctuation in ERα levels was associated with a reciprocal elevation in IGF1/2 expression at D7 and reduction at D30. ERα degradation by RAD1901 treatment enhanced IGF1 levels, suggesting ERα inhibits IGF1 expression. E2 treatment alleviates the prematurity-induced maldevelopment of DG and cognitive dysfunctions by regulating ERα and IGF1 levels.

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.

Molecules at Play in Cancer.

Despite its wide range of incidence, cancer can spontaneously occur in any part of the body and invade regions other than the originally affected tissue [...].

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.

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.

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.

PPAR-γ activation enhances myelination and neurological recovery in premature rabbits with intraventricular hemorrhage.

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.

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.

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.

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.

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.

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.

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.

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.

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.