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

Meat quality differences and protein molecular mechanisms affecting meat flavor in different breeds of Tibetan sheep analyzed by 4D label-free quantitative proteomics

To evaluate the meat quality of the new breed of Panou sheep, the longissimus dorsi (LD) muscles of 1.5-year-old Panou sheep and the local breed of Oula sheep were selected for comparative analysis in terms of meat quality, and the molecular mechanisms influencing flavor were investigated using 4D label-free proteomics technology. The results revealed that the fiber density, tenderness, and brightness of the Panou sheep meat were lower than those of the Oula sheep, and the composition of amino acids and flavor substances made it possible to determine that the Panou sheep meat has a high-quality and distinctive flavor. Proteomic analysis indicated that the metabolic pathways that may be associated with meat flavor are amino acid catabolism and sugar metabolism. This study explored the role of proteins in the regulation of meat flavor in Tibetan sheep, which provides a reference for the identification of meat products and subsequent breed improvement.

Detecting the footprint of selection on the genome of Girgentana goat, a popular ancient breed

Girgentana goats are an ancient breed with distinctive morphological, adaptive, and production traits, making this population an interesting model for studying the genetic architecture underlying these traits. These special features result from natural and human-mediated selection. In this study, we aimed to detect potential signatures of selection in the Girgentana genome by combining the following statistical methods: the integrated haplotype score (iHS), the standardised log-ratio of the integrated site-specific extended haplotype homozygosity test between pairs of populations (Rsb), the runs of homozygosity (ROH) islands and the population differentiation index (FST). A composite dataset of 206 Girgentana and 334 animals from 13 goat populations across Northern and Southern Italy was analysed. All animals were genotyped using an Illumina Goat medium-density BeadChip. Multidimensional scaling and neighbour-joining analyses revealed a clear separation of the three major clades, coinciding with Girgentana, Northern, and Southern Italian goats. Twelve regions putatively under selection were detected using iHS and Rsb, whereas 16 hotspot regions were identified using FST and ROH. Notably, a candidate region on chromosome 01 was consistently identified in all four tests. This region, along with other candidate regions, includes several genes associated with adaptive immunity, reproduction, and body size traits. The Girgentana breed showed signals of ongoing selection in a region of chromosome 6 encompassing several milk quality genes, such as caseins (CSN2, CSN1S2, and CSN3). Our study provides a glimpse into the genomic regions harbouring genes that presumably affect the desired features of Girgentana. This highlights the importance of ancient breeds in providing essential genetic traits for adapting livestock to increasing climate change challenges.

Genetic advancements and future directions in ruminant livestock breeding: from reference genomes to multiomics innovations

Ruminant livestock provide a rich source of products, such as meat, milk, and wool, and play a critical role in global food security and nutrition. Over the past few decades, genomic studies of ruminant livestock have provided valuable insights into their domestication and the genetic basis of economically important traits, facilitating the breeding of elite varieties. In this review, we summarize the main advancements for domestic ruminants in reference genome assemblies, population genomics, and the identification of functional genes or variants for phenotypic traits. These traits include meat and carcass quality, reproduction, milk production, feed efficiency, wool and cashmere yield, horn development, tail type, coat color, environmental adaptation, and disease resistance. Functional genomic research is entering a new era with the advancements of graphical pangenomics and telomere-to-telomere (T2T) gap-free genome assembly. These advancements promise to improve our understanding of domestication and the molecular mechanisms underlying economically important traits in ruminant livestock. Finally, we provide new perspectives and future directions for genomic research on ruminant genomes. We suggest how ever-increasing multiomics datasets will facilitate future studies and molecular breeding in livestock, including the potential to uncover novel genetic mechanisms underlying phenotypic traits, to enable more accurate genomic prediction models, and to accelerate genetic improvement programs.

Potential Genetic Markers Associated with Environmental Adaptability in Herbivorous Livestock

Herbivorous livestock, such as cattle, sheep, goats, horses, and donkeys, play a crucial role in agricultural production and possess remarkable resilience to extreme environmental conditions, driven by complex genetic mechanisms. Recent advancements in high-throughput sequencing, genome assembly, and environmental data integration have enabled a deeper understanding of the genetic basis of their environmental adaptation. This review identifies key genes associated with high-altitude, heat, cold, and drought adaptation, providing insights into the molecular mechanisms underlying these traits. By elucidating these genetic adaptations, our study aims to support conservation efforts, inform selective breeding programs, and enhance agricultural productivity, ultimately contributing to sustainable livestock farming and economic benefits for farmers.

Telomere-to-telomere sheep genome assembly identifies variants associated with wool fineness

Ongoing efforts to improve sheep reference genome assemblies still leave many gaps and incomplete regions, resulting in a few common failures and errors in genomic studies. Here, we report a 2.85-Gb gap-free telomere-to-telomere genome of a ram (T2T-sheep1.0), including all autosomes and the X and Y chromosomes. This genome adds 220.05 Mb of previously unresolved regions and 754 new genes to the most updated reference assembly ARS-UI_Ramb_v3.0; it contains four types of repeat units (SatI, SatII, SatIII and CenY) in centromeric regions. T2T-sheep1.0 has a base accuracy of more than 99.999%, corrects several structural errors in previous reference assemblies and improves structural variant detection in repetitive sequences. Alignment of whole-genome short-read sequences of global domestic and wild sheep against T2T-sheep1.0 identifies 2,664,979 new single-nucleotide polymorphisms in previously unresolved regions, which improves the population genetic analyses and detection of selective signals for domestication (for example, ABCC4) and wool fineness (for example, FOXQ1).

Whole-genome resequencing to investigate the genetic diversity and mechanisms of plateau adaptation in Tibetan sheep

INTRODUCTION

Tibetan sheep, economically important animals on the Qinghai-Tibet Plateau, have diversified into numerous local breeds with unique characteristics through prolonged environmental adaptation and selective breeding. However, most current research focuses on one or two breeds, and lacks a comprehensive representation of the genetic diversity across multiple Tibetan sheep breeds. This study aims to fill this gap by investigating the genetic structure, diversity and high-altitude adaptation of 6 Tibetan sheep breeds using whole-genome resequencing data.

RESULTS

Six Tibetan sheep breeds were investigated in this study, and whole-genome resequencing data were used to investigate their genetic structure and population diversity. The results showed that the 6 Tibetan sheep breeds exhibited distinct separation in the phylogenetic tree; however, the levels of differentiation among the breeds were minimal, with extensive gene flow observed. Population structure analysis broadly categorized the 6 breeds into 3 distinct ecological types: plateau-type, valley-type and Euler-type. Analysis of unique single-nucleotide polymorphisms (SNPs) and selective sweeps between Argali and Tibetan sheep revealed that Tibetan sheep domestication was associated primarily with sensory and signal transduction, nutrient absorption and metabolism, and growth and reproductive characteristics. Finally, comprehensive analysis of selective sweep and transcriptome data suggested that Tibetan sheep breeds inhabiting different altitudes on the Qinghai-Tibet Plateau adapt by enhancing cardiopulmonary function, regulating body fluid balance through renal reabsorption, and modifying nutrient digestion and absorption pathways.

CONCLUSION

In this study, we investigated the genetic diversity and population structure of 6 Tibetan sheep breeds in Qinghai Province, China. Additionally, we analyzed the domestication traits and investigated the unique adaptation mechanisms residing varying altitudes in the plateau region of Tibetan sheep. This study provides valuable insights into the evolutionary processes of Tibetan sheep in extreme environments. These findings will also contribute to the preservation of genetic diversity and offer a foundation for Tibetan sheep diversity preservation and plateau animal environmental adaptation mechanisms.

A novel codominant plumage color pattern of white breast patches in WugangTong geese was controlled by EDNRB2

Two basic plumage color patterns are observed in adult geese: solid grey (G) or colorless white (W). However, a Chinese indigenous breed, the Wugangtong goose (WGT), continues to be subject to selective breeding efforts as it displays segregation of plumage colors, including G, W, and a novel color pattern designated Wb (G with white breast circles). The underlying genetic mechanisms responsible for the Wb phenotype are yet to be determined. The current study employed the population differentiation index (FST) to analyze 90 geese exhibiting diverse plumage colors, identifying the fifth intron of EDNRB2 as a particularly noteworthy region with the highest FST values. Sanger sequencing of the region surrounding the EDNRB2 gene identified a 14-bp insertion within exon 3 as the causal mutation. The heterozygosity of this 14-bp insertion and wild-type alleles was completely associated with the Wb phenotype, thereby substantiating the codominant nature of the G and W phenotypes. An inter-species corroborated this finding cross between the graylag (no 14-bp insertion) and the swan goose (homozygous for the 14-bp insertion) breeds, as hybrids from this cross exhibited the Wb phenotype. Transcriptomes from white breast patches and gray dorsal skins of 4 Wb geese were compared. A significant downregulation of genes involved in melanin synthesis and melanocyte development was observed, including EDRNB2 and MLANA. The downregulation of MLANA indicated that the mutated EDNRB2 resulted in melanocyte loss in specific body regions, as MLANA is a marker gene for melanocytes. The findings were corroborated by melanin staining using the Mansson-Fontana method, which revealed no melanin particles deposited in the white breast patches. In summary, the gray plumage color was codominant to the white color in WGT geese, and plumage color variations were controlled by EDNRB2. The findings of our study offer valuable and practical guidance for the purification of plumage colors among WGT, whether through traditional phenotype selection or molecular breeding methods.

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.

Integrating dynamic high-throughput phenotyping and genetic analysis to monitor growth variation in foxtail millet

BACKGROUND

Foxtail millet [Setaria italica (L.) Beauv] is a C4 graminoid crop cultivated mainly in the arid and semiarid regions of China for more than 7000 years. Its grain highly nutritious and is rich in starch, protein, essential vitamins such as carotenoids, folate, and minerals. To expand the utilisation of foxtail millet, efficient and precise methods for dynamic phenotyping of its growth stages are needed. Traditional foxtail millet monitoring methods have high labour costs and are inefficient and inaccurate, impeding the precise evaluation of foxtail millet genotypic variation.

RESULTS

This study introduces a high-throughput imaging system (HIS) with advanced image processing techniques to enhance monitoring efficiency and data quality. The HIS can accurately extract a range of key growth feature parameters, such as plant height (PH), convex hull area (CHA), side projected area (SPA) and colour distribution, from foxtail millet images. Compared with traditional manual measurements, this HIS improved data quality and phenotyping of the key foxtail millet growth traits. High-throughput phenotyping combined with a genome-wide association study (GWAS) revealed genetic loci associated with dynamic growth traits, particularly plant height (PH), in foxtail millet. The loci were linked to genes involved in the gibberellic acid (GA) synthesis pathway related to PH.

CONCLUSION

The HIS developed in this study enables the efficient and dynamic monitoring of foxtail millet phenotypic traits. It significantly improves the quality of data obtained for phenotyping key growth traits. The integration of high-throughput phenotyping with GWAS provides new insights into the genetic underpinnings of dynamic growth traits, particularly plant height, by identifying associated genetic loci in the GA synthesis pathway. This methodological advancement opens new avenues for the precise phenotyping and exploration of genetic resources in foxtail millet, potentially enhancing its utilisation.

A novel 22-bp InDel within FGF7 gene is significantly associated with growth traits in goat

Fibroblast growth factor 7 (FGF7) is involved in lipid metabolism, which is considered as a candidate gene with close relation with muscle development by eGWAs and RNA-Seq analyses. To date, limited research has been conducted on the relationship between FGF7 gene and growth traits. The main objective of this work was to further investigate the association between novel InDel within FGF7 gene and growth traits in goat. Herein, FGF7 mRNA expression levels were investigated in various Fuqing goat tissues. We found that FGF7 gene was expressed in six adult goat tissues with the highest mRNA levels in adipose tissue. This result suggested that FGF7 gene might play a critical role in fat deposition. We also detected potential polymorphisms in Fuqing, Nubian and Jianyang Daer breeds. A 22-bp InDel polymorphism in FGF7 gene was detected in 396 goats and the three genotypes were designated as II, ID, and DD. Correlation analysis revealed that InDel polymorphism was significantly associated with growth traits (P < 0.05). Goats with genotypes ID and/or II had superior growth traits compared to those with genotype DD. In summary, our findings suggested that the 22-bp InDel within FGF7 gene could act as a molecular marker to improve the growth traits of goats in breeding programs.

Whole-genome sequencing identifies functional genes for environmental adaptation in Chinese sheep

Sheep (Ovis aries), among the first domesticated species, are now globally widespread and exhibit remarkable adaptability to diverse environments. In this study, we perform whole-genome sequencing of 266 animals from 18 distinct Chinese sheep populations, each displaying unique phenotypes indicative of adaptation to varying environmental conditions. Integrating 131 environmental factors with single nucleotide polymorphism variations, we conduct a comprehensive genetic-environmental association analysis. This analysis identifies 35 key genes likely integral to the environmental adaptation of sheep. The functions of these genes include fat tail formation (HOXA10, HOXA11, JAZF1), wool characteristics (FER, FGF5, MITF, PDE4B), horn phenotypes (RXFP2), reproduction (HIBADH, TRIM71, C6H4orf22), and growth traits (ADGRL3, TRHDE). Notably, we observe a significant correlation between the frequency of missense mutations in the PAPSS2 and RXFP2 genes and variations in altitude. Our study reveals candidate genes for adaptive variation in sheep and demonstrates the diversity in how sheep adapt to their environment.

Genome-Wide Association Study Reveals the Genetic Architecture of Growth and Meat Production Traits in a Chicken F2 Resource Population

BACKGROUND/OBJECTIVES

For genomic selection to enhance the efficiency of broiler production, finding SNPs and candidate genes that define the manifestation of main selected traits is essential. We conducted a genome-wide association study (GWAS) for growth and meat productivity traits of roosters from a chicken F2 resource population (n = 152).

METHODS

The population was obtained by crossing two breeds with contrasting phenotypes for performance indicators, i.e., Russian White (slow-growing) and Cornish White (fast-growing). The birds were genotyped using the Illumina Chicken 60K SNP iSelect BeadChip. After LD filtering of the data, 54,188 SNPs were employed for the GWAS analysis that allowed us to reveal significant specific associations for phenotypic traits of interest and economic importance.

RESULTS

At the threshold value of p < 9.2 × 10-7, 83 SNPs associated with body weight at the age of 28, 42, and 63 days were identified, as well as 171 SNPs associated with meat qualities (average daily gain, slaughter yield, and dressed carcass weight and its components). Moreover, 34 SNPs were associated with a group of three or more traits, including 15 SNPs significant for a group of growth traits and 5 SNPs for a group of meat productivity indicators. Relevant to these detected SNPs, nine prioritized candidate genes associated with the studied traits were revealed, including WNT2, DEPTOR, PPA2, UNC80, DDX51, PAPPA, SSC4D, PTPRU, and TLK2.

CONCLUSIONS

The found SNPs and candidate genes can serve as genetic markers for growth and meat performance characteristics in chicken breeding in order to achieve genetic improvement in broiler production.

Genome-wide discovery of selection signatures in four Anatolian sheep breeds revealed by ddRADseq

High-density genomic data analyzed by accurate statistical methods are of potential to enlighten past breeding practices such as selection by unraveling fixed regions. In this study, four native Turkish sheep breeds (80 samples) were genotyped via 296.097 single nucleotide polymorphisms (SNPs) detected by double-digest restriction site-associated DNA (ddRADseq) library preparation combined with the Illumina HiSeq X Ten instrument in order to identify genes under selection pressure. A total of 32, 136, 133, and 119 protein-coding genes were detected under selection pressure by runs of homozygosity (ROH), integrated haplotype score (iHS), the ratio of extended haplotype homozygosity (Rsb), and fixation index (FST) approaches, respectively. Of these, a total of 129 genes were identified by at least two statistical models which overlapped with a total of 52 quantitative trait loci (QTL)-associated SNPs, known to be related to fiber diameter, milk content, body weight, carcass traits, some blood parameters, and entropion. A total of six genes under selection pressure were validated by three statistical approaches five of which are of potential to be integrated into animal breeding since they were associated with wool fiber diameter (ZNF208B), behaviors related to neurocognitive development (CBX1 and NFE2L1), adaptation to high-altitude (SDK1), and anxiety causing internal stress (GSG1L). The sixth gene (COPZ1) turned out to play an important role in coping with different types of cancer in mammals. In particular, ROH analysis uncovered significant findings that the Güney Karaman (GKR) had experienced different selection practices than the Akkaraman (AKR) breed. Moreover, some genes specifically under selection in the GKR breed turned out to be associated with olfaction (OR6K6, OR6N1, OR6N2, and OR4C16), survival during the gestation period (PRR15L), and heat stress (CDK5RAP9). The results of this study imply that GKR may become genetically different from the AKR breed at the genome level due to most probably experiencing different adaptation processes occurring in raised climatic conditions. These differences should be conserved to face future challenges, while other native Turkish sheep breeds could be monitored via genome-wide high-density SNP data to obtain deeper knowledge about the effects of natural selection.

Genetic structure and selective sweeps in Kirghiz sheep using SNP50K bead chip

The objective of this study is to analyze environmental genetic selection signals in large-scale sheep populations with conflicting environmental adaptations, aiming to identify and isolate genes associated with environmental adaptations in sheep populations. Kirghiz sheep, which inhabit high-altitude environments year-round, demonstrate the ability to adapt to extreme conditions. In this study, 42 Kirghiz sheep, 24 Tien-Shan in Kyrgyzstan sheep, 189 Qira black sheep, and 160 Chinese Merino sheep were genotyped using Illumina Ovine SNP50K chip. Regions exhibiting a selection signal threshold of 5%, as well as PI analysis and haplotype statistical scanning gene data were annotated, and intersecting genes were identified as candidate genes. Through Fst and haplotype statistical analysis revealed the key gene PDGFD and its vicinity's impact on fat deposition in sheep tails. Additionally, Fst and PI analysis uncovered genes related to high-altitude adaptation as well as those linked to animal growth and reproduction.Further GO and KEGG enrichment pathway analyses unveiled pathways associated with high-altitude adaptation such as negative regulation of peptidyl-tyrosine phosphorylation and xenobiotic metabolism processes.This investigation into the adaptability of Kirghiz sheep provides theoretical support and practical guidance for the conservation and genetic enhancement of Kirghiz sheep germplasm resources.

The genetic characterization of germplasm and identification of the litter size trait associated candidate genes in Dexin mutton and fine-wool sheep

Xinjiang is a major province of sheep breeding in China, which plays an important role in meeting people's needs for meat products, increasing farmers' income and sustainable development of animal husbandry. However, the genetic differentiation relationship between breeds was not clear, and most sheep had low fecundity, which seriously restricted the efficient development of sheep industry. Therefore, this study used the whole genome resequencing to detect the genetic variation of Dexin mutton and fine-wool sheep, explored the selected regions and important genes of the litter size traits, analyzed the genetic mechanism of reproductive traits, and provided new insights for the high fecundity breeding of sheep. A total of 5,236.338 G genome data and 35,884,037 SNPs were obtained. Furthermore, we identified 39 selection signals spanning candidate genes, 99 genes were significantly associated related to growth, reproduction and immunity, among which, BRIP1, BMPR1B, BMP4, NGF, etc. genes, and MAKP signaling pathway, Fanconi anemia pathway and Thyroid hormone signaling pathway and other signaling pathways were significantly correlated with litter size trait. Among them, we identified NGF, TrKA and BRIP1 genes was the important genes for sheep litter size traits and the mutation frequencies of 9 SNPs in BRIP1 gene were significantly different in domestic sheep in the world. The research provided new insights for the breeding of self-cultivated meat fine-wool sheep.

Genome-wide comparative analysis reveals selection signatures for reproduction traits in prolific Suffolk sheep

The identification of genome-wide selection signatures can reveal the potential genetic mechanisms involved in the generation of new breeds through natural or artificial selection. In this study, we screened the genome-wide selection signatures of prolific Suffolk sheep, a new strain of multiparous mutton sheep, to identify candidate genes for reproduction traits and unravel the germplasm characteristics and population genetic evolution of this new strain of Suffolk sheep. Whole-genome resequencing was performed at an effective sequencing depth of 20× for genomic diversity and population structure analysis. Additionally, selection signatures were investigated in prolific Suffolk sheep, Suffolk sheep, and Hu sheep using fixation index (F ST) and heterozygosity H) analysis. A total of 5,236.338 Gb of high-quality genomic data and 28,767,952 SNPs were obtained for prolific Suffolk sheep. Moreover, 99 selection signals spanning candidate genes were identified. Twenty-three genes were significantly associated with KEGG pathway and Gene Ontology terms related to reproduction, growth, immunity, and metabolism. Through selective signal analysis, genes such as ARHGEF4, CATIP, and CCDC115 were found to be significantly correlated with reproductive traits in prolific Suffolk sheep and were highly associated with the mTOR signaling pathway, the melanogenic pathway, and the Hippo signaling pathways, among others. These results contribute to the understanding of the evolution of artificial selection in prolific Suffolk sheep and provide candidate reproduction-related genes that may be beneficial for the establishment of new sheep breeds.

Comparative whole-genome resequencing to uncover selection signatures linked to litter size in Hu Sheep and five other breeds

Hu sheep (HS), a breed of sheep carrying the FecB mutation gene, is known for its "year-round estrus and multiple births" and is an ideal model for studying the high fecundity mechanisms of livestock. Through analyzing and comparing the genomic selection features of Hu sheep and other sheep breeds, we identified a series of candidate genes that may play a role in Hu sheep's high fecundity mechanisms. In this study, we conducted whole-genome resequencing on six breeds and screened key mutations significantly correlated with high reproductive traits in sheep. Notably, the CC2D1B gene was selected by the fixation index (FST) and the cross-population composite likelihood ratio (XP-CLR) methods in HS and other five breeds. It was worth noting that the CC2D1B gene in HS was different from that in other sheep breeds, and seven missense mutations have been identified. Furthermore, the linkage disequilibrium (LD) analysis revealed a strong linkage disequilibrium in this specific gene region. Subsequently, by performing different grouping based on FecB genotypes in Hu sheep, genome-wide selective signal analysis screened several genes related to reproduction, such as BMPR1B and PPM1K. Besides, FST analysis identified functional genes related to reproductive traits, including RHEB, HSPA2, PPP1CC, HVCN1, and CCDC63. Additionally, a missense mutation was found in the CCDC63 gene and the haplotype was different between the high reproduction (HR) group and low reproduction (LR) group in HS. In summary, we discovered genetic differentiation among six distinct breeding sheep breeds at the whole genome level. Additionally, we identified a set of genes which were associated with reproductive performance in Hu sheep and visualized how these genes differed in different breeds. These findings laid a theoretical foundation for understanding genetic mechanisms behind high prolific traits in sheep.

FecB mutation and litter size are associated with a 90-base pair deletion in BMPR1B in East Friesian and Hu crossbred sheep

Litter size is a critical economic trait in livestock, but only a few studies have focused on associated indel mutations in BMPR1B, a key regulator of ovulation and litter size in sheep. We evaluated the effects of BMPR1B mutations on the reproductive performance of sheep. We used Hu, East Friesian, and East Friesian/Hu crossbred sheep as experimental subjects and identified a novel 90 bp deletion in BMPR1B, which coincides with the c.746A > G (FecB mutation) genotype. The correlation between the two loci and litter size was then evaluated. We identified three genotypes for the Del-90bp locus, namely, II, ID, and DD, and three genotypes for the c.746A > G locus, namely ++, B+, and BB. Both Del-90bp and c.746A > G significantly affected the litter size of Hu and East Friesian/Hu crossbred sheep. Linkage disequilibrium analysis revealed a strong linkage disequilibrium between these loci in Hu sheep and the F1 population (r2 > 0.33), which suggests that detecting this 90 bp deletion might be a simple method to identify the likely carriers of c.746A > G. However, the function of this 90-bp deletion still needs further exploration. We provide genetic data that can be used as a reference for the breeding of improved prolific traits in sheep.

Genome-Wide Association Study of Fiber Diameter in Alpacas

The aim of this study was the identification of candidate genomic regions associated with fiber diameter in alpacas. DNA samples were collected from 1011 female Huacaya alpacas from two geographical Andean regions in Peru (Pasco and Puno), and three alpaca farms within each region. The samples were genotyped using an Affymetrix Custom Alpaca genotyping array containing 76,508 SNPs. After the quality controls, 960 samples and 51,742 SNPs were retained. Three association study methodologies were performed. The GWAS based on a linear model allowed us to identify 11 and 35 SNPs (-log10(p-values) > 4) using information on all alpacas and alpacas with extreme values of fiber diameter, respectively. The haplotype and marker analysis method allowed us to identify nine haplotypes with standardized haplotype heritability higher than six standard deviations. The selection signatures based on cross-population extended haplotype homozygosity (XP-EHH) allowed us to identify 180 SNPs with XP-EHH values greater than |3|. Four candidate regions with adjacent SNPs identified via two association methods of analysis are located on VPA6, VPA9, VPA29 and one chromosomally unassigned scaffold. This study represents the first analysis of alpaca whole genome association with fiber diameter, using a recently assembled alpaca SNP microarray.

Identification of copy number variation in Tibetan sheep using whole genome resequencing reveals evidence of genomic selection

BACKGROUND

Copy number variation (CNV) is an important source of structural variation in the mammalian genome. CNV assays present a new method to explore the genomic diversity of environmental adaptations in animals and plants and genes associated with complex traits. In this study, the genome-wide CNV distribution characteristics of 20 Tibetan sheep from two breeds (10 Oula sheep and 10 Panou sheep) were analysed using whole-genome resequencing to investigate the variation in the genomic structure of Tibetan sheep during breeding.

RESULTS

CNVs were detected using CNVnator, and the overlapping regions of CNVs between individual sheep were combined. Among them, a total of 60,429 CNV events were detected between the indigenous sheep breed (Oula) and the synthetic sheep breed (Panou). After merging the overlapping CNVs, 4927 CNV regions (CNVRs) were finally obtained. Of these, 4559 CNVRs were shared by two breeds, and there were 368 differential CNVRs. Deletion events have a higher percentage of occurrences than duplication events. Functional enrichment analysis showed that the shared CNVRs were significantly enriched in 163 GO terms and 62 KEGG pathways, which were mainly associated with organ development, neural regulation, immune regulation, digestion and metabolism. In addition, 140 QTLs overlapped with some of the CNVRs at more than 1 kb, such as average daily gain QTL, body weight QTL, and total lambs born QTL. Many of the CNV-overlapping genes such as PPP3CA, SSTR1 and FASN, overlap with the average daily weight gain and carcass weight QTL regions. Moreover, VST analysis showed that XIRP2, ABCB1, CA1, ASPA and EEF2 differed significantly between the synthetic breed and local sheep breed. The duplication of the ABCB1 gene may be closely related to adaptation to the plateau environment in Panou sheep, which deserves further study. Additionally, cluster analysis, based on all individuals, showed that the CNV clustering could be divided into two origins, indicating that some Tibetan sheep CNVs are likely to arise independently in different populations and contribute to population differences.

CONCLUSIONS

Collectively, we demonstrated the genome-wide distribution characteristics of CNVs in Panou sheep by whole genome resequencing. The results provides a valuable genetic variation resource and help to understand the genetic characteristics of Tibetan sheep. This study also provides useful information for the improvement and breeding of Tibetan sheep in the future.

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.

Position of Hungarian Merino among other Merinos, within-breed genetic similarity network and markers associated with daily weight gain

OBJECTIVE

In this study, we aimed to position the Hungarian Merino among other Merino-derived sheep breeds, explore the characteristics of our sampled animals' genetic similarity network within the breed, and highlight single nucleotide polymorphisms (SNPs) associated with daily weight-gain.

METHODS

Hungarian Merino (n = 138) was genotyped on Ovine SNP50 Bead Chip (Illumina, San Diego, CA, USA) and positioned among 30 Merino and Merino-derived breeds (n = 555). Population characteristics were obtained via PLINK, SVS, Admixture, and Treemix software, within-breed network was analysed with python networkx 2.3 library. Daily weight gain of Hungarian Merino was standardised to 60 days and was collected from the database of the Association of Hungarian Sheep and Goat Breeders. For the identification of loci associated with daily weight gain, a multi-locus mixed-model was used.

RESULTS

Supporting the breed's written history, the closest breeds to Hungarian Merino were Estremadura and Rambouillet (pairwise FST values are 0.035 and 0.036, respectively). Among Hungarian Merino, a highly centralised connectedness has been revealed by network analysis of pairwise values of identity-by-state, where the animal in the central node had a betweenness centrality value equal to 0.936. Probing of daily weight gain against the SNP data of Hungarian Merinos revealed five associated loci. Two of them, OAR8_17854216.1 and s42441.1 on chromosome 8 and 9 (-log10P>22, false discovery rate<5.5e-20) and one locus on chromosome 20, s28948.1 (-log10P = 13.46, false discovery rate = 4.1e-11), were close to the markers reported in other breeds concerning daily weight gain, six-month weight, and post-weaning gain.

CONCLUSION

The position of Hungarian Merino among other Merino breeds has been determined. We have described the similarity network of the individuals to be applied in breeding practices and highlighted several markers useful for elevating the daily weight gain of Hungarian Merino.

High-Altitude Stress Orchestrates mRNA Expression and Alternative Splicing of Ovarian Follicle Development Genes in Tibetan Sheep

High-altitude stress threatens the survival rate of Tibetan sheep and reduces their fertility. However, the molecular basis of this phenomenon remains elusive. Here, we used RNA-seq to elucidate the transcriptome dynamics of high-altitude stress in Tibetan sheep ovaries. In total, 104 genes were characterized as high-altitude stress-related differentially expressed genes (DEGs). In addition, 36 DEGs contributed to ovarian follicle development, and 28 of them were downregulated under high-altitude stress. In particular, high-altitude stress significantly suppressed the expression of two ovarian lymphatic system marker genes: LYVE1 and ADAMTS-1. Network analysis revealed that luteinizing hormone (LH)/follicle-stimulating hormone (FSH) signaling-related genes, such as EGR1, FKBP5, DUSP1, and FOS, were central regulators in the DEG network, and these genes were also suppressed under high-altitude stress. As a post-transcriptional regulation mechanism, alternative splicing (AS) is ubiquitous in Tibetan sheep. High-altitude stress induced 917 differentially alternative splicing (DAS) events. High-altitude stress modulated DAS in an AS-type-specific manner: suppressing skipped exon events but increasing retained intron events. C2H2-type zinc finger transcription factors and RNA processing factors were mainly enriched in DAS. These findings revealed high-altitude stress repressed ovarian development by suppressing the gene expression of LH/FSH hormone signaling genes and inducing intron retention of C2H2-type zinc finger transcription factors.

Whole-genome resequencing reveals domestication and signatures of selection in Ujimqin, Sunit, and Wu Ranke Mongolian sheep breeds

Objective: The current study aimed to perform whole-genome resequencing of Chinese indigenous Mongolian sheep breeds including Ujimqin, Sunit, and Wu Ranke sheep breeds (UJMQ, SNT, WRK) and deeply analyze genetic variation, population structure, domestication, and selection for domestication traits among these Mongolian sheep breeds.Methods: Blood samples were collected from a total of 60 individuals comprising 20 WRK, 20 UJMQ, and 20 SNT. For genome sequencing, about 1.5 μg of genomic DNA was used for library construction with an insert size of about 350 bp. Pair-end sequencing were performed on Illumina NovaSeq platform, with the read length of 150 bp at each end. We then investigated the domestication and signatures of selection in these sheep breeds.Results: According to the population and demographic analyses, WRK and SNT populations were very similar, which were different from UJMQ populations. Genome wide association study identified 468 and 779 significant loci from SNT vs UJMQ, and UJMQ vs WRK, respectively. However, only 3 loci were identified from SNT vs WRK. Genomic comparison and selective sweep analysis among these sheep breeds suggested that genes associated with regulation of secretion, metabolic pathways including estrogen metabolism and amino acid metabolism, and neuron development have undergone strong selection during domestication.Conclusion: Our findings will facilitate the understanding of Chinese indigenous Mongolian sheep breeds domestication and selection for complex traits and provide a valuable genomic resource for future studies of sheep and other domestic animal breeding.

lncRNA-miRNA-mRNA ceRNA Network Involved in Sheep Prolificacy: An Integrated Approach

Understanding the molecular pattern of fertility is considered as an important step in breeding of different species, and despite the high importance of the fertility, little success has been achieved in dissecting the interactome basis of sheep fertility. However, the complex mechanisms associated with prolificacy in sheep have not been fully understood. Therefore, this study aimed to use competitive endogenous RNA (ceRNA) networks to evaluate this trait to better understand the molecular mechanisms responsible for fertility. A competitive endogenous RNA (ceRNA) network of the corpus luteum was constructed between Romanov and Baluchi sheep breeds with either good or poor genetic merit for prolificacy using whole-transcriptome analysis. First, the main list of lncRNAs, miRNAs, and mRNA related to the corpus luteum that alter with the breed were extracted, then miRNA−mRNA and lncRNA−mRNA interactions were predicted, and the ceRNA network was constructed by integrating these interactions with the other gene regulatory networks and the protein−protein interaction (PPI). A total of 264 mRNAs, 14 lncRNAs, and 34 miRNAs were identified by combining the GO and KEGG enrichment analyses. In total, 44, 7, 7, and 6 mRNAs, lncRNAs, miRNAs, and crucial modules, respectively, were disclosed through clustering for the corpus luteum ceRNA network. All these RNAs involved in biological processes, namely proteolysis, actin cytoskeleton organization, immune system process, cell adhesion, cell differentiation, and lipid metabolic process, have an overexpression pattern (Padj < 0.01). This study increases our understanding of the contribution of different breed transcriptomes to phenotypic fertility differences and constructed a ceRNA network in sheep (Ovis aries) to provide insights into further research on the molecular mechanism and identify new biomarkers for genetic improvement.

Selection Signature Analyses Revealed Genes Associated With Adaptation, Production, and Reproduction in Selected Goat Breeds in Kenya

Artificial and natural selection in livestock is expected to leave unique footprints on their genomes. Goat breeds in Kenya have evolved for survival, breeding, and production in various harsh ecological areas, and their genomes are likely to have acquired unique alleles for adaptation to such diverse production environments and other traits of economic importance. To investigate signals of selection for some selected goat breeds in Kenya, Alpine (n = 29), Galla (n = 12), Saanen (n = 24), and Toggenburg (n = 31) were considered. A total of 53,347 single-nucleotide polymorphisms (SNPs) generated using the Illumina GoatSNP50 BeadChip were analyzed. After quality control, 47,663 autosomal single-nucleotide polymorphisms remained for downstream analyses. Several complementary approaches were applied for the following analyses: integrated Haplotype Score (iHS), cross-population-extended haplotype homozygosity (XP-EHH), hapFLK, and FLK. A total of 404 top genomic regions were identified across all the four breeds, based on the four complementary analyses. Out of the 16 identified putative selection signature regions by the intersection of multiple-selective signal analyses, most of the putative regions were found to overlap significantly with the iHS and XP-EHH analyses on chromosomes 3, 4, 10, 15, 22, and 26. These regions were enriched with some genes involved in pathways associated directly or indirectly with environmental adaptation regulating immune responses (e.g., HYAL1 and HYAL3), milk production (e.g., LEPR and PDE4B), and adaptability (e.g., MST1 and PCK). The results revealed few intersect between breeds in genomic selection signature regions. In general, this did not present the typical classic selection signatures as predicted due to the complex nature of the traits. The results support that some various selection pressures (e.g., environmental challenges, artificial selection, and genome admixture challenges) have molded the genome of goat breeds in Kenya. Therefore, the research provides new knowledge on the conservation and utilization of these goat genetic resources in Kenya. In-depth research is needed to detect precise genes connected with adaptation and production in goat breeds in Kenya.

Evaluation of Genetic Parameters of Growth of Pelibuey and Blackbelly Sheep through Pedigree in Mexico

Simple Summary

The importance of sheep production in Mexico has been increasing in recent years. Animal performance can be improved through continuous selection and cumulative genetic management. Studies on non-genetic factors that influence the growth characteristics of sheep, due to their positive genetic correlation with other live weights, can present data of great interest to small producers with limited resources. In the present study, genealogical and functional information from the historical archive of the Mexican Association of Sheep Breeders (AMCO) was used. The objective was to estimate the heritability and genetic correlations of the growth of ewes born and weaned at different times of the same year from different herds with pedigree registration. The expected results could be incorporated into genetic selection programs with a favorable impact on the local economy of small sheep producers in Mexico.

Abstract

Birth weight (BW) and weaning weight (WW) data from Pelibuey and Blackbelly lambs belonging to the Asociación Mexicana de Criadores de Ovinos (AMCO) were used with the objective of estimating genetic parameters (heritability and genetic correlations) and analyzing the growth characteristics of ewes born and weaned at different times of the same year from different herds with pedigree registration. In the case of Pelibuey lambs, the animal model included the weaning weight at 75 days of age, considering the direct additive genetic effect, maternal additive genetic effect, covariance between direct and maternal effects, as well as the permanent environmental effect of the mother. The direct estimators of heritability for Pelibuey were BW = 0.01 ± 0.021 and WW = 0.31 ± 0.074 and for Blackbelly they were BW = 0.05 ± 0.042 and WW = 0.41 ± 0.146. In the case of the maternal heritability estimators in Pelibuey they were BW = 0.02 ± 0.040 and WW = 0.21 ± 0.121 and for Blackbelly they were BW = 0.12 ± 0.054 and WW = 0.28 ± 0.121. The magnitude of the estimates of genetic correlations between direct and maternal effects for adjusted weaning weight at 75 days of age indicate that genetic progress may be slow in a breeding program. However, these selection results could be included in the short term in the breeding programs for the Pelibuey and Blackbelly breeds in Mexico, for livestock development in low-income rural areas.

Variations in HIF-1α Contributed to High Altitude Hypoxia Adaptation via Affected Oxygen Metabolism in Tibetan Sheep

Simple Summary

Hypoxia-inducible factors (HIFs) play an important role in the adaptation of animals to high-altitude hypoxia. In high-altitude indigenous species, variation in the hypoxia-inducible factor-1 alpha (HIF-1α) gene has been reported in Tibetans, yaks and Tibetan horses, but has not been investigated in Tibetan sheep, and is not known if it might affect high-altitude hypoxia adaptation in these sheep. In this study, Kompetitive Allele-Specific PCR (KASP) was used for genotyping of ovine HIF-1α and investigated the effect of variation in HIF-1α on the high-altitude hypoxia adaptation of Tibetan sheep. The results suggest that ovine HIF-1α variants may promote the ability of oxygen utilization in Tibetan sheep and it may serve as a genetic marker for improving high-altitude hypoxia adaptability.

Abstract

The Tibetan sheep is an indigenous species of the Tibetan plateau and has been well adapted to high-altitude hypoxia. In comparison to lowland sheep breeds, the blood gas indicators have changed and the HIFs signaling pathway is activated in Tibetan sheep. These phenotypic and genetic alterations in Tibetan sheep are thought to be an important basis for adaptation to high-altitude hypoxia and variation in genes encoding the subunits that make up HIFs, such as HIF-1α can affect blood gas indicators. In this study, exons 9, 10, 12 of the HIF-1α gene were sequenced to find variations and 3 SNPs were detected, and these 3 SNPs were genotyped by KASP in 341 Hu sheep and 391 Tibetan sheep. In addition, 197 Hu sheep, 160 Tibetan sheep and 12 Gansu alpine merino sheep were used for blood gas indicators analysis. The results showed significant differences between the blood gas indicators of high-altitude breeds (Tibetan sheep and Gansu alpine merino sheep) and low-altitude breeds (Hu sheep), implying that the differences in blood gas indicators are mainly caused by differences in altitude. The haplotype combinations H2H3 and H1H3 were more frequent in the Tibetan sheep population, H2H3 increases O₂ carrying capacity by increasing hematocrit and hemoglobin concentrations; H1H3 makes O₂ dissociate more readily from oxyhemoglobin by decreasing partial pressure of oxygen and oxygen saturation. These results suggest that variants at the HIF-1α promote the ability of oxygen utilization in Tibetan sheep, which may underpin the survival and reproduction of Tibetan sheep on the Tibetan plateau.

The Quest for Genes Involved in Adaptation to Climate Change in Ruminant Livestock

Livestock radiated out from domestication centres to most regions of the world, gradually adapting to diverse environments, from very hot to sub-zero temperatures and from wet and humid conditions to deserts. The climate is changing; generally global temperature is increasing, although there are also more extreme cold periods, storms, and higher solar radiation. These changes impact livestock welfare and productivity. This review describes advances in the methodology for studying livestock genomes and the impact of the environment on animal production, giving examples of discoveries made. Sequencing livestock genomes has facilitated genome-wide association studies to localize genes controlling many traits, and population genetics has identified genomic regions under selection or introgressed from one breed into another to improve production or facilitate adaptation. Landscape genomics, which combines global positioning and genomics, has identified genomic features that enable animals to adapt to local environments. Combining the advances in genomics and methods for predicting changes in climate is generating an explosion of data which calls for innovations in the way big data sets are treated. Artificial intelligence and machine learning are now being used to study the interactions between the genome and the environment to identify historic effects on the genome and to model future scenarios.

British Sheep Breeds as a Part of World Sheep Gene Pool Landscape: Looking into Genomic Applications

Sheep farming has been an important sector of the UK's economy and rural life for many centuries. It is the favored source of wool, meat and milk products. In the era of exponential progress in genomic technologies, we can now address the questions of what is special about UK sheep breed genotypes and how they differ genetically form one another and from other countries. We can reflect how their natural history has been determined at the level of their genetic code and what traces have been left in their genomes because of selection for phenotypic traits. These include adaptability to certain environmental conditions and management, as well as resistance to disease. Application of these advancements in genetics and genomics to study sheep breeds of British domestic selection has begun and will continue in order to facilitate conservation solutions and production improvement.