Ground tit genome reveals avian adaptation to living at high altitudes in the Tibetan plateau

Cardiovascular Plasticity and Adaptation of High-Altitude Birds and Mammals

Exposure to a hypoxic environment at high altitudes imposes severe pressure on animals living there, which utilize substantial cardiovascular and respiratory responses to meet the physiological challenge of oxygen requirement. These responses may result from phenotypic plasticity through short-term exposure (i.e., within a generation) to a new environment or shaped by adaptation (i.e., many generations) through long-term evolution. For example, plasticity triggers a sympathetic-mediated adrenergic response, resulting in an elevation of heart rate and hypoxia-induced pulmonary vasoconstriction that eventually contributes to pulmonary hypertension in some animals. Adaptation to high altitudes can drive an increase in muscular capillarization and adaptive cardiac growth, which promote oxygen diffusion and transportation. Exposure to a high-altitude hypoxic environment stimulates excessive erythropoiesis, which has maladaptive effects and contributes to chronic mountain sickness. Maladaptation caused by plasticity at early stages can be reversed during adaptation. Despite extensive research on high-altitude adaptation, the phenotypic changes and genetic variations in cardiovascular systems responding to high-altitude hypoxia remain insufficiently integrated across taxa. While genomic and transcriptomic studies have advanced our understanding, a cross-taxa comparison of cardiovascular adaptations is still incomplete. We here review recent literature on phenotypic plasticity, adaptations, and genetic and transcriptional basis of cardiovascular systems of mammals and birds living in high altitudes with respect to their duration of exposure at high altitudes. By integrating and comparing data across mammalian and avian species, we aim to provide a framework for understanding the plasticity and adaptation of the cardiovascular system in high-altitude environments.

Chromosome-level assembly of Pseudopodoces humilis genome: A resource for avian evolutionary studies

Pseudopodoces humilis is a small passerine bird predominantly found in the mid-latitude regions of the Tibetan Plateau in Asia. A chromosome-level reference genome assembly for P. humilis was generated using PacBio CLR with Hi-C. The final genome assembly spans approximately 1.096 Gb, consisting of 1,968 contigs with a Contig N50 of 32.246 Mb, and was evaluated to be 95.60% complete using BUSCO. Hi-C chromosome mapping resulted in 33 chromosome sequences, which enabled the ordering and orientation of 329 contigs, with chromosome lengths ranging from 2.08 Mb to 152.13 Mb, covering 95.85% of the total genome sequence. Repetitive sequences comprised 144.91 Mb of the genome. A total of 381 tRNA, 507 non-coding RNA (ncRNA), and 205 rRNA were identified. In addition, we identified 17,108 protein-coding genes and 29,473 proteins, comprising a total of 17,236,726 amino acids. This high-quality genome assembly provides a strong genomic foundation for exploring critical questions in evolutionary genetics, phylogenomics, and the molecular mechanisms of adaptation - key areas for understanding biodiversity and species resilience amidst changing environments.

The special adaptation to hypoxia facilitated the expansion of the Asian house rat (Rattus tanezumi) into Tibet but not other Rattus species

Rattus species are thought to live only at altitudes less than 2500 m, but the Asian house rat (R. tanezumi) (RT) has recently expanded to altitudes greater than 3500 m in China. Other Rattus species, especially brown rats (R. norvegicus) (RN), still reach only low altitudes on the Tibetan Plateau. Comparative genomics revealed the positive selection of hypoxia-inducible transcription factors 1 and 2 (HIFs) in RT, with the rapid evolution of HIF pathway genes in RT and Mus musculus (MM) but not RN or R. rattus. Population genomics revealed that genes associated with energy metabolism and oxygen transport were positively selected in RT compared with the other four Rattus species, and two specific substitutions (arginine 31 serine and leucine 33 methionine) were identified in the hemoglobin subunit beta (HBB) in RT. The above results suggested that RT possesses unique genetic adaptations to hypoxia, which was further confirmed by behavioral experiments on RT and RN. Normobaric hypoxia significantly reduced locomotion in RN but not in RT. Moreover, through intraspecific transcriptome analysis, the expression of Hbb and genes related to angiogenesis, oxygen transport, and glycolysis was upregulated, and the expression of genes associated with immunological functions in the liver, lungs, and/or sperm was downregulated in RT compared to those in RN. Interspecific transcriptome analysis further revealed that HIF-1α plays a role in modulating the hypoxic adaptation of RT rather than RN. Our work provides genomic, behavioral, and physiological insights into why RT, but not other Rattus species, could invade the Tibetan Plateau.

Whole-genome relaxed selection and molecular constraints in Triplophysa under adapted Qinghai-Tibetan Plateau

High-altitude environments are inhospitable, but Triplophysa, the largest taxon among the three major fish groups in the Qinghai-Tibetan Plateau (QTP), is an exception. However, the evolutionary profiling of the common ancestor and its contribution to the adaptation of existing QTP native species is unclear.We researched the comparative genomics of Triplophysa species and found that the genome-wide genes of Triplophysa and its ancestry have the characteristics of rapid evolution.Moreover, the rapid evolution of the ancestral genes was caused by relaxed selection. Natural selection analysis showed that more ancestral relaxed selection genes were under strongly purifying selection and showed higher expression in QTP endemic Triplophysa species.The change in natural selection might be associated with the adaptation to QTP. It should be noted that SPT5 homolog, DSIF elongation factor subunit (supt5h) experienced relaxed selection in common ancestral populations of Triplophysa but under purifying selection in extant species, which might be related to hypoxia adaptation of QTP. In summary, the extant species in different environments were used to infer the evolutionary profile of the common ancestor and to identify candidate genes based on changes in natural selection. Our work might provide new clues for understanding adaptation to extreme environments.

Selection on the vascular-remodeling BMPER gene is associated with altitudinal adaptation in an insular lizard

High altitude imposes several extreme constraints on life, such as low oxygen pressure and high levels of ultraviolet radiation, which require specialized adaptations. Many studies have focused on how endothermic vertebrates respond to these challenging environments, but there is still uncertainty on how ectotherms adapt to these conditions. Here, we used whole-genome sequencing of low-altitude (100-600 m) and high-altitude (3,550 m) populations of the wide-ranging Tenerife lizard Gallotia galloti to uncover signatures of selection for altitudinal adaptation. The studied populations show reduced differentiation, sharing similar patterns of genetic variation. Selective sweep mapping suggests that signatures of adaptation to high altitude are not widespread across the genome, clustering in a relatively small number of genomic regions. One of these regions contains BMPER, a gene involved with vascular remodeling, and that has been associated with hypoxia-induced angiogenic response. By genotyping samples across 2 altitudinal transects, we show that allele frequency changes at this locus are not gradual, but rather show a well-defined shift above ca. 1,900 m. Transcript and protein structure analyses on this gene suggest that putative selection likely acts on noncoding variation. These results underline how low oxygen pressure generates the most consistent selective constraint in high-altitude environments, to which vertebrates with vastly contrasting physiological profiles need to adapt in the context of ongoing climate change.

Whole-genome resequencing landscape of adaptive evolution in Relict gull (Larus relictus)

BACKGROUND

The relict gull (Larus relictus, Charadriiformes, Laridae) classified as vulnerable in the IUCN Red List is defined as a first-class national protected bird in China. However, our knowledge of the evolutionary history of L. relictus is limited. Here, we performed whole-genome resequencing of L. relictus (n = 14) and L. brunnicephalus (n = 3) to explore the genetic relationships and population structures and understand their adaptive evolution.

RESULTS

The whole genome resequencing generated 667.55 Gb clean reads with an average sequencing depth of ~ 29×. The genomic variant analysis identified 13,717,267 heterozygous SNPs in the samples obtained from 17 individuals. Population genetic diversity analysis revealed that low nucleotide diversity (0.00029) and no obvious population structure in L. relictus. Demographic history revealed that from 180 to 5 kya (thousand years ago), the effective population size (Ne) of L. relictus exhibited declines (24,000 to 5,000), with a very low range population size (2,200 to 5,000). In contrast, from 100 to 80 kya, L. brunnicephalus peaked in ancestral Ne, followed by distinct declines at ~ 70 kya (100,000 to 16,000). The findings identified several genes associated with the correlated changed life-history traits of L. relictus, including BMP4 involved in beak adaptation; HAND2, NEUROG1, COL11A2, and EDNRB involved in the evolution of the palate, soft palate, and tongue; PIGR and PLCB2 involved in an enhanced response to bitter taste by sensing chemical secretions released by staple food substrate insects to activate protective mechanisms. Furthermore, thirty-four genes related to sperm development and activity, including KLHL10 and TEKT3, were identified in the expanded gene family. In addition, MED1, CNOT9, NR5A1, and PATZ1 were involved in enhanced male hormone secretion and a high density of candidate genes associated with embryonic development were identified. The findings indicated that the L. relictus population was in a male-biased diffusion mode; the function of the TEKT3 gene showed that males played a dominant role in brooding, which enhanced their attraction to females. Our study revealed that significant enrichment of olfactory signaling pathway genes, including OR14C36, OR14J1, OR14I1, and OR14A16; inner ear development-related, including PTN, PTPN11, GATA2, ATP8B1, and MYO15A; and those related to hypoxic adaptation to high-altitude breeding and iris colour.

CONCLUSIONS

Based on the results and the knowledge of this organism biology and habitat use, we infer that less adaptive evolutionary pressure on vision in L. relictus were related with their feeding behaviour and adaptation. In summary, this comprehensive analysis provides insights into the evolutionary features of L. relictus and a new perspective for scientific research on L. relictus to effectively determine its future survival viability.

A chromosome-scale reference genome assembly for Triplophysa lixianensis

In this study, we constructed a chromosome-scale reference genome assembly for Lixian plateau loach, Triplophysa lixianensis, by integration of MGI short-read, PacBio HiFi long-read and Hi-C sequencing technologies. A 668-Mb haplotypic genome assembly was obtained for a female T. lixianensis, and 98.91% of the assembled sequences were anchored into 25 chromosomes. This assembly owned a moderate repeat content (35.63%) and an annotation of 23,774 protein-coding genes, among them 94.15% were predicted with functions. The assembled genome of T. lixianensis shared a good syntenic relationship with previously published data of its relative T. dalaica. Taken together, our genome data presented here provide a valuable genetic resource for in-depth evolutionary and functional studies, as well as molecular breeding and conservation of this valuable fish species to elevate its ecological and economical values.

Genomic Insights Into Genetic Basis of Evolutionary Conservatism and Innovation in Frogs

Examining closely related species evolving in similar environments offers valuable insights into the mechanisms driving phylogenetic conservatism and evolutionary lability. This can elucidate the intricate relationship between inheritance and environmental factors. Nonetheless, the precise genomic dynamics and molecular underpinnings of this process remain enigmatic. This study explores the evolutionary conservatism and adaptation exhibited by two closely related high-altitude frog species: Nanorana parkeri and N. pleskei. We assembled a high-quality genome for Tibetan N. pleskei and compared it to the genomes of N. parkeri and their lowland relatives. Our findings reveal that these two Tibetan frog species diverged approximately 16.6 million years ago, pointing to a possible ancestral colonization of high-elevation habitats. Following this colonization, significant adaptive evolution occurred in both coding and non-coding regions of the ancestral lineage. This evolution led to notable phenotypic alterations, as evidenced by the reduced body size. Also, due to purifying selection, most ancestral adaptive features persisted in descendant species, indicating a strong element of evolutionary conservatism. However, descendant species evolved novel adaptations to exacerbated environmental challenges in the Tibet Plateau, mainly related to hypoxia response. Furthermore, our analysis underscores the critical role of regulatory variations in descendant adaptive evolution. Notably, hub genes in networks, such as EGLN3, accumulated more variations in regulatory regions as they were transmitted from ancestors to descendants. In sum, our study sheds light on the profound and lasting impact of genetic heritage on species' adaptive evolution.

Selection signatures for high altitude adaptation in livestock: A review

High altitude adapted livestock species (cattle, yak, goat, sheep, and horse) has critical role in the human socioeconomic sphere and acts as good source of animal source products including milk, meat, and leather, among other things. These species sustain production and reproduction even in harsh environments on account of adaptation resulting from continued evolution of beneficial traits. Selection pressure leads to various adaptive strategies in livestock whose footprints are evident at the different genomic sites as the "Selection Signature". Scrutiny of these signatures provides us crucial insight into the evolutionary process and domestication of livestock adapted to diverse climatic conditions. These signatures have the potential to change the sphere of animal breeding and further usher the selection programmes in right direction. Technological revolution and recent strides made in genomic studies has opened the routes for the identification of selection signatures. Numerous statistical approaches and bioinformatics tools have been developed to detect the selection signature. Consequently, studies across years have identified candidate genes under selection region found associated with numerous traits which have a say in adaptation to high-altitude environment. This makes it pertinent to have a better understanding about the selection signature, the ways to identify and how to utilize them for betterment of livestock populations as well as farmers. This review takes a closer look into the general concept, various methodologies, and bioinformatics tools commonly employed in selection signature studies and summarize the results of recent selection signature studies related to high-altitude adaptation in various livestock species. This review will serve as an informative and useful insight for researchers and students in the field of animal breeding and evolutionary biology.

Genome-Wide Selection Signals Reveal Candidate Genes Associated with Plateau Adaptation in Tibetan Sheep

Tibetan sheep have developed unique adaptations for survival in the Qinghai-Tibet Plateau environment. However, the functional genes and molecular mechanisms that regulate hypoxia adaptation have not been fully characterized. In this study, based on the whole-genome resequencing data for Tibetan sheep at different altitudes, the population differentiation index (FST) and nucleotide diversity ratio (θπ ratio) were evaluated in populations of 20 Oula sheep (3501 m altitude, OL), 20 Zashijia sheep (4369 m altitude, ZSJ), and 20 Awang sheep (4643 m altitude, AW) to reveal candidate loci related to high-altitude hypoxia. We screened 728 and 524 candidate genes in the AW vs. OL and ZSJ vs. OL groups, respectively, of which 134 genes were jointly screened. Candidate genes were mainly enriched in the Ras, melanoma, melanogenesis, VEGF, and PPAR signaling pathways. HIF1AN, PDGFA, PDGFD, ANXA2, SOCS2, NOXA1, WNT7B, MMP14, GNG2, ATF6, PGAM2, PPP3R1, GSTCD, and PPARA may play important roles in the high-altitude adaptation of Tibetan sheep. In conclusion, this study provides valuable insights into the genes and molecular mechanisms that underlie high-altitude hypoxia adaptation in Tibetan sheep.

Trait Variation and Spatiotemporal Dynamics across Avian Secondary Contact Zones

A secondary contact zone (SCZ) is an area where incipient species or divergent populations may meet, mate, and hybridize. Due to the diverse patterns of interspecific hybridization, SCZs function as field labs for illuminating the on-going evolutionary processes of speciation and the establishment of reproductive isolation. Interspecific hybridization is widely present in avian populations, making them an ideal system for SCZ studies. This review exhaustively summarizes the variations in unique traits within avian SCZs (vocalization, plumage, beak, and migratory traits) and the various movement patterns of SCZs observed in previous publications. It also highlights several potential future research directions in the genomic era, such as the relationship between phenotypic and genomic differentiation in SCZs, the genomic basis of trait differentiation, SCZs shared by multiple species, and accurate predictive models for forecasting future movements under climate change and human disturbances. This review aims to provide a more comprehensive understanding of speciation processes and offers a theoretical foundation for species conservation.

Comparative skin histological and transcriptomic analysis of Rana kukunoris with two different skin colors

This study compares the skin structures of Rana kukunoris with two different skin colors living in the same area of Haibei in the Northeastern Qinghai-Tibet Plateau. The skin thickness of the khaki R. kukunoris was significantly greater than that of the brown R. kukunoris (P < 0.01), and significantly more mucous and granular glands were present on the dorsal skin of the khaki frog (P < 0.05). Meanwhile, the melanocytes on the dorsal skin of the brown frog were significantly larger than those on the khaki one (P < 0.05). Morphological changes in the expansion and aggregation of melanocytes seemed to deepen the skin color of R. kukunoris. Moreover, transcriptome sequencing identified tyrosine metabolism, melanogenesis, and riboflavin metabolism as the main pathways involved in melanin formation and metabolism in brown R. kukunoris. TYR, MC1R was upregulated as the skin color of R. kukunoris was deepened and contributed to melanin production and metabolism. In contrast, the khaki frog had significantly more upregulated genes and metabolic pathways related to autoimmunity. The khaki frog appeared to defend against ultraviolet (UV) radiation-induced damage by secreting mucus and small molecular peptides, whereas the brown frog protected itself by distributing a large amount of melanin. Hence, the different skin colors of R. kukunoris might represent different adaptation strategies for survival in the intense UV radiation environment of the Qinghai-Tibet Plateau.

Coping with alpine habitats: genomic insights into the adaptation strategies of Triplostegia glandulifera (Caprifoliaceae)

How plants find a way to thrive in alpine habitats remains largely unknown. Here we present a chromosome-level genome assembly for an alpine medicinal herb, Triplostegia glandulifera (Caprifoliaceae), and 13 transcriptomes from other species of Dipsacales. We detected a whole-genome duplication event in T. glandulifera that occurred prior to the diversification of Dipsacales. Preferential gene retention after whole-genome duplication was found to contribute to increasing cold-related genes in T. glandulifera. A series of genes putatively associated with alpine adaptation (e.g. CBFs, ERF-VIIs, and RAD51C) exhibited higher expression levels in T. glandulifera than in its low-elevation relative, Lonicera japonica. Comparative genomic analysis among five pairs of high- vs low-elevation species, including a comparison of T. glandulifera and L. japonica, indicated that the gene families related to disease resistance experienced a significantly convergent contraction in alpine plants compared with their lowland relatives. The reduction in gene repertory size was largely concentrated in clades of genes for pathogen recognition (e.g. CNLs, prRLPs, and XII RLKs), while the clades for signal transduction and development remained nearly unchanged. This finding reflects an energy-saving strategy for survival in hostile alpine areas, where there is a tradeoff with less challenge from pathogens and limited resources for growth. We also identified candidate genes for alpine adaptation (e.g. RAD1, DMC1, and MSH3) that were under convergent positive selection or that exhibited a convergent acceleration in evolutionary rate in the investigated alpine plants. Overall, our study provides novel insights into the high-elevation adaptation strategies of this and other alpine plants.

Comparative mitogenomic analysis of subterranean and surface amphipods (Crustacea, Amphipoda) with special reference to the family Crangonyctidae

Mitochondrial genomes play important roles in studying genome evolution, phylogenetic analyses, and species identification. Amphipods (Class Malacostraca, Order Amphipoda) are one of the most ecologically diverse crustacean groups occurring in a diverse array of aquatic and terrestrial environments globally, from freshwater streams and lakes to groundwater aquifers and the deep sea, but we have a limited understanding of how habitat influences the molecular evolution of mitochondrial energy metabolism. Subterranean amphipods likely experience different evolutionary pressures on energy management compared to surface-dwelling taxa that generally encounter higher levels of predation and energy resources and live in more variable environments. In this study, we compared the mitogenomes, including the 13 protein-coding genes involved in the oxidative phosphorylation (OXPHOS) pathway, of surface and subterranean amphipods to uncover potentially different molecular signals of energy metabolism between surface and subterranean environments in this diverse crustacean group. We compared base composition, codon usage, gene order rearrangement, conducted comparative mitogenomic and phylogenomic analyses, and examined evolutionary signals of 35 amphipod mitogenomes representing 13 families, with an emphasis on Crangonyctidae. Mitogenome size, AT content, GC-skew, gene order, uncommon start codons, location of putative control region (CR), length of rrnL and intergenic spacers differed between surface and subterranean amphipods. Among crangonyctid amphipods, the spring-dwelling Crangonyx forbesi exhibited a unique gene order, a long nad5 locus, longer rrnL and rrnS loci, and unconventional start codons. Evidence of directional selection was detected in several protein-encoding genes of the OXPHOS pathway in the mitogenomes of surface amphipods, while a signal of purifying selection was more prominent in subterranean species, which is consistent with the hypothesis that the mitogenome of surface-adapted species has evolved in response to a more energy demanding environment compared to subterranean amphipods. Overall, gene order, locations of non-coding regions, and base-substitution rates points to habitat as an important factor influencing the evolution of amphipod mitogenomes.

Metabolomics and transcriptomics of embryonic livers reveal hypoxia adaptation of Tibetan chickens

BACKGROUND

Exploring the hypoxia adaptation mechanism of Tibetan chicken is of great significance for revealing the survival law of Tibetan chicken and plateau animal husbandry production. To investigate the hypoxia adaptation of Tibetan chickens (TBCs), an integrative metabolomic-transcriptomic analysis of the liver on day 18 of embryonic development was performed. Dwarf laying chickens (DLCs), a lowland breed, were used as a control.

RESULTS

A total of 1,908 metabolites were identified in both TBCs and DLCs. Energy metabolism and amino acid metabolism related differentially regulated metabolites (DRMs) were significantly enriched under hypoxia. Important metabolic pathways including the TCA cycle and arginine and proline metabolism were screened; PCK1, SUCLA2, and CPS1 were found to be altered under hypoxic conditions. In addition, integrated analysis suggested potential differences in mitochondrial function, which may play a crucial role in the study of chicken oxygen adaptation.

CONCLUSIONS

These results suggest that hypoxia changed the gene expression and metabolic patterns of embryonic liver of TBCs compared to DLCs. Our study provides a basis for uncovering the molecular regulation mechanisms of hypoxia adaptation in TBCs with the potential application of hypoxia adaptation research for other animals living on the Qinghai-Tibet plateau, and may even contribute to the study of diseases caused by hypoxia.

Coping with extremes: Alternations in diet, gut microbiota, and hepatic metabolic functions in a highland passerine

In wild animals, diet and gut microbiota interactions are critical moderators of metabolic functions and are highly contingent on habitat conditions. Challenged by the extreme conditions of high-altitude environments, the strategies implemented by highland animals to adjust their diet and gut microbial composition and modulate their metabolic substrates remain largely unexplored. By employing a typical human commensal species, the Eurasian tree sparrow (Passer montanus, ETS), as a model species, we studied the differences in diet, digestive tract morphology and enzyme activity, gut microbiota, and metabolic energy profiling between highland (the Qinghai-Tibet Plateau, QTP; 3230 m) and lowland (Shijiazhuang, Hebei; 80 m) populations. Our results showed that highland ETSs had enlarged digestive organs and longer small intestinal villi, while no differences in key digestive enzyme activities were observed between the two populations. The 18S rRNA sequencing results revealed that the dietary composition of highland ETSs were more animal-based and less plant-based than those of the lowland ones. Furthermore, 16S rRNA sequencing results suggested that the intestinal microbial communities were structurally segregated between populations. PICRUSt metagenome predictions further indicated that the expression patterns of microbial genes involved in material and energy metabolism, immune system and infection, and xenobiotic biodegradation were strikingly different between the two populations. Analysis of liver metabolomics revealed significant metabolic differences between highland and lowland ETSs in terms of substrate utilization, as well as distinct sex-specific alterations in glycerophospholipids. Furthermore, the interplay between diet, liver metabolism, and gut microbiota suggests a dietary shift resulting in corresponding changes in gut microbiota and metabolic functions. Our findings indicate that highland ETSs have evolved to optimize digestion and absorption, rely on more protein-rich foods, and possess gut microbiota tailored to their dietary composition, likely adaptive physiological and ecological strategies adopted to cope with extreme highland environments.

The draft genome of the Temminck's tragopan (Tragopan temminckii) with evolutionary implications

BACKGROUND

High-quality genome data of birds play a significant role in the systematic study of their origin and adaptive evolution. The Temminck's tragopan (Tragopan temminckii) (Galliformes, Phasianidae), a larger pheasant, is one of the most abundant and widely distributed species of the genus Tragopan, and was defined as class II of the list of national key protected wild animals in China. The absence of a sequenced genome has restricted previous evolutionary trait studies of this taxa.

RESULTS

The whole genome of the Temminck's tragopan was sequenced using Illumina and PacBio platform, and then de novo assembled and annotated. The genome size was 1.06 Gb, with a contig N50 of 4.17 Mb. A total of 117.22 Mb (11.00%) repeat sequences were identified. 16,414 genes were predicted using three methods, with 16,099 (98.08%) annotated as functional genes based on five databases. In addition, comparative genome analyses were conducted across 12 Galliformes species. The results indicated that T. temminckii was the first species to branch off from the clade containing Lophura nycthemera, Phasianus colchicus, Chrysolophus pictus, Syrmaticus mikado, Perdix hodgsoniae, and Meleagris gallopavo, with a corresponding divergence time of 31.43 million years ago (MYA). Expanded gene families associated with immune response and energy metabolism were identified. Genes and pathways associated with plumage color and feather development, immune response, and energy metabolism were found in the list of positively selected genes (PSGs).

CONCLUSIONS

A genome draft of the Temminck's tragopan was reported, genome feature and comparative genome analysis were described, and genes and pathways related to plumage color and feather development, immune response, and energy metabolism were identified. The genomic data of the Temminck's tragopan considerably contribute to the genome evolution and phylogeny of the genus Tragopan and the whole Galliformes species underlying ecological adaptation strategies.

Gene family expansion analysis and identification of the histone family in Spodoptera frugiperda

Spodoptera frugiperda, a major invasive pest, causes severe damage to various economically important crops. Previous comparative genomics studies have revealed a close association between the invasiveness of S. frugiperda and its genome. In recent years, a vast amount of genome from lepidopteran species has become available, offering an opportunity for a more detailed and comprehensive understanding of the biological characteristics of S. frugiperda. In this study, we conducted a comprehensive comparative genomics analysis of S. frugiperda using genome from 46 lepidopteran species. We found the highest number of gene family expansion events in S. frugiperda, indicating that gene family expansion is a crucial mechanism in its adaptive evolution. The expanded gene families are enriched in various biological processes, including nutrient metabolism, development, stress response, reproduction, and immune processes, suggesting that the expansion of these gene families likely contributes to the strong environmental adaptability of S. frugiperda. Furthermore, we identified the expansion of histone gene families in S. frugiperda which resulted from chromosome segmental duplications after the divergence from closely related species. Expression analysis of histone genes indicated that certain members might exert an influence on the growth and reproduction processes of S. frugiperda. Overall, our study deepens our understanding of the biological characteristics of S. frugiperda, providing a theoretical basis for the comprehensive management and sustained control of S. frugiperda and other lepidopteran pests in the future.

Comparative analyses reveal potential genetic mechanisms for high-altitude adaptation of Schizopygopsis fishes based on chromosome-level genomes

The schizothoracine fishes, widely distributed in the Qinghai-Tibetan Plateau and its adjacent areas, are considered as ideal models for investigation of high-altitude adaptation. Schizophygopsis are one group of the highly specialized schizothoracine fishes, and the genetic basis for their high-altitude adaptation is poorly understood. In this study, we performed comparative genomics analyses to investigate the potential genetic mechanisms for high-altitude adaptation of Schizopygopsis malacanthus and Schizopygopsis pylzovi based on the chromosome-level genomes. Functional enrichment analysis revealed that many expanded gene families in Schizopygopsis were associated with immune response while many contracted gene families were functionally associated with olfaction. Among the 123 positively selected genes (PSGs), angpt2a was detected in HIF-1 signaling pathway and possibly related to the hypoxia adaptation of Schizopygopsis. Furthermore, two PSGs cox15 and ndufb10 were distributed in thermogenesis, and there was a Schizopygopsis-specific missense mutation in cox15 (Gln115Glu), which possibly contributed to the cold temperature adaptation of the Schizopygopsis. Kyoto Encyclopedia of Genes and Genomes enrichment of the PSGs revealed three significant pathways including metabolic pathways, cell cycle, and homologous recombination and Gene Ontology enrichment analysis of the PSGs revealed several categories associated with DNA repair, cellular response to DNA damage stimulus, and metabolic process. Chromosome-scale characterization of olfactory receptor (OR) repertoires indicated that Schizopygopsis had the least number of OR genes, and the OR gene contraction was possibly caused by the limited food variety and the environmental factors such as lower air pressure, lower humidity, and lower temperature. Our study will help expand our understanding of the potential adaptive mechanism of Schizopygopsis to cope with the high-altitude conditions.

Genome and transcriptome of Ips nitidus provide insights into high-altitude hypoxia adaptation and symbiosis

Ips nitidus is a well-known conifer pest that has contributed significantly to spruce forest disturbance in the Qinghai-Tibet Plateau and seriously threatens the ecological balance of these areas. We report a chromosome-level genome of I. nitidus determined by PacBio and Hi-C technology. Phylogenetic inference showed that it diverged from the common ancestor of I. typographus ∼2.27 mya. Gene family expansion in I. nitidus was characterized by DNA damage repair and energy metabolism, which may facilitate adaptation to high-altitude hypoxia. Interestingly, differential gene expression analysis revealed upregulated genes associated with high-altitude hypoxia adaptation and downregulated genes associated with detoxification after feeding and tunneling in fungal symbiont Ophiostoma bicolor-colonized substrates. Our findings provide evidence of the potential adaptability of I. nitidus to conifer host, high-altitude hypoxia and insight into how fungal symbiont assist in this process. This study enhances our understanding of insect adaptation, symbiosis, and pest management.

Genomic analysis of Tibetan ground tits identifies molecular adaptations associated with cooperative breeding

Cooperative breeding is a sophisticated altruistic social behavior that helps social animals to adapt to harsh environments. The Tibetan ground tit, Pseudopodoces humilis, is a high-altitude bird endemic to the Tibetan plateau. Recently, it has become an exciting system for studying the evolution of facultative cooperative breeding. To test for molecular adaptations associated with cooperative breeding, we resequenced the whole genome of ground tits from 6 wild populations that display remarkable variation in the frequency of cooperative breeding. Population structure analyses showed that the 6 populations were divided into 4 lineages, which is congruent with the major geographical distribution of the sampling sites. Using genome-wide selective sweep analysis, we identified putative positively selected genes (PSGs) in groups of tits that displayed high and low cooperative breeding rates. The total number of PSGs varied from 146 to 722 in high cooperative breeding rate populations, and from 272 to 752 in low cooperative breeding rate populations. Functional enrichment analysis of these PSGs identified several significantly enriched ontologies related to oxytocin signaling, estrogen signaling, and insulin secretion. PSGs involved in these functional ontologies suggest that molecular adaptations in hormonal regulation may have played important roles in shaping the evolution of cooperative breeding in the ground tit. Taken together, our study provides candidate genes and functional ontologies involved in molecular adaptations associated with cooperative breeding in Tibetan ground tits, and calls for a better understanding of the genetic roles in the evolution of cooperative breeding.

Genomic Basis of Adaptive Divergence in Leg Length between Ground- and Tree-Dwelling Species within a Bird Family

Hind limbs of tetrapods vary greatly in length and the variability can be associated with locomotor adaptation. Although the phenotypic evolution has been well documented, the underlying genetic basis remains poorly understood. We address this issue by integrating comparative genomics and functional prediction with a study system consisting of ground-dwelling, long-legged and tree-dwelling, short-legged species within the avian family Paridae. Genome-wide divergence and phenotypic correlation analyses jointly identified five highly divergent genomic regions that are significantly related with the difference in leg length between these two groups. Gene annotation for these regions detected three genes involved in skeletal development, that is, PTPA, BRINP1, and MIGA2, with the first one being under the strongest selection. Furthermore, four single nucleotide polymorphisms (SNPs) in the coding region of PTPA can well distinguish the two groups with distinct leg length. Among the four SNPs, one is non-synonymous mutation, and according to the prediction for protein structure and function, it can affect the 3D structure of the encoded protein by altering the corresponding amino acid's position. The alleles of PTPA were found in all sequenced species of the orders Palaeognathae and Psittaciformes, which typically take a ground locomotion style. A whole-genome scanning across bird species uncovered that the four SNPs are more likely to be present in resident passerines with increased leg length/wing length ratios (a proxy of leg-dependent locomotion efficiency). Our findings provide insight into the molecular evolution of locomotion performance based on leg morphology in birds.

Autumn migration of black-necked crane (Grus nigricollis) on the Qinghai-Tibetan and Yunnan-Guizhou plateaus

Despite previous research efforts, the majority migration routes of the black-necked cranes (Grus nigricollis) have remained veiled. In this study, we utilized satellite telemetry data from 45 cranes between 2015 and 2021 to unveil critical insights. Our results revealed 11 distinct autumn migration routes and one sedentary flock, of which eight routes and the sedentary flock were previously undocumented. Our findings highlighted the remarkable diversity in the migration routes of black-necked cranes, especially in terms of migration orientations, spatial-temporal patterns, and altitudinal movement patterns. Cranes breeding on the eastern, northern, and central Qinghai-Tibetan Plateau migrated southward, while those on the northern slopes of the Himalayas migrated eastward, westward, northward, or opted to remain sedentary. Moreover, we expanded the known range of migration distances to 84-1520 km at both ends (excluding sedentary individuals) and identified two long-term (Da Qaidam and Chaka) and one short-term (Gyatong grassland) stopover sites. Furthermore, our study revealed that the breeding colonies in the Qilian Mountains on the northeastern Qinghai-Tibetan Plateau utilized long-term stopover sites before embarking on significant altitude ascent, while other flocks displayed more urgent migration patterns, preferring to roost only at night. By unveiling the near-complete autumn migration routes of black-necked cranes, our research has contributed to discovering the critical habitats and connectivity among various breeding colonies, which is instrumental in developing effective seasonal conservation plans.

Avian malaria parasite infections do not affect personality in the chestnut thrush (Turdus rubrocanus) on the Qinghai-Tibet Plateau

Personality traits, the consistent individual behavioral differences, are currently gaining much attention in studies of natural bird populations. However, associations between personality traits and parasite infections are not often investigated. Even less attention has been given to studies of birds in the high-elevation region such as the Tibetan plateau. This research aims to examine the relationship between avian malaria parasites and two personality traits in a population of the Chestnut Thrush (Turdus rubrocanus) breed in the Tibetan plateau. Our results revealed no evidence of sex bias in malaria parasite prevalence. Furthermore, we found no effect of infection status on two personality scores: activity and boldness. Additionally, no effects on the activity level or boldness were observed for different parasite lineages of Haemoproteus, Leucocytozoon, the sex of the birds, or their interactions. Similarly, we did not find any relationship between activity level and boldness with nestling numbers, sex, or their interactions. Notably, individuals with a larger number of offspring tended to display greater boldness. Our findings indicate that blood parasite infections are common in this population but do not significantly impact the personality of the birds.

Population genomic analysis reveals distinct demographics and recent adaptation in the black flying fox (Pteropus alecto)

As the only mammalian group capable of powered flight, bats have many unique biological traits. Previous comparative genomic studies in bats have focused on long-term evolution. However, the micro-evolutionary processes driving recent evolution are largely under-explored. Using resequencing data from 50 black flying foxes (Pteropus alecto), one of the model species for bats, we find that black flying fox has much higher genetic diversity and lower levels of linkage disequilibrium than most of the mammalian species. Demographic inference reveals strong population fluctuations (>100 fold) coinciding with multiple historical events including the last glacial change and Toba super eruption, suggesting that the black flying fox is a very resilient species with strong recovery abilities. While long-term adaptation in the black flying fox is enriched in metabolic genes, recent adaptation in the black flying fox has a unique landscape where recently selected genes are not strongly enriched in any functional category. The demographic history and mode of adaptation suggest that black flying fox might be a well-adapted species with strong evolutionary resilience. Taken together, this study unravels a vibrant landscape of recent evolution for the black flying fox and sheds light on several unique evolutionary processes for bats comparing to other mammalian groups.

Drivers of plateau adaptability in cashmere goats revealed by genomic and transcriptomic analyses

BACKGROUND

The adaptive evolution of plateau indigenous animals is a current research focus. However, phenotypic adaptation is complex and may involve the interactions between multiple genes or pathways, many of which remain unclear. As a kind of livestock with important economic value, cashmere goat has a high ability of plateau adaptation, which provides us with good materials for studying the molecular regulation mechanism of animal plateau adaptation.

RESULTS

In this study, 32 Jiangnan (J) and 32 Tibetan (T) cashmere goats were sequenced at an average of 10. Phylogenetic, population structure, and linkage disequilibrium analyses showed that natural selection or domestication has resulted in obvious differences in genome structure between the two breeds. Subsequently, 553 J vs. T and 608 T vs. J potential selected genes (PSGs) were screened. These PSGs showed potential relationships with various phenotypes, including myocardial development and activity (LOC106502520, ATP2A2, LOC102181869, LOC106502520, MYL2, ISL1, and LOC102181869 genes), pigmentation (MITF and KITLG genes), hair follicles/hair growth (YAP1, POGLUT1, AAK1, HES1, WNT1, PRKAA1, TNKS, WNT5A, VAX2, RSPO4, CSNK1G1, PHLPP2, CHRM2, PDGFRB, PRKAA1, MAP2K1, IRS1, LPAR1, PTEN, PRLR, IBSP, CCNE2, CHAD, ITGB7, TEK, JAK2, and FGF21 genes), and carcinogenesis (UBE2R2, PIGU, DIABLO, NOL4L, STK3, MAP4, ADGRG1, CDC25A, DSG3, LEPR, PRKAA1, IKBKB, and ABCG2 genes). Phenotypic analysis showed that Tibetan cashmere goats has finer cashmere than Jiangnan cashmere goats, which may allow cashmere goats to better adapt to the cold environment in the Tibetan plateau. Meanwhile, KRTs and KAPs expression in Jiangnan cashmere goat skin was significantly lower than in Tibetan cashmere goat.

CONCLUSIONS

The mutations in these PSGs maybe closely related to the plateau adaptation ability of cashmere goats. In addition, the expression differences of KRTs and KAPs may directly determine phenotypic differences in cashmere fineness between the two breeds. In conclusion, this study provide a reference for further studying plateau adaptive mechanism in animals and goat breeding.

Divergent contributions of coding and noncoding sequences to initial high-altitude adaptation in passerine birds endemic to the Qinghai-Tibet Plateau

Early events in the evolution of an ancestral lineage can shape the adaptive patterns of descendant species, but the evolutionary mechanisms driving initial adaptation from an ancestor remain largely unexplored. High-altitude adaptations have been extensively explored from the viewpoint of protein-coding genes; however, the contribution of noncoding regions remains relatively neglected. Here, we integrate genomic and transcriptomic data to investigate adaptive evolution in the ancestor of three high-altitude snowfinch species endemic to the Qinghai-Tibet Plateau. Our genome-wide scan for adaptation in the snowfinch ancestor identifies strong adaptation signals in functions of development and metabolism for the coding genes, but in functions of the nervous system development for noncoding regions. This pattern is exclusive to the snowfinch ancestor compared to a control ancestral lineage subject to weak selection. Changes in noncoding regions in the snowfinch ancestor, especially those nearest to coding genes, may be disproportionately associated with the differential expression of genes in the brain tissue compared to other tissues. Extensive gene expression in the brain tissue can be further altered via genetic regulatory networks of transcription factors harbouring potential accelerated regulatory regions (e.g., the development-related transcription factor YEATS4). Altogether, our study provides new evidence concerning how coding and noncoding sequences work through decoupled pathways in initial adaptation to the selective pressure of high-altitude environments. The analysis highlights the idea that noncoding sequences may be promising elements in facilitating the rapid evolution and adaptation to high altitudes.

Convergent molecular evolution of thermogenesis and circadian rhythm in Arctic ruminants

The muskox and reindeer are the only ruminants that have evolved to survive in harsh Arctic environments. However, the genetic basis of this Arctic adaptation remains largely unclear. Here, we compared a de novo assembled muskox genome with reindeer and other ruminant genomes to identify convergent amino acid substitutions, rapidly evolving genes and positively selected genes among the two Arctic ruminants. We found these candidate genes were mainly involved in brown adipose tissue (BAT) thermogenesis and circadian rhythm. Furthermore, by integrating transcriptomic data from goat adipose tissues (white and brown), we demonstrated that muskox and reindeer may have evolved modulating mitochondrion, lipid metabolism and angiogenesis pathways to enhance BAT thermogenesis. In addition, results from co-immunoprecipitation experiments prove that convergent amino acid substitution of the angiogenesis-related gene hypoxia-inducible factor 2alpha (HIF2A), resulting in weakening of its interaction with prolyl hydroxylase domain-containing protein 2 (PHD2), may increase angiogenesis of BAT. Altogether, our work provides new insights into the molecular mechanisms involved in Arctic adaptation.

Gut microbiome signatures of extreme environment adaption in Tibetan pig

Tibetan pigs (TPs) can adapt to the extreme environments in the Tibetan plateau implicated by their self-genome signals, but little is known about roles of the gut microbiota in the host adaption. Here, we reconstructed 8210 metagenome-assembled genomes from TPs (n = 65) living in high-altitude and low-altitude captive pigs (87 from China-CPs and 200 from Europe-EPs) that were clustered into 1050 species-level genome bins (SGBs) at the threshold of 95% average nucleotide identity. 73.47% of SGBs represented new species. The gut microbial community structure analysis based on 1,048 SGBs showed that TPs was significantly different from low-altitude captive pigs. TP-associated SGBs enabled to digest multiple complex polysaccharides, including cellulose, hemicellulose, chitin and pectin. Especially, we found TPs showed the most common enrichment of phyla Fibrobacterota and Elusimicrobia, which were involved in the productions of short- and medium-chain fatty acids (acetic acid, butanoate and propanoate; octanomic, decanoic and dodecanoic acids), as well as in the biosynthesis of lactate, 20 essential amino acids, multiple B vitamins (B1, B2, B3, B5, B7 and B9) and cofactors. Unexpectedly, Fibrobacterota solely showed powerful metabolic capacity, including the synthesis of acetic acid, alanine, histidine, arginine, tryptophan, serine, threonine, valine, B2, B5, B9, heme and tetrahydrofolate. These metabolites might contribute to host adaptation to high-altitude, such as energy harvesting and resistance against hypoxia and ultraviolet radiation. This study provides insights into understanding the role of gut microbiome played in mammalian high-altitude adaptation and discovers some potential microbes as probiotics for improving animal health.

Weak gene-gene interaction facilitates the evolution of gene expression plasticity

BACKGROUND

Individual organisms may exhibit phenotypic plasticity when they acclimate to different conditions. Such plastic responses may facilitate or constrain the adaptation of their descendant populations to new environments, complicating their evolutionary trajectories beyond the genetic blueprint. Intriguingly, phenotypic plasticity itself can evolve in terms of its direction and magnitude during adaptation. However, we know little about what determines the evolution of phenotypic plasticity, including gene expression plasticity. Recent laboratory-based studies suggest dominance of reversing gene expression plasticity-plastic responses that move the levels of gene expression away from the new optima. Nevertheless, evidence from natural populations is still limited.

RESULTS

Here, we studied gene expression plasticity and its evolution in the montane and lowland populations of an elevationally widespread songbird-the Rufous-capped Babbler (Cyanoderma ruficeps)-with reciprocal transplant experiments and transcriptomic analyses; we set common gardens at altitudes close to these populations' native ranges. We confirmed the prevalence of reversing plasticity in genes associated with altitudinal adaptation. Interestingly, we found a positive relationship between magnitude and degree of evolution in gene expression plasticity, which was pertinent to not only adaptation-associated genes but also the whole transcriptomes from multiple tissues. Furthermore, we revealed that genes with weaker expressional interactions with other genes tended to exhibit stronger plasticity and higher degree of plasticity evolution, which explains the positive magnitude-evolution relationship.

CONCLUSIONS

Our experimental evidence demonstrates that species may initiate their adaptation to new habitats with genes exhibiting strong expression plasticity. We also highlight the role of expression interdependence among genes in regulating the magnitude and evolution of expression plasticity. This study illuminates how the evolution of phenotypic plasticity in gene expression facilitates the adaptation of species to challenging environments in nature.

A de novo assembled genome of the Tibetan Partridge (Perdix hodgsoniae) and its high-altitude adaptation

The Tibetan Partridge (Perdix hodgsoniae) is an endemic species distributed in high-altitude areas of 3600-5600 m on the Qinghai-Tibet Plateau. To explore how the species is adapted to the high elevation environment, we assembled a draft genome based on both the Illumina and PacBio sequencing platforms with its population genetics and genomics analysis. In total, 134.74 Gb short reads and 30.81 Gb long reads raw data were generated. The 1.05-Gb assembled genome had a contig N50 of 4.56 Mb, with 91.94% complete BUSCOs. The 17 457 genes were annotated, and 11.35% of the genome was composed of repeat sequences. The phylogenetic tree showed that P. hodgsoniae was located at the basal position of the clade, including Golden Pheasant (Chrysolophus pictus), Common Pheasant (Phasianus colchicus), and Mikado Pheasant (Syrmaticus mikado). We found that 1014, 2595, and 2732 of the 6641 one-to-one orthologous genes were under positive selection in P. hodgsoniae, detected using PAML, BUSTED, and aBSREL programs, respectively, of which 965 genes were common under positive selection with 3 different programs. Several positively selected genes and immunity pathways relevant to high-altitude adaptation were detected. Gene family evolution showed that 99 gene families experienced significant expansion events, while 6 gene families were under contraction. The total number of olfactory receptor genes was relatively low in P. hodgsoniae. Genomic data provide an important resource for a further study on the evolutionary history of P. hodgsoniae, which provides a new insight into its high-altitude adaptation mechanisms.

Enlarged fins of Tibetan catfish provide new evidence of adaptation to high plateau

The uplift of the Tibetan Plateau significantly altered the geomorphology and climate of the Euroasia by creating large mountains and rivers. Fishes are more likely to be affected relative to other organisms, as they are largely restricted to river systems. Faced with the rapidly flowing water in the Tibetan Plateau, a group of catfish has evolved greatly enlarged pectoral fins with more numbers of fin-rays to form an adhesive apparatus. However, the genetic basis of these adaptations in Tibetan catfishes remains elusive. In this study, we performed comparative genomic analyses based on the chromosome-level genome of Glyptosternum maculatum in family Sisoridae and detected some proteins with conspicuously high evolutionary rates in particular in genes involved in skeleton development, energy metabolism, and hypoxia response. We found that the hoxd12a gene evolved faster and a loss-of-function assay of hoxd12a supports a potential role for this gene in shaping the enlarged fins of these Tibetan catfishes. Other genes with amino acid replacements and signatures of positive selection included proteins involved in low temperature (TRMU) and hypoxia (VHL) responses. Functional assays reveal that the G. maculatumTRMU allele generates more mitochondrial ATP than the ancestral allele found in low-altitude fishes. Functional assays of VHL alleles suggest that the G. maculatum allele has lower transactivation activity than the low-altitude forms. These findings provide a window into the genomic underpinnings of physiological adaptations that permit G. maculatum to survive in the harsh environment of the Tibetan Himalayas that mirror those that are convergently found in other vertebrates such as humans.

Genomic architecture underlying morphological and physiological adaptation to high elevation in a songbird

Organisms often acquire physiological and morphological modifications to conquer ecological challenges when colonizing new environments which lead to their adaptive evolution. However, deciphering the genomic mechanism of ecological adaptation is difficult because ecological environments are often too complex for straightforward interpretation. Thus, we examined the adaptation of a widespread songbird-the rufous-capped babbler (Cyanoderma ruficeps)-to a relatively simple system: distinct environments across elevational gradients on the mountainous island of Taiwan. We focused on the genomic sequences of 43 birds from five populations to show that the Taiwan group split from its sister group in mainland China around 1-2 million years ago (Ma) and colonized the montane habitats of Taiwan at least twice around 0.03-0.22 Ma. The montane and lowland Taiwan populations diverged with gene flow between them, suggesting strong selection associated with different elevations. We found that the montane babblers had smaller beaks than the lowland ones, consistent with Allen's rule, and identified candidate genes-COL9A1 and SOX11-underlying the beak size changes. We also found that altitudinally divergent mutations were mostly located in noncoding regions and tended to accumulate in chromosomal inversions and autosomes. The altitudinally divergent mutations might regulate genes related to haematopoietic, metabolic, immune, auditory and vision functions, as well as cerebrum morphology and plumage development. The results reveal the genomic bases of morphological and physiological adaptation in this species to the low temperature, hypoxia and high UV light environment at high elevation. These findings improve our understanding of how ecological adaptation drives population divergence from the perspective of genomic architecture.

The metabolic adaptation in wild vertebrates via omics approaches

Metabolism is the basis for sustaining life and essential to the adaptive evolution of organisms. With the development of high-throughput sequencing technology, genetic mechanisms of adaptive evolution, including metabolic adaptation, have been extensively resolved by omics approaches, but a deep understanding of genetic and epigenetic metabolic adaptation is still lacking. Exploring metabolic adaptations from genetic and epigenetic perspectives in wild vertebrates is vital to understanding species evolution, especially for the early stages of adaptative evolution. Herein, we summarize the advances in our understanding of metabolic adaptations via omics approaches in wild vertebrates based on three types of cases: extreme environment, periodically changing environment, and changes of species characteristics. We conclude that the understanding of the formation of metabolic adaptations at the genetic level alone can well identify the adaptive genetic variation that has developed during evolution, but cannot resolve the potential impact of metabolic adaptations on the adaptative evolution in the future. Thus, it seems imperative to include epigenomics and metabolomics in the study of adaptation, and that in the future genomic and epigenetic data should be integrated to understand the formation of metabolic adaptation of wild vertebrate organisms.

Transcriptome analysis of pika heart tissue reveals mechanisms underlying the adaptation of a keystone species on the roof of the world

High-altitude environments impose intense stresses on living organisms and drive striking phenotypic and genetic adaptations, such as hypoxia resistance, cold tolerance, and increases in metabolic capacity and body mass. As one of the most successful and dominant mammals on the Qinghai-Tibetan Plateau (QHTP), the plateau pika (Ochotona curzoniae) has adapted to the extreme environments of the highest altitudes of this region and exhibits tolerance to cold and hypoxia, in contrast to closely related species that inhabit the peripheral alpine bush or forests. To explore the potential genetic mechanisms underlying the adaptation of O. curzoniae to a high-altitude environment, we sequenced the heart tissue transcriptomes of adult plateau pikas (comparing specimens from sites at two different altitudes) and Gansu pikas (O. cansus). Differential expression analysis and weighted gene co-expression network analysis (WGCNA) were used to identify differentially expressed genes (DEGs) and their primary functions. Key genes and pathways related to high-altitude adaptation were identified. In addition to the biological processes of signal transduction, energy metabolism and material transport, the identified plateau pika genes were mainly enriched in biological pathways such as the negative regulation of smooth muscle cell proliferation, the apoptosis signalling pathway, the cellular response to DNA damage stimulus, and ossification involved in bone maturation and heart development. Our results showed that the plateau pika has adapted to the extreme environments of the QHTP via protection against cardiomyopathy, tissue structure alterations and improvements in the blood circulation system and energy metabolism. These adaptations shed light on how pikas thrive on the roof of the world.

Latitudinal but not elevational variation in blood glucose level is linked to life history across passerine birds

Macrophysiological research is vital to our understanding of mechanisms underpinning global life history variation and adaptation to diverse environments. Here, we examined latitudinal and elevational variation in a key substrate of energy metabolism and an emerging physiological component of pace-of-life syndromes, blood glucose concentration. Our data, collected from 61 European temperate and 99 Afrotropical passerine species, revealed that baseline blood glucose increases with both latitude and elevation, whereas blood glucose stress response shows divergent directions, being stronger at low latitudes and high elevations. Low baseline glucose in tropical birds, compared to their temperate counterparts, was mainly explained by their low fecundity, consistent with the slow pace-of-life syndrome in the tropics. In contrast, elevational variation in this trait was decoupled from fecundity, implying a unique montane pace-of-life syndrome combining slow-paced life histories with fast-paced physiology. The observed patterns suggest that pace-of-life syndromes do not evolve along the single fast-slow axis.

Chromosome-level Genome Assembly of the High-altitude Leopard (Panthera pardus) Sheds Light on Its Environmental Adaptation

The leopard (Panthera pardus) has the largest natural distribution from low- to high-altitude areas of any wild felid species, but recent studies have revealed that leopards have disappeared from large areas, probably owing to poaching, a decline of prey species, and habitat degradation. Here, we reported the chromosome-scale genome assembly of the high-altitude leopard (HL) based on nanopore sequencing and high-throughput chromatin conformation capture (Hi-C) technology. Panthera genomes revealed similar repeat composition, and there was an appreciably conserved synteny between HL and the other two Panthera genomes. Divergence time analysis based on the whole genomes revealed that the HL and the low-altitude leopard differentiate from a common ancestor ∼2.2 Ma. Through comparative genomics analyses, we found molecular genetic signatures that may reflect high-altitude adaptation of the HL. Three HL-specific missense mutations were detected in two positively selected genes, that is, ITGA7 (Ala112Gly, Asp113Val, and Gln115Pro) and NOTCH2 (Ala2398Ser), which are likely to be associated with hypoxia adaptation. The chromosome-level genome of the HL provides valuable resources for the investigation of high-altitude adaptation and protection management of the vulnerable leopard.

Chromosome-level assembly of Gymnocypris eckloni genome

Gymnocypris eckloni is widely distributed in isolated lakes and the upper reaches of the Yellow River and play significant roles in the trophic web of freshwater communities. In this study, we generated a chromosome-level genome of G. eckloni using PacBio, Illumina and Hi-C sequencing data. The genome consists of 23 pseudo-chromosomes that contain 918.68 Mb of sequence, with a scaffold N50 length of 43.54 Mb. In total, 23,157 genes were annotated, representing 94.80% of the total predicted protein-coding genes. The phylogenetic analysis showed that G. eckloni was most closely related to C. carpio with an estimated divergence time of ~34.8 million years ago. For G. eckloni, we identified a high-quality genome at the chromosome level. This genome will serve as a valuable genomic resource for future research on the evolution and ecology of the schizothoracine fish in the Qinghai-Tibetan Plateau.

Temperature acclimation in hot-spring snakes and the convergence of cold response

Animals have evolved sophisticated temperature-sensing systems and mechanisms to detect and respond to ambient temperature changes. As a relict species endemic to the Qinghai-Tibet Plateau, hot-spring snake (Thermophis baileyi) survived the dramatic changes in climate that occurred during plateau uplift and ice ages, providing an excellent opportunity to explore the evolution of temperature sensation in ectotherms. Based on distributional information and behavioral experiments, we found that T. baileyi prefer hot-spring habitats and respond more quickly to warmth than other two snakes, suggesting that T. baileyi may evolve an efficient thermal-sensing system. Using high-quality chromosome-level assembly and comparative genomic analysis, we identified cold acclimation genes experiencing convergent acceleration in high-altitude lineages. We also discovered significant evolutionary changes in thermosensation- and thermoregulation-related genes, including the transient receptor potential (TRP) channels. Among these genes, TRPA1 exhibited three species-specific amino acid replacements, which differed from those found in infrared imaging snakes, implying different temperature-sensing molecular strategies. Based on laser-heating experiments, the T. baileyi-specific mutations in TRPA1 resulted in an increase in heat-induced opening probability and thermal sensitivity of the ion channels under the same degree of temperature stimulation, which may help the organism respond to temperature changes more quickly. These results provide insight into the genetic mechanisms underpinning the evolution of temperature-sensing strategies in ectotherms as well as genetic evidence of temperature acclimation in this group.

•

Hot-spring snakes prefer hot-spring habitats on the Qinghai-Tibet Plateau

•

Genetic variation in the snakes contribute to the temperature acclimation

•

Unique mutations in TRPA1 increase thermal sensitivity of the ion channel

•

Different temperature-sensing strategies existed across snakes

Genetic variation in HIF-2α attenuates ventilatory sensitivity and carotid body growth in chronic hypoxia in high-altitude deer mice

KEY POINTS

High-altitude natives of many species have experienced natural selection on the gene encoding HIF-2α, Epas1, including high-altitude populations of deer mice. HIF-2α regulates ventilation and carotid body growth in hypoxia, so the genetic variants in Epas1 in high-altitude natives may underlie evolved changes in control of breathing. Deer mice from controlled crosses between high- and low-altitude populations were used to examine the effects of Epas1 genotype on an admixed genomic background. The high-altitude variant was associated with reduced ventilatory chemosensitivity and carotid body growth in chronic hypoxia, but had no effects on haematology. The results help us better understand the genetic basis for the unique physiological phenotype of high-altitude natives.

ABSTRACT

The gene encoding HIF-2α, Epas1, has experienced a history of natural selection in many high-altitude taxa, but the functional role of mutations in this gene are still poorly understood. We investigated the influence of the high-altitude variant of Epas1 in North American deer mice (Peromyscus maniculatus) on control of breathing and carotid body growth during chronic hypoxia. We created hybrids between high- and low-altitude populations of deer mice to disrupt linkages between genetic loci so physiological effects of Epas1 alleles (Epas1^H and Epas1^L , respectively) could be examined on an admixed genomic background. In general, chronic hypoxia (4 weeks at 12 kPa O₂ ) enhanced ventilatory chemosensitivity (assessed as the acute ventilatory response to hypoxia), increased total ventilation and arterial O₂ saturation during progressive poikilocapnic hypoxia, and increased haematocrit and blood haemoglobin content across genotypes. However, effects of chronic hypoxia on ventilatory chemosensitivity were attenuated in mice that were homozygous for the high-altitude Epas1 allele (Epas1^H/H ). Carotid body growth and glomus cell hyperplasia, which was strongly induced in Epas1^L/L mice in chronic hypoxia, was not observed in Epas1^H/H mice. Epas1 genotype also modulated the effects of chronic hypoxia on metabolism and body temperature depression in hypoxia, but had no effects on haematological traits. These findings confirm the important role of HIF-2α in modulating ventilatory sensitivity and carotid body growth in chronic hypoxia, and show that genetic variation in Epas1 is responsible for evolved changes in the control of breathing and metabolism in high-altitude deer mice. Abstract figure legend ventilation and carotid body growth in hypoxia, so we investigated the role genetic variants in Epas1 in highaltitude deer mice on the control of breathing. In the lab, hybrids between high- and lowaltitude populations of deer mice were created to disrupt linkages between genetic loci so physiological effects of Epas1 alleles (Epas1H and Epas1L, respectively) could be examined on an admixed genomic background. The high-altitude variant was associated with reduced ventilatory chemosensitivity and carotid body growth after 4 weeks of chronic hypoxia, compared to mice homozygous for the low-altitude allele (Epas1LL). These results help us better understand the genetic basis for the unique physiological phenotype of high-altitude natives. This article is protected by copyright. All rights reserved.

Chromosome-level genome assemblies of four wild peach species provide insights into genome evolution and genetic basis of stress resistance

Background

Peach (Prunus persica) is an economically important stone fruit crop in Rosaceae and widely cultivated in temperate and subtropical regions, emerging as an excellent material to study the interaction between plant and environment. During its genus, there are four wild species of peach, all living in harsh environments. For example, one of the wild species, P. mira, originates from the Qinghai-Tibet Plateau (QTP) and exhibits strong cold/ultraviolet ray environmental adaptations. Although remarkable progresses in the gene discovery of fruit quality-related traits in peach using previous assembled genome were obtained, genomic basis of the response of these wild species to different geographical environments remains unclear.

Results

To uncover key genes regulating adaptability in different species and analyze the role of genetic variations in resistance formation, we performed de novo genome assembling of four wild relatives of peach (P. persica), P. mira, P. davidiana, P. kansuensis, and P. ferganensis and resequenced 175 peach varieties. The phylogenetic tree showed that the divergence time of P. mira and other wild relatives of peach was 11.5 million years ago, which was consistent with the drastic crustal movement of QTP. Abundant genetic variations were identified in four wild species when compared to P. persica, and the results showed that plant-pathogen interaction pathways were enriched in genes containing small insertions and deletions and copy number variations in all four wild relatives of peach. Then, the data were used to identify new genes and variations regulating resistance. For example, presence/absence variations which result from a hybridization event that occurred between P. mira and P. dulcis enhanced the resistance of their putative hybrid, P. davidiana. Using bulked segregant analysis, we located the nematode resistance locus of P. kansuensis in chromosome 2. Within the mapping region, a deletion in the promoter of one NBS-LRR gene was found to involve the resistance by regulating gene expression. Furthermore, combined with RNA-seq and selective sweeps analysis, we proposed that a deletion in the promoter of one CBF gene was essential for high-altitude adaptation of P. mira through increasing its resistance to low temperature.

Conclusions

In general, the reference genomes assembled in the study facilitate our understanding of resistance mechanism of perennial fruit crops, and provide valuable resources for future breeding and improvement.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01342-y.

Incomplete Concordance Between Host Phylogeny and Gut Microbial Community in Tibetan Wetland Birds

Gut microbial communities of animals play key roles in host evolution, while the relationship between gut microbiota and host evolution in Tibetan birds remains unknown. Herein, we sequenced the gut microbiota of 67 wild birds of seven species dwelling in the Tibetan wetlands. We found an obvious species-specific structure of gut microbiota among these plateau birds whose habitats were overlapped. Different from plateau mammals, there was no strict synergy between the hierarchical tree of gut microbial community and species phylogeny. In brown-headed gulls (Larus brunnicephalus) as an example, the structure of gut microbiota differed in different habitats, and the relative abundance of bacteria, such as Lactobacillus, Streptococcus, Paracoccus, Lachnospiraceae, and Vibrio, significantly correlated with altitude. Finally, we found various pathogenic bacteria in the birds of these plateau wetlands, and the interspecific differences were related to their diet and living environments.

Comparative genomic analysis of high-altitude adaptation for Mongolia Mastiff, Tibetan Mastiff, and Canis Lupus

Tibetan Mastiff has adapted to the extreme environment of the Qinghai-Tibetan Plateau. Yet, the underlying mechanisms of its high-altitude-adaptation and origin remains elusive. Here, we generated the draft genomes of Mongolia Mastiff, Tibetan Mastiff, and Canis Lupus. The phylogenetic tree uncovered that Tibetan Mastiff and Mongolia Mastiff were derived from Canis Lupus species. The comparative genomic analyses identified that the expansion of gene families related to DNA repair and damage response, and contraction related to ATPase activity revealed the genetic adaptations of Tibetan Mastiff and Canis Lupus to high altitude. In addition, the Tibetan Mastiff and Canis Lupus had signals of positive selection for genes involved in fatty-acid α/β- oxidation for highland adaptation. Notably, the positively selected TERT of Tibetan Mastiff should be an adaptive trait for correcting DNA damage. These findings suggested that the Tibetan Mastiff and Canis Lupus evolves basic strategies for adaptation to high altitude.

Characterization of olfactory receptor repertoires provides insights into the high-altitude adaptation of the yak based on the chromosome-level genome

Olfaction in vertebrates plays pivotal parts in many aspects, such as localizing prey or food, mating behavior, avoiding predators, and social communication. Yak (Bos grunniens) is the only Bos species that can thrive in high-altitude areas. In view of the critical role of olfactory receptors (ORs) in the specific recognition of diverse stimuli, investigating the evolutionary dynamics of ORs in the yak means a lot. In this study, we used the chromosome-level genome of the yak to identify the ORs genes and discussed the effects of high altitude on the yak's olfaction by comparing the yak with other low-altitude living Bos species (Bos frontalis (gayal), Bos gaurus (gaur), Bos indicus (zebu) and Bos taurus (cattle)). The yak had 400 OR genes, including 264 functional genes, 16 partial genes and 120 OR pseudo genes. There were 387 OR genes mapped to yak 31 chromosomes, and chromosomes 13 and 8 had the most OR genes and functional OR genes. Among these five Bos species, yak had the least number of OR gene subfamilies, OR genes and functional OR genes, while the total number of OR genes in gayal (n = 784) was almost twice as many as that of yak, indicating that the olfaction of yak may be less developed. In addition, the phylogenetic relationships of the functional Bos OR genes were illustrated, which comprised 79 families and 466 subfamilies distributed in two classes (Class I and Class II). There were 76 OR gene subfamilies shared by these five Bos species and 17 OR gene subfamilies were unique to the yak. The potential odor specificity of 44 yak OR genes was identified through the similarity to human OR protein sequences. Remarkably, yak lacks β-ionone and Isovaleric acid(IVA)-related ORs, which may be related to the decline of high-altitude herbaceous plant diversity and underdeveloped yak sweat glands. The conserved motifs of OR genes were highly conserved in Bos species. These results provided a solid foundation for further studies on the molecular mechanisms of the yak's adaptation to the high-altitude environment in olfaction.

A Multilevel Assessment of Plasticity in Response to High-Altitude Environment for Agama Lizards

Upslope range shifting has been documented in diverse species in response to global warming. Plasticity, which refers to the ability of organisms to alter their phenotypes in changing environments, is crucial for the survival of those that newly migrated to a high-altitude environment. The scope and mechanisms of plasticity across biological levels, however, have rarely been examined. We used two agama lizards (genus Phrynocephalus) as model systems and a transplant experiment to comprehensively assess their plasticity on multiple organization levels. Two low-altitude (934 m) agama species, Phrynocephalus axillaris (oviparous) and P. forsythii (viviparous), were transplanted to a high-altitude site (3,400 m). After acclimation for 6 weeks in seminatural enclosures, plasticity was measured from bite force, tail display behavior, gene expression, and metabolome. Both lizards were capable of acclimating to the high-altitude environment without sacrificing their performance in bite force, but they also showed high plasticity in tail display behavior by either decreasing the intensity of a specific display component (P. forsythii) or by the trade-off between display components (P. axillaris). Genes and metabolites associated with lipids, especially fatty acid metabolism, exhibited significant differentiation in expression, compared to individuals from their native habitats. Improved fatty acid storage and metabolism appeared to be a common response among animals at high altitudes. Despite distinct reproductive modes that may differ in response to physiological pressure, the two lizards demonstrated high concordance in plasticity when they faced a novel environment at high altitudes. Taken together, lizards likely acclimate to high-altitude environments by reducing behavioral activity and increasing energy efficiency after range shifting. Our results provide new insights into our understanding of phenotypic plasticity and its importance in today’s changing climate.

The Discovery of Chicken Foxp3 Demands Redefinition of Avian Regulatory T Cells

Since the publication of the first chicken genome sequence, we have encountered genes playing key roles in mammalian immunology, but being seemingly absent in birds. One of those was, until recently, Foxp3, the master transcription factor of regulatory T cells in mammals. Therefore, avian regulatory T cell research is still poorly standardized. In this study we identify a chicken ortholog of Foxp3 We prove sequence homology with known mammalian and sauropsid sequences, but also reveal differences in major domains. Expression profiling shows an association of Foxp3 and CD25 expression levels in CD4⁺CD25⁺ peripheral T cells and identifies a CD4^-CD25⁺Foxp3^high subset of thymic lymphocytes that likely represents yet undescribed avian regulatory T precursor cells. We conclude that Foxp3 is existent in chickens and that it shares certain functional characteristics with its mammalian ortholog. Nevertheless, pathways for regulatory T cell development and Foxp3 function are likely to differ between mammals and birds. The identification and characterization of chicken Foxp3 will help to define avian regulatory T cells and to analyze their functional properties and thereby advance the field of avian immunology.

Comparative Genomics and Evolution of Avian Specialized Traits

Genomic data are important for understanding the origin and evolution of traits. Under the context of rapidly developing of sequencing technologies and more widely available genome sequences, researchers are able to study evolutionary mechanisms of traits via comparative genomic methods. Compared with other vertebrates, bird genomes are relatively small and exhibit conserved synteny with few repetitive elements, which makes them suitable for evolutionary studies. Increasing genomic progress has been reported on the evolution of powered flight, body size variation, beak morphology, plumage colouration, high-elevation colonization, migration, and vocalization. By summarizing previous studies, we demonstrate the genetic bases of trait evolution, highlighting the roles of small-scale sequence variation, genomic structural variation, and changes in gene interaction networks. We suggest that future studies should focus on improving the quality of reference genomes, exploring the evolution of regulatory elements and networks, and combining genomic data with morphological, ecological, behavioural, and developmental biology data.

Genome sequencing and transcriptomic analysis of the Andean killifish Orestias ascotanensis reveals adaptation to high-altitude aquatic life

Orestias ascotanensis (Cyprinodontidae) is a teleost pupfish endemic to springs feeding into the Ascotan saltpan in the Chilean Altiplano (3,700 m.a.s.l.) and represents an opportunity to study adaptations to high-altitude aquatic environments. We have de novo assembled the genome of O. ascotanensis at high coverage. Comparative analysis of the O. ascotanensis genome showed an overall process of contraction, including loss of genes related to G-protein signaling, chemotaxis and signal transduction, while there was expansion of gene families associated with microtubule-based movement and protein ubiquitination. We identified 818 genes under positive selection, many of which are involved in DNA repair. Additionally, we identified novel and conserved microRNAs expressed in O. ascotanensis and its closely-related species, Orestias gloriae. Our analysis suggests that positive selection and expansion of genes that preserve genome stability are a potential adaptive mechanism to cope with the increased solar UV radiation to which high-altitude animals are exposed to.

Parallel genomic responses to historical climate change and high elevation in East Asian songbirds

Parallel evolution can be expected among closely related taxa exposed to similar selective pressures. However, parallelism is typically stronger at the phenotypic level, while genetic solutions to achieve these phenotypic similarities may differ. For polygenic traits, the availability of standing genetic variation (i.e., heterozygosity) may influence such genetic nonparallelism. Here, we examine the extent to which high-elevation adaptation is parallel-and whether the level of parallelism is affected by heterozygosity-by analyzing genomes of 19 Paridae species distributed across East Asia with a dramatic east-west elevation gradient. We find that western highlands endemic parids have consistently lower levels of heterozygosity-likely the result of late-Pleistocene demographic contraction-than do parids found exclusively in eastern lowlands, which remained unglaciated during the late Pleistocene. Three widespread species (east to west) have high levels of heterozygosity similar to that observed in eastern species, although their western populations are less variable than eastern ones. Comparing genomic responses to extreme environments of the Qinghai-Tibet Plateau, we find that the most differentiated genomic regions between each high-elevation taxon and its low-elevation relative are significantly enriched for genes potentially related to the oxygen transport cascade and/or thermogenesis. Despite no parallelism at particular genes, high similarity in gene function is found among comparisons. Furthermore, parallelism is not higher in more heterozygous widespread parids than in highland endemics. Thus, in East Asian parids, parallel functional response to extreme elevation appears to rely on different genes, with differences in heterozygosity having no effect on the degree of genetic parallelism.

The De Novo Genome Sequencing of Silver Pheasant (Lophura nycthemera)

Silver pheasant (Lophura nycthemera) belongs to Phasianidae, Galliformes, which exhibits high subspecific differentiation. In this study, we assembled a novel genome based on 98.42 Gb of Illumina sequencing data and 30.20 Gb of PacBio sequencing data. The size of the final assembled genome was 1.01 Gb, with a contig N50 of 6.96 Mb. Illumina paired-end reads (94.96%) were remapped to the contigs. The assemble genome shows high completeness, with a complete BUSCO score of 92.35% using the avian data set. A total of 16,747 genes were predicted from the generated assembly, and 16,486 (98.44%) of the genes were annotated. The average length of genes, exons, and introns were 19,827.53, 233.69, and 1841.19 bp, respectively. Noncoding RNAs included 208 miRNAs, 40 rRNAs, and 264 tRNAs, and a total of 189 pseudogenes were identified; 116.31 Mb (11.47%) of the genome consisted of repeat sequences, with the greatest proportion of LINEs. This assembled genome provides a valuable reference genome for further studies on the evolutionary history and conversion genetics of L. nycthemera and the phylogenomics of the Galliformes lineage.