Ancient and modern genomes unravel the evolutionary history of the rhinoceros family.

Only five species of the once-diverse Rhinocerotidae remain, making the reconstruction of their evolutionary history a challenge to biologists since Darwin. We sequenced genomes from five rhinoceros species (three extinct and two living), which we compared to existing data from the remaining three living species and a range of outgroups. We identify an early divergence between extant African and Eurasian lineages, resolving a key debate regarding the phylogeny of extant rhinoceroses. This early Miocene (∼16 million years ago [mya]) split post-dates the land bridge formation between the Afro-Arabian and Eurasian landmasses. Our analyses also show that while rhinoceros genomes in general exhibit low levels of genome-wide diversity, heterozygosity is lowest and inbreeding is highest in the modern species. These results suggest that while low genetic diversity is a long-term feature of the family, it has been particularly exacerbated recently, likely reflecting recent anthropogenic-driven population declines.

Climate-driven flyway changes and memory-based long-distance migration.

Millions of migratory birds occupy seasonally favourable breeding grounds in the Arctic1, but we know little about the formation, maintenance and future of the migration routes of Arctic birds and the genetic determinants of migratory distance. Here we established a continental-scale migration system that used satellite tracking to follow 56 peregrine falcons (Falco peregrinus) from 6 populations that breed in the Eurasian Arctic, and resequenced 35 genomes from 4 of these populations. The breeding populations used five migration routes across Eurasia, which were probably formed by longitudinal and latitudinal shifts in their breeding grounds during the transition from the Last Glacial Maximum to the Holocene epoch. Contemporary environmental divergence between the routes appears to maintain their distinctiveness. We found that the gene ADCY8 is associated with population-level differences in migratory distance. We investigated the regulatory mechanism of this gene, and found that long-term memory was the most likely selective agent for divergence in ADCY8 among the peregrine populations. Global warming is predicted to influence migration strategies and diminish the breeding ranges of peregrine populations of the Eurasian Arctic. Harnessing ecological interactions and evolutionary processes to study climate-driven changes in migration can facilitate the conservation of migratory birds.

Historic Sampling of a Vanishing Beast: Population Structure and Diversity in the Black Rhinoceros.

The black rhinoceros (Diceros bicornis L.) is a critically endangered species historically distributed across sub-Saharan Africa. Hunting and habitat disturbance have diminished both its numbers and distribution since the 19th century, but a poaching crisis in the late 20th century drove them to the brink of extinction. Genetic and genomic assessments can greatly increase our knowledge of the species and inform management strategies. However, when a species has been severely reduced, with the extirpation and artificial admixture of several populations, it is extremely challenging to obtain an accurate understanding of historic population structure and evolutionary history from extant samples. Therefore, we generated and analyzed whole genomes from 63 black rhinoceros museum specimens collected between 1775 and 1981. Results showed that the black rhinoceros could be genetically structured into six major historic populations (Central Africa, East Africa, Northwestern Africa, Northeastern Africa, Ruvuma, and Southern Africa) within which were nested four further subpopulations (Maasailand, southwestern, eastern rift, and northern rift), largely mirroring geography, with a punctuated north-south cline. However, we detected varying degrees of admixture among groups and found that several geographical barriers, most prominently the Zambezi River, drove population discontinuities. Genomic diversity was high in the middle of the range and decayed toward the periphery. This comprehensive historic portrait also allowed us to ascertain the ancestry of 20 resequenced genomes from extant populations. Lastly, using insights gained from this unique temporal data set, we suggest management strategies, some of which require urgent implementation, for the conservation of the remaining black rhinoceros diversity.

Multiple genomic solutions for local adaptation in two closely related species (sheep and goats) facing the same climatic constraints.

The question of how local adaptation takes place remains a fundamental question in evolutionary biology. The variation of allele frequencies in genes under selection over environmental gradients remains mainly theoretical and its empirical assessment would help understanding how adaptation happens over environmental clines. To bring new insights to this issue we set up a broad framework which aimed to compare the adaptive trajectories over environmental clines in two domesticated mammal species co-distributed in diversified landscapes. We sequenced the genomes of 160 sheep and 161 goats extensively managed along environmental gradients, including temperature, rainfall, seasonality and altitude, to identify genes and biological processes shaping local adaptation. Allele frequencies at putatively adaptive loci were rarely found to vary gradually along environmental gradients, but rather displayed a discontinuous shift at the extremities of environmental clines. Of the 430 candidate adaptive genes identified, only 6 were orthologous between sheep and goats and those responded differently to environmental pressures, suggesting different putative mechanisms involved in local adaptation in these two closely related species. Interestingly, the genomes of the 2 species were impacted differently by the environment, genes related to signatures of selection were most related to altitude, slope and rainfall seasonality for sheep, and summer temperature and spring rainfall for goats. The diversity of candidate adaptive pathways may result from a high number of biological functions involved in the adaptations to multiple eco-climatic gradients, and a differential role of climatic drivers on the two species, despite their co-distribution along the same environmental gradients. This study describes empirical examples of clinal variation in putatively adaptive alleles with different patterns in allele frequency distributions over continuous environmental gradients, thus showing the diversity of genetic responses in adaptive landscapes and opening new horizons for understanding genomics of adaptation in mammalian species and beyond.

The importance of well protected forests for the conservation genetics of West African colobine monkeys.

In tropical forests, anthropogenic activities are major drivers of the destruction and degradation of natural habitats, causing severe biodiversity loss. African colobine monkeys (Colobinae) are mainly folivore and strictly arboreal primates that require large forests to subsist, being among the most vulnerable of all nonhuman primates. The Western red colobus Piliocolobus badius and the King colobus Colobus polykomos inhabit highly fragmented West African forests, including the Cantanhez Forests National Park (CFNP) in Guinea-Bissau. Both species are also found in the largest and best-preserved West African forest-the Taï National Park (TNP) in Ivory Coast. Colobine monkeys are hunted for bushmeat in both protected areas, but these exhibit contrasting levels of forest fragmentation, thus offering an excellent opportunity to investigate the importance of well-preserved forests for the maintenance of evolutionary potential in these arboreal primates. We estimated genetic diversity, population structure, and demographic history by using microsatellite loci and mitochondrial DNA. We then compared the genetic patterns of the colobines from TNP with the ones previously obtained for CFNP and found contrasting genetic patterns. Contrary to the colobines from CFNP that showed very low genetic diversity and a strong population decline, the populations in TNP still maintain high levels of genetic diversity and we found no clear signal of population decrease in Western red colobus and a limited decrease in King colobus. These results suggest larger and historically more stable populations in TNP compared to CFNP. We cannot exclude the possibility that the demographic effects resulting from the recent increase of bushmeat hunting are not yet detectable in TNP using genetic data. Nevertheless, the fact that the TNP colobus populations are highly genetically diverse and maintain large effective population sizes suggests that well-preserved forests are crucial for the maintenance of populations, species, and probably for the evolutionary potential in colobines.

Genetic diversity goals and targets have improved, but remain insufficient for clear implementation of the post-2020 global biodiversity framework.

Genetic diversity among and within populations of all species is necessary for people and nature to survive and thrive in a changing world. Over the past three years, commitments for conserving genetic diversity have become more ambitious and specific under the Convention on Biological Diversity's (CBD) draft post-2020 global biodiversity framework (GBF). This Perspective article comments on how goals and targets of the GBF have evolved, the improvements that are still needed, lessons learned from this process, and connections between goals and targets and the actions and reporting that will be needed to maintain, protect, manage and monitor genetic diversity. It is possible and necessary that the GBF strives to maintain genetic diversity within and among populations of all species, to restore genetic connectivity, and to develop national genetic conservation strategies, and to report on these using proposed, feasible indicators.

Historical Introgression from Wild Relatives Enhanced Climatic Adaptation and Resistance to Pneumonia in Sheep.

How animals, particularly livestock, adapt to various climates and environments over short evolutionary time is of fundamental biological interest. Further, understanding the genetic mechanisms of adaptation in indigenous livestock populations is important for designing appropriate breeding programs to cope with the impacts of changing climate. Here, we conducted a comprehensive genomic analysis of diversity, interspecies introgression, and climate-mediated selective signatures in a global sample of sheep and their wild relatives. By examining 600K and 50K genome-wide single nucleotide polymorphism data from 3,447 samples representing 111 domestic sheep populations and 403 samples from all their seven wild relatives (argali, Asiatic mouflon, European mouflon, urial, snow sheep, bighorn, and thinhorn sheep), coupled with 88 whole-genome sequences, we detected clear signals of common introgression from wild relatives into sympatric domestic populations, thereby increasing their genomic diversities. The introgressions provided beneficial genetic variants in native populations, which were significantly associated with local climatic adaptation. We observed common introgression signals of alleles in olfactory-related genes (e.g., ADCY3 and TRPV1) and the PADI gene family including in particular PADI2, which is associated with antibacterial innate immunity. Further analyses of whole-genome sequences showed that the introgressed alleles in a specific region of PADI2 (chr2: 248,302,667-248,306,614) correlate with resistance to pneumonia. We conclude that wild introgression enhanced climatic adaptation and resistance to pneumonia in sheep. This has enabled them to adapt to varying climatic and environmental conditions after domestication.

Genomic erosion in a demographically recovered bird species during conservation rescue.

The pink pigeon (Nesoenas mayeri) is an endemic species of Mauritius that has made a remarkable recovery after a severe population bottleneck in the 1970s to early 1990s. Prior to this bottleneck, an ex situ population was established from which captive-bred individuals were released into free-living subpopulations to increase population size and genetic variation. This conservation rescue led to rapid population recovery to 400-480 individuals, and the species was twice downlisted on the International Union for the Conservation of Nature (IUCN) Red List. We analyzed the impacts of the bottleneck and genetic rescue on neutral genetic variation during and after population recovery (1993-2008) with restriction site-associated sequencing, microsatellite analyses, and quantitative genetic analysis of studbook data of 1112 birds from zoos in Europe and the United States. We used computer simulations to study the predicted changes in genetic variation and population viability from the past into the future. Genetic variation declined rapidly, despite the population rebound, and the effective population size was approximately an order of magnitude smaller than census size. The species carried a high genetic load of circa 15 lethal equivalents for longevity. Our computer simulations predicted continued inbreeding will likely result in increased expression of deleterious mutations (i.e., a high realized load) and severe inbreeding depression. Without continued conservation actions, it is likely that the pink pigeon will go extinct in the wild within 100 years. Conservation rescue of the pink pigeon has been instrumental in the recovery of the free-living population. However, further genetic rescue with captive-bred birds from zoos is required to recover lost variation, reduce expression of harmful deleterious variation, and prevent extinction. The use of genomics and modeling data can inform IUCN assessments of the viability and extinction risk of species, and it helps in assessments of the conservation dependency of populations.

La paloma rosada (Nesoenas mayeri) es una especie endémica de Mauricio que se ha recuperado impresionantemente después de un grave cuello de botella poblacional a principios de la década de 1970 que duró hasta inicios de la década de 1990. Antes de este cuello de botella se había establecido una población ex situ de la cual se liberaban individuos reproducidos en cautiverio a las subpoblaciones en libertad para incrementar la variación genética y el tamaño poblacional. Este rescate de conservación derivó en una recuperación rápida de la población (400-480 individuos) y la especie cambió positivamente de categoría dos veces en la Lista Roja de la Unión Internacional para la Conservación de la Naturaleza (UICN). Analizamos los impactos del cuello de botella y el rescate genético sobre la variación genética neutral durante y después de la recuperación poblacional (de 1993 a 2008) mediante secuenciación RAD, análisis de microsatélites y análisis genéticos cuantitativos de los datos del libro genealógico de 1112 aves ubicadas en zoológicos de Europa y los Estados Unidos. Usamos simulaciones por computadora para estudiar los cambios pronosticados en la variación genética y en la viabilidad poblacional del pasado hacia el futuro. La variación genética declinó rápidamente, a pesar de la recuperación poblacional, y el tamaño efectivo de la población fue aproximadamente un orden de magnitud más pequeño que el tamaño del censo. La especie contó con una carga genética elevada de casi 15 equivalentes letales para la longevidad. Nuestras simulaciones pronostican que la endogamia continua probablemente resultará en un incremento en la expresión de mutaciones deletéreas (es decir, una carga realizada elevada) y en una depresión endogámica severa. Sin acciones continuas para la conservación, es probable que la paloma rosada esté extinta en vida libre dentro de cien años. El rescate de conservación de la paloma rosada ha sido fundamental en la recuperación de la población silvestre; sin embargo, se requiere de un rescate genético adicional con las aves de reproducción en cautiverio de los zoológicos para recuperar la variación perdida, reducir la expresión de la variación deletérea dañina y prevenir la extinción. El uso de la genómica y los datos modelados puede orientar las valoraciones de la UICN sobre la viabilidad y el riesgo de extinción de las especies, además de que ayuda en la evaluación de la dependencia que tienen las poblaciones de la conservación.

Broad maternal geographic origin of domestic sheep in Anatolia and the Zagros.

We investigated the controversial origin of domestic sheep (Ovis aries) using large samples of contemporary and ancient domestic individuals and their closest wild relatives: the Asiatic mouflon (Ovis gmelini), the urial (Ovis vignei) and the argali (Ovis ammon). A phylogeny based on mitochondrial DNA, including 213 new cytochrome-b sequences of wild Ovism confirmed that O. gmelini is the maternal ancestor of sheep and precluded mtDNA contributions from O. vignei (and O. gmelini × O. vignei hybrids) to domestic lineages. We also produced 54 new control region sequences showing shared haplogroups (A, B, C and E) between domestic sheep and wild O. gmelini which localized the domestication center in eastern Anatolia and central Zagros, excluding regions further east where exclusively wild haplogroups were found. This overlaps with the geographic distribution of O. gmelini gmelini, further suggesting that the maternal origin of domestic sheep derives from this subspecies. Additionally, we produced 57 new CR sequences of Neolithic sheep remains from a large area covering Anatolia to Europe, showing the early presence of at least three mitochondrial haplogroups (A, B and D) in Western colonization routes. This confirmed that sheep domestication was a large-scale process that captured diverse maternal lineages (haplogroups).

Draft genome of a biparental beetle species, Lethrus apterus.

BACKGROUND: The lack of an understanding about the genomic architecture underpinning parental behaviour in subsocial insects displaying simple parental behaviours prevents the development of a full understanding about the evolutionary origin of sociality. Lethrus apterus is one of the few insect species that has biparental care. Division of labour can be observed between parents during the reproductive period in order to provide food and protection for their offspring. RESULTS: Here, we report the draft genome of L. apterus, the first genome in the family Geotrupidae. The final assembly consisted of 286.93 Mbp in 66,933 scaffolds. Completeness analysis found the assembly contained 93.5% of the Endopterygota core BUSCO gene set. Ab initio gene prediction resulted in 25,385 coding genes, whereas homology-based analyses predicted 22,551 protein coding genes. After merging, 20,734 were found during functional annotation. Compared to other publicly available beetle genomes, 23,528 genes among the predicted genes were assigned to orthogroups of which 1664 were in species-specific groups. Additionally, reproduction related genes were found among the predicted genes based on which a reduction in the number of odorant- and pheromone-binding proteins was detected. CONCLUSIONS: These genes can be used in further comparative and functional genomic researches which can advance our understanding of the genetic basis and hence the evolution of parental behaviour.

Interspecific Gene Flow and the Evolution of Specialization in Black and White Rhinoceros.

Africa's black (Diceros bicornis) and white (Ceratotherium simum) rhinoceros are closely related sister-taxa that evolved highly divergent obligate browsing and grazing feeding strategies. Although their precursor species Diceros praecox and Ceratotherium mauritanicum appear in the fossil record ∼5.2 Ma, by 4 Ma both were still mixed feeders, and were even spatiotemporally sympatric at several Pliocene sites in what is today Africa's Rift Valley. Here, we ask whether or not D. praecox and C. mauritanicum were reproductively isolated when they came into Pliocene secondary contact. We sequenced and de novo assembled the first annotated black rhinoceros reference genome and compared it with available genomes of other black and white rhinoceros. We show that ancestral gene flow between D. praecox and C. mauritanicum ceased sometime between 3.3 and 4.1 Ma, despite conventional methods for the detection of gene flow from whole genome data returning false positive signatures of recent interspecific migration due to incomplete lineage sorting. We propose that ongoing Pliocene genetic exchange, for up to 2 My after initial divergence, could have potentially hindered the development of obligate feeding strategies until both species were fully reproductively isolated, but that the more severe and shifting paleoclimate of the early Pleistocene was likely the ultimate driver of ecological specialization in African rhinoceros.

Global Commitments to Conserving and Monitoring Genetic Diversity Are Now Necessary and Feasible.

Global conservation policy and action have largely neglected protecting and monitoring genetic diversity-one of the three main pillars of biodiversity. Genetic diversity (diversity within species) underlies species' adaptation and survival, ecosystem resilience, and societal innovation. The low priority given to genetic diversity has largely been due to knowledge gaps in key areas, including the importance of genetic diversity and the trends in genetic diversity change; the perceived high expense and low availability and the scattered nature of genetic data; and complicated concepts and information that are inaccessible to policymakers. However, numerous recent advances in knowledge, technology, databases, practice, and capacity have now set the stage for better integration of genetic diversity in policy instruments and conservation efforts. We review these developments and explore how they can support improved consideration of genetic diversity in global conservation policy commitments and enable countries to monitor, report on, and take action to maintain or restore genetic diversity.

Dispersal and genetic structure in a tropical small mammal, the Bornean tree shrew (Tupaia longipes), in a fragmented landscape along the Kinabatangan River, Sabah, Malaysia.

BACKGROUND: Constraints in migratory capabilities, such as the disruption of gene flow and genetic connectivity caused by habitat fragmentation, are known to affect genetic diversity and the long-term persistence of populations. Although negative population trends due to ongoing forest loss are widespread, the consequence of habitat fragmentation on genetic diversity, gene flow and genetic structure has rarely been investigated in Bornean small mammals. To fill this gap in knowledge, we used nuclear and mitochondrial DNA markers to assess genetic diversity, gene flow and the genetic structure in the Bornean tree shrew, Tupaia longipes, that inhabits forest fragments of the Lower Kinabatangan Wildlife Sanctuary, Sabah. Furthermore, we used these markers to assess dispersal regimes in male and female T. longipes. RESULTS: In addition to the Kinabatangan River, a known barrier for dispersal in tree shrews, the heterogeneous landscape along the riverbanks affected the genetic structure in this species. Specifically, while in larger connected forest fragments along the northern riverbank genetic connectivity was relatively undisturbed, patterns of genetic differentiation and the distribution of mitochondrial haplotypes in a local scale indicated reduced migration on the strongly fragmented southern riverside. Especially, oil palm plantations seem to negatively affect dispersal in T. longipes. Clear sex-biased dispersal was not detected based on relatedness, assignment tests, and haplotype diversity. CONCLUSION: This study revealed the importance of landscape connectivity to maintain migration and gene flow between fragmented populations, and to ensure the long-term persistence of species in anthropogenically disturbed landscapes.

Transcription-Associated Mutation Promotes RNA Complexity in Highly Expressed Genes-A Major New Source of Selectable Variation.

Alternatively spliced transcript isoforms are thought to play a critical role for functional diversity. However, the mechanism generating the enormous diversity of spliced transcript isoforms remains unknown, and its biological significance remains unclear. We analyzed transcriptomes in saker falcons, chickens, and mice to show that alternative splicing occurs more frequently, yielding more isoforms, in highly expressed genes. We focused on hemoglobin in the falcon, the most abundantly expressed genes in blood, finding that alternative splicing produces 10-fold more isoforms than expected from the number of splice junctions in the genome. These isoforms were produced mainly by alternative use of de novo splice sites generated by transcription-associated mutation (TAM), not by the RNA editing mechanism normally invoked. We found that high expression of globin genes increases mutation frequencies during transcription, especially on nontranscribed DNA strands. After DNA replication, transcribed strands inherit these somatic mutations, creating de novo splice sites, and generating multiple distinct isoforms in the cell clone. Bisulfate sequencing revealed that DNA methylation may counteract this process by suppressing TAM, suggesting DNA methylation can spatially regulate RNA complexity. RNA profiling showed that falcons living on the high Qinghai-Tibetan Plateau possess greater global gene expression levels and higher diversity of mean to high abundance isoforms (reads per kilobases per million mapped reads ≥18) than their low-altitude counterparts, and we speculate that this may enhance their oxygen transport capacity under low-oxygen environments. Thus, TAM-induced RNA diversity may be physiologically significant, providing an alternative strategy in lifestyle evolution.

Genomic selection strategies for breeding adaptation and production in dairy cattle under climate change.

Livestock production both contributes to and is affected by global climate change, and substantial modifications will be required to increase its climate resilience. In this context, reliance on dominant commercial livestock breeds, featuring small effective population sizes, makes current production strategies vulnerable if their production is restricted to environments, which may be too costly to support under future climate scenarios. The adaptability of animal populations to future environments will therefore become important. To help evaluate the role of genetics in climate adaptation, we compared selection strategies in dairy cattle using breeding simulations, where genomic selection was used on two negatively correlated traits for production (assumed to be moderately heritable) and adaptation (assumed to have low heritability). Compared with within-population breeding, genomic introgression produced a more positive genetic change for both production and adaptation traits. Genomic introgression from highly adapted but low production value populations into highly productive but low adaptation populations was most successful when the adaptation trait was given a lower selection weight than the production trait. Genomic introgression from highly productive population to highly adapted population was most successful when the adaptation trait was given a higher selection weight than the production trait. Both these genomic introgression schemes had the lowest risk of inbreeding. Our results suggest that both adaptation and production can potentially be improved simultaneously by genomic introgression.

Contrasting evolutionary history, anthropogenic declines and genetic contact in the northern and southern white rhinoceros (Ceratotherium simum).

The white rhinoceros (Ceratotherium simum) has a discontinuous African distribution, which is limited by the extent of sub-Saharan grasslands. The southern population (SWR) declined to its lowest number around the turn of the nineteenth century, but recovered to become the world's most numerous rhinoceros. In contrast, the northern population (NWR) was common during much of the twentieth century, declining rapidly since the 1970s, and now only two post-reproductive individuals remain. Despite this species's conservation status, it lacks a genetic assessment of its demographic history. We therefore sampled 232 individuals from extant and museum sources and analysed ten microsatellite loci and the mtDNA control region. Both marker types reliably partitioned the species into SWR and NWR, with moderate nuclear genetic diversity and only three mtDNA haplotypes for the species, including historical samples. We detected ancient interglacial demographic declines in both populations. Both populations may also have been affected by recent declines associated with the colonial expansion for the SWR, and with the much earlier Bantu migrations for the NWR. Finally, we detected post-divergence secondary contact between NWR and SWR, possibly occurring as recently as the last glacial maximum. These results suggest the species was subjected to regular periods of fragmentation and low genetic diversity, which may have been replenished upon secondary contact during glacial periods. The species's current situation thus reflects prehistoric declines that were exacerbated by anthropogenic pressure associated with the rise of late Holocene technological advancement in Africa. Importantly, secondary contact suggests a potentially positive outcome for a hybrid rescue conservation strategy, although further genome-wide data are desirable to corroborate these results.

Population Genomics Reveals Low Genetic Diversity and Adaptation to Hypoxia in Snub-Nosed Monkeys.

Snub-nosed monkeys (genus Rhinopithecus) are a group of endangered colobines endemic to South Asia. Here, we re-sequenced the whole genomes of 38 snub-nosed monkeys representing four species within this genus. By conducting population genomic analyses, we observed a similar load of deleterious variation in snub-nosed monkeys living in both smaller and larger populations and found that genomic diversity was lower than that reported in other primates. Reconstruction of Rhinopithecus evolutionary history suggested that episodes of climatic variation over the past 2 million years, associated with glacial advances and retreats and population isolation, have shaped snub-nosed monkey demography and evolution. We further identified several hypoxia-related genes under selection in R. bieti (black snub-nosed monkey), a species that exploits habitats higher than any other nonhuman primate. These results provide the first detailed and comprehensive genomic insights into genetic diversity, demography, genetic burden, and adaptation in this radiation of endangered primates.

Population transcriptomes reveal synergistic responses of DNA polymorphism and RNA expression to extreme environments on the Qinghai-Tibetan Plateau in a predatory bird.

Low oxygen and temperature pose key physiological challenges for endotherms living on the Qinghai-Tibetan Plateau (QTP). Molecular adaptations to high-altitude living have been detected in the genomes of Tibetans, their domesticated animals and a few wild species, but the contribution of transcriptional variation to altitudinal adaptation remains to be determined. Here we studied a top QTP predator, the saker falcon, and analysed how the transcriptome has become modified to cope with the stresses of hypoxia and hypothermia. Using a hierarchical design to study saker populations inhabiting grassland, steppe/desert and highland across Eurasia, we found that the QTP population is already distinct despite having colonized the Plateau <2000 years ago. Selection signals are limited at the cDNA level, but of only seventeen genes identified, three function in hypoxia and four in immune response. Our results show a significant role for RNA transcription: 50% of upregulated transcription factors were related to hypoxia responses, differentiated modules were significantly enriched for oxygen transport, and importantly, divergent EPAS1 functional variants with a refined co-expression network were identified. Conservative gene expression and relaxed immune gene variation may further reflect adaptation to hypothermia. Our results exemplify synergistic responses between DNA polymorphism and RNA expression diversity in coping with common stresses, underpinning the successful rapid colonization of a top predator onto the QTP. Importantly, molecular mechanisms underpinning highland adaptation involve relatively few genes, but are nonetheless more complex than previously thought and involve fine-tuned transcriptional responses and genomic adaptation.

Mitogenomic Meta-Analysis Identifies Two Phases of Migration in the History of Eastern Eurasian Sheep.

Despite much attention, history of sheep (Ovis aries) evolution, including its dating, demographic trajectory and geographic spread, remains controversial. To address these questions, we generated 45 complete and 875 partial mitogenomic sequences, and performed a meta-analysis of these and published ovine mitochondrial DNA sequences (n = 3,229) across Eurasia. We inferred that O. orientalis and O. musimon share the most recent female ancestor with O. aries at approximately 0.790 Ma (95% CI: 0.637-0.934 Ma) during the Middle Pleistocene, substantially predating the domestication event (∼8-11 ka). By reconstructing historical variations in effective population size, we found evidence of a rapid population increase approximately 20-60 ka, immediately before the Last Glacial Maximum. Analyses of lineage expansions showed two sheep migratory waves at approximately 4.5-6.8 ka (lineages A and B: ∼6.4-6.8 ka; C: ∼4.5 ka) across eastern Eurasia, which could have been influenced by prehistoric West-East commercial trade and deliberate mating of domestic and wild sheep, respectively. A continent-scale examination of lineage diversity and approximate Bayesian computation analyses indicated that the Mongolian Plateau region was a secondary center of dispersal, acting as a "transportation hub" in eastern Eurasia: Sheep from the Middle Eastern domestication center were inferred to have migrated through the Caucasus and Central Asia, and arrived in North and Southwest China (lineages A, B, and C) and the Indian subcontinent (lineages B and C) through this region. Our results provide new insights into sheep domestication, particularly with respect to origins and migrations to and from eastern Eurasia.