Recent geological events and intrinsic behavior influence the population genetic structure of the chiru and tibetan gazelle on the tibetan plateau

Uncovering rearrangements in the Tibetan antelope via population-derived genome refinement and comparative analysis with homologous species

Introduction: The Tibetan antelope (Pantholops hodgsonii) is a remarkable mammal thriving in the extreme Qinghai-Tibet Plateau conditions. Despite the availability of its genome sequence, limitations in the scaffold-level assembly have hindered a comprehensive understanding of its genomics. Moreover, comparative analyses with other Bovidae species are lacking, along with insights into genome rearrangements in the Tibetan antelope. Methods: Addressing these gaps, we present a multifaceted approach by refining the Tibetan Antelope genome through linkage disequilibrium analysis with data from 15 newly sequenced samples. Results: The scaffold N50 of the refined reference is 3.2 Mbp, surpassing the previous version by 1.15-fold. Our annotation analysis resulted in 50,750 genes, encompassing 29,324 novel genes not previously study. Comparative analyses reveal 182 unique rearrangements within the scaffolds, contributing to our understanding of evolutionary dynamics and species-specific adaptations. Furthermore, by conducting detailed genomic comparisons and reconstructing rearrangements, we have successfully pioneered the reconstruction of the X-chromosome in the Tibetan antelope. Discussion: This effort enhances our comprehension of the genomic landscape of this species.

Vigilance behaviour during the calving season in female Tibetan antelopes (Pantholopshodgsonii)

Tibetan antelopes (Pantholopshodgsonii) migrate great distances to specific delivery and calving areas. In the current study, we investigated calving site selection and vigilance behaviour during delivery and nursing in migratory female Tibetan antelopes at Zonag Lake. According to observations and analysis, the females were distributed south of Zonag Lake, where vegetation was abundant. We determined their dates of migration (crossing the Qinghai-Tibet Highway observation site), showing a shift of one month during the period from June in 2008 to May 2021. Results also showed that 81.4% of females expressed high vigilance behaviour during calving and nursing compared to those without calves (7.1%). From delivery until calf standing, females were highly vigilant and spent considerable time scanning, with 96% of females showing vigilance behaviour. Females with calves (average 9.94 ± 0.62 s) spent more time on vigilance behaviour than females without calves (average 6.25 ± 1.38 s). Females with newborns spent the greatest amount of time being vigilant (average 51.63 ± 4.24 s). These results not only identify basic Tibetan antelope calving behaviour, but also provide scientific analysis and evidence for further ethological research on female Tibetan antelopes.

Phylogeography of the Plateau Pika (Ochotona curzoniae) in Response to the Uplift of the Qinghai-Tibet Plateau

The evolution and current distribution of species on the Qinghai-Tibet Plateau have been significantly impacted by historical occurrences, including the uplift of the plateau and the Quaternary climate upheaval. As a remnant species, the plateau pika (Ochotona curzoniae) is a great model for researching historical events. In this study, 302 samples from 42 sample sites were utilized to analyze the impact of historical events on the evolution and distribution pattern of plateau pikas. The genetic diversity, patterns of differentiation, and historical dynamics of the plateau pika were investigated using molecular markers that included four mitochondrial genes (COI, D-loop, Cytb, and 12S rRNA) and three nuclear genes (GHR, IRBP, and RAG1). The results showed that: (1) The genetic diversity of the plateau pika was high in the Tibetan Plateau (Hd = 0.9997, π = 0.01205), and the plateau pika evolved into five lineages that occupied different geographical areas, with lineage 1 (Group 1) in the south of the Yarlung Zangbo River, lineage 2 (Group 2) in the hinterland of the plateau, lineage 3 (Group 3) in the northeastern part of the plateau, lineage 4 (Group 4) in the Hengduan Mountains, and lineage 5 (Group 5) in the eastern part of the plateau. (2) The gene flow among the five lineages was low, and the differentiation level was high (Nm < 0.25; Fst > 0.25), indicating that the geographical barriers between the five lineages, such as the Yarlung Zangbo River, the Qaidam-Ghuong-Guide Basin, and the Lancang River, effectively promoted the population differentiation of the plateau pika. (3) The plateau pika first spread from the Hengduan Mountains to the entire Qinghai-Tibet Plateau and then conducted small-scale migration and dispersal in several refuges across the plateau in response to climate changes during the glacial and interglacial periods. (4) Except for Group 1 and Group 4, all the other populations exhibited a rapid expansion between 0.06 and 0.01 Mya, but the expansion was considerably delayed or halted by the effects of climate change during the last glacial maximum (0.02 Mya). Overall, the plateau pika on the Qinghai-Tibet Plateau exhibits high genetic diversity, and topographic obstacles, including mountains, valleys, and basins, created by the uplift of the plateau and climatic changes since the Quaternary period have played an important role in the differentiation and historical dynamics of the plateau pika population.

Revisiting the Woolly wolf (Canis lupus chanco) phylogeny in Himalaya: Addressing taxonomy, spatial extent and distribution of an ancient lineage in Asia

Of the sub-species of Holarctic wolf, the Woolly wolf (Canis lupus chanco) is uniquely adapted to atmospheric hypoxia and widely distributed across the Himalaya, Qinghai Tibetan Plateau (QTP) and Mongolia. Taxonomic ambiguity still exists for this sub-species because of complex evolutionary history anduse of limited wild samples across its range in Himalaya. We document for the first time population genetic structure and taxonomic affinity of the wolves across western and eastern Himalayan regions from samples collected from the wild (n = 19) using mitochondrial control region (225bp). We found two haplotypes in our data, one widely distributed in the Himalaya that was shared with QTP and the other confined to Himachal Pradesh and Uttarakhand in the western Himalaya, India. After combining our data withpublished sequences (n = 83), we observed 15 haplotypes. Some of these were shared among different locations from India to QTP and a few were private to geographic locations. A phylogenetic tree indicated that Woolly wolves from India, Nepal, QTP and Mongolia are basal to other wolves with shallow divergence (K2P; 0.000-0.044) and high bootstrap values. Demographic analyses based on mismatch distribution and Bayesian skyline plots (BSP) suggested a stable population over a long time (~million years) with signs of recent declines. Regional dominance of private haplotypes across its distribution range may indicate allopatric divergence. This may be due to differences in habitat characteristics, availability of different wild prey species and differential deglaciation within the range of the Woolly wolf during historic time. Presence of basal and shallow divergence within-clade along with unique ecological requirements and adaptation to hypoxia, the Woolly wolf of Himalaya, QTP, and Mongolian regions may be considered as a distinct an Evolutionary Significant Unit (ESU). Identifying management units (MUs) is needed within its distribution range using harmonized multiple genetic data for effective conservation planning.

Diversity and Evolution of Sensor Histidine Kinases in Eukaryotes

Histidine kinases (HKs) are primary sensor proteins that act in cell signaling pathways generically referred to as "two-component systems" (TCSs). TCSs are among the most widely distributed transduction systems used by both prokaryotic and eukaryotic organisms to detect and respond to a broad range of environmental cues. The structure and distribution of HK proteins are now well documented in prokaryotes, but information is still fragmentary for eukaryotes. Here, we have taken advantage of recent genomic resources to explore the structural diversity and the phylogenetic distribution of HKs in the prominent eukaryotic supergroups. Searches of the genomes of 67 eukaryotic species spread evenly throughout the phylogenetic tree of life identified 748 predicted HK proteins. Independent phylogenetic analyses of predicted HK proteins were carried out for each of the major eukaryotic supergroups. This allowed most of the compiled sequences to be categorized into previously described HK groups. Beyond the phylogenetic analysis of eukaryotic HKs, this study revealed some interesting findings: 1) characterization of some previously undescribed eukaryotic HK groups with predicted functions putatively related to physiological traits; 2) discovery of HK groups that were previously believed to be restricted to a single kingdom in additional supergroups, and 3) indications that some evolutionary paths have led to the appearance, transfer, duplication, and loss of HK genes in some phylogenetic lineages. This study provides an unprecedented overview of the structure and distribution of HKs in the Eukaryota and represents a first step toward deciphering the evolution of TCS signaling in living organisms.

The roles of calving migration and climate change in the formation of the weak genetic structure in the Tibetan antelope (Pantholops hodgsonii)

Geographical barriers and distance can reduce gene exchange among animals, resulting in genetic divergence of geographically isolated populations. The Tibetan antelope (Pantholops hodgsonii) has a geographical range of approximately 1600 km across the Qinghai-Tibet Plateau, which comprises a series of tall mountains and big rivers. However, previous studies indicate that there is little genetic differentiation among their geographically delineated populations. To better understand the genetic structure of P. hodgsonii populations, we collected 145 samples from the 3 major calving regions, taking into consideration their various calving grounds and migration routes. We used a combination of mitochondrial sequences (Cyt b, ATPase, D-loop and COX I) to investigate the genetic structure and the evolutionary divergence of the populations. Significant, albeit weak, genetic differentiation was detected among the 3 geographical populations. Analysis of the genetic divergence process revealed that the animals gradually entered a period of rapid genetic differentiation approximately 60 000 years ago. The calving migration of P. hodgsonii cannot be the main cause of their weak genetic structure because this cannot fully homogenize the genetic pool. Instead, the geological and climatic events as well as the coupling vegetation succession process during this period have been suggested to greatly contribute to the genetic structure and the expansion of genetic diversity.

Palghat gap reveals presence of two diverged populations of Nilgiri tahr (Nilgiritragus hylocrius) in Western Ghats, India

Genetic analysis is an important tool in understanding population structure, genetic diversity, and phylogenetics of endangered species likely to be affected by microevolution and anthropogenic factors. Western Ghats landscape is one of the identified biodiversity hotspots in India, and micro-evolutionary processes are observed in this landscape due to the presence of the gaps in the mountain ranges. Nilgiri tahr is endemic to and distributed in this landscape while very little is known about genetic characteristics, population structure and impact of these gaps on the species. In the present study, two different populations of Nilgiri tahr from the north (NPG) and south (SPG) of Palghat gap (PG) were studied using the cytochrome b gene (Cyt b; 310 bp) of mtDNA genome in the Western Ghats, India. Two variable sites were observed in the Cyt b fragment while the mean pairwise genetic distance between these two populations was 0.007. All the samples phylogenetically clustered in either north or south of PG. The presence of shallow divergence indicates the presence of suitable habitat in past which may have facilitated movement between NPG and SPG. A subsequent change in Paleo-climatic conditions and gradual formation of PG may have resulted in population diversification during the Pleistocene. Besides, Forensically Informative Nucleotide Sequence (FINS) observed would help in geo-assigning any individual from NPG or SPG to understand the likely influences on population demography due to poaching.

Genetic diversity of the Tibetan antelope (Pantholops hodgsonii) population of Ladakh, India, its relationship with other populations and conservation implications

BACKGROUND

The Tibetan antelope (Pantholops hodgsonii), or chiru, is an endangered antelope, distributed in China [Xinjiang, Xizang, Qinghai, Zhuolaihu Lake (Breeding habitat)], and India (Aksai Chin and Ladakh). There is a global demand for the species prized wool, which is used in weaving shahtoosh shawls. Over the years, the population of the Tibetan antelope has drastically declined from more than a million to a few thousand individuals, mainly due to poaching. Field studies undertaken in Ladakh, India also indicated winter migration of the population to Tibet. Migration between winter and calving habitats is well established to be female-biased across the Qinghai Tibetan Plateau (QTP). For effective conservation planning, genetic characterization is considered the best way to understand the likely impact of threats for ensuring the long-term viability of the population. In this regard, genetic characteristics of all Chinese populations are well-studied using mitochondrial and microsatellite markers, but information is lacking for the Indian population. Therefore, using the control region marker, we document for the first time the genetic variation of the Indian population of the Tibetan antelope, the extent of migration and its relationships with other populations of China.

RESULTS

The partial fragment of control region (259 bp) marker was successfully amplified in 30 Tibetan antelope samples that were collected from the Chang Chenmo Valley in eastern Ladakh, India. We also retrieved control region sequences (n = 88) available in the public domain from GenBank of different Chinese populations. Low levels of nucleotide (π; 0.004) and haplotype (hd; 0.543) diversity were observed in the Indian population when compared to Chinese populations (π = 0.01357-0.02048 and hd = 0.889-0.986). Commonly used indices (Tajima's D and Fu's Fs) were analyzed for inferring the demographic history of the Indian populations, and all values were negative indicating population expansion or demographic equilibrium, though nonsignificant. We observed five haplotypes in the Indian population, and these were not reported in previously studied populations of QTP. Bayesian-based phylogenetic analysis indicates the presence of four clades, however, the posterior probability support for three of these clades is weak (<0.5). Of these, the Indian population formed a distinct clade, whereas the Chinese populations exhibited shared haplotypes, and no geographic structure was observed. Median-joining network analysis was conducted for 46 haplotypes in the overall population, except the samples from India which showed a star-like topology. The Indian population is separated by one median vector from the Chinese population.

CONCLUSIONS

The present study revealed the presence of different sub-clades in the Bayesian phylogenetic tree and five new haplotypes only in the Indian population or sampling location. Furthermore, in the phylogenetic tree, Indian haplotypes of Tibetan antelopes were clustered with the haplotype reported in the Chinese population of the Xinjiang region. Median-joining network analysis showed shared haplotypes pattern in all populations of QTP except the samples from India which showed new haplotypes. Given the presence of low nucleotide and haplotype diversity in eastern Ladakh populations and limited information available for populations of the western side in its range, we suggest to include genetic studies of Tibetan antelope populations around Aksai Chin (Fig. 1) under the proposed transboundary agenda between India and China and assess relationships with other populations. Such understanding would enable the planning of conservation strategies for ensuring long-term survival of westernmost populations in its range, and if required, it would establish connectivity with the other populations.

Morphological and Genetic Variation along a North-to-South Transect in Stipa purpurea, a Dominant Grass on the Qinghai-Tibetan Plateau: Implications for Response to Climate Change

Estimating the potential of species to cope with rapid environmental climatic modifications is of vital importance for determining their future viability and conservation. The variation between existing populations along a climatic gradient may predict how a species will respond to future climate change. Stipa purpurea is a dominant grass species in the alpine steppe and meadow of the Qinghai-Tibetan Plateau (QTP). Ecological niche modelling was applied to S. purpurea, and its distribution was found to be most strongly correlated with the annual precipitation and the mean temperature of the warmest quarter. We established a north-to-south transect over 2000 km long on the QTP reflecting the gradients of temperature and precipitation, and then we estimated the morphological by sampling fruited tussocks and genetic divergence by using 11 microsatellite markers between 20 populations along the transect. Reproductive traits (the number of seeds and reproductive shoots), the reproductive-vegetative growth ratio and the length of roots in the S. purpurea populations varied significantly with climate variables. S. purpurea has high genetic diversity (He = 0.585), a large effective population size (Ne >1,000), and a considerable level of gene flow between populations. The S. purpurea populations have a mosaic genetic structure: some distant populations (over 1000 km apart) clustered genetically, whereas closer populations (< 100 km apart) had diverged significantly, suggesting local adaptation. Asymmetrical long-distance inter-population gene flow occurs along the sampling transect and might be mediated by seed dispersal via migratory herbivores, such as the chiru (Pantholops hodgsonii). These findings suggest that population performance variation and gene flow both facilitate the response of S. purpurea to climate change.

Niche divergence accelerates evolution in Asian endemic Procapra gazelles

Ecological niche divergence and adaptation to new environments are thought to play important roles in driving speciation. Whether recently evolved species show evidence for niche divergence or conservation is vital towards understanding the role of ecology in the process of speciation. The genus Procapra is an ancient, monophyletic lineage endemic to Asia that contains three extant species (P. gutturosa, P. przewalskii and P. picticaudata). These species mainly inhabit the Qinghai-Tibetan and Mongolian Plateaus, and today have primarily allopatric distributions. We applied a series of geographic information system-based analyses to test for environmental variation and niche divergence among these three species. We found substantial evidence for niche divergence in species' bioclimatic preferences, which supports the hypothesis that niche divergence accelerates diversification in Procapra. Our results provide important insight into the evolutionary history of ungulates in Asia and help to elucidate how environmental changes accelerate lineage diversification.