Post-transcriptional regulation through alternative splicing in the lungs of Tibetan pigs under hypoxia

Abstract

In multicellular organisms, alternative splicing (AS) is central to the regulation of multiple biological processes. To further elucidate the adaptive strategy of AS in the lungs of Tibetan pigs in response to hypoxia, we identified and analyzed five basic AS types and 59,930 AS events in 18,179 genes. We found that approximately 65.10% of the total expressed genes underwent AS in the lungs of Tibetan pigs at a high altitude (TH). The frequencies of AS events were similar among the different groups (5.06–5.30 events in each gene on average). Skipped exons (SEs) were the predominant type of AS event, followed by mutually exclusive exons (MXEs), alternative 3′ splice sites (A3SSs) and alternative 5′ splice sites (A5SSs). Retained introns (RIs), the remaining type of AS event, showed lower frequencies. Further comparison analysis of differentially expressed genes (DEGs) and differentially spliced genes (DSGs) identified 2,209 differential splicing events in the above 18,000 expressed genes, including 918 increased and 1,291 decreased splicing events between the TH and Tibetan pigs at a low altitude (TL) groups. We identified 227 hypoxia-related genes involved in lung development that were differentially regulated through AS. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis clearly identified many DEGs and DSGs at high or low altitude. Seven pathways in the top 20 enriched KEGG terms overlapped for the DEGs and DSGs, including the chemokine signaling pathway, B cell receptor signaling pathway, and cytokine–cytokine receptor interaction, which exert many immunoregulatory and inflammatory actions critical to the lung under hypoxia. Twelve pathways overlapped in hypoxic DEGs and DSGs and included antigen processing, presentation and biosynthesis. ​GO analysis of the DEGs and DSGs among the four groups showed that numerous GO terms were enriched in the biological category, and the proportion of genes with downregulated expression was greater among 227 hypoxic genes than that of all genes. The results suggest that AS plays an essential role in the regulation of gene expression during hypoxia and that numerous genes involved in lung development are differentially regulated through AS.