Molecular phylogenetic relationship of snow finch complex (genera Montifringilla, Pyrgilauda, and Onychostruthus) from the Tibetan plateau

Museomics help resolving the phylogeny of snowfinches (Aves, Passeridae, Montifringilla and allies)

Historical specimens from museum collections provide a valuable source of material also from remote areas or regions of conflict that are not easily accessible to scientists today. With this study, we are providing a taxon-complete phylogeny of snowfinches using historical DNA from whole skins of an endemic species from Afghanistan, the Afghan snowfinch, Pyrgilauda theresae. To resolve the strong conflict between previous phylogenetic hypotheses, we generated novel mitogenome sequences for selected taxa and genome-wide SNP data using ddRAD sequencing for all extant snowfinch species endemic to the Qinghai-Tibet Plateau (QTP) and for an extended intraspecific sampling of the sole Central and Western Palearctic snowfinch species (Montifringilla nivalis). Our phylogenetic reconstructions unanimously refuted the previously suggested paraphyly of genus Pyrgilauda. Misplacement of one species-level taxon (Onychostruthus tazcanowskii) in previous snowfinch phylogenies was undoubtedly inferred from chimeric mitogenomes that included heterospecific sequence information. Furthermore, comparison of novel and previously generated sequence data showed that the presumed sister-group relationship between M. nivalis and the QTP endemic M. henrici was suggested based on flawed taxonomy. Our phylogenetic reconstructions based on genome-wide SNP data and on mitogenomes were largely congruent and supported reciprocal monophyly of genera Montifringilla and Pyrgilauda with monotypic Onychostruthus being sister to the latter. The Afghan endemic P. theresae likely originated from a rather ancient Pliocene out-of-Tibet dispersal probably from a common ancestor with P. ruficollis. Our extended trans-Palearctic sampling for the white-winged snowfinch, M. nivalis, confirmed strong lineage divergence between an Asian and a European clade dated to 1.5 - 2.7 million years ago (mya). Genome-wide SNP data suggested subtle divergence among European samples from the Alps and from the Cantabrian mountains.

Shennongjia-Wushan Mountains-One cryptic glacial refugium introduced by the phylogeographical study of the Geometridae moth Ourapteryx szechuana Wehrli

The origin and evolution of biodiversity in the Shennongjia and Wushan Mountains, located in central China, are little known. In this study, we used Ourapteryx szechuana, which is widely distributed in China and northern Nepal, to explore whether these mountains acted as glacial refugia during climate oscillations of the Quaternary. In total, 192 samples of O. szechuana were collected throughout much of the distribution range. Phylogenetic analysis, molecular dating, demographic history reconstructions, and MAXENT were used to investigate the evolutionary history and differentiation mechanisms and predict the potential species distributions during four different periods. The phylogenetic tree and the star-like median-joining network strongly supported two reciprocally monophyletic and allopatric lineages. Lineage I was restricted to the Shennongjia and Wushan Mountains. The divergence time of O. szechuana from its sister species O. thibetaria was approximately 1.94 Ma. The differentiation processes of the two intraspecific lineages occurred at approximately 0.47 Ma. The demographic history reconstruction and the ecological niche model suggested that Lineage II experienced an expansion after the LGM (Last Glacial Maximum), whereas Lineage I did not experience any expansion. Our results suggested the Naynayxungla glaciation promoted the divergence of the two lineages by restricting them to different refugia. The valleys of the Shennongjia-Wushan Mountains may have kept stable and warm (thus ice-free) environments during Quaternary glaciations, allowing this region to act as a glacial refugia. Our studies show that the Shennongjia and Wushan Mountains are likely to be important but little studied glacial refugia for the insect and thus worthy of more attention.

The evolution of ancestral and species-specific adaptations in snowfinches at the Qinghai-Tibet Plateau

Species in a shared environment tend to evolve similar adaptations under the influence of their phylogenetic context. Using snowfinches, a monophyletic group of passerine birds (Passeridae), we study the relative roles of ancestral and species-specific adaptations to an extreme high-elevation environment, the Qinghai-Tibet Plateau. Our ancestral trait reconstruction shows that the ancestral snowfinch occupied high elevations and had a larger body mass than most nonsnowfinches in Passeridae. Subsequently, this phenotypic adaptation diversified in the descendant species. By comparing high-quality genomes from representatives of the three phylogenetic lineages, we find that about 95% of genes under positive selection in the descendant species are different from those in the ancestor. Consistently, the biological functions enriched for these species differ from those of their ancestor to various degrees (semantic similarity values ranging from 0.27 to 0.5), suggesting that the three descendant species have evolved divergently from the initial adaptation in their common ancestor. Using a functional assay to a highly selective gene, DTL, we demonstrate that the nonsynonymous substitutions in the ancestor and descendant species have improved the repair capacity of ultraviolet-induced DNA damage. The repair kinetics of the DTL gene shows a twofold to fourfold variation across the ancestor and the descendants. Collectively, this study reveals an exceptional case of adaptive evolution to high-elevation environments, an evolutionary process with an initial adaptation in the common ancestor followed by adaptive diversification of the descendant species.

Coping with extremes: convergences of habitat use, territoriality, and diet in summer but divergences in winter between two sympatric snow finches on the Qinghai-Tibet Plateau

On the Qinghai-Tibet Plateau, extreme environmental conditions have imposed intense selective pressure on the evolution of phenotypic traits of wild animals. To date, limited information is available on behavioral and ecological traits concerning niche differentiation among sympatric animals on the Qinghai-Tibet Plateau, especially during winter when the environments are most severe. Here, we studied the seasonal variations in habitat occurrence, territorial behavior, and diet in two sympatric snow finches (the white-rumped snow finch, Onychostruthus taczanowskii, WRSF; and the rufous-necked snow finch, Pyrgilauda ruficollis, RNSF) to determine convergence and divergence of ecological traits in such severe climatic conditions. Our results showed that: (i) WRSF occupied rural areas as a dominant species throughout the annual cycle while RNSF occupied the rural areas in summer and then shifted to human-occupied areas in winter and spring; (ii) WRSFs exhibited robust aggressive behavior and territoriality during winter relative to RNSFs; (iii) the diets of both species varied with the season but did not vary between species except that WRSF ate significantly more seeds but RNSF consumed more starchy material derived from human food waste during winter. Therefore, the separations in the spatial niche and territoriality between WRSF and RNSF, especially in winter, may contribute to alleviating the pressure of interspecific competition, and promoting the coexistence of the two sympatric snow finches in the extreme environments on the Qinghai-Tibet Plateau.

Coping with extremes: Remarkably blunt adrenocortical responses to acute stress in two sympatric snow finches on the Qinghai-Tibet Plateau during winter relative to other seasons

The extreme climatic conditions (ECCs) of the Qinghai-Tibet Plateau impose strong selective pressures on the evolution of phenotypic traits in free-living animals. It is not well understood how animals on the Qinghai-Tibet Plateau modify their adrenocortical functions in response to both predictable and unpredictable events of ECCs, especially when the available resources are lowest during the wintering life-history stage. To uncover potential physiological mechanisms, we studied the life history stage dependent features of morphology, the plasma corticosterone response to acute stress and brain glucocorticoid receptor (GR) and mineralocorticoid receptor (MR) mRNA expression in two sympatric snow finches: the white-rumped snow finch (Onychostruthus taczanowskii, WRSF); and the rufous-necked snow finch, Pyrgilauda ruficollis, RNSF) in Qinghai Province, China. Our results showed that (a) baseline corticosterone and stressor-induced corticosterone levels significantly varied with life history stage, but not between the species; (b) in WRSF, GR mRNA expression in the paraventricular nucleus was higher in the wintering stage compared to the pre-basic molt stage. There were no differences in hippocampus MR mRNA expression between stages in either species; (c) in the wintering stage, the suppression of corticosterone secretion in both species was an unexpected strategy in free-living animals. Both convergent and divergent phenotypic traits of adrenocortical responses to acute stress in two sympatric snow finches contribute to our understanding of the coping mechanisms of closely related species in the severe winter on the Qinghai-Tibet Plateau.

The geography and timing of genetic divergence in the lizard Phrynocephalus theobaldi on the Qinghai-Tibetan plateau

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.

Opposite companion effect on flight initiation distance in sympatric species: plateau pika (Ochotona curzoniae) and White-rumped Snowfinch (Onychostruthus taczanowskii)

There are many factors influencing prey’s risk perception and escape decision during predator–prey encounters. The distance at which animals move away from perceived danger (often quantified as flight initiation distance or FID) has been used by behavioral ecologists to understand the economics of antipredator behavior. Using general linear models, we investigated escape decision-making processes in plateau pika (Ochotona curzoniae (Hodgson, 1858)) and White-rumped Snowfinch (Onychostruthus taczanowskii (Prjevalsky, 1876), formerly known as Montifringilla taczanowskii Przewalski, 1876) together and we found that (i) there are significant positive correlations between starting distance and FID in both species; (ii) pika escapes at a longer distance from an approaching intruder when it is far from its burrow; (iii) foraging animals tolerate closer distances than watching ones, both in pikas and in Snowfinches; (iv) conspecifics seem to have no effect on pikas’ escaping behavior, while the appearance of Snowfinches dramatically decreases the FID of pikas. On the contrary, conspecifics significantly decrease the FID of Snowfinches, while the presence of pikas has no effects. These findings provide new evidence, which are consistent with optimal escape theory. These novel results in multispecies interactions prompt us to pay attention to the potential relationship between pikas and Snowfinches, as well as the “bird–pika in one hole” phenomenon.

An Update of Wallace’s Zoogeographic Regions of the World

Modern attempts to produce biogeographic maps focus on the distribution of species, and the maps are typically drawn without phylogenetic considerations. Here, we generate a global map of zoogeographic regions by combining data on the distributions and phylogenetic relationships of 21,037 species of amphibians, birds, and mammals. We identify 20 distinct zoogeographic regions, which are grouped into 11 larger realms. We document the lack of support for several regions previously defined based on distributional data and show that spatial turnover in the phylogenetic composition of vertebrate assemblages is higher in the Southern than in the Northern Hemisphere. We further show that the integration of phylogenetic information provides valuable insight on historical relationships among regions, permitting the identification of evolutionarily unique regions of the world.

Pitfalls in comparisons of genetic distances: a case study of the avian family Acrocephalidae

Genetic distances are increasingly being used for identification and species delimitation, especially since the introduction of "barcoding". While for phylogenetic inferences great care is generally taken to choose the best-fit evolutionary model, this is usually neglected in calculating genetic distances. Moreover, distances obtained from others than best-fit models, different lengths of sequences, and even different loci are often freely compared. We examined the influence of different methods on calculating genetic distances using mitochondrial cytochrome b sequences for the passerine family Acrocephalidae. We found substantial differences between: (1) corrected distances based on the best-fit model (TrN+Γ) vs. uncorrected p-distances; (2) distances calculated based on different parts of the same gene; and (3) distances calculated using the methods of "complete deletion" vs. "pairwise deletion" for sequences that included uncertain nucleotides. All these methodological differences affected comparisons between species and potential taxonomical conclusions. We suggest that (1) different loci are incomparable. (2) Only perfectly homologous regions (same length, same part of locus) should be compared. (3) In the case of sequences with some uncertain nucleotides, only distances calculated by the method of "complete deletion" are fully comparable. (4) Only distances based on the optimal substitution model should be used. (5) Even within the same locus, corrected genetic distances are unique to the study in which they are calculated, as they are conditional on the particular dataset and model selected for that dataset.

Phylogeography of Phrynocephalus erythrurus from the Qiangtang Plateau of the Tibetan Plateau

Phrynocephalus erythrurus of the Qiangtang Plateau occupies the highest regions of any reptile on earth. Here, we report mitochondrial DNA haplotypes sampled throughout the distribution of P. erythrurus and analyze patterns of genetic divergence among populations. The species diverged into two major lineages/subspecies at 3.67mya corresponding to the Northern and Southern Qiangtang Plateau. The Northern Qiangtang lineage diverged into two subpopulations at 2.76mya separated by the Beilu River Region and Wulanwula Mountains. Haplotypes from the southern Qiangtang lineage diverged 0.98mya as a star-shaped pattern. Analyses of molecular variance indicated that most of the observed genetic variation occurred among populations/regions implying long-term interruptions to gene flow. There was no evidence of sudden recent range expansions within any of the clades/lineages. NCPA infers allopatric fragmentation and restricted gene flow as the most likely mechanisms of population differentiation. Our results also indicate the presence of at least three refugia since the Hongya glaciation. Mountain movement and glaciations since mid-Pliocene are considered to have shaped phylogenetic patterns of P. erythrurus. P. erythrurus parva is suggested as a valid subspecies of P. erythrurus. Using four calibration points, we estimate an evolutionary rate of 0.762% divergence per lineage per million years for a mitochondrial genomic segment comprising the genes encoding ND2, tRNA(Trp) and tRNA(Ala).

The effect of landscape features on population genetic structure in Yunnan snub-nosed monkeys (Rhinopithecus bieti) implies an anthropogenic genetic discontinuity

The Tibetan Plateau is one of the top 10 biodiversity hotspots in the world and acts as a modern harbour for many rare species because of its relatively pristine state. In this article, we report a landscape genetic study on the Yunnan snub-nosed monkey (Rhinopithecus bieti), a primate endemic to the Tibetan Plateau. DNA was extracted from blood, tissue and fecal samples of 135 wild individuals representing 11 out of 15 extant monkey groups. Ten microsatellite loci were used to characterize patterns of genetic diversity. The most striking feature of the population structure is the presence of five subpopulations with distinct genetic backgrounds and unique spatial regions. The population structure of R. bieti appears to be shaped by anthropogenic landscape features as gene flow between subpopulations is strongly impeded by arable land, highways and human habitation. A partial Mantel test showed that 36.23% (r = 0.51, P = 0.01) of the genetic distance was explained by habitat gaps after controlling for the effect of geographical distance. Only 4.92% of the genetic distance was explained by geographical distance in the partial Mantel test, and no significant correlation was found. Estimation of population structure history indicates that environmental change during the last glacial maximum and human impacts since the Holocene, or a combination of both, have shaped the observed population structure of R. bieti. Increasing human activity on the Plateau, especially that resulting in habitat fragmentation, is becoming an important factor in shaping the genetic structure and evolutionary potential of species inhabiting this key ecosystem.

A molecular phylogeny of the Sylvia cantillans complex: cryptic species within the Mediterranean basin

The subalpine warbler Sylvia cantillans is formally considered a polytypic species, with four subspecies, European S. c. cantillans, albistriata, moltonii (recently resumed name: subalpina) and North African S. c.inornata. They are very similar in external morphology but clearly differ in their vocalizations. We evaluated their uncertain taxonomic status reconstructing the phylogenetic and phylogeographic relationships among populations sampled across major biogeographical areas in the European species' range, using nucleotide sequences of the mitochondrial cytochrome b gene (mtDNA cyt b). A variety of phylogenetic analyses concordantly led to identify four major groups, only partially corresponding to the three European nominal subspecies. Phylogenetic trees showed a monophyletic group including all moltonii individuals, well diverged from all other taxa. Populations taxonomically assigned to cantillans were polyphyletic being split into two distinct clades (western and southern cantillans), with monophyletic albistriata closely related to southern cantillans. Individuals of moltonii and southern cantillans sampled in sites of sympatry in central Italy were assigned to their respective groups, with perfect concordance between phenotypic and genetic identifications. All findings indicate that moltonii should be ranked as a distinct species. Former subspecies cantillans is polyphyletic, but additional data are needed to define the taxonomic status of its two clades. Albistriata is phylogenetically related to southern cantillans and should be provisionally kept as a subspecies of S. cantillans. The cantillans complex thus provides an interesting case-study illustrating geographical structuring across small geographical ranges, and it exemplifies speciation through differentiation in allopatry leading to reproductive isolation after a secondary contact.