Whole-Genome Sequencing of Native Sheep Provides Insights into Rapid Adaptations to Extreme Environments

A time-resolved multi-omics atlas of transcriptional regulation in response to high-altitude hypoxia across whole-body tissues

High-altitude hypoxia acclimatization requires whole-body physiological regulation in highland immigrants, but the underlying genetic mechanism has not been clarified. Here we use sheep as an animal model for low-to-high altitude translocation. We generate multi-omics data including whole-genome sequences, time-resolved bulk RNA-Seq, ATAC-Seq and single-cell RNA-Seq from multiple tissues as well as phenotypic data from 20 bio-indicators. We characterize transcriptional changes of all genes in each tissue, and examine multi-tissue temporal dynamics and transcriptional interactions among genes. Particularly, we identify critical functional genes regulating the short response to hypoxia in each tissue (e.g., PARG in the cerebellum and HMOX1 in the colon). We further identify TAD-constrained cis-regulatory elements, which suppress the transcriptional activity of most genes under hypoxia. Phenotypic and transcriptional evidence indicate that antenatal hypoxia could improve hypoxia tolerance in offspring. Furthermore, we provide time-series expression data of candidate genes associated with human mountain sickness (e.g., BMPR2) and high-altitude adaptation (e.g., HIF1A). Our study provides valuable resources and insights for future hypoxia-related studies in mammals.

Elucidating the Role of circTIAM1 in Guangling Large-Tailed Sheep Adipocyte Proliferation and Differentiation via the miR-485-3p/PLCB1 Pathway

Fat tissue-a vital energy storage organ-is intricately regulated by various factors, including circular RNA, which plays a significant role in modulating fat development and lipid metabolism. Therefore, this study aims to clarify the regulatory mechanism of sheep adipocyte proliferation and differentiation by investigating the involvement of circTIAM1, miR-485-3p, and its target gene PLCB1. Through previous sequencing data, circTIAM1 was identified in sheep adipocytes, with its circularization mechanism elucidated, confirming its cytoplasmic localization. Experimental evidence from RNase R treatment and transcription inhibitors highlighted that circTIAM1 is more stable than linear RNA. Additionally, circTIAM1 promoted sheep adipocyte proliferation and differentiation. Furthermore, bioinformatic analysis demonstrated a robust interaction between miR-485-3p and circTIAM1. Further experiments revealed that miR-485-3p inhibits fat cell proliferation and differentiation by inhibiting PLCB1, with circTIAM1 alleviating the inhibitory effect via competitive binding. In summary, our findings elucidate the mechanism through which circTIAM1 regulates Guangling Large-Tailed sheep adipocyte proliferation and differentiation via the miR-485-3p-PLCB1 pathway, offering a novel perspective for further exploring fat metabolism regulation.

Ancient mitogenomes from Pre-Pottery Neolithic Central Anatolia and the effects of a Late Neolithic bottleneck in sheep (Ovis aries)

Occupied between ~10,300 and 9300 years ago, the Pre-Pottery Neolithic site of Aşıklı Höyük in Central Anatolia went through early phases of sheep domestication. Analysis of 629 mitochondrial genomes from this and numerous sites in Anatolia, southwest Asia, Europe, and Africa produced a phylogenetic tree with excessive coalescences (nodes) around the Neolithic, a potential signature of a domestication bottleneck. This is consistent with archeological evidence of sheep management at Aşıklı Höyük which transitioned from residential stabling to open pasturing over a millennium of site occupation. However, unexpectedly, we detected high genetic diversity throughout Aşıklı Höyük's occupation rather than a bottleneck. Instead, we detected a tenfold demographic bottleneck later in the Neolithic, which caused the fixation of mitochondrial haplogroup B in southwestern Anatolia. The mitochondrial genetic makeup that emerged was carried from the core region of early Neolithic sheep management into Europe and dominates the matrilineal diversity of both its ancient and the billion-strong modern sheep populations.

Whole-genome resequencing of Chinese indigenous sheep provides insight into the genetic basis underlying climate adaptation

BACKGROUND

Chinese indigenous sheep are valuable resources with unique features and characteristics. They are distributed across regions with different climates in mainland China; however, few reports have analyzed the environmental adaptability of sheep based on their genome. We examined the variants and signatures of selection involved in adaptation to extreme humidity, altitude, and temperature conditions in 173 sheep genomes from 41 phenotypically and geographically representative Chinese indigenous sheep breeds to characterize the genetic basis underlying environmental adaptation in these populations.

RESULTS

Based on the analysis of population structure, we inferred that Chinese indigenous sheep are divided into four groups: Kazakh (KAZ), Mongolian (MON), Tibetan (TIB), and Yunnan (YUN). We also detected a set of candidate genes that are relevant to adaptation to extreme environmental conditions, such as drought-prone regions (TBXT, TG, and HOXA1), high-altitude regions (DYSF, EPAS1, JAZF1, PDGFD, and NF1) and warm-temperature regions (TSHR, ABCD4, and TEX11). Among all these candidate genes, eight ABCD4, CNTN4, DOCK10, LOC105608545, LOC121816479, SEM3A, SVIL, and TSHR overlap between extreme environmental conditions. The TSHR gene shows a strong signature for positive selection in the warm-temperature group and harbors a single nucleotide polymorphism (SNP) missense mutation located between positions 90,600,001 and 90,650,001 on chromosome 7, which leads to a change in the protein structure of TSHR and influences its stability.

CONCLUSIONS

Analysis of the signatures of selection uncovered genes that are likely related to environmental adaptation and a SNP missense mutation in the TSHR gene that affects the protein structure and stability. It also provides information on the evolution of the phylogeographic structure of Chinese indigenous sheep populations. These results provide important genetic resources for future breeding studies and new perspectives on how animals can adapt to climate change.

Multi-omics data elucidate parasite-host-microbiota interactions and resistance to Haemonchus contortus in sheep

BACKGROUND

The integration of molecular data from hosts, parasites, and microbiota can enhance our understanding of the complex biological interactions underlying the resistance of hosts to parasites. Haemonchus contortus, the predominant sheep gastrointestinal parasite species in the tropics, causes significant production and economic losses, which are further compounded by the diminishing efficiency of chemical control owing to anthelmintic resistance. Knowledge of how the host responds to infection and how the parasite, in combination with microbiota, modulates host immunity can guide selection decisions to breed animals with improved parasite resistance. This understanding will help refine management practices and advance the development of new therapeutics for long-term helminth control.

METHODS

Eggs per gram (EPG) of feces were obtained from Morada Nova sheep subjected to two artificial infections with H. contortus and used as a proxy to select animals with high resistance or susceptibility for transcriptome sequencing (RNA-seq) of the abomasum and 50 K single-nucleotide genotyping. Additionally, RNA-seq data for H. contortus were generated, and amplicon sequence variants (ASV) were obtained using polymerase chain reaction amplification and sequencing of bacterial and archaeal 16S ribosomal RNA genes from sheep feces and rumen content.

RESULTS

The heritability estimate for EPG was 0.12. GAST, GNLY, IL13, MGRN1, FGF14, and RORC genes and transcripts were differentially expressed between resistant and susceptible animals. A genome-wide association study identified regions on chromosomes 2 and 11 that harbor candidate genes for resistance, immune response, body weight, and adaptation. Trans-expression quantitative trait loci were found between significant variants and differentially expressed transcripts. Functional co-expression modules based on sheep genes and ASVs correlated with resistance to H. contortus, showing enrichment in pathways of response to bacteria, immune and inflammatory responses, and hub features of the Christensenellaceae, Bacteroides, and Methanobrevibacter genera; Prevotellaceae family; and Verrucomicrobiota phylum. In H. contortus, some mitochondrial, collagen-, and cuticle-related genes were expressed only in parasites isolated from susceptible sheep.

CONCLUSIONS

The present study identified chromosome regions, genes, transcripts, and pathways involved in the elaborate interactions between the sheep host, its gastrointestinal microbiota, and the H. contortus parasite. These findings will assist in the development of animal selection strategies for parasite resistance and interdisciplinary approaches to control H. contortus infection in sheep.

Multiomics Analyses Provide New Insight into Genetic Variation of Reproductive Adaptability in Tibetan Sheep

Domestication and artificial selection during production-oriented breeding have greatly shaped the level of genomic variability in sheep. However, the genetic variation associated with increased reproduction remains elusive. Here, two groups of samples from consecutively monotocous and polytocous sheep were collected for genome-wide association, transcriptomic, proteomic, and metabolomic analyses to explore the genetic variation in fecundity in Tibetan sheep. Genome-wide association study revealed strong associations between BMPR1B (p.Q249R) and litter size, as well as between PAPPA and lambing interval; these findings were validated in 1,130 individuals. Furthermore, we constructed the first single-cell atlas of Tibetan sheep ovary tissues and identified a specific mural granulosa cell subtype with PAPPA-specific expression and differential expression of BMPR1B between the two groups. Bulk RNA-seq indicated that BMPR1B and PAPPA expressions were similar between the two groups of sheep. 3D protein structure prediction and coimmunoprecipitation analysis indicated that mutation and mutually exclusive exons of BMPR1B are the main mechanisms for prolific Tibetan sheep. We propose that PAPPA is a key gene for stimulating ovarian follicular growth and development, and steroidogenesis. Our work reveals the genetic variation in reproductive performance in Tibetan sheep, providing insights and valuable genetic resources for the discovery of genes and regulatory mechanisms that improve reproductive success.

Skin transcriptomic analysis reveals candidate genes and pathways associated with thermotolerance in hair sheep

The skin plays an important role in thermoregulation. Identification of genes on the skin that contribute to increased heat tolerance can be used to select animals with the best performance in warm environments. Our objective was to identify candidate genes associated with the heat stress response in the skin of Santa Ines sheep. A group of 80 sheep assessed for thermotolerance was kept in a climatic chamber for 8 days at a stress level temperature of 36 °C (10 am to 04 pm) and a maintenance temperature of 28 °C (04 pm to 10 am). Two divergent groups, with seven animals each, were formed after ranking them by thermotolerance using rectal temperature. From skin biopsy samples, total RNA was extracted, quantified, and used for RNA-seq analysis. 15,989 genes were expressed in sheep skin samples, of which 4 genes were differentially expressed (DE; FDR < 0.05) and 11 DE (FDR 0.05-0.177) between the two divergent groups. These genes are involved in cellular protection against stress (HSPA1A and HSPA6), ribosome assembly (28S, 18S, and 5S ribosomal RNA), and immune response (IGHG4, GNLY, CXCL1, CAPN14, and SAA-4). The candidate genes and main pathways related to heat tolerance in Santa Ines sheep require further investigation to understand their response to heat stress in different climatic conditions and under solar radiation. It is essential to verify whether these genes and pathways are present in different breeds and to understand the relationship between heat stress and other genes identified in this study.

Characterization of Complete Mitochondrial Genome of Badri Breed of Bos indicus (Bovidae: Bovinae): Selection Pressure and Comparative Analysis

High-altitude mammals are often subject to specific environmental obstacles, which exert selective pressure on their physiological and morphological traits, hence driving their evolutionary processes. It is anticipated that these circumstances will lead to the adaptive evolution of protein-coding genes (PCGs) in the mitochondrial genome, which play a crucial role in the oxidative phosphorylation system. In this study, we have generated the complete mitochondrial genome of the Badri breed of Bos indicus inhabiting a high-altitude environment to test the signatures of adaptive evolution on PCGs and their phylogenetic relationships. The complete mitogenome of the Badri breed is 16,339 bp and most tRNAs showed typical clover-leaf secondary structure with a few exceptions, like trnS1 and trnS2 without DHU arm and trnK without DHU loop. Comparative analysis of PCGs indicated that cox1 is the most conserved, while atp6 is the most variable gene. Moreover, the ratios of non-synonymous to synonymous substitution rates indicated the purifying selection (Ka/Ks < 1) in the protein-coding genes that shape the diversity in mitogenome of Bos indicus. Furthermore, Branch-site model (BSM) suggested that cox1, cox2, nad3, nad4L, and nad6 underwent stronger purifying selection (ω < 1) than other PCGs in 15 breeds of 4 species, including Badri. BSM also detected 10 positive sites in PCGs and one in 13 PCGs concatenated dataset. Additional analyses in Datamonkey indicated 11 positive sites and 23 purifying sites in the concatenated dataset, a relaxation of selection strength in nad3, and no evidence of episodic diversifying selection in any PCGs. Phylogeny revealed the sister relationship of the Badri with other breeds of Bos indicus as well as Bos frontalis (Gayal-2). The mitogenome of the Badri breed is an important genomic resource for conservation genetics of this species and also contributes to the understanding of the adaptive evolution of mitochondrial protein coding genes.

Characterization of runs of homozygosity islands in American mink using whole-genome sequencing data

The genome-wide analysis of runs of homozygosity (ROH) islands can be an effective strategy for identifying shared variants within a population and uncovering important genomic regions related to complex traits. The current study performed ROH analysis to characterize the genome-wide patterns of homozygosity, identify ROH islands and annotated genes within these candidate regions using whole-genome sequencing data from 100 American mink (Neogale vison). After sequence processing, variants were called using GATK and Samtools pipelines. Subsequent to quality control, 8,373,854 bi-allelic variants identified by both pipelines remained for further analysis. A total of 34,652 ROH segments were identified in all individuals, among which shorter segments (0.3-1 Mb) were abundant throughout the genome, approximately accounting for 84.39% of all ROH. Within these segments, we identified 63 ROH islands housing 156 annotated genes. The genes located in ROH islands were associated with fur quality (EDNRA, FGF2, FOXA2 and SLC24A4), body size/weight (MYLK4, PRIM2, FABP2, EYS and PHF3), immune capacity (IL2, IL21, PTP4A1, SEMA4C, JAK2, CCNA2 and TNIP3) and reproduction (ADAD1, KHDRBS2, INSL6, PGRMC2 and HSPA4L). Furthermore, Gene Ontology and KEGG pathway enrichment analyses revealed 56 and 9 significant terms (FDR-corrected p-value < 0.05), respectively, among which cGMP-PKG signalling pathway, regulation of actin cytoskeleton, and calcium signalling pathway were highlighted due to their functional roles in growth and fur characteristics. This is the first study to present ROH islands in American mink. The candidate genes from ROH islands and functional enrichment analysis suggest possible signatures of selection in response to the mink breeding targets, such as increased body length, reproductive performance and fur quality. These findings contribute to our understanding of genetic characteristics, and provide complementary information to assist with implementation of breeding strategies for genetic improvement in American mink.

Whole-genome resequencing provides insights into the diversity and adaptation to desert environment in Xinjiang Mongolian cattle

BACKGROUND

Xinjiang Mongolian cattle is an indigenous breed that inhabits the Taklimakan Desert and is characterized by its small body size. However, the genomic diversity, origin, and genetic basis underlying the adaptation to the desert environment have been poorly studied.

RESULTS

We analyzed patterns of Xinjiang Mongolian cattle genetic variation by sequencing 20 genomes together with seven previously sequenced genomes and comparing them to the 134 genomes of nine representative breeds worldwide. Among the breeds of Bos taurus, we found the highest nucleotide diversity (0.0024) associated with the lower inbreeding coefficient (2.0110-6), the lowest linkage disequilibrium (r2 = 0.3889 at distance of 10 kb), and the highest effective population size (181 at 20 generations ago) in Xinjiang Mongolian cattle. The genomic diversity pattern could be explained by a limited introgression of Bos indicus genes. More importantly, similarly to desert-adapted camel and same-habitat sheep, we also identified signatures of selection including genes, GO terms, and/or KEGG pathways controlling water reabsorption and osmoregulation, metabolic regulation and energy balance, as well as small body size in Xinjiang Mongolian cattle.

CONCLUSIONS

Our results imply that Xinjiang Mongolian cattle might have acquired distinct genomic diversity by virtue of the introgression of Bos indicus, which helps understand the demographic history. The identification of selection signatures can provide novel insights into the genomic basis underlying the adaptation of Xinjiang Mongolian cattle to the desert environment.

The genomics of adaptation to climate in European great tit (Parus major) populations

The recognition that climate change is occurring at an unprecedented rate means that there is increased urgency in understanding how organisms can adapt to a changing environment. Wild great tit (Parus major) populations represent an attractive ecological model system to understand the genomics of climate adaptation. They are widely distributed across Eurasia and they have been documented to respond to climate change. We performed a Bayesian genome-environment analysis, by combining local climate data with single nucleotide polymorphisms genotype data from 20 European populations (broadly spanning the species' continental range). We found 36 genes putatively linked to adaptation to climate. Following an enrichment analysis of biological process Gene Ontology (GO) terms, we identified over-represented terms and pathways among the candidate genes. Because many different genes and GO terms are associated with climate variables, it seems likely that climate adaptation is polygenic and genetically complex. Our findings also suggest that geographical climate adaptation has been occurring since great tits left their Southern European refugia at the end of the last ice age. Finally, we show that substantial climate-associated genetic variation remains, which will be essential for adaptation to future changes.

Revelation of genetic diversity and genomic footprints of adaptation in Indian pig breeds

In the present study, the genetic diversity measures among four Indian domestic breeds of pig namely Agonda Goan, Ghurrah, Ghungroo, and Nicobari, of different agro-climatic regions of country were explored and compared with European commercial breeds, European wild boar and Chinese domestic breeds. The double digest restriction site-associated DNA sequencing (ddRADseq) data of Indian pigs (102) and Landrace (10 animals) were generated and whole genome sequencing data of exotic pigs (60 animals) from public data repository were used in the study. The principal component analysis (PCA), admixture analysis and phylogenetic analysis revealed that Indian breeds were closer in ancestry to Chinese breeds than European breeds. European breeds exhibited highest genetic diversity measures among all the considered breeds. Among Indian breeds, Agonda Goan and Ghurrah were found to be more genetically diverse than Nicobari and Ghungroo. The selection signature regions in Indian pigs were explored using iHS and XP-EHH, and during iHS analysis, it was observed that genes related to growth, reproduction, health, meat quality, sensory perception and behavior were found to be under selection pressure in Indian pig breeds. Strong selection signatures were recorded in 24.25-25.25 Mb region of SSC18, 123.25-124 Mb region of SSC15 and 118.75-119.5 Mb region of SSC2 in most of the Indian breeds upon pairwise comparison with European commercial breeds using XP-EHH. These regions were harboring some important genes such as EPHA4 for thermotolerance, TAS2R16, FEZF1, CADPS2 and PTPRZ1 for adaptability to scavenging system of rearing, TRIM36 and PGGT1B for disease resistance and CCDC112, PIAS1, FEM1B and ITGA11 for reproduction.

Whole-Genome Resequencing Reveals Genetic Diversity and Wool Trait-Related Genes in Liangshan Semi-Fine-Wool Sheep

Understanding the genetic makeup of local sheep breeds is essential for their scientific conservation and sustainable utilization. The Liangshan semi-fine-wool sheep (LSS), a Chinese semi-fine-wool breed renowned for its soft wool, was analyzed using whole-genome sequencing data including 35 LSS, 84 sheep from other domestic breeds, and 20 Asiatic mouflons. We investigated the genetic composition of LSS by conducting analyses of the population structure, runs of homozygosity, genomic inbreeding coefficients, and selection signature. Our findings indicated that LSS shares greater genetic similarity with Border Leicester and Romney sheep than with Tibetan (TIB), Yunnan (YNS), and Chinese Merino sheep. Genomic analysis indicated low to moderate inbreeding coefficients, ranging from 0.014 to 0.154. In identifying selection signals across the LSS genome, we pinpointed 195 candidate regions housing 74 annotated genes (e.g., IRF2BP2, BVES, and ALOX5). We also found the overlaps between the candidate regions and several known quantitative trait loci related to wool traits, such as the wool staple length and wool fiber diameter. A selective sweep region, marked by the highest value of cross-population extended haplotype homozygosity, encompassed IRF2BP2-an influential candidate gene affecting fleece fiber traits. Furthermore, notable differences in genotype frequency at a mutation site (c.1051 + 46T > C, Chr25: 6,784,190 bp) within IRF2BP2 were observed between LSS and TIB and YNS sheep (Fisher's exact test, p < 2.2 × 10-16). Taken together, these findings offer insights crucial for the conservation and breeding enhancement of LSS.

When the tap runs dry: The multi-tissue gene expression and physiological responses of water deprived Peromyscus eremicus

The harsh and dry conditions of desert environments have resulted in genomic adaptations, allowing for desert organisms to withstand prolonged drought, extreme temperatures, and limited food resources. Here, we present a comprehensive exploration of gene expression across five tissues (kidney, liver, lung, gastrointestinal tract, and hypothalamus) and 19 phenotypic measurements to explore the whole-organism physiological and genomic response to water deprivation in the desert-adapted cactus mouse (Peromyscus eremicus). The findings encompass the identification of differentially expressed genes and correlative analysis between phenotypes and gene expression patterns across multiple tissues. Specifically, we found robust activation of the vasopressin renin-angiotensin-aldosterone system (RAAS) pathways, whose primary function is to manage water and solute balance. Animals reduce food intake during water deprivation, and upregulation of PCK1 highlights the adaptive response to reduced oral intake via its actions aimed at maintained serum glucose levels. Even with such responses to maintain water balance, hemoconcentration still occurred, prompting a protective downregulation of genes responsible for the production of clotting factors while simultaneously enhancing angiogenesis which is thought to maintains tissue perfusion. In this study, we elucidate the complex mechanisms involved in water balance in the desert-adapted cactus mouse, P. eremicus. By prioritizing a comprehensive analysis of whole-organism physiology and multi-tissue gene expression in a simulated desert environment, we describe the complex and successful response of regulatory processes.

Genetic Diversity, Population Structure and Selection Signature in Begait Goats Revealed by Whole-Genome Sequencing

Goats belong to a group of animals called small ruminants and are critical sources of livelihood for rural people. Genomic sequencing can provide information ranging from basic knowledge about goat diversity and evolutionary processes that shape genomes to functional information about genes/genomic regions. In this study, we exploited a whole-genome sequencing data set to analyze the genetic diversity, population structure and selection signatures of 44 individuals belonging to 5 Ethiopian goat populations: 12 Aberegalle (AB), 5 Afar (AF), 11 Begait (BG), 12 Central highlands (CH) and 5 Meafure (MR) goats. Our results revealed the highest genetic diversity in the BG goat population compared to the other goat populations. The pairwise genetic differentiation (FST) among the populations varied and ranged from 0.011 to 0.182, with the closest pairwise value (0.003) observed between the AB and CH goats and a distant correlation (FST = 0.182) between the BG and AB goats, indicating low to moderate genetic differentiation. Phylogenetic tree, ADMIXTURE and principal component analyses revealed a classification of the five Ethiopian goat breeds in accordance with their geographic distribution. We also found three top genomic regions that were detected under selection on chromosomes 2, 5 and 13. Moreover, this study identified different candidate genes related to milk characteristics (GLYCAM1 and SRC), carcass (ZNF385B, BMP-7, PDE1B, PPP1R1A, FTO and MYOT) and adaptive and immune response genes (MAPK13, MAPK14, SCN7A, IL12A, EST1 DEFB116 and DEFB119). In conclusion, this information could be helpful for understanding the genetic diversity and population structure and selection scanning of these important indigenous goats for future genetic improvement and/or as an intervention mechanism.

Pan-Omics in Sheep: Unveiling Genetic Landscapes

Multi-omics-integrated analysis, known as panomics, represents an advanced methodology that harnesses various high-throughput technologies encompassing genomics, epigenomics, transcriptomics, proteomics, and metabolomics. Sheep, playing a pivotal role in agricultural sectors due to their substantial economic importance, have witnessed remarkable advancements in genetic breeding through the amalgamation of multiomics analyses, particularly with the evolution of high-throughput technologies. This integrative approach has established a robust theoretical foundation, enabling a deeper understanding of sheep genetics and fostering improvements in breeding strategies. The comprehensive insights obtained through this approach shed light on diverse facets of sheep development, including growth, reproduction, disease resistance, and the quality of livestock products. This review primarily focuses on the application of principal omics analysis technologies in sheep, emphasizing correlation studies between multiomics data and specific traits such as meat quality, wool characteristics, and reproductive features. Additionally, this paper anticipates forthcoming trends and potential developments in this field.

IBD sharing patterns as intra-breed admixture indicators in small ruminants

In this study, we investigated how IBD patterns shared between individuals of the same breed could be informative of its admixture level, with the underlying assumption that the most admixed breeds, i.e. the least genetically isolated, should have a much more fragmented genome. We considered 111 goat breeds (i.e. 2501 individuals) and 156 sheep breeds (i.e. 3304 individuals) from Europe, Africa and Asia, for which beadchip SNP genotypes had been performed. We inferred the breed's level of admixture from: (i) the proportion of the genome shared by breed's members (i.e. "genetic integrity level" assessed from ADMIXTURE software analyses), and (ii) the "AV index" (calculated from Reynolds' genetic distances), used as a proxy for the "genetic distinctiveness". In both goat and sheep datasets, the statistical analyses (comparison of means, Spearman correlations, LM and GAM models) revealed that the most genetically isolated breeds, also showed IBD profiles made up of more shared IBD segments, which were also longer. These results pave the way for further research that could lead to the development of admixture indicators, based on the characterization of intra-breed shared IBD segments, particularly effective as they would be independent of the knowledge of the whole genetic landscape in which the breeds evolve. Finally, by highlighting the fragmentation experienced by the genomes subjected to crossbreeding carried out over the last few generations, the study reminds us of the need to preserve local breeds and the integrity of their adaptive architectures that have been shaped over the centuries.

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.

Landscape genomics reveals adaptive divergence of indigenous sheep in different ecological environments of Xinjiang, China

Sheep are important livestock animals that have evolved under various ecological pressures. Xinjiang is a region with diverse and harsh environments that have shaped many local sheep breeds with unique characteristics and environmental adaptability. However, these breeds are losing ecological flexibility due to the promotion of intensive farming practices. Here we sequenced 14 local sheep breeds from Xinjiang and analyzed their genetic structure and gene flow with other sheep breeds from neighboring regions. The Tibetan Plateau was the geographic origin of Xinjiang native sheep evolution. We performed genome-environment association analysis and identified Bio9: Mean Temperature of Driest Quarter and Bio15: Precipitation Seasonality as the key environmental factors affecting Xinjiang local sheep and the key genes involved in their survival and adaptation. We classified Xinjiang native sheep breeds into six groups based on their differential genes by pairwise selective sweep analysis and Community Network Analysis. We analyzed transcriptome expression data of 832 sheep tissues and detected tissue-specific enrichment of six group-specific genes in different biological systems. Our results revealed the genetic basis of year-round estrus, drought tolerance, hypoxia resistance, and cold tolerance traits of Xinjiang sheep breeds. Moreover, we proposed conservation strategies for Xinjiang local sheep breeds and provided theoretical guidance for breeding new sheep breeds under global extreme environments.

The impacts of thermal heterogeneity across microhabitats on post-settlement selection of intertidal mussels

Rapid genetic selection is critical for allowing natural populations to adapt to different thermal environments such as those that occur across intertidal microhabitats with high degrees of thermal heterogeneity. To address the question of how thermal regimes influence selection and adaptation in the intertidal black mussel Mytilisepta virgata, we continuously recorded environmental temperatures in both tidal pools and emergent rock microhabitats and then assessed genetic differentiation, gene expression patterns, RNA editing level, and cardiac performance. Our results showed that the subpopulations in the tidal pool and on emergent rocks had different genetic structures and exhibited different physiological and molecular responses to high-temperature stress. These results indicate that environmental heterogeneity across microhabitats is important for driving genetic differentiation and shed light on the importance of post-settlement selection for adaptively modifying the genetic composition and thermal responses of these intertidal mussels.

Morphophysiological Responses of the Goat Mammary Gland to Water Scarcity in Arid and Semi-Arid Environments: Are They Enough to Generate Adaptation to New Climatic Challenges?

Due to climate change, diverse territories of the planet will suffer from water restrictions. Goats are perceived as the most resilient ruminants in this scenario. So, various studies have focused on describing how a lower water intake influences milk production, especially in breeds adapted to desert environments. In water-stress situations, goats lose up to 32% of their body weight (BW), the rate of passage is reduced, and the digestibility of the feed increases. When goats consume water again, the rumen prevents hemolysis and osmotic shock from occurring. Regarding milk production, the response varies depending on the breed and the level of water restriction, maintaining the milk volume or reducing it by up to 41%. Systemically, it decreases the urinary volume and glomerular filtration rate, increasing blood osmolality and the vasopressin (ADH) concentration. Studies are scarce regarding changes in blood flow to the mammary gland, but there would be a reduction in blood flow velocity of up to 40% without changing blood pressure. New studies must be undertaken to determine which breeds or crosses are the best adapted to changing environmental conditions and to improve our understanding of the changes that occur at the morphophysiological level of the caprine mammary gland.

Genome-wide association study for reproduction traits in Colombian Creole Blanco Orejinegro cattle

The profitability of the beef cattle production system relies heavily on reproductive traits. Unfortunately, certain traits, such as age at first calving (AFC), calving interval (CI), and gestation length (GL), can pose challenges in traditional breeding programs because of their low heritability (0.01-0.12) and sex-limited characteristics. Another important aspect is the conservation of the genetic resources of animals adapted to the Colombian regions, which implies the preservation and rational use of the creole breeds in the country market. Therefore, this study aimed to identify genomic regions in the creole cattle breed Blanco Orejinegro (BON) that influence the reproductive traits in females. The dataset comprised 439 animals and 118,116 single-nucleotide polymorphisms' (SNPs) markers. The GS3 program was used to identify the SNP effects employing the BAYES Cπ methodology. The number of SNPs with effect for AFC was 25, 1527 for CI, and 23 for GL. Some of the genes found associated with reproductive and growth traits as well as immune response and environmental adaptation ECE1, EPH, EPHB2, SMARCAL1, IGFBP5, IGFBP2, FCGRT, EGFR, MUL1, PINK1, STPG1, CNGB1, TGFB1, OXTR, IL22RA1, MYOM3, OXTR, CNR2, HIVEP3, CTPS1, CXCL8, FCGRT, MREG, TMEM169, PECR, and MC1R. Our results evidenced a high contribution of the genetic architecture of the Colombian creole cattle breed Blanco Orejinegro that may impact should be included in implementing genetic improvement and conservation programs.

A whole genome scan reveals distinct features of selection in Zhaotong cattle of Yunnan province

Over the years, indigenous cattle have not only played an essential role in securing primary food sources but have also been utilized for labor by humans, making them invaluable genetic resources. The Zhaotong cattle, a native Chinese breed from the Yunnan province, possess excellent meat quality and resistance to heat and humidity. Here we used whole genome sequencing data of 104 animals to delve into the population structure, genomic diversity and potential positive selection signals in Zhaotong cattle. The findings of this study demonstrate that the genetic composition of Zhaotong cattle was primarily derived from Chinese indicine cattle and East Asian cattle. The nucleotide diversity of Zhaotong cattle was only lower than that of Chinese indicine cattle, which was much higher than that of other taurine cattle. Genome-wide selection scans detected a series of positive candidate regions containing multiple key genes related to bone development and metabolism (CA10, GABRG3, GLDN and NOTUM), meat quality traits (ALG8, LINGO2, MYO5B, PRKG1 and GABRB1), immune response (ADA2, BMF, LEF1 and PAK6) and heat resistance (EIF2AK4 and LEF1). In summary, this study supplies essential genetic insights into the genome diversity within Zhaotong cattle and provides a foundational framework for comprehending the genetic basis of indigenous cattle breeds.

Whole-genome resequencing of Dorper and Hu sheep to reveal selection signatures associated with important traits

Dorper and Hu sheep exhibit different characteristics in terms of reproduction, growth, and meat quality. Comparison of the genomes of two breeds help to reveal important genomic information. In this study, whole genome resequencing of 30 individuals (Dorper, DB and Hu sheep, HY) identified 15,108,125 single nucleotide polymorphisms (SNPs). Population differentiation (Fst) and cross population extended haplotype homozygosity (XP-EHH) were performed for selective signal analysis. In total, 106 and 515 overlapped genes were present in both the Fst results and XP-EHH results in HY vs DB and in DB vs HY, respectively. In HY vs DB, 106 genes were enriched in 12 GO terms and 83 KEGG pathways, such as ATP binding (GO:0005524) and PI3K-Akt signaling pathway (oas04151). In DB vs HY, 515 genes were enriched in 109 GO terms and 215 KEGG pathways, such as skeletal muscle cell differentiation (GO:0035914) and MAPK signaling pathway (oas04010). According to the annotation results, we identified a series of candidate genes associated with reproduction (UNC5C, BMPR1B, and GLIS1), meat quality (MECOM, MEF2C, and MYF6), and immunity (GMDS, GALK1, and ITGB4). Our investigation has uncovered genomic information for important traits in sheep and provided a basis for subsequent studies of related traits.

Recent advances in the genomic resources for sheep

Sheep (Ovis aries) provide a vital source of protein and fibre to human populations. In coming decades, as the pressures associated with rapidly changing climates increase, breeding sheep sustainably as well as producing enough protein to feed a growing human population will pose a considerable challenge for sheep production across the globe. High quality reference genomes and other genomic resources can help to meet these challenges by: (1) informing breeding programmes by adding a priori information about the genome, (2) providing tools such as pangenomes for characterising and conserving global genetic diversity, and (3) improving our understanding of fundamental biology using the power of genomic information to link cell, tissue and whole animal scale knowledge. In this review we describe recent advances in the genomic resources available for sheep, discuss how these might help to meet future challenges for sheep production, and provide some insight into what the future might hold.

When the tap runs dry: the physiological effects of acute experimental dehydration in Peromyscus eremicus

Desert organisms have evolved physiological, biochemical and genomic mechanisms to survive the extreme aridity of desert environments. Studying desert-adapted species provides a unique opportunity to investigate the survival strategies employed by organisms in some of the harshest habitats on Earth. Two of the primary challenges faced in desert environments are maintaining water balance and thermoregulation. We collected data in a simulated desert environment and a captive colony of cactus mice (Peromyscus eremicus) and used lab-based experiments with real time physiological measurements; energy expenditure, water loss rate and respiratory exchange rate, to characterize the response to water deprivation. Mice without access to water had significantly lower energy expenditures and in turn, reduced water loss compared to mice with access to water after the first 24 h of the experiment. Additionally, we observed significant mass loss that is probably due to dehydration-associated anorexia a response to limit fluid loss by reducing waste and the solute load as well as allowing water reabsorption from the kidneys and gastrointestinal tract. Finally, we observed body temperature correlated with sex, with males without access to water maintaining body temperature when compared with hydrated males, whereas body temperature decreased for females without access to water, suggesting daily metabolic depression in females.

Global genetic diversity, introgression, and evolutionary adaptation of indicine cattle revealed by whole genome sequencing

Indicine cattle, also referred to as zebu (Bos taurus indicus), play a central role in pastoral communities across a wide range of agro-ecosystems, from extremely hot semiarid regions to hot humid tropical regions. However, their adaptive genetic changes following their dispersal into East Asia from the Indian subcontinent have remained poorly documented. Here, we characterize their global genetic diversity using high-quality whole-genome sequencing data from 354 indicine cattle of 57 breeds/populations, including major indicine phylogeographic groups worldwide. We reveal their probable migration into East Asia was along a coastal route rather than inland routes and we detected introgression from other bovine species. Genomic regions carrying morphology-, immune-, and heat-tolerance-related genes underwent divergent selection according to Asian agro-ecologies. We identify distinct sets of loci that contain promising candidate variants for adaptation to hot semi-arid and hot humid tropical ecosystems. Our results indicate that the rapid and successful adaptation of East Asian indicine cattle to hot humid environments was promoted by localized introgression from banteng and/or gaur. Our findings provide insights into the history and environmental adaptation of indicine cattle.

Whole Genome Resequencing Reveals Selection Signals Related to Wool Color in Sheep

Wool color is controlled by a variety of genes. Although the gene regulation of some wool colors has been studied in relative depth, there may still be unknown genetic variants and control genes for some colors or different breeds of wool that need to be identified and recognized by whole genome resequencing. Therefore, we used whole genome resequencing data to compare and analyze sheep populations of different breeds by population differentiation index and nucleotide diversity ratios (Fst and θπ ratio) as well as extended haplotype purity between populations (XP-EHH) to reveal selection signals related to wool coloration in sheep. Screening in the non-white wool color group (G1 vs. G2) yielded 365 candidate genes, among which PDE4B, GMDS, GATA1, RCOR1, MAPK4, SLC36A1, and PPP3CA were associated with the formation of non-white wool; an enrichment analysis of the candidate genes yielded 21 significant GO terms and 49 significant KEGG pathways (p < 0.05), among which 17 GO terms and 21 KEGG pathways were associated with the formation of non-white wool. Screening in the white wool color group (G2 vs. G1) yielded 214 candidate genes, including ABCD4, VSX2, ITCH, NNT, POLA1, IGF1R, HOXA10, and DAO, which were associated with the formation of white wool; an enrichment analysis of the candidate genes revealed 9 significant GO-enriched pathways and 19 significant KEGG pathways (p < 0.05), including 5 GO terms and 12 KEGG pathways associated with the formation of white wool. In addition to furthering our understanding of wool color genetics, this research is important for breeding purposes.

Genome-wide identification of candidate copy number polymorphism genes associated with complex traits of Tibetan-sheep

Copy number variation (CNV) is a genetic structural polymorphism important for phenotypic diversity and important economic traits of livestock breeds, and it plays an important role in the desired genetic variation. This study used whole genome sequencing to detect the CNV variation in the genome of 6 local Tibetan sheep groups. We detected 69,166 CNV events and 7230 copy number variable regions (CNVRs) after merging the overlapping CNVs, accounting for 2.72% of the reference genome. The CNVR length detected ranged from 1.1 to 1693.5 Kb, with a total length of 118.69 Mb and an average length of 16.42 Kb per CNVR. Functional GO cluster analysis showed that the CNVR genes were mainly involved in sensory perception systems, response to stimulus, and signal transduction. Through CNVR-based Vst analysis, we found that the CACNA2D3 and CTBP1 genes related to hypoxia adaptation, the HTR1A gene related to coat color, and the TRNAS-GGA and PIK3C3 genes related to body weight were all strongly selected. The findings of our study will contribute novel insights into the genetic structural variation underlying hypoxia adaptation and economically important traits in Tibetan sheep.

Whole-Genome Resequencing Reveals Selection Signal Related to Sheep Wool Fineness

Wool fineness affects the quality of wool, and some studies have identified about forty candidate genes that affect sheep wool fineness, but these genes often reveal only a certain proportion of the variation in wool thickness. We further explore additional genes associated with the fineness of sheep wool. Whole-genome resequencing of eight sheep breeds was performed to reveal selection signals associated with wool fineness, including four coarse wool and four fine/semi-fine wool sheep breeds. Multiple methods to reveal selection signals (Fst and θπ Ratio and XP-EHH) were applied for sheep wool fineness traits. In total, 269 and 319 genes were annotated in the fine wool (F vs. C) group and the coarse wool (C vs. F) group, such as LGR4, PIK3CA, and SEMA3C and NFIB, OPHN1, and THADA. In F vs. C, 269 genes were enriched in 15 significant GO Terms (p < 0.05) and 38 significant KEGG Pathways (p < 0.05), such as protein localization to plasma membrane (GO: 0072659) and Inositol phosphate metabolism (oas 00562). In C vs. F, 319 genes were enriched in 21 GO Terms (p < 0.05) and 16 KEGG Pathways (p < 0.05), such as negative regulation of focal adhesion assembly (GO: 0051895) and Axon guidance (oas 04360). Our study has uncovered genomic information pertaining to significant traits in sheep and has identified valuable candidate genes. This will pave the way for subsequent investigations into related traits.

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.

Whole-genome sequencing reveals adaptations of hairy-footed jerboas (Dipus, Dipodidae) to diverse desert environments

BACKGROUND

Environmental conditions vary among deserts across the world, spanning from hyper-arid to high-elevation deserts. However, prior genomic studies on desert adaptation have focused on desert and non-desert comparisons overlooking the complexity of conditions within deserts. Focusing on the adaptation mechanisms to diverse desert environments will advance our understanding of how species adapt to extreme desert environments. The hairy-footed jerboas are well adapted to diverse desert environments, inhabiting high-altitude arid regions, hyper-arid deserts, and semi-deserts, but the genetic basis of their adaptation to different deserts remains unknown.

RESULTS

Here, we sequenced the whole genome of 83 hairy-footed jerboas from distinct desert zones in China to assess how they responded under contrasting conditions. Population genomics analyses reveal the existence of three species in hairy-footed jerboas distributed in China: Dipus deasyi, Dipus sagitta, and Dipus sowerbyi. Analyses of selection between high-altitude desert (elevation ≥ 3000m) and low-altitude desert (< 500m) populations identified two strongly selected genes, ATR and HIF1AN, associated with intense UV radiation and hypoxia in high-altitude environments. A number of candidate genes involved in energy and water homeostasis were detected in the comparative genomic analyses of hyper-arid desert (average annual precipitation < 70mm) and arid desert (< 200mm) populations versus semi-desert (> 360mm) populations. Hyper-arid desert animals also exhibited stronger adaptive selection in energy homeostasis, suggesting water and resource scarcity may be the main drivers of desert adaptation in hairy-footed jerboas.

CONCLUSIONS

Our study challenges the view of deserts as homogeneous environments and shows that distinct genomic adaptations can be found among desert animals depending on their habitats.

Whole-genome resequencing of the native sheep provides insights into the microevolution and identifies genes associated with reproduction traits

BACKGROUND

Sheep genomes undergo numerous genes losses, gains and mutation that generates genome variability among breeds of the same species after long time natural and artificial selection. However, the microevolution of native sheep in northwest China remains elusive. Our aim was to compare the genomes and relevant reproductive traits of four sheep breeds from different climatic environments, to unveil the selection challenges that this species cope with, and the microevolutionary differences in sheep genomes. Here, we resequenced the genomes of 4 representative sheep breeds in northwest China, including Kazakh sheep and Duolang sheep of native breeds, and Hu sheep and Suffolk sheep of exotic breeds with different reproductive characteristics.

RESULTS

We found that these four breeds had a similar expansion experience from ~ 10,000 to 1,000,000 years ago. In the past 10,000 years, the selection intensity of the four breeds was inconsistent, resulting in differences in reproductive traits. We explored the sheep variome and selection signatures by F_ST and θπ. The genomic regions containing genes associated with different reproductive traits that may be potential targets for breeding and selection were detected. Furthermore, non-synonymous mutations in a set of plausible candidate genes and significant differences in their allele frequency distributions across breeds with different reproductive characteristics were found. We identified PAK1, CYP19A1 and PER1 as a likely causal gene for seasonal reproduction in native sheep through qPCR, Western blot and ELISA analyses. Also, the haplotype frequencies of 3 tested gene regions related to reproduction were significantly different among four sheep breeds.

CONCLUSIONS

Our results provide insights into the microevolution of native sheep and valuable genomic information for identifying genes associated with important reproductive traits in sheep.

When the tap runs dry: The physiological effects of acute experimental dehydration in Peromyscus eremicus

Desert organisms have evolved physiological, biochemical, and genomic mechanisms to survive the extreme aridity of desert environments. Studying desert-adapted species provides a unique opportunity to investigate the survival strategies employed by organisms in some of the harshest habitats on Earth. Two of the primary challenges faced in desert environments are maintaining water balance and thermoregulation. We collected data in a simulated desert environment and a captive colony of cactus mice (Peromyscus eremicus) and used lab-based experiments with real time physiological measurements to characterize the response to water-deprivation. Mice without access to water had significantly lower energy expenditures and in turn, reduced water loss compared to mice with access to water after the first 24 hours of the experiment. Additionally, we observed significant weight loss likely related to dehydration-associated anorexia a response to limit fluid loss by reducing waste and the solute load as well as allowing water reabsorption from the kidneys and gastrointestinal tract. Finally, we observed body temperature correlated with sex, with males without access to water maintaining body temperature when compared to hydrated males while body temperature decreased for females without access to water compared to hydrated, suggesting daily torpor in females.

Genomic Landscape of Copy Number Variations and Their Associations with Climatic Variables in the World's Sheep

Sheep show characteristics of phenotypic diversity and adaptation to diverse climatic regions. Previous studies indicated associations between copy number variations (CNVs) and climate-driven adaptive evolution in humans and other domestic animals. Here, we constructed a genomic landscape of CNVs (n = 39,145) in 47 old autochthonous populations genotyped at a set of high-density (600 K) SNPs to detect environment-driven signatures of CNVs using a multivariate regression model. We found 136 deletions and 52 duplications that were significantly (Padj. < 0.05) associated with climatic variables. These climate-mediated selective CNVs are involved in functional candidate genes for heat stress and cold climate adaptation (e.g., B3GNTL1, UBE2L3, and TRAF2), coat and wool-related traits (e.g., TMEM9, STRA6, RASGRP2, and PLA2G3), repairing damaged DNA (e.g., HTT), GTPase activity (e.g., COPG), fast metabolism (e.g., LMF2 and LPIN3), fertility and reproduction (e.g., SLC19A1 and CCDC155), growth-related traits (e.g., ADRM1 and IGFALS), and immune response (e.g., BEGAIN and RNF121) in sheep. In particular, we identified significant (Padj. < 0.05) associations between probes in deleted/duplicated CNVs and solar radiation. Enrichment analysis of the gene sets among all the CNVs revealed significant (Padj. < 0.05) enriched gene ontology terms and pathways related to functions such as nucleotide, protein complex, and GTPase activity. Additionally, we observed overlapping between the CNVs and 140 known sheep QTLs. Our findings imply that CNVs can serve as genomic markers for the selection of sheep adapted to specific climatic conditions.

Comparative genomic analysis uncovers candidate genes related with milk production and adaptive traits in goat breeds

During the process of animal domestication, both natural and artificial selection cause variation in allele frequencies among populations. Identifying genomic areas of selection in domestic animals may aid in the detection of genomic areas linked to ecological and economic traits. We studied genomic variation in 140 worldwide goat individuals, including 75 Asian, 30 African and 35 European goats. We further carried out comparative population genomics to detect genomic regions under selection for adaptability to harsh conditions in local Asian ecotypes and also milk production traits in European commercial breeds. In addition, we estimated the genetic distances among 140 goat individuals. The results showed that among all studied goat groups, local breeds from West and South Asia emerged as an independent group. Our search for selection signatures in local goats from West and South Asia revealed candidate genes related to adaptation to hot climate (HSPB6, HSF4, VPS13A and NBEA genes) and immune response (IL7, IL5, IL23A and LRFN5) traits. Furthermore, selection signatures in European commercial goats involved several milk production related genes, such as VPS13C, NCAM2, TMPRSS15, CSN3 and ABCG2. The identified candidate genes could be the fundamental genetic resource for enhancement of goat production and environmental-adaptive traits, and as such they should be used in goat breeding programs to select more efficient breeds.

Resequencing of Rosa rugosa accessions revealed the history of population dynamics, breed origin, and domestication pathways

BACKGROUND

Rosa rugosa is a shrub that originated in China and has economic and ecological value. However, during the development of R. rugosa, the genetic background was chaotic, and the genetic structure among different wild populations was unclear, as well as wild and cultivated accessions. Here, we report whole-genome resequencing of wild and cultivated R. rugosa accessions.

RESULTS

A total of 19,041,284 SNPs were identified in 188 R. rugosa accessions and 3 R. chinensis accessions by resequencing. Population genetic analysis revealed that cultivated and wild groups were separated very early. All R. rugosa accessions were divided into 8 categories based on genetic structure: (1) Weihai, Yantai, and Liaoning category, (2) Jilin category, and (3) Hammonasset category (above three are wild); (4) traditional varieties, (5) hybrids between R. rugosa and R. chinensis, (6) Zizhi Rose, (7) Kushui Rose, (8) hybrids between R. rugosa and R. multiflora. We found that the heterozygosity and genetic diversity of wild accessions were generally lower than those of cultivated individuals. The genes that were selected during cultivation were identified, and it was found that these genes were mainly related to environmental adaptation and growth.

CONCLUSIONS

The Jilin population was the oldest population and later migrated to Liaoning and then migrated to Yantai and Weihai by sea regression in the Bohai Basin. The Hammonasset naturalized population probably originated from the Jilin population and then experienced separate differentiation. The long-term asexual reproduction pattern of R. rugosa decreased genetic diversity in the wild population. During R. rugosa cultivation, the ancestors of the Jilin population were involved in breeding traditional varieties, after which almost no wild individuals were engaged in breeding. However, in recent decades, cross breeding of R. rugosa started the utilization of wild germplasms. In comparison, some other species play important roles in variety formation. Few genes related to economic traits were selected, suggesting no directional domestication in the R. rugosa cultivation process.

Whole genome sequencing of simmental cattle for SNP and CNV discovery

BACKGROUD

The single nucleotide polymorphisms (SNPs) and copy number variations (CNVs) are two major genomic variants, which play crucial roles in evolutionary and phenotypic diversity.

RESULTS

In this study, we performed a comprehensive analysis to explore the genetic variations (SNPs and CNVs) of high sperm motility (HSM) and poor sperm motility (PSM) Simmental bulls using the high-coverage (25×) short-read next generation sequencing and single-molecule long reads sequencing data. A total of ~ 15 million SNPs and 2,944 CNV regions (CNVRs) were detected in Simmental bulls, and a set of positive selected genes (PSGs) and CNVRs were found to be overlapped with quantitative trait loci (QTLs) involving immunity, muscle development, reproduction, etc. In addition, we detected two new variants in LEPR, which may be related to the artificial breeding to improve important economic traits. Moreover, a set of genes and pathways functionally related to male fertility were identified. Remarkably, a CNV on SPAG16 (chr2:101,427,468 - 101,429,883) was completely deleted in all poor sperm motility (PSM) bulls and half of the bulls in high sperm motility (HSM), which may play a crucial role in the bull-fertility.

CONCLUSIONS

In conclusion, this study provides a valuable genetic variation resource for the cattle breeding and selection programs.

Genome-wide analysis emancipates genomic diversity and signature of selection in Altay white-headed cattle of Xinjiang, China

Altay white-headed cattle have not received enough attention for several reasons. Due to irrational breeding and selection practices, the number of pure Altay white-headed cattle has decreased significantly and the breed is now on the eve of extinction. The genomic characterization will be a crucial step towards understanding the genetic basis of productivity and adaptability to survival under native Chinese agropastoral systems; nevertheless, no attempt has been made in Altay white-headed cattle. In the current study, we compared the genomes of 20 Altay white-headed cattle to the genomes of 144 individuals in representative breeds. Population genetic diversity revealed that the nucleotide diversity of Altay white-headed cattle was less than that of indicine breeds and comparable to that of Chinese taurus cattle. Using population structure analysis, we also found that Altay white-headed cattle carried the ancestry of the European and East Asian cattle lineage. In addition, we used three different methods (FST, θπ ratio and XP-EHH) to investigate the adaptability and white-headed phenotype of Altay white-headed cattle and compared it with Bohai black cattle. We found EPB41L5, SCG5 and KIT genes on the list of the top one percent genes, these genes might have an association with environmental adaptability and the white-headed phenotype for this breed. Our research reveals the distinctive genomic features of Altay white-headed cattle at the genome-wide level.

Effects of New Mutations in BMPRIB, GDF9, BMP15, LEPR, and B4GALNT2 Genes on Litter Size in Sheep

Prolificacy is a crucial characteristic of livestock, particularly for species such as sheep that have many births. The objectives of this study were as follows: (1) to investigate the genetic diversity of the 13 new and 7 known variants in the BMPRIB, GDF9, BMP15, LEPR, and B4GALNT2 genes in Ujimqin (UM), the F1 population of Dorper × Ujimqin crossbred (DPU), the F1 population of Suffolk × Ujimqin crossbred (SFKU), Sonid sheep (SN), Tan sheep (Tan), Hu sheep (Hu), and Small-tailed Han sheep (STH) sheep breeds/populations; (2) to perform an association analysis of the above 20 variants with litter size in 325 UM, 304 DPU, and 66 SFKU sheep populations; (3) to compare the frequencies of the litter-size-related alleles of these 20 variants among 8 sheep breeds/populations (the above seven sheep breeds + Mongolia sheep breed). With the use of the Sequenom MassARRAY®SNP assay technology, these 20 mutations were genotyped. The association analysis results showed that the c.746A>G (FecB) mutation in BMPR1B was significantly associated with the litter size of UM and DPU, the c.994A>G (FecGA) in GDF9 was significantly associated with the litter size of SFKU, and the c.31_33CTTinsdel (B1) in BMP15 was significantly associated with the litter size of UM. Our findings might provide valuable genetic markers for expanding sheep litter sizes.

Deciphering local adaptation of native Indian cattle (Bos indicus) breeds using landscape genomics and in-silico prediction of deleterious SNP effects on protein structure and function

India has 50 registered breeds of native cattle (Bos indicus) which are locally adapted to diverse environmental conditions. This study aimed to investigate the genomic basis of adaptation of native Indian cattle and to predict the impact of key SNPs on the amino acid changes that affect protein function. The Illumina 777 K BovineHD BeadChip was used to genotype 178 native cattle belonging to contrasting landscapes and agro-climatic conditions. The genotype-environment association was investigated with R. SamBada, using 5,74,382 QC passed SNPs and 11 predictor variables (10 multi-collinearity controlled environmental variables and 1 variable as "score of PCA" on ancestry coefficients of individuals). In total, 1,12,780 models were selected as significant (q < 0.05) based on G score. The pathway ontology of the annotated genes revealed many important pathways and genes having a direct and indirect role in cold and hot adaptation. Only ten SNP variants had a SIFT score of < 0.05 (deleterious), and only two of them, each lying in the genes CRYBA1 and USP18, were predicted to be deleterious with high confidence. RaptorX predicted the tertiary structures of proteins encoded by wild and mutant variants of these genes. The quality of the models was determined using Ramachandran plots and RaptorX parameters, indicating that they are accurate. RaptorX and I-Mutant 2.0 softwares revealed significant differences among wild and mutant proteins. Adaptive alleles identified in the present investigation might be responsible for the local adaptation of these cattle breeds.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03493-3.

Long divergent haplotypes introgressed from wild sheep are associated with distinct morphological and adaptive characteristics in domestic sheep

The worldwide sheep population comprises more than 1000 breeds. Together, these exhibit a considerable morphological diversity, which has not been extensively investigated at the molecular level. Here, we analyze whole-genome sequencing individuals of 1,098 domestic sheep from 154 breeds, and 69 wild sheep from seven Ovis species. On average, we detected 6.8%, 1.0% and 0.2% introgressed sequence in domestic sheep originating from Iranian mouflon, urial and argali, respectively, with rare introgressions from other wild species. Interestingly, several introgressed haplotypes contributed to the morphological differentiations across sheep breeds, such as a RXFP2 haplotype from Iranian mouflon conferring the spiral horn trait, a MSRB3 haplotype from argali strongly associated with ear morphology, and a VPS13B haplotype probably originating from urial and mouflon possibly associated with facial traits. Our results reveal that introgression events from wild Ovis species contributed to the high rate of morphological differentiation in sheep breeds, but also to individual variation within breeds. We propose that long divergent haplotypes are a ubiquitous source of phenotypic variation that allows adaptation to a variable environment, and that these remain intact in the receiving population probably due to reduced recombination.

Selection signature analysis reveals RDH5 performed key function in vision during sheep domestication process

Abstract. As one of the most successful domesticated animals in the Neolithic age, sheep gradually migrated all over the world with human activities. During the domestication process, remarkable changes have taken place in morphology, physiology, and behavior, resulting in different breeds with different characters via artificial and natural selection. However, the genetic background responsible for these phenotypic variations remains largely unclear. Here, we used whole genome resequencing technology to compare and analyze the genome differences between Asiatic mouflon wild sheep (Ovis orientalis) and Hu sheep (Ovis aries). A total of 755 genes were positively selected in the process of domestication and selection, and the genes related to sensory perception had directional evolution in the autosomal region, such as OPRL1, LEF1, TAS1R3, ATF6, VSX2, MYO1A, RDH5, and some novel genes. A missense mutation of c.T722C/p.M241T in exon 4 of RDH5 existing in sheep were found, and the T allele was completely fixed in Hu sheep. In addition, the mutation with the C allele reduced the retinol dehydrogenase activity encoding by RDH5, which can impair retinoic acid metabolism and further influenced the visual cycle. Overall, our results showed significant enrichment for positively selected genes involved in sensory perception development during sheep domestication; RDH5 and its variants may be related to the retinal degeneration in sheep. We infer that the wild sheep ancestors with weaker visual sensitivity were weeded out by humans, and the mutation was selective, swept by the dual pressures of natural and artificial selection.

Chromosome-level genome assembly and population genomics of Mongolian racerunner (Eremias argus) provide insights into high-altitude adaptation in lizards

BACKGROUND

Although the extreme environmental adaptation of organisms is a hot topic in evolutionary biology, genetic adaptation to high-altitude environment remains poorly characterized in ectothermic animals. Squamates are among the most diverse terrestrial vertebrates, with tremendous ecological plasticity and karyotype diversity, and are a unique model system to investigate the genetic footprints of adaptation.

RESULTS

We report the first chromosome-level assembly of the Mongolian racerunner (Eremias argus) and our comparative genomics analyses found that multiple chromosome fissions/fusions events are unique to lizards. We further sequenced the genomes of 61 Mongolian racerunner individuals that were collected from altitudes ranging from ~ 80 to ~ 2600 m above sea level (m.a.s.l.). Population genomic analyses revealed many novel genomic regions under strong selective sweeps in populations endemic to high altitudes. Genes embedded in those genomic regions are mainly associated with energy metabolism and DNA damage repair pathways. Moreover, we identified and validated two substitutions of PHF14 that may enhance the lizards' tolerance to hypoxia at high altitudes.

CONCLUSIONS

Our study reveals the molecular mechanism of high-altitude adaptation in ectothermic animal using lizard as a research subject and provides a high-quality lizard genomic resource for future research.

Whole-body adipose tissue multi-omic analyses in sheep reveal molecular mechanisms underlying local adaptation to extreme environments

The fat tail of sheep is an important organ that has evolved to adapt to extreme environments. However, the genetic mechanisms underlying the fat tail phenotype remain poorly understood. Here, we characterize transcriptome and lipidome profiles and morphological changes in 250 adipose tissues from two thin-tailed and three fat-tailed sheep populations in summer and winter. We implement whole-genome selective sweep tests to identify genetic variants related to fat-tails. We identify a set of functional genes that show differential expression in the tail fat of fat-tailed and thin-tailed sheep in summer and winter. These genes are significantly enriched in pathways, such as lipid metabolism, extracellular matrix (ECM) remodeling, molecular transport, and inflammatory response. In contrast to thin-tailed sheep, tail fat from fat-tailed sheep show slighter changes in adipocyte size, ECM remodeling, and lipid metabolism, and had less inflammation in response to seasonal changes, indicating improved homeostasis. Whole-genome selective sweep tests identify genes involved in preadipocyte commitment (e.g., BMP2, PDGFD) and terminal adipogenic differentiation (e.g., VEGFA), which could contribute to enhanced adipocyte hyperplasia. Altogether, we establish a model of regulatory networks regulating adipose homeostasis in sheep tails. These findings improve our understanding of how adipose homeostasis is maintained, in response to extreme environments in animals.

Identification of candidate genes related to highland adaptation from multiple Chinese local chicken breeds by whole genome sequencing analysis

Understanding the genetic mechanism of highland adaptation is of great importance for breeding improvement of Tibetan chickens (TBC). Some studies of TBC have identified some candidate genes and pathways from multiple subgroups, but the related genetic mechanisms remain largely unknown. Different genetic backgrounds and the independent genetic basis of highland adaptation make it difficult to identity the selective region of highland adaptation with all TBC samples. In this study, we conducted pre-analysis in a large-scale population to select a TBC subgroup with the purest and highest level of highland-specific lineage for the further analysis. Finally, the 37 samples from a TBC subgroup and 19 Lahsa White chickens were used to represent the highland group for further analysis with 80 samples from five Chinese local lowland breeds as controls. Population structure analysis revealed that highland adaptation significantly affected population stratification in Chinese local chicken breeds. Genome-wide selection signal analysis identified 201 candidate genes associated with highland adaptation of TBC, and these genes were significantly enriched in calcium signaling, vascular smooth muscle contraction and the cellular response to oxidative stress pathways. Additionally, we identified a narrow 1.76 kb region containing an overlapping region between HBZ and an active enhancer, and our identified region showed a highly significant signal. The highland group selected the haplotype with high activity to improve the oxygen-carrying capacity, thus being adapted to a hypoxic environment. We also found that STX2 was significantly selected in the highland group, thus potentially reducing the oxidative stress caused by hypoxia, and that STX2 exhibited the opposite effects on highland adaptation and reproductive traits. Our findings advance our understanding of extreme environment adaptation of highland chickens, and provide some variants and genes beneficial to TBC genetic breeding improvement.

Molecular markers for thermo-tolerance are associated with reproductive and physiological traits in Pelibuey ewes raised in a semiarid environment

Pelibuey sheep exhibit reproductive activity through the year, but warm weather lowers their fertility and demonstrates physiological limitations of environmental heat stress. Single nucleotide polymorphisms (SNPs) associated with heat stress tolerance in sheep have been reported previously. The objective was to validate the association of seven thermo-tolerance SNP markers with reproductive and physiological traits in Pelibuey ewes raised in a semiarid region. Pelibuey ewes were assigned to a cool (January 1st^.- March 31st^.; n = 101) or warm (April 1st^.- August 31st^.; n = 104) experimental group. All ewes were exposed to fertile rams and assessed for pregnancy diagnosis 90 days later; lambing day was reported at birth. These data served to calculate the reproductive traits of services per conception, prolificacy, days to estrus, days to conception, conception rate and lambing rate. Rectal temperature, rump/leg skin temperature and respiratory rate were measured and reported as physiological traits. Blood samples were collected and processed to extract DNA, which was genotyped using the TaqMan allelic discrimination method and qPCR. A mixed effects statistical model was used to validate associations between SNP genotypes and phenotypic traits. The SNPs rs421873172, rs417581105 and rs407804467 were confirmed as markers associated with reproductive and physiological traits (P < 0.05), and these SNPs were in the genes PAM, STAT1 and FBXO11, respectively. Interestingly, these SNP markers resulted as predictors for the evaluated traits but only in ewes from the warm group, which indicated their association with heat-stress tolerance. An additive SNP effect was confirmed with the highest contribution (P < 0.01) of the SNP rs417581105 for the evaluated traits. Reproductive performance improved (P < 0.05) and physiological parameters decreased in ewes carrying favorable SNP genotypes. In conclusion, three thermo-tolerance SNP markers were associated with improved reproductive and physiological traits in a prospective population of heat-stressed ewes raised in a semiarid environment.

Whole genome sequencing revealed genetic diversity, population structure, and selective signature of Panou Tibetan sheep

BACKGROUND

The detection of selective traits in different populations can not only reveal current mechanisms of artificial selection for breeding, but also provide new insights into phenotypic variation in new varieties and the search for genes associated with important traits. Panou sheep is a cultivated breed of Tibetan sheep in China with stable genetic performance, consistent appearance and fast growth and development after decades of artificial selection and cultivation. Due to long-term adaptation to the high altitude, cold and hypoxic environment in the plateau area, they may have formed a unique gene pool that is different from other Tibetan sheep breeds. To explore the genetic resources of Panou sheep, we used next-generation sequencing technology for the first time to investigate the genome-wide population structure, genetic diversity, and candidate signatures of positive selection in Panou sheep.

RESULTS

Comparative genomic analysis with the closely related species Oula sheep (a native breed of Tibetan sheep in China) was used to screen the population selection signal of Panou sheep. Principal component analysis and neighbor joining tree showed that Panou sheep and Oula sheep had differences in population differentiation. Furthermore, analyses of population structure, they came from the same ancestor, and when K = 2, the two populations could be distinguished. Panou sheep exhibit genetic diversity comparable to Oula sheep, as shown by observed heterozygosity, expected heterozygosity and runs of homozygosity. Genome-wide scanning using the Fst and π ratio methods revealed a list of potentially selected related genes in Panou sheep compared to Oula sheep, including histone deacetylase 9 (HDAC9), protein tyrosine kinase 2 (PTK2), microphthalmia-related transcription factor (MITF), vesicular amine transporter 1 (VAT1), trichohyalin-like 1 (TCHHL1), amine oxidase, copper containing 3 (AOC3), interferon-inducible protein 35 (IFI35).

CONCLUSIONS

The results suggest that traits related to growth and development and plateau adaptation may be selection targets for the domestication and breeding improvement of Tibetan sheep. This study provides the fundamental footprints for Panou sheep breeding and management.

Similar adaptative mechanism but divergent demographic history of four sympatric desert rodents in Eurasian inland

Phenotypes associated with metabolism and water retention are thought to be key to the adaptation of desert species. However, knowledge on the genetic changes and selective regimes on the similar and divergent ways to desert adaptation in sympatric and phylogenetically close desert organisms remains limited. Here, we generate a chromosome level genome assembly for Northern three-toed jerboa (Dipus sagitta) and three other high-quality genome assemblies for Siberian jerboa (Orientallactaga sibirica), Midday jird (Meriones meridianus), and Desert hamster (Phodopus roborovskii). Genomic analyses unveil that desert adaptation of the four species mainly result from similar metabolic pathways, such as arachidonic acid metabolism, thermogenesis, oxidative phosphorylation, insulin related pathway, DNA repair and protein synthesis and degradation. However, the specific evolved genes in the same adaptative molecular pathway often differ in the four species. We also reveal similar niche selection but different demographic histories and sensitivity to climate changes, which may be related to the diversified genomic adaptative features. In addition, our study suggests that nocturnal rodents have evolved some specific adaptative mechanism to desert environments compared to large desert animals. Our genomic resources will provide an important foundation for further research on desert genetic adaptations.