Identification of the Key Pathways and Genes in Hypoxia Pulmonary Arterial Hypertension Following Intrauterine Growth Retardation

High-throughput sequencing and weighted gene co-expression network analysis (WGCNA) were used to identify susceptibility modules and genes in liver tissue for the hypoxic pulmonary arterial hypertension (PAH) animal model following intrauterine growth retardation (IUGR). A total of 5,000 genes were clustered into eight co-expression modules via WGCNA. Module blue was mostly significantly correlated with the IUGR–hypoxia group. Gene Ontology analysis showed that genes in the module blue were mainly enriched in the fatty acid metabolic process, lipid modification, and fatty acid catabolic process. The Kyoto Encyclopedia of Genes and Genomes enrichment analyses showed that the genes in module blue were mainly associated with fatty acid metabolism, PPAR signaling pathway, and biosynthesis of unsaturated fatty acids. In addition, the maximal clique centrality method was used to identify the hub genes in the subnetworks, and the obtained results were verified using real-time quantitative PCR. Finally, we identified that four genes including Cyp2f4, Lipc, Acadl, and Hacl1 were significantly associated with IUGR-hypoxia. Our study identified a module and several key genes that acted as essential components in the etiology of the long-term metabolic consequences in hypoxia PAH following IUGR.