Genomic analysis identified a potential novel molecular mechanism for high-altitude adaptation in sheep at the Himalayas

Selection signatures for high altitude adaptation in livestock: A review

High altitude adapted livestock species (cattle, yak, goat, sheep, and horse) has critical role in the human socioeconomic sphere and acts as good source of animal source products including milk, meat, and leather, among other things. These species sustain production and reproduction even in harsh environments on account of adaptation resulting from continued evolution of beneficial traits. Selection pressure leads to various adaptive strategies in livestock whose footprints are evident at the different genomic sites as the "Selection Signature". Scrutiny of these signatures provides us crucial insight into the evolutionary process and domestication of livestock adapted to diverse climatic conditions. These signatures have the potential to change the sphere of animal breeding and further usher the selection programmes in right direction. Technological revolution and recent strides made in genomic studies has opened the routes for the identification of selection signatures. Numerous statistical approaches and bioinformatics tools have been developed to detect the selection signature. Consequently, studies across years have identified candidate genes under selection region found associated with numerous traits which have a say in adaptation to high-altitude environment. This makes it pertinent to have a better understanding about the selection signature, the ways to identify and how to utilize them for betterment of livestock populations as well as farmers. This review takes a closer look into the general concept, various methodologies, and bioinformatics tools commonly employed in selection signature studies and summarize the results of recent selection signature studies related to high-altitude adaptation in various livestock species. This review will serve as an informative and useful insight for researchers and students in the field of animal breeding and evolutionary biology.

Genomic structure analysis and construction of DNA fingerprint for four sheep populations

The Yongdeng Qishan sheep (QS) is a sheep population found locally in China. To gain in-depth knowledge of its population characteristics, three control groups were chosen, comprising the Lanzhou fat-tailed sheep (LFT), TAN sheep (TAN), and Minxian black fur sheep (MBF), inhabiting the nearby environments. This study genotyped a total of 120 individuals from four sheep populations: QS, LFT, TAN, and MBF. Using Specific-Locus Amplified Fragment Sequencing, we conducted genetic diversity, population structure, and selective sweep analysis, and constructed the fingerprint of each population. In total, there were 782 535 single nucleotide polymorphism (SNP) variations identified, with most being situated within regions that are intergenic or intronic. The genetic diversity analysis revealed that the QS population exhibited lower genetic diversity compared to the other three populations. Consistent results were obtained from the principal component, phylogenetic tree, and population structure analysis, indicating significant genetic differences between QS and the other three populations. However, a certain degree of differentiation was observed within the QS population. The linkage disequilibrium (LD) patterns among the four populations showed clear distinctions, with the QS group demonstrating the most rapid LD decline. Kinship analysis supported the findings of population structure, dividing the 90 QS individuals into two subgroups consisting of 23 and 67 individuals. Selective sweep analysis identified a range of genes associated with reproduction, immunity, and adaptation to high-altitude hypoxia. These genes hold potential as candidate genes for marker-assisted selection breeding. Additionally, a total of 86 523 runs of homozygosity (ROHs) were detected, showing non-uniform distribution across chromosomes, with chromosome 1 having the highest coverage percentage and chromosome 26 the lowest. In the high-frequency ROH islands, 79 candidate genes were associated with biological processes such as reproduction and fat digestion and absorption. Furthermore, a DNA fingerprint was constructed for the four populations using 349 highly polymorphic SNPs. In summary, our research delves into the genetic diversity and population structure of QS population. The construction of DNA fingerprint profiles for each population can provide valuable references for the identification of sheep breeds both domestically and internationally.

Refining the genomic profiles of North African sheep breeds through meta-analysis of worldwide genomic SNP data

Introduction

The development of reproducible tools for the rapid genotyping of thousands of genetic markers (SNPs) has promoted cross border collaboration in the study of sheep genetic diversity on a global scale.

Methods

In this study, we collected a comprehensive dataset of 239 African and Eurasian sheep breeds genotyped at 37,638 filtered SNP markers, with the aim of understanding the genetic structure of 22 North African (NA) sheep breeds within a global context.

Results and discussion

We revealed asubstantial enrichment of the gene pool between the north and south shores of the Mediterranean Sea, which corroborates the importance of the maritime route in the history of livestock. The genetic structure of North African breeds mirrors the differential composition of genetic backgrounds following the breed history. Indeed, Maghrebin sheep stocks constitute a geographically and historically coherent unit with any breed-level genetic distinctness among them due to considerable gene flow. We detected a broad east-west pattern describing the most important trend in NA fat-tailed populations, exhibited by the genetic closeness of Egyptian and Libyan fat-tailed sheep to Middle Eastern breeds rather than Maghrebin ones. A Bayesian FST scan analysis revealed a set of genes with potentially key adaptive roles in lipid metabolism (BMP2, PDGFD VEGFA, TBX15, and WARS2), coat pigmentation (SOX10, PICK1, PDGFRA, MC1R, and MTIF) and horn morphology RXFP2) in Tunisian sheep. The local ancestry method detected a Merino signature in Tunisian Noire de Thibar sheep near the SULF1gene introgressed by Merino's European breeds. This study will contribute to the general picture of worldwide sheep genetic diversity.

Chromosome-level genome assembly of the kiang (Equus kiang) illuminates genomic basis for its high-altitude adaptation

The kiang (Equus kiang) can only be observed in the Qinghai-Tibet Plateau (QTP). The kiang displayed excellent athletic performance in the high-altitude environment, which attracted wide interest in the investigation of the potential adaptive mechanisms to the extreme environment. Here, we assembled a chromosome-level genome of the kiang based on Hi-C sequencing technology. A total of 324.14 Gb clean data were generated, and the chromosome-level genome with 26 chromosomes (25 + X) and scaffold N50 of 101.77 Mb was obtained for the kiang. The genomic synteny analysis revealed large-scale chromosomal rearrangement during the evolution process of Equus species. Phylogenetic and divergence analyses revealed that the kiang was the sister branch to the ass and diverged from a common ancestor at approximately 13.5 Mya. The expanded gene families were mainly related to the hypoxia response, metabolism, and immunity. The kiang suffered a significant loss of olfaction-related genes, which might indicate decreased olfactory sensibility. Positively selected genes (PSGs) detected in the kiang were mainly associated with hypoxia response. Especially, there were two species-specific missense amino acid mutations in the PSG STAT3 annotated in the hypoxia-inducible factor 1 signal pathway, which may play an important role in the high-altitude adaptation of the kiang. Moreover, structure variations in the kiang genome were also identified, which possibly contributed to the high-altitude adaptation of the kiang. Comparative analysis revealed a lot of species-specific insertions and deletions in the kiang genome, such as PIK3CB and AKT with 3258 and 189 bp insertions in the intron region, respectively, possibly affecting the expression and regulation of hypoxia-related downstream pathways. This study provided valuable genomic resources, and our findings help a better understanding of the underlying adaptive strategies to the high-altitude environment in the kiang.

The mitochondrial genome of the mountain wooly tapir, Tapirus pinchaque and a formal test of the effect of altitude on the adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates

BACKGROUND

The harsh conditions of high-altitude environments are known to drive the evolution of physiological and morphological traits in endothermic animals. These conditions are expected to result in the adaptive evolution of protein coding genes encoded in mitochondrial genomes that are vital for the oxidative phosphorylation pathway. In this study, we formally tested for signatures of adaptive evolution on mitochondrial protein coding genes in Tapirus pinchaque and other odd-toed ungulates inhabiting high-elevation environments.

RESULTS

The AT-rich mitochondrial genome of T. pinchaque is 16,750 bp long. A phylomitogenomic analysis supports the monophyly of the genus Tapirus and families in the Perissodactyla. The ratio of non-synonymous to synonymous substitutions demonstrated that all mitochondrial genes undergo purifying selection in T. pinchaque and other odd ungulates living at high elevations. Over this negative background selection, Branch Models suggested that cox3 and nad6 might be undergoing stronger purifying selection than other mitochondrial protein coding genes. Furthermore, Site Models suggested that one and four sites in nad2 and nad5, respectively, could be experiencing positive selection. However, these results were supported by Likelihood Ratio Tests but not Bayesian Empirical Bayes posterior probabilities. Additional analyses (in DataMonkey) indicated a relaxation of selection strength in nad6, evidence of episodic diversifying selection in cob, and revealed episodic positive/diversifying selection signatures for two sites in nad1, and one site each in nad2 and nad4.

CONCLUSION

The mitochondrial genome of T. pinchaque is an important genomic resource for conservation of this species and this study contributes to the understanding of adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates inhabiting high-altitude environments.

New developments in the biology of fibroblast growth factors

The fibroblast growth factor (FGF) family is composed of 18 secreted signaling proteins consisting of canonical FGFs and endocrine FGFs that activate four receptor tyrosine kinases (FGFRs 1-4) and four intracellular proteins (intracellular FGFs or iFGFs) that primarily function to regulate the activity of voltage-gated sodium channels and other molecules. The canonical FGFs, endocrine FGFs, and iFGFs have been reviewed extensively by us and others. In this review, we briefly summarize past reviews and then focus on new developments in the FGF field since our last review in 2015. Some of the highlights in the past 6 years include the use of optogenetic tools, viral vectors, and inducible transgenes to experimentally modulate FGF signaling, the clinical use of small molecule FGFR inhibitors, an expanded understanding of endocrine FGF signaling, functions for FGF signaling in stem cell pluripotency and differentiation, roles for FGF signaling in tissue homeostasis and regeneration, a continuing elaboration of mechanisms of FGF signaling in development, and an expanding appreciation of roles for FGF signaling in neuropsychiatric diseases. This article is categorized under: Cardiovascular Diseases > Molecular and Cellular Physiology Neurological Diseases > Molecular and Cellular Physiology Congenital Diseases > Stem Cells and Development Cancer > Stem Cells and Development.

The genomic basis of high-elevation adaptation in wild house mice (Mus musculus domesticus) from South America

Understanding the genetic basis of environmental adaptation in natural populations is a central goal in evolutionary biology. The conditions at high elevation, particularly the low oxygen available in the ambient air, impose a significant and chronic environmental challenge to metabolically active animals with lowland ancestry. To understand the process of adaptation to these novel conditions and to assess the repeatability of evolution over short timescales, we examined the signature of selection from complete exome sequences of house mice (Mus musculus domesticus) sampled across two elevational transects in the Andes of South America. Using phylogenetic analysis, we show that house mice colonized high elevations independently in Ecuador and Bolivia. Overall, we found distinct responses to selection in each transect and largely nonoverlapping sets of candidate genes, consistent with the complex nature of traits that underlie adaptation to low oxygen availability (hypoxia) in other species. Nonetheless, we also identified a small subset of the genome that appears to be under parallel selection at the gene and SNP levels. In particular, three genes (Col22a1, Fgf14, and srGAP1) bore strong signatures of selection in both transects. Finally, we observed several patterns that were common to both transects, including an excess of derived alleles at high elevation, and a number of hypoxia-associated genes exhibiting a threshold effect, with a large allele frequency change only at the highest elevations. This threshold effect suggests that selection pressures may increase disproportionately at high elevations in mammals, consistent with observations of some high-elevation diseases in humans.

The newborn sheep translational model for pulmonary arterial hypertension of the neonate at high altitude

Chronic hypoxia during gestation induces greater occurrence of perinatal complications such as intrauterine growth restriction, fetal hypoxia, newborn asphyxia, and respiratory distress, among others. This condition may also cause a failure in the transition of the fetal to neonatal circulation, inducing pulmonary arterial hypertension of the neonate (PAHN), a syndrome that involves pulmonary vascular dysfunction, increased vasoconstrictor tone and pathological remodeling. As this syndrome has a relatively low prevalence in lowlands (~7 per 1000 live births) and very little is known about its prevalence and clinical evolution in highlands (above 2500 meters), our understanding is very limited. Therefore, studies on appropriate animal models have been crucial to comprehend the mechanisms underlying this pathology. Considering the strengths and weaknesses of any animal model of human disease is fundamental to achieve an effective and meaningful translation to clinical practice. The sheep model has been used to study the normal and abnormal cardiovascular development of the fetus and the neonate for almost a century. The aim of this review is to highlight the advances in our knowledge on the programming of cardiopulmonary function with the use of high-altitude newborn sheep as a translational model of PAHN.

Analysis of selection signatures reveals important insights into the adaptability of high-altitude Indian sheep breed Changthangi

Changthangi is a high-altitude sheep breed of India that is adapted to cold and hypoxic climate of Himalayas. In the present study, we analysed population structure of Changthangi and contrasted it with selected Indian and European commercial sheep breeds to detect genomic regions under positive selection. The Illumina OvineSNP50v1 genotype data on 292 animals from seven different sheep breeds i.e., Changthangi (n = 29), Garole (n = 26), Deccani (n = 24), Tibetan (n = 37), Rambouillet (n = 102) and Australian Merino (n = 50) was used. European Mouflon (n = 24) was used as an out-group for studying the stratification and phylogenetic lineage. While the principal component analysis (PCA) revealed Changthangi to cluster with Tibetan sheep; TREEMIX and ADMIXTURE results also detected the introgression of lowland Indian sheep inheritance in Changthangi. Changthangi sheep were compared with other breed groups as reference i.e., commercial (Australian Merino and Rambouillet), Indian (Deccani, Garole and Tibetan) and breeds inhabiting plains (Australian Merino, Rambouillet, Deccani and Garole). Genomic comparisons of Changthangi using cross population extended haplotype homozygosity (XP-EHH) showed multiple functional regions present on Ovis aries (Oar) chromosomes 2, 3, 6 and 18 to be under selection in Changthangi sheep. These regions were related with adaptation to climatic and hypoxic stressors, fleece characteristics and functioning of immune and reproductive systems. UCP genes, associated with adaptation to cold and hypoxic conditions, were the main loci under positive selection in Changthangi sheep population. The selection signals in Indian and European commercial sheep breeds were mainly associated with body weight and carcass traits. Furthermore, selection signals found in different comparisons were found to be part of different quantitative trait loci (QTLs) associated with important traits in different breed classes. The genes present in these regions are suitable candidates for future studies on the genetic mechanisms underlying high-altitude adaptation.

Genomics of Adaptations in Ungulates

Ungulates are a group of hoofed animals that have long interacted with humans as essential sources of food, labor, clothing, and transportation. These consist of domesticated, feral, and wild species raised in a wide range of habitats and biomes. Given the diverse and extreme environments inhabited by ungulates, unique adaptive traits are fundamental for fitness. The documentation of genes that underlie their genomic signatures of selection is crucial in this regard. The increasing availability of advanced sequencing technologies has seen the rapid growth of ungulate genomic resources, which offers an exceptional opportunity to understand their adaptive evolution. Here, we summarize the current knowledge on evolutionary genetic signatures underlying the adaptations of ungulates to different habitats.

Multiple mechanisms drive genomic adaptation to extreme O2 levels in Drosophila melanogaster

To detect the genomic mechanisms underlying evolutionary dynamics of adaptation in sexually reproducing organisms, we analyze multigenerational whole genome sequences of Drosophila melanogaster adapting to extreme O₂ conditions over an experiment conducted for nearly two decades. We develop methods to analyze time-series genomics data and predict adaptive mechanisms. Here, we report a remarkable level of synchronicity in both hard and soft selective sweeps in replicate populations as well as the arrival of favorable de novo mutations that constitute a few asynchronized sweeps. We additionally make direct experimental observations of rare recombination events that combine multiple alleles on to a single, better-adapted haplotype. Based on the analyses of the genes in genomic intervals, we provide a deeper insight into the mechanisms of genome adaptation that allow complex organisms to survive harsh environments.

The genomic details of adaptation to extreme environments remain challenging to characterize. Using new methods to analyze flies experimentally evolved to survive extreme O₂ conditions, the authors find a surprising level of synchronicity in selective sweeps, de novo mutations and adaptive recombination events.

Genome-wide comparative analyses reveal selection signatures underlying adaptation and production in Tibetan and Poll Dorset sheep

The identification of genome-wide selection signatures can provide insights on the mechanisms of natural and/or artificial selection and uncover genes related to biological functions and/or phenotypes. Tibetan sheep are an important livestock in Tibet, providing meat and wool for Tibetans who are renown for breeding livestock that adapt well to high altitudes. Using whole-genome sequences with an effective sequencing depth of 5×, we investigated the genomic diversity and structure and, identified selection signatures of White Tibetan, Oula and Poll Dorset sheep. We obtained 30,163,679 Single Nucleotide Polymorphisms (SNPs) and 5,388,372 indels benchmarked against the ovine Oar_v4.0 genome assembly. Next, using F_ST, ZHp and XP-EHH approaches, we identified selection signatures spanning a set of candidate genes, including HIF1A, CAPN3, PRKAA1, RXFP2, TRHR and HOXA10 that are associated with pathways and GO categories putatively related to hypoxia responses, meat traits and disease resistance. Candidate genes and GO terms associated with coat color were also identified. Finally, quantification of blood physiological parameters, revealed higher levels of mean corpuscular hemoglobin measurement and mean corpuscular hemoglobin concentration in Tibetan sheep compared with Poll Dorset, suggesting a greater oxygen-carrying capacity in the Tibetan sheep and thus better adaptation to high-altitude hypoxia. In conclusion, this study provides a greater understanding of genome diversity and variations associated with adaptive and production traits in sheep.

Current scenario of and future perspective for scientific research in Nepal

Scientific research can act as the 'bedrock for development' which can provide a solid foundation for the overall socioeconomic transformation in a country through invention and innovation. The inclusion of seven groups of Nepalese monuments in the United Nations Educational, Scientific and Cultural Organization (UNESCO)'s world heritage sites provides evidence of the richness in architecture and civil engineering in the country, which also show how Nepal was influenced by arts and religious philosophy. The government of Nepal (GoN) has established different scientific departments, universities, and research institutions, and has tried to emphasize the application of science and research for the development of the nation. These institutions, however, have inadequate resources, exist in a poor academic and research environment and are subject to overt political influence. Despite these various problems, a variety of pioneer research and development activities have been conducted, which show positive rays of hope. This review presents history, current situation, progress and future perspective for scientific research and development in Nepal.

Microsatellite and mitochondrial DNA analyses unveil the genetic structure of native sheep breeds from three major agro-ecological regions of India

Sheep farming has been fundamental to many civilizations in the world and is practiced in India since antiquity. Several thousand years of adaptation to local environmental conditions and selective breeding have evolved 44 sheep breeds in India. They are paramount in terms of economic, scientific, and cultural heritage. Genetic characterization information is imperative for sustainable utilization and conservation of ovine heritage. In this study, the genetic diversity, differentiation, and structure of 11 indigenous sheep breeds from three different agro-ecological zones of India were explored with genomic microsatellite loci and mitochondrial DNA (D loop). The estimated diversity parameters indicated that populations retained high levels of genetic diversity (Na = 8.27 ± 0.17; Ho = 0.65 ± 0.01), which provides an optimistic viewpoint for their survival. However, significant inbreeding was also observed in the nine populations. Moderate genetic differentiation existed among the groups (F_ST = 0.129 ± 0.012), and most likely clusters existing in the dataset are seven. Phylogenetic clustering was in line with the geographical locations of sheep populations. Mitochondrial sequences revealed high haplotype diversity with the existence of maternal haplogroups A, B, and C, and signals of population expansion. Decreased genetic diversity and unique maternal lineage (C) in endangered Tibetan and Bonpala sheep breed, warrant their immediate scientific management. Overall, the quantitative data reported here on the extant variability, and genetic relationships among the Indian sheep breeds, provide critically important inputs that will be valuable for the decision-making process on their management, both for the conservation of endangered breeds, and formulation of breeding programs to check genetic erosion.

Genomic analysis of worldwide sheep breeds reveals PDGFD as a major target of fat-tail selection in sheep

Background

Fat tail is a unique trait in sheep acquired during domestication. Several genomic analyses have been conducted in sheep breeds from limited geographic origins to identify the genetic factors underlying this trait. Nevertheless, these studies obtained different candidates. The results of these regional studies were easily biased by the breed structures.

Results

To minimize the bias and distinguish the true candidates, we used an extended data set of 968 sheep representing 18 fat-tailed breeds and 14 thin-tailed breeds from around the world, and integrated two statistical tests to detect selection signatures, including Genetic Fixation Index (F_ST) and difference of derived allele frequency (ΔDAF). The results showed that platelet derived growth factor D (PDGFD) exhibited the highest genetic differentiation between fat- and thin-tailed sheep breeds. Analysis of sequence variation identified that a 6.8-kb region within the first intron of PDGFD is likely the target of positive selection and contains regulatory mutation(s) in fat-tailed sheep. Histological and gene expression analyses demonstrated that PDGFD expression is associated with maturation and hemostasis of adipocytes. Further retrospective analysis of public transcriptomic datasets revealed that PDGFD expression is down-regulated during adipogenesis in both human and mouse, and is higher in fat tissues of obese individuals than that in lean individuals.

Conclusions

These results reveal that PDGFD is the predominant factor for the fat tail phenotype in sheep by contributing to adiopogenesis and maintaining the hemostasis of mature adipocytes. This study provides insights into the selection of fat-tailed sheep and has important application to animal breeding, as well as obesity-related human diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07210-9.

The draft genome of the blood pheasant (Ithaginis cruentus): Phylogeny and high-altitude adaptation

The blood pheasant (Ithaginis cruentus), the only species in the genus Ithaginis, lives in an extremely inhospitable high-altitude environment, coping with hypoxia and ultraviolet (UV) radiation. To further investigate the phylogeny of Phasianidae species based on complete genomes and understand the molecular genetic mechanisms of the high-altitude adaptation of the blood pheasant, we de novo assembled and annotated the complete genome of the blood pheasant. The blood pheasant genome size is 1.04 Gb with scaffold N50 of 10.88 Mb. We identified 109.92 Mb (10.62%) repetitive elements, 279,037 perfect microsatellites, and 17,209 protein-coding genes. The phylogenetic tree of Phasianidae based on whole genomes revealed three highly supported major clades with the blood pheasant included in the "erectile clade." Comparative genomics analysis showed that many genes were positively selected in the blood pheasant, which was associated with response to hypoxia and/or UV radiation. More importantly, among these positively selected genes (PSGs) which were related to high-altitude adaptation, sixteen PSGs had blood pheasant-specific missense mutations. Our data and analysis lay solid foundation to the study of Phasianidae phylogeny and provided new insights into the potential adaptation mechanisms to the high altitude employed by the blood pheasant.

Genetic diversity and genetic origin of Lanping black-boned sheep investigated by genome-wide single-nucleotide polymorphisms (SNPs)

Lanping black-boned sheep was first discovered in the 1950s in Lanping county of China and characterized by black pigmentation on skin and internal organs. Due to the novel and unique trait, the genetic background of Lanping black-boned sheep is of great interest. Here, we genotyped genome-wide SNPs (single nucleotide polymorphisms) of Lanping black-boned sheep and Lanping normal sheep using Illumina OvineSNP50 BeadChip to investigate the genetic diversity and genetic origin of Lanping black-boned sheep. We also downloaded a subset SNP dataset of two Tibet-lineage sheep breeds and four other sheep breeds from the International Sheep Genomics Consortium (ISGC) as a reference for interpreting. Lanping black-boned sheep had a lower genetic diversity level when compared to seven other sheep breeds. Principal component analysis (PCA) showed that Lanping black-boned sheep and Lanping normal sheep were clustered into the Asian group, but there was no clear separation between the two breeds. Structure analysis demonstrated a high ancestry coefficient in Lanping black-boned sheep and Lanping normal sheep. However, the two populations were separated into two distinct branches in a neighbor-joining (NJ) tree. We further evaluated the genetic divergence using population F ST , which showed that the genetic differentiation that existed between Lanping black-boned sheep and Lanping normal sheep was higher than that between Tibet sheep and Changthangi sheep, which revealed that Lanping black-boned sheep is a different breed from Lanping normal sheep on the genetic level. In addition, structure analysis and NJ tree showed that Lanping black-boned sheep had a relatively close relation with Tibet sheep. The results reported herein are a first step toward understanding the genetic background of Lanping black-boned sheep, and it will provide informative knowledge on the unique genetic resource conservation and mechanism of novel breed formation.

On the origin of European sheep as revealed by the diversity of the Balkan breeds and by optimizing population-genetic analysis tools

Background

In the Neolithic, domestic sheep migrated into Europe and subsequently spread in westerly and northwesterly directions. Reconstruction of these migrations and subsequent genetic events requires a more detailed characterization of the current phylogeographic differentiation.

Results

We collected 50 K single nucleotide polymorphism (SNP) profiles of Balkan sheep that are currently found near the major Neolithic point of entry into Europe, and combined these data with published genotypes from southwest-Asian, Mediterranean, central-European and north-European sheep and from Asian and European mouflons. We detected clines, ancestral components and admixture by using variants of common analysis tools: geography-informative supervised principal component analysis (PCA), breed-specific admixture analysis, across-breed \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{4}$$\end{document}f4 profiles and phylogenetic analysis of regional pools of breeds. The regional Balkan sheep populations exhibit considerable genetic overlap, but are clearly distinct from the breeds in surrounding regions. The Asian mouflon did not influence the differentiation of the European domestic sheep and is only distantly related to present-day sheep, including those from Iran where the mouflons were sampled. We demonstrate the occurrence, from southeast to northwest Europe, of a continuously increasing ancestral component of up to 20% contributed by the European mouflon, which is assumed to descend from the original Neolithic domesticates. The overall patterns indicate that the Balkan region and Italy served as post-domestication migration hubs, from which wool sheep reached Spain and north Italy with subsequent migrations northwards. The documented dispersal of Tarentine wool sheep during the Roman period may have been part of this process. Our results also reproduce the documented 18th century admixture of Spanish Merino sheep into several central-European breeds.

Conclusions

Our results contribute to a better understanding of the events that have created the present diversity pattern, which is relevant for the management of the genetic resources represented by the European sheep population.

The fine-scale genetic structure and selection signals of Chinese indigenous pigs

Genome-wide SNP profiling has yielded insights into the genetic structure of China indigenous pigs, but has focused on a limited number of populations. Here, we present an analysis of population structure and signals of positive selection in 42 Chinese pig populations that represent the most extensive pig phenotypic diversity in China, using genotype data of 1.1 million SNPs on customized Beadchips. This unravels the fine-scale genetic diversity, phylogenic relationships, and population structure of these populations, which shows remarkably concordance between genetic clusters and geography with few exceptions. We also reveal the genetic contribution to North Chinese pigs from European modern pigs. Furthermore, we identify possible targets of selection in the Tibetan pig, including the well-characterized hypoxia gene (EPAS1) and several previously unrecognized candidates. Intriguingly, the selected haplotype in the EPAS1 gene is associated with higher hemoglobin contents in Tibetan pigs, which is different from the protective role of EPAS1 in the high-altitude adaptation in Tibetan dogs and their owners. Additionally, we present evidence for the causality between EDNRB variants and the two-end-black (TEB) coat color phenotype in all Chinese pig populations except the Jinhua pig. We hypothesize that distinct targets have been independently selected for the formation of the TEB phenotype in Chinese pigs of different geographic origins. This highlights the importance of characterizing population-specific genetic determinants for heritable phenotype in diverse pig populations.

Selection signatures for high-altitude adaptation in ruminants

High-altitude areas are important socio-economical habitats with ruminants serving as a major source of food and commodities for humans. Living at high altitude, however, is extremely challenging, predominantly due to the exposure to hypoxic conditions, but also because of cold temperatures and limited feed for livestock. To survive in high-altitude environments over the long term, ruminants have evolved adaptation strategies, e.g. physiological and morphological modifications, which allow them to cope with these harsh conditions. Identification of such selection signatures in ruminants may contribute to more informed breeding decisions, and thus improved productivity. Moreover, studying the genetic background of altitude adaptation in ruminants provides insights into a common molecular basis across species and thus a better understanding of the physiological basis of this adaptation. In this paper, we review the major effects of high altitude on the mammalian body and highlight some of the most important short-term (coping) and genetically evolved (adaptation) physiological modifications. We then discuss the genetic architecture of altitude adaptation and target genes that show evidence of being under selection based on recent studies in various species, with a focus on ruminants. The yak is presented as an interesting native species that has adapted to the high-altitude regions of Tibet. Finally, we conclude with implications and challenges of selection signature studies on altitude adaptation in general. We found that the number of studies on genetic mechanisms that enable altitude adaptation in ruminants is growing, with a strong focus on identifying selection signatures, and hypothesise that the investigation of genetic data from multiple species and regions will contribute greatly to the understanding of the genetic basis of altitude adaptation.

Convergent evolution in human and domesticate adaptation to high-altitude environments

Humans and their domestic animals have lived and thrived in high-altitude environments worldwide for thousands of years. These populations have developed a number of adaptations to survive in a hypoxic environment, and several genomic studies have been conducted to identify the genes that drive these adaptations. Here, we discuss the various adaptations and genetic variants that have been identified as adaptive in human and domestic animal populations and the ways in which convergent evolution has occurred as these populations have adapted to high-altitude environments. We found that human and domesticate populations have adapted to hypoxic environments in similar ways. Specific genes and biological pathways have been involved in high-altitude adaptation for multiple populations, although the specific variants differ between populations. Additionally, we found that the gene EPAS1 is often a target of selection in hypoxic environments and has been involved in multiple adaptive introgression events. High-altitude environments exert strong selective pressures, and human and animal populations have evolved in convergent ways to cope with a chronic lack of oxygen.

This article is part of the theme issue ‘Convergent evolution in the genomics era: new insights and directions'.

High-density genotyping reveals signatures of selection related to acclimation and economically important traits in 15 local sheep breeds from Russia

Background

Domestication and centuries of selective breeding have changed genomes of sheep breeds to respond to environmental challenges and human needs. The genomes of local breeds, therefore, are valuable sources of genomic variants to be used to understand mechanisms of response to adaptation and artificial selection. As a step toward this we performed a high-density genotyping and comprehensive scans for signatures of selection in the genomes from 15 local sheep breeds reared across Russia.

Results

Results demonstrated that the genomes of Russian sheep breeds contain multiple regions under putative selection. More than 50% of these regions matched with intervals identified in previous scans for selective sweeps in sheep genomes. These regions contain well-known candidate genes related to morphology, adaptation, and domestication (e.g., KITLG, KIT, MITF, and MC1R), wool quality and quantity (e.g., DSG@, DSC@, and KRT@), growth and feed intake (e.g., HOXA@, HOXC@, LCORL, NCAPG, LAP3, and CCSER1), reproduction (e.g., CMTM6, HTRA1, GNAQ, UBQLN1, and IFT88), and milk-related traits (e.g., ABCG2, SPP1, ACSS1, and ACSS2). In addition, multiple genes that are putatively related to environmental adaptations were top-ranked in selected intervals (e.g., EGFR, HSPH1, NMUR1, EDNRB, PRL, TSHR, and ADAMTS5). Moreover, we observed that multiple key genes involved in human hereditary sensory and autonomic neuropathies, and genetic disorders accompanied with an inability to feel pain and environmental temperatures, were top-ranked in multiple or individual sheep breeds from Russia pointing to a possible mechanism of adaptation to harsh climatic conditions.

Conclusions

Our work represents the first comprehensive scan for signatures of selection in genomes of local sheep breeds from the Russian Federation of both European and Asian origins. We confirmed that the genomes of Russian sheep contain previously identified signatures of selection, demonstrating the robustness of our integrative approach. Multiple novel signatures of selection were found near genes which could be related to adaptation to the harsh environments of Russia. Our study forms a basis for future work on using Russian sheep genomes to spot specific genetic variants or haplotypes to be used in efforts on developing next-generation highly productive breeds, better suited to diverse Eurasian environments.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5537-0) contains supplementary material, which is available to authorized users.

A comprehensive comparison of four species of Onchidiidae provides insights on the morphological and molecular adaptations of invertebrates from shallow seas to wetlands

The Onchidiidae family is ideal for studying the evolution of marine invertebrate species from sea to wetland environments. However, comparative studies of Onchidiidae species are rare. A total of 40 samples were collected from four species (10 specimens per onchidiid), and their histological and molecular differences were systematically evaluated to elucidate the morphological foundations underlying the adaptations of these species. A histological analysis was performed to compare the structures of respiratory organs (gill, lung sac, dorsal skin) among onchidiids, and transcriptome sequencing of four representative onchidiids was performed to investigate the molecular mechanisms associated with their respective habitats. Twenty-six SNP markers of Onchidium reevesii revealed some DNA polymorphisms determining visible traits. Non-muscle myosin heavy chain II (NMHC II) and myosin heavy chain (MyHC), which play essential roles in amphibian developmental processes, were found to be differentially expressed in different onchidiids and tissues. The species with higher terrestrial ability and increased integrated expression of Os-MHC (NMHC II gene) and the MyHC gene, illustrating that the expression levels of these genes were associated with the evolutionary degree. This study provides a comprehensive analysis of the adaptions of a diverse and widespread group of invertebrates, the Onchidiidae. Some onchidiids can breathe well through gills and skin when under seawater, and some can breathe well through lung sacs and skin when in wetlands. A histological comparison of respiratory organs and the relative expression levels of two genes provided insights into the adaptions of onchidiids that allowed their transition from shallow seas to wetlands. This work provides a valuable reference and might encourage further study.

Genomic footprints of dryland stress adaptation in Egyptian fat-tail sheep and their divergence from East African and western Asia cohorts

African indigenous sheep are classified as fat-tail, thin-tail and fat-rump hair sheep. The fat-tail are well adapted to dryland environments, but little is known on their genome profiles. We analyzed patterns of genomic variation by genotyping, with the Ovine SNP50K microarray, 394 individuals from five populations of fat-tail sheep from a desert environment in Egypt. Comparative inferences with other East African and western Asia fat-tail and European sheep, reveal at least two phylogeographically distinct genepools of fat-tail sheep in Africa that differ from the European genepool, suggesting separate evolutionary and breeding history. We identified 24 candidate selection sweep regions, spanning 172 potentially novel and known genes, which are enriched with genes underpinning dryland adaptation physiology. In particular, we found selection sweeps spanning genes and/or pathways associated with metabolism; response to stress, ultraviolet radiation, oxidative stress and DNA damage repair; activation of immune response; regulation of reproduction, organ function and development, body size and morphology, skin and hair pigmentation, and keratinization. Our findings provide insights on the complexity of genome architecture regarding dryland stress adaptation in the fat-tail sheep and showcase the indigenous stocks as appropriate genotypes for adaptation planning to sustain livestock production and human livelihoods, under future climates.