Elevation as a selective force on mitochondrial respiratory chain complexes of the Phrynocephalus lizards in the Tibetan plateau

Gut Virome of the World's Highest-Elevation Lizard Species (Phrynocephalus erythrurus and Phrynocephalus theobaldi) Reveals Versatile Commensal Viruses

The gut virome is a reservoir of diverse symbiotic and pathogenic viruses coevolving with their hosts, and yet limited research has explored the gut viromes of highland-dwelling rare species. Using viral metagenomic analysis, the viral communities of the Phrynocephalus lizards living in the Qinghai-Tibet Plateau were investigated. Phage-encoded functional genes and antibiotic resistance genes (ARGs) were analyzed. The viral communities of different lizard species were all predominated by bacteriophages, especially the Caudovirales order. The virome of Phrynocephalus erythrurus living around the Namtso Lake possessed a unique structure, with the greatest abundance of the Parvoviridae family and the highest number of exclusive viral species. Several vertebrate-infecting viruses were discovered, including caliciviruses, astroviruses, and parvoviruses. Phylogenetic analyses demonstrated that the virus hallmark genes of bacteriophages possessed high genetic diversity. After functional annotation, the majority of phage-associated functional genes were classified in the energy metabolism category. In addition, plenty of ARGs belonging to the multidrug category were discovered, and five ARGs were exclusive to the virome from Phrynocephalus theobaldi. This study provided the first insight into the structure and function of the virome in highland lizards, contributing to the protection of threatened lizard species. Also, our research is of exemplary significance for the gut virome research of lizard species and other cold-blooded and highland-dwelling animals, prompting a better understanding of the interspecific differences and transmission of commensal viruses. IMPORTANCE The Phrynocephalus lizards inhabiting the Qinghai-Tibet Plateau (QTP) are considered to be the highest-altitude lizard species in the world, and they have been added to the IUCN list of threatened species. Living in the QTP with hypoxic, arid, natural conditions, the lizards presented a unique pattern of gut virome, which could provide both positive and negative effects, such as the enrichment of functional genes and the dissemination of antibiotic resistance genes (ARGs). This work provides the foundation for further research on the gut virome in these endangered lizard species and other cold-blooded and highland-dwelling animals, contributing to the maintenance of ecological balance on the plateau.

De novo Assembly, Annotation, and Analysis of Transcriptome Data of the Ladakh Ground Skink Provide Genetic Information on High-Altitude Adaptation

The Himalayan Arc is recognized as a global biodiversity hotspot. Among its numerous cryptic and undiscovered organisms, this composite high-mountain ecosystem harbors many taxa with adaptations to life in high elevations. However, evolutionary patterns and genomic features have been relatively rarely studied in Himalayan vertebrates. Here, we provide the first well-annotated transcriptome of a Greater Himalayan reptile species, the Ladakh Ground skink Asymblepharus ladacensis (Squamata: Scincidae). Based on tissues from the brain, an embryonic disc, and pooled organ material, using pair-end Illumina NextSeq 500 RNAseq, we assembled ~77,000 transcripts, which were annotated using seven functional databases. We tested ~1600 genes, known to be under positive selection in anurans and reptiles adapted to high elevations, and potentially detected positive selection for 114 of these genes in Asymblepharus. Even though the strength of these results is limited due to the single-animal approach, our transcriptome resource may be valuable data for further studies on squamate reptile evolution in the Himalayas as a hotspot of biodiversity.

A Re-Assessment of Positive Selection on Mitochondrial Genomes of High-Elevation Phrynocephalus Lizards

Due to their integral roles in oxidative phosphorylation, mitochondrially encoded proteins represent common targets of selection in response to altitudinal hypoxia across high-altitude taxa. While previous studies revealed evidence of positive selection on mitochondrial genomes of high-altitude Phrynocephalus lizards, their conclusions were restricted by out-of-date phylogenies and limited taxonomic sampling. Using topologies derived from both nuclear and mitochondrial DNA phylogenies, we re-assessed the evidence of positive selection on the mitochondrial genomes of high-altitude Phrynocephalus. We sampled representative species from all four main lineages and sequenced the mitochondrial genome of P. maculatus, a putative sister taxon to the high-altitude group. Positive selection was assessed through two widely used branch-site tests: the branch-site model in PAML and BUSTED in HyPhy. No evidence of positive selection on mitochondrial genes was detected on branches leading to two most recent common ancestors of high-altitude species; however, we recovered evidence of positive selection on COX1 on the P. forsythii branch, which represents a reversal from high- to low-elevation environments. A positively selected site therein marked a threonine to valine substitution at position 419. We suggest this bout of selection occurred as the ancestors of P. forsythii re-colonized lower altitude environments north of the Tibetan Plateau. Despite their role in oxidative phosphorylation, we posit that mitochondrial genes are unlikely to have represented historical targets of selection for high-altitude adaptation in Phrynocephalus. Consequently, future studies should address the roles of nuclear genes and differential gene expression.