The hypoxic environment imposes severe selective pressure on species living at high altitude. To understand the genetic bases of adaptation to high altitude in dogs, we performed whole-genome sequencing of 60 dogs including five breeds living at continuous altitudes along the Tibetan Plateau from 800 to 5100 m as well as one European breed. More than 150× sequencing coverage for each breed provides us with a comprehensive assessment of the genetic polymorphisms of the dogs, including Tibetan Mastiffs. Comparison of the breeds from different altitudes reveals strong signals of population differentiation at the locus of hypoxia-related genes including endothelial Per-Arnt-Sim (PAS) domain protein 1 (EPAS1) and beta hemoglobin cluster. Notably, four novel nonsynonymous mutations specific to high-altitude dogs are identified atEPAS1, one of which occurred at a quite conserved site in the PAS domain. The association testing betweenEPAS1genotypes and blood-related phenotypes on additional high-altitude dogs reveals that the homozygous mutation is associated with decreased blood flow resistance, which may help to improve hemorheologic fitness. Interestingly,EPAS1was also identified as a selective target in Tibetan highlanders, though no amino acid changes were found. Thus, our results not only indicate parallel evolution of humans and dogs in adaptation to high-altitude hypoxia, but also provide a new opportunity to study the role ofEPAS1in the adaptive processes.

Tibetan chicken, unlike their lowland counterparts, exhibit specific adaptations to high-altitude conditions. The genetic mechanisms of such adaptations in highland chickens were determined by resequencing the genomes of four highland (Tibetan and Lhasa White) and four lowland (White Leghorn, Lindian, and Chahua) chicken populations. Our results showed an evident genetic admixture in Tibetan chickens, suggesting a history of introgression from lowland gene pools. Genes showing positive selection in highland populations were related to cardiovascular and respiratory system development, DNA repair, response to radiation, inflammation, and immune responses, indicating a strong adaptation to oxygen scarcity and high-intensity solar radiation. The distribution of allele frequencies of nonsynonymous single nucleotide polymorphisms between highland and lowland populations was analyzed using chi-square test, which showed that several differentially distributed genes with missense mutations were enriched in several functional categories, especially in blood vessel development and adaptations to hypoxia and intense radiation. RNA sequencing revealed that several differentially expressed genes were enriched in gene ontology terms related to blood vessel and respiratory system development. Several candidate genes involved in the development of cardiorespiratory system (FGFR1, CTGF, ADAM9, JPH2, SATB1, BMP4, LOX, LPR, ANGPTL4, and HYAL1), inflammation and immune responses (AIRE, MYO1F, ZAP70, DDX60, CCL19, CD47, JSC, and FAS), DNA repair, and responses to radiation (VCP, ASH2L, and FANCG) were identified to play key roles in the adaptation to high-altitude conditions. Our data provide new insights into the unique adaptations of highland animals to extreme environments.

Transcriptome profiles established using high-throughput sequencing can be effectively used for screening genome-wide differentially expressed genes (DEGs). RNA sequences (from RNA-seq) and microRNA sequences (from miRNA-seq) from the tissues of longissimus dorsi muscle of two indigenous Chinese pig breeds (Diannan Small-ear pig [DSP] and Tibetan pig [TP]) and two introduced pig breeds (Landrace [LL] and Yorkshire [YY]) were examined using HiSeq 2000 to identify and compare the differential expression of functional genes related to muscle growth and lipid deposition. We obtained 27.18 G clean data through the RNA-seq and detected that 18,208 genes were positively expressed and 14,633 of them were co-expressed in the muscle tissues of the four samples. In all, 315 DEGs were found between the Chinese pig group and the introduced pig group, 240 of which were enriched with functional annotations from the David database and significantly enriched in 27 Gene Ontology (GO) terms that were mainly associated with muscle fiber contraction, cadmium ion binding, response to organic substance and contractile fiber part. Based on functional annotation, we identified 85 DEGs related to growth traits that were mainly involved in muscle tissue development, muscle system process, regulation of cell development, and growth factor binding, and 27 DEGs related to lipid deposition that were mainly involved in lipid metabolic process and fatty acid biosynthetic process. With miRNA-seq, we obtained 23.78 M reads and 320 positively expressed miRNAs from muscle tissues, including 271 known pig miRNAs and 49 novel miRNAs. In those 271 known miRNAs, 20 were higher and 10 lower expressed in DSP-TP than in LL-YY. The target genes of the 30 miRNAs were mainly participated in MAPK, GnRH, insulin and Calcium signaling pathway and others involved cell development, growth and proliferation, etc. Combining the DEGs and the differentially expressed (DE) miRNAs, we drafted a network of 46 genes and 18 miRNAs for regulating muscle growth and a network of 15 genes and 16 miRNAs for regulating lipid deposition. We identified that CAV2, MYOZ2, FRZB, miR-29b, miR-122, miR-145-5p and miR-let-7c, etc, were key genes or miRNAs regulating muscle growth, and FASN, SCD, ADORA1, miR-4332, miR-182, miR-92b-3p, miR-let-7a and miR-let-7e, etc, were key genes or miRNAs regulating lipid deposition. The quantitative expressions of eight DEGs and seven DE miRNAs measured with real-time PCR certified that the results of differential expression genes or miRNAs were reliable. Thus, 18,208 genes and 320 miRNAs were positively expressed in porcine longissimus dorsi muscle. We obtained 85 genes and 18 miRNAs related to muscle growth and 27 genes and 16 miRNAs related to lipid deposition, which provided new insights into molecular mechanism of the economical traits in pig.

The domestic pig currently indigenous to the Tibetan highlands is supposed to have been introduced during a continuous period of colonization by the ancestors of modern Tibetans. However, there is no direct genetic evidence of either the local origin or exotic migration of the Tibetan pig. We analyzed mtDNA hypervariable segment I (HVI) variation of 218 individuals from seven Tibetan pig populations and 1,737 reported mtDNA sequences from domestic pigs and wild boars across Asia. The Bayesian consensus tree revealed a main haplogroup M and twelve minor haplogroups, which suggested a large number of small scale in situ domestication episodes. In particular, haplogroups D1 and D6 represented two highly divergent lineages in the Tibetan highlands and Island Southeastern Asia, respectively. Network analysis of haplogroup M further revealed one main subhaplogroup M1 and two minor subhaplogroups M2 and M3. Intriguingly, M2 was mainly distributed in Southeastern Asia, suggesting for a local origin. Similar with haplogroup D6, M3 was mainly restricted in Island Southeastern Asia. This pattern suggested that Island Southeastern Asia, but not Southeastern Asia, might be the center of domestication of the so-called Pacific clade (M3 and D6 here) described in previous studies. Diversity gradient analysis of major subhaplogroup M1 suggested three local origins in Southeastern Asia, the middle and downstream regions of the Yangtze River, and the Tibetan highlands, respectively. We identified two new origin centers for domestic pigs in the Tibetan highlands and in the Island Southeastern Asian region.

Growth rate and meat quality, two economically important traits in pigs, are controlled by multiple genes and biological pathways. In the present study, we performed a proteomic analysis of longissimus dorsi muscle from six-month-old pigs from two Chinese native mini-type breeds (TP and DSP) and two introduced western breeds (YY and LL) using isobaric tag for relative and absolute quantification (iTRAQ). In total, 4,815 peptides corresponding to 969 proteins were detected. Comparison of expression patterns between TP-DSP and YY-LL revealed 288 differentially expressed proteins (DEPs), of which 169 were up-regulated and 119 were down-regulated. Functional annotation suggested that 28 DEPs were related to muscle growth and 15 to lipid deposition. Protein interaction network predictions indicated that differences in muscle growth and muscle fibre between TP-DSP and YY-LL groups were regulated by ALDOC, ENO3, PGK1, PGK2, TNNT1, TNNT3, TPM1, TPM2, TPM3, MYL3, MYH4, and TNNC2, whereas differences in lipid deposition ability were regulated by LPL, APOA1, APOC3, ACADM, FABP3, ACADVL, ACAA2, ACAT1, HADH, and PECI. Twelve DEPs were analysed using parallel reaction monitoring to confirm the reliability of the iTRAQ analysis. Our findings provide new insights into key proteins involved in muscle growth and lipid deposition in the pig.

Tibetan pigs that inhabit the Tibetan Plateau exhibit striking phenotypic and physiological differences from lowland pigs, and have adapted well to extreme conditions. However, the mechanisms involved in regulating gene expression at high altitude in these animals are not fully understood. In this study, we obtained transcriptomic and proteomic data from the heart tissues of Tibetan and Yorkshire pigs raised in the highlands (TH and YH) and lowlands (TL and YL) via RNA-seq and iTRAQ (isobaric tags for relative and absolute quantitation) analyses, respectively. Comparative analyses of TH vs. YH, TH vs.TL, TL vs. YL, and YH vs. YL yielded 299, 169, 242, and 368 differentially expressed genes (DEGs), and 473, 297, 394, and 297 differentially expressed proteins (DEPs), respectively. By functional annotation of these DEGs and DEPs, genes that were enriched in the HIF-1 signaling pathway (NPPA, ERK2, ENO3, and EGLN3), VEGF signaling pathway (ERK2, A2M, FGF1, CTGF, and DPP4), and hypoxia-related processes (CRYAB, EGLN3, TGFB2, DPP4, and ACE) were identified as important candidate genes for high-altitude adaptation in the Tibetan pig. This study enhances our understanding of the molecular mechanisms involved in hypoxic adaptation in pigs, and furthers our understanding of human hypoxic diseases.