Search

Identification of the role of the MC1R gene has provided major insights into variation in skin pigmentation in several organisms, including humans, but the evolutionary genetics of this variation is less well established. Variation in this gene and its relationship with degree of melanism was analyzed in one of the world's highest-elevation lizards, Phrynocephalus theobaldi from the Qinghai-Tibetan Plateau. Individuals from the low-elevation group were shown to have darker dorsal pigmentation than individuals from a high-elevation group. The existence of climatic variation across these elevations was quantified, with lower elevations exhibiting higher air pressure, temperatures, and humidity, but less wind and insolation. Analysis of the MC1R gene in 214 individuals revealed amino acid differences at five sites between intraspecific sister lineages from different elevations, with two sites showing distinct fixed residues at low elevations. Three of the four single-nucleotide polymorphisms that underpinned these amino acid differences were highly significant outliers, relative to the generalized MC1R population structuring, suggestive of selection. Transfection of cells with an MC1R allele from a lighter high-elevation population caused a 43% reduction in agonist-induced cyclic AMP accumulation, and hence lowered melanin synthesis, relative to transfection with an allele from a darker low-elevation population. The high-elevation allele led to less efficient integration of the MC1R protein into melanocyte membranes. Our study identifies variation in the degree of melanism that can be explained by four or fewer MC1R substitutions. We establish a functional link between these substitutions and melanin synthesis and demonstrate elevation-associated shifts in their frequencies., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Many species have been established on the basis of morphology, with markers such as mtDNA used to confirm the existence of independent historical lineages. Discordance between morphology and gene trees makes this less straightforward. Genotyping by sequencing (GBS) was used to analyse general genomic divergence across two recognized high altitude lizard species found in the eastern Qinghai-Tibetan Plateau. One of the species (Phyrnocephalus guinanensis) is found on a large area of sand dune habitat and distinguished from the other (P. putjatia) by morphology. We found that the primary pattern of genomic divergence is discordant with these morphological species: northern P. putjatia populations from around the large saline Qinghai lake are genomically distinct from P. putjatia and P. guinanensis populations located south of the Qinghai South and Riyue Mountains. Two competing historical scenarios were assessed using approximate Bayesian computation which unequivocally favoured a split between populations separated by the Qinghai South and Riyue mountains over a split between morphological species. The findings indicate that historical vicariance due to geographical features underpins the phylogenetic split rather than ecology-mediated divergence between sand dune and non-sand areas which i) is predicted by the mtDNA tree (showing the utility of this marker in species delimitation) and ii) demonstrates the unsuitability of the morphology-based taxonomy (indicating that large morphological differences do not always reflect historical lineages). In addition, we found a clear signal of isolation-by-distance around the periphery of Qinghai lake which suggests: i) a high level of resolution by GBS for detecting local divergence and ii) restricted gene flow over relatively short geographic distances. Overall, we show how morphological variation can mislead taxonomic conclusions and the utility of GBS for resolving these issues., (Copyright © 2019 Elsevier Inc. All rights reserved.)

The Qinghai-Tibetan Plateau (QTP) represents one of the earth's most significant physical features and there is increasing interest in the historical generation of biodiversity within this region. We hypothesized that there should be clear geographically coherent genetic structuring within one of the world's highest altitude lizards, Phrynocephalus theobaldi, due to considerable historical population fragmentation in this environment. This was tested using a major mitochondrial DNA (mtDNA) survey and sequencing of two nuclear markers (AME and RAG-1) from P. theobaldi, from across the southern QTP. A Bayesian method (BPEC) was used to detect four geographically structured mtDNA clusters. A Bayesian phylogenetic tree, together with associated dating analyses, supported four corresponding evolutionary lineages with a timing of 3.74-7.03 Ma for the most basal P. theobaldi split and Pliocene splits of 2.97-5.79 Ma and 2.40-5.39 Ma in the two daughter lineages. Himalayan uplift and changes in the Jilong basin may have contributed to these divergences, but uplift of the Gangdese mountains is rejected due to its timing. The nuclear markers appeared to be sorted between the four mtDNA groups, and species delimitation analyses supported the four phylogeographical groups as candidate species. The study contributes to our understanding of biodiversity on the QTP.

The Qinghai-Tibetan Plateau (QTP) is an important biogeographical area and has recently become a focus for biodiversity studies. Phyrnocephalus lizards form a widespread Eurasian group with oviparous and viviparous reproductive modes, but two previous mtDNA studies of species from around the QTP have provided different phylogenetic hypotheses. We analysed three loci (mtDNA, RAG-1, AME) from all recognised Chinese Phrynocephalus species to reconstruct the speciation history of the group and to estimate species divergence times. The effects of mtDNA partitioning strategy on phylogenetic inference were examined. Bayes factor comparisons of marginal likelihoods (mLs) estimated using stepping-stone sampling revealed that partitioning strategy had a major impact on mL. Nevertheless, it had a negligible effect on the inferred tree topology. The impact of hard-bound uniform or equivalent soft-bound gamma speciation time calibration priors as well as the use of a fixed topology (as opposed to integration over all possible species histories) on divergence time estimation were also assessed, and found to have little impact on posterior estimates. All three gene trees and the species tree supported the hypothesis that the Chinese species form oviparous and viviparous sister clades. This was in agreement with an early mtDNA study but differed from a subsequent reanalysis of the mtDNA data. Inclusion of mtDNA from more widely distributed Phrynocephalus (from previous studies) indicates that the oviparous P. interscapularis from Central Asia lies outside the clade of Chinese viviparous and oviparous species, but that other Asian oviparous species lie within the Chinese oviparous clade. The median of the posterior on the divergence time of Chinese oviparous and viviparous species was 9.7 Ma ago (95% interval: 7.2-13.0 Ma ago), which coincides with major uplifting of the QTP and indicates that viviparity evolved when this clade became restricted to regions of high elevation. We also found that cladogenesis within the viviparous clade began around 5 Ma ago whereas those in the oviparous clade began around 8.6 Ma ago. We establish more robust estimates of divergence times and relationships within this important group and so provide improved insights into the origins of Phrynocephalus diversity across the QTP., (Copyright © 2013 Elsevier Inc. All rights reserved.)

Phrynocephalus vlangalii is restricted to dry sand or Gobi desert highlands between major mountain ranges in the Qinghai (Tibetan) Plateau. Mitochondrial DNA (mtDNA) sequence (partial ND2, tRNA(Trp) and partial tRNA(Ala)) was obtained from 293 Phrynocephalus sampled from 34 sites across the plateau. Partitioned Bayesian and maximum parsimony phylogenetic analyses revealed that P. vlangalii and two other proposed species (P. erythrus and P. putjatia) together form a monophyletic mtDNA clade which, in contrast with previous studies, does not include P. theobaldi and P. zetangensis. The main P. vlangalli clade comprises seven well-supported lineages that correspond to distinct geographical areas with little or no overlap, and share a most recent common ancestor at 5.06 +/- 0.68 million years ago (mya). This is much older than intraspecific lineages in other Tibetan animal groups. Analyses of molecular variance indicated that most of the observed genetic variation occurred among populations/regions implying long-term interruption of maternal gene flow. A combined approach based on tests of population expansion, estimation of node dates, and significance tests on clade areas indicated that phylogeographical structuring has been primarily shaped by three main periods of plateau uplift during the Pliocene and Pleistocene, specifically 3.4 mya, 2.5 mya and 1.7 mya. These periods corresponded to the appearance of several mountain ranges that formed physical barriers between lineages. Populations from the Qaidam Basin are shown to have undergone major demographic and range expansions in the early Pleistocene, consistent with colonization of areas previously covered by the huge Qaidam palaeolake, which desiccated at this time. The study represents one of the most detailed phylogeographical analyses of the Qinghai Plateau to date and shows how geological events have shaped current patterns of diversity.