Generation and characterization of the blood transcriptome of Macaca thibetana and comparative analysis with M. mulatta

Physiology and Transcriptomics Analysis Reveal the Contribution of Lungs on High-Altitude Hypoxia Adaptation in Tibetan Sheep

The Tibetan sheep is an indigenous species on the Tibetan plateau with excellent adaptability to high-altitude hypoxia and is distributed at altitudes of 2500–5000 m. The high-altitude hypoxia adaptation of Tibetan sheep requires adaptive reshaping of multiple tissues and organs, especially the lungs. To reveal the mechanisms of adaptation at the tissue and molecular levels in the lungs of Tibetan sheep under hypoxic conditions at different altitudes, we performed light and electron microscopic observations, transcriptomic sequencing, and enzyme-linked immunosorbent assay studies on the lungs of Tibetan sheep from three altitudes (2500, 3500, and 4500 m). The results showed that in addition to continuous increase in pulmonary artery volume, thickness, and elastic fiber content with altitude, Tibetan sheep increase the hemoglobin concentration at an altitude of 3500 m, while they decrease the Hb concentration and increase the surface area of gas exchange and capacity of the blood at an altitude of 4500 m. Other than that, some important differentially expressed genes related to angiogenesis (FNDC1, HPSE, and E2F8), vasomotion and fibrogenesis (GJA4, FAP, COL1A1, COL1A2, COL3A1, and COL14A1), and gas transport (HBB, HBA1, APOLD1, and CHL1) were also identified; these discoveries at the molecular level explain to some extent the physiological findings. In conclusion, the lungs of Tibetan sheep adopt different strategies when adapting to different altitudes, and these findings are valuable for understanding the basis of survival of indigenous species on the Tibetan plateau.

Genomic Copy Number Variation Study of Nine Macaca Species Provides New Insights into Their Genetic Divergence, Adaptation, and Biomedical Application

Copy number variation (CNV) can promote phenotypic diversification and adaptive evolution. However, the genomic architecture of CNVs among Macaca species remains scarcely reported, and the roles of CNVs in adaptation and evolution of macaques have not been well addressed. Here, we identified and characterized 1,479 genome-wide hetero-specific CNVs across nine Macaca species with bioinformatic methods, along with 26 CNV-dense regions and dozens of lineage-specific CNVs. The genes intersecting CNVs were overrepresented in nutritional metabolism, xenobiotics/drug metabolism, and immune-related pathways. Population-level transcriptome data showed that nearly 46% of CNV genes were differentially expressed across populations and also mainly consisted of metabolic and immune-related genes, which implied the role of CNVs in environmental adaptation of Macaca. Several CNVs overlapping drug metabolism genes were verified with genomic quantitative polymerase chain reaction, suggesting that these macaques may have different drug metabolism features. The CNV-dense regions, including 15 first reported here, represent unstable genomic segments in macaques where biological innovation may evolve. Twelve gains and 40 losses specific to the Barbary macaque contain genes with essential roles in energy homeostasis and immunity defense, inferring the genetic basis of its unique distribution in North Africa. Our study not only elucidated the genetic diversity across Macaca species from the perspective of structural variation but also provided suggestive evidence for the role of CNVs in adaptation and genome evolution. Additionally, our findings provide new insights into the application of diverse macaques to drug study.

The Transcriptomic Landscape of Yaks Reveals Molecular Pathways for High Altitude Adaptation

Yak is one of the largest native mammalian species at the Himalayas, the highest plateau area in the world with an average elevation of >4,000 m above the sea level. Yak is well adapted to high altitude environment with a set of physiological features for a more efficient blood flow for oxygen delivery under hypobaric hypoxia. Yet, the genetic mechanism underlying its adaptation remains elusive. We conducted a cross-tissue, cross-altitude, and cross-species study to characterize the transcriptomic landscape of domestic yaks. The generated multi-tissue transcriptomic data greatly improved the current yak genome annotation by identifying tens of thousands novel transcripts. We found that among the eight tested tissues (lung, heart, kidney, liver, spleen, muscle, testis, and brain), lung and heart are two key organs showing adaptive transcriptional changes and >90% of the cross-altitude differentially expressed genes in lung display a nonlinear regulation. Pathways related to cell survival and proliferation are enriched, including PI3K-Akt, HIF-1, focal adhesion, and ECM–receptor interaction. These findings, in combination with the comprehensive transcriptome data set, are valuable to understanding the genetic mechanism of hypoxic adaptation in yak.

Comparative genome-wide survey of single nucleotide variation uncovers the genetic diversity and potential biomedical applications among six Macaca species

Macaca is of great importance in evolutionary and biomedical research. Aiming at elucidating genetic diversity patterns and potential biomedical applications of macaques, we characterized single nucleotide variations (SNVs) of six Macaca species based on the reference genome of Macaca mulatta. Using eight whole-genome sequences, representing the most comprehensive genomic SNV study in Macaca to date, we focused on discovery and comparison of nonsynonymous SNVs (nsSNVs) with bioinformatic tools. We observed that SNV distribution patterns were generally congruent among the eight individuals. Outlier tests of nsSNV distribution patterns detected 319 bins with significantly distinct genetic divergence among macaques, including differences in genes associated with taste transduction, homologous recombination, and fat and protein digestion. Genes with specific nsSNVs in various macaques were differentially enriched for metabolism pathways, such as glycolysis, protein digestion and absorption. On average, 24.95% and 11.67% specific nsSNVs were putatively deleterious according to PolyPhen2 and SIFT4G, respectively, among which the shared deleterious SNVs were located in 564–1981 genes. These genes displayed enrichment signals in the ‘obesity-related traits’ disease category for all surveyed macaques, confirming that they were suitable models for obesity related studies. Additional enriched disease categories were observed in some macaques, exhibiting promising potential for biomedical application. Positively selected genes identified by PAML in most tested Macaca species played roles in immune and nervous system, growth and development, and fat metabolism. We propose that metabolism and body size play important roles in the evolutionary adaptation of macaques.