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.

Population structure and selection signal analysis of indigenous sheep from the southern edge of the Taklamakan Desert.

Analyzing the genetic diversity and selection characteristics of sheep (Ovis aries) holds significant value in understanding their environmental adaptability, enhancing breeding efficiency, and achieving effective conservation and rational utilization of genetic resources. In this study, we utilized Illumina Ovine SNP 50 K BeadChip data from four indigenous sheep breeds from the southern margin of the Taklamakan Desert (Duolang sheep: n = 36, Hetian sheep: n = 74, Kunlun sheep: n = 27, Qira black sheep: n = 178) and three foreign meat sheep breeds (Poll Dorset sheep: n = 105, Suffolk sheep: n = 153, Texel sheep: n = 150) to investigate the population structure, genetic diversity, and genomic signals of positive selection within the indigenous sheep. According to the Principal component analysis (PCA), the Neighbor-Joining tree (NJ tree), and Admixture, we revealed distinct clustering patterns of these seven sheep breeds based on their geographical distribution. Then used Cross Population Extended Haplotype Homozygosity (XP-EHH), Fixation Index (FST), and Integrated Haplotype Score (iHS), we identified a collective set of 32 overlapping genes under positive selection across four indigenous sheep breeds. These genes are associated with wool follicle development and wool traits, desert environmental adaptability, disease resistance, reproduction, and high-altitude adaptability. This study reveals the population structure and genomic selection characteristics in the extreme desert environments of native sheep breeds from the southern edge of the Taklimakan Desert, providing new insights into the conservation and sustainable use of indigenous sheep genetic resources in extreme environments. Additionally, these findings offer valuable genetic resources for sheep and other mammals to adapt to global climate change.

Whole-genome analysis reveals distinct adaptation signatures to diverse environments in Chinese domestic pigs.

BACKGROUND: Long-term natural and artificial selection has resulted in many genetic footprints within the genomes of pig breeds across distinct agroecological zones. Nevertheless, the mechanisms by which these signatures contribute to phenotypic diversity and facilitate environmental adaptation remain unclear. RESULTS: Here, we leveraged whole-genome sequencing data from 82 individuals from 6 domestic pig breeds originating in tropical, high-altitude, and frigid regions. Population genetic analysis suggested that habitat isolation significantly shaped the genetic diversity and contributed to population stratification in local Chinese pig breeds. Analysis of selection signals revealed regions under selection for adaptation in tropical (55.5 Mb), high-altitude (43.6 Mb), and frigid (17.72 Mb) regions. The potential functions of the selective sweep regions were linked to certain complex traits that might play critical roles in different geographic environments, including fat coverage in frigid environments and blood indicators in tropical and high-altitude environments. Candidate genes under selection were significantly enriched in biological pathways involved in environmental adaptation. These pathways included blood circulation, protein degradation, and inflammation for adaptation to tropical environments; heart and lung development, hypoxia response, and DNA damage repair for high-altitude adaptation; and thermogenesis, cold-induced vasodilation (CIVD), and the cell cycle for adaptation to frigid environments. By examining the chromatin state of the selection signatures, we identified the lung and ileum as two candidate functional tissues for environmental adaptation. Finally, we identified a mutation (chr1: G246,175,129A) in the cis-regulatory region of ABCA1 as a plausible promising variant for adaptation to tropical environments. CONCLUSIONS: In this study, we conducted a genome-wide exploration of the genetic mechanisms underlying the adaptability of local Chinese pig breeds to tropical, high-altitude, and frigid environments. Our findings shed light on the prominent role of cis-regulatory elements in environmental adaptation in pigs and may serve as a valuable biological model of human plateau-related disorders and cardiovascular diseases.

Recombination and selection trajectory of the monkeypox virus during its adaptation in the human population.

Monkeypox, caused by the monkeypox virus (MPXV), was historically confined to West and Central Africa but has now spread globally. Recombination and selection play crucial roles in the evolutionary adaptation of MPXV; however, the evolution of MPXV and its relationship with the recent, ground-breaking monkeypox epidemic remains poorly understood. To gain insights into the evolutionary dynamics of MPXV, comprehensive in silico recombination and selection analyses were conducted based on MPXV whole genome sequence data. Three types of recombination were identified: five ancestor-sharing interspecies recombination events, six specific interspecies recombination events and four intraspecies recombination events. The results highlight the prevalent occurrence of recombination in MPXV, with 73.3% occurring in variable regions of the genome. Selection analysis was performed from three dimensions: proteins around recombination regions, proteins from recombinant ancestors and MPXV branches, and whole-genome gene analysis. Results revealed 2 and 7 proteins under positive selection in the first two dimensions, respectively. These proteins are mainly involved in infection immunity, apoptosis regulation and viral virulence. Whole-genome analysis detected 25 genes under positive selection, mainly associated with immune response and viral regulation. Understanding their evolutionary patterns will help predict and prevent cross-species transmission, zoonotic outbreaks and potential human epidemics.

Genome-wide selection reveals candidate genes associated with multiple teats in Hu sheep.

Increasing the number of teats in sheep helps to improve the survival rate of sheep lambs after birth. In order to analyze the candidate genes related to the formation of multiple teats in Hu sheep, the present study was conducted to investigate the genetic pattern of multiple teats in Hu sheep. In this study, based on genome-wide data from 157 Hu sheep, Fst, xp-EHH, Pi and iHS signaling were performed, and the top 5% signal regions of each analyzed result were annotated based on the Oar_v4.0 for sheep. The results show that a total of 142 SNP loci were selected. We found that PTPRG, TMEM117 and LRP1B genes were closely associated with polypodium formation in Hu sheep, in addition, among the candidate genes related to polypodium we found genes such as TMEM117, SLC25A21 and NCKAP5 related to milk traits. The present study screened out candidate genes for the formation of multiple teats at the genomic level in Hu sheep.

Whole-Genome Resequencing Revealed Selective Signatures for Growth Traits in Hu and Gangba Sheep.

A genomic study was conducted to uncover the selection signatures in sheep that show extremely significant differences in growth traits under the same breed, age in months, nutrition level, and management practices. Hu sheep from Gansu Province and Gangba sheep from the Tibet Autonomous Region in China were selected. We collected whole-genome data from 40 sheep individuals (24 Hu sheep and 16 Gangba sheep), through whole-genome sequencing. Selection signals were analyzed using parameters such as FST, π ratio, and Tajima's D. We have identified several candidate genes that have undergone strong selection, particularly those associated with growth traits. Specifically, five growth-related genes were identified in both the Hu sheep group (HDAC1, MYH7B, LCK, ACVR1, GNAI2) and the Gangba sheep group (RBBP8, ACSL3, FBXW11, PLAT, CRB1). Additionally, in a genomic region strongly selected in both the Hu and Gangba sheep groups (Chr 22: 51,425,001-51,500,000), the growth-associated gene CYP2E1 was identified, further highlighting the genetic factors influencing growth characteristics in these breeds. This study analyzes the genetic basis for significant differences in sheep phenotypes, identifies candidate genes related to sheep growth traits, lays the foundation for molecular genetic breeding in sheep, and accelerates the genetic improvement in livestock.

Whole genome resequencing-based analysis of plateau adaptation in Meiren yak (Bos grunniens).

The Meiren yak is an important genetic resource in Gansu Province, China. In this study, we aimed to explore the evolutionary history and population structure of the genetic resource of Meiren yak and to mine the characteristic genes of Meiren yak. We analysed a total of 93 yaks of eight yak breeds based on whole genome resequencing combined with population genomics and used θπ ratio and Fst method to screen the selected sites in the genome region. The results proved that Meiren yak can be used as a potential genetic resource in Gansu Province. The genes in Meiren yak with positive selection in selection signal analysis were subjected to the Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses, which indicated that the genes were related to the adaptability to high altitude and hypoxic environment. By analysing the genetic variation of Meiren yak at the genome-wide level, this study provided a theoretical basis for genetic improvement of Meiren yak and for the development of high-quality yak resources.

Genome-wide screening for selection signatures in native and cosmopolitan cattle breeds reared in Türkiye.

Via long-term natural and artificial selection pressure, homozygosity may extend across the genome, leaving genomic patterns called selection signatures. This study is the first attempt to assess genome-wide selection signatures in six native Turkish and two cosmopolitan cattle breeds by 211.119 bi-allelic SNPs recovered using the double digest restriction associated DNA sequencing method. The integrated haplotype score (iHS) statistic was utilised to reveal selection signatures within populations, whereas the cross-population extended haplotype homozygosity (XP-EHH) and fixation index (FST ) approaches were preferred to reveal differently fixed genomic regions between native Turkish and cosmopolitan cattle breeds. Selection signatures in 142 genomic regions containing 305 genes were detected within eight cattle breeds by iHS statistics. The XP-EHH and FST approaches revealed that 197 and 114 SNPs were under selection pressure, respectively, which overlapped with 144 and 190 genes, respectively. A total of 18 genes were detected by at least two approaches. Six genes related to disease resistance (TTP2), meat yield (DIAPH3 and METTL21C), meat quality (ZNF24 and ZNF397) and first calving interval (ZSCAN30) turned out to be differently fixed between native Turkish and cosmopolitan cattle breeds, as they were identified by both XP-EHH and FST approaches. In addition, the iHS approach revealed that eight genes associated with visual modality (LSGN), olfaction (MOXD2, OR4C1F and OR4C1F), and immune response (TRBV3-1 and CLDN10) were under selection pressure in both native and cosmopolitan cattle breeds. Owing to their being significantly related to survival traits, these regions may have played a key role in cattle genome evolution. Future studies utilising denser genetic data are required to obtain deeper knowledge on effects of natural and artificial selection in Anatolian cattle breeds. © 2023 Stichting International Foundation for Animal Genetics.

Identification of Selection Signals on the X-Chromosome in East Adriatic Sheep: A New Complementary Approach.

Sheep are one of the most important livestock species in Croatia, found mainly in the Mediterranean coastal and mountainous regions along the East Adriatic coast, well adapted to the environment and mostly kept extensively. Our main objective was therefore to map the positive selection of the X-chromosome (18,983 SNPs that passed quality control), since nothing is known about the adaptation genes on this chromosome for any of the breeds from the Balkan cluster. Analyses were performed on a sample of eight native Croatian breeds (101 females and 100 males) representing the East Adriatic metapopulation and on 10 mouflons (five females and males), all sampled in Croatia. Three classical within-population approaches (extreme Runs of Homozygosity islands, integrated Haplotype Score, and number of Segregating Sites by Length) were applied along with our new approach called Haplotype Richness Drop (HRiD), which uses only the information contained in male haplotypes. We have also shown that phylogenetic analyses, such as the Median-joining network, can provide additional information when performed with the selection signals identified by HRiD. Our new approach identifies positive selection signals by searching for genomic regions that exhibit a sudden decline in haplotype richness. In total, we identified 14 positive selection signals, 11 using the classical approach and three using the HRiD approach, all together containing 34 annotated genes. The most reliable selection signal was mapped by all four approaches in the same region, overlapping between 13.17 and 13.60 Mb, and assigned to the CA5B, ZRSR2, AP1S2, and GRPR genes. High repeatability (86%) of results was observed, as 12 identified selection signals were also confirmed in other studies with sheep. HRiD offers an interesting possibility to be used complementary to other approaches or when only males are genotyped, which is often the case in genomic breeding value estimations. These results highlight the importance of the X-chromosome in the adaptive architecture of domestic ruminants, while our novel HRiD approach opens new possibilities for research.

Investigation of the Genetic Architecture of Pigs Subjected to Breeding Intensification.

Pigs are strategically important animals for the agricultural industry. An assessment of genetic differentiation between pigs, undergone and not undergone to selection intensification, is of particular interest. Our research was conducted on two groups of Large White pigs grown on the same farm but in different years. A total of 165 samples were selected with 78 LW_А (n = 78, the Russian selection) and LW_B (n = 87, a commercial livestock). For genotyping, we used GeneSeek® GGP Porcine HD Genomic Profiler v1 (Illumina Inc, San Diego, CA, USA). To define breeding characteristics of selection, we used smoothing FST and segment identification of HBD (Homozygous-by-Descent). The results of smoothing FST showed 20 areas of a genome with strong ejection regions of the genome located on all chromosomes except SSC2, SSC3, and SSC8. The average realized autozygosity in Large White pigs of native selection was in (LW_A)-0.21, in LW_В-0.29. LW_А showed 13,338 HBD segments, 171 per one animal, and LW_B-15,747 HBD segments, 181 per one animal. The ejections found by the smoothing FST method were partially localized in the HBD regions. In these areas, the genes ((NCBP1, PLPPR1, GRIN3A, NBEA, TRPC4, HS6ST3, NALCN, SMG6, TTC3, KCNJ6, IKZF2, OBSL1, CARD10, ETV6, VWF, CCND2, TSPAN9, CDH13, CEP128, SERPINA11, PIK3CG, COG5, BCAP29, SLC26A4) were defined. The revealed genes can be of special interest for further studying their influence on an organism of an animal since they can act as candidate genes for selection-significant traits.

Genome-wide analysis of genetic diversity and artificial selection in Large White pigs in Russia.

Breeding practices adopted at different farms are aimed at maximizing the profitability of pig farming. In this work, we have analyzed the genetic diversity of Large White pigs in Russia. We compared genomes of historic and modern Large White Russian breeds using 271 pig samples. We have identified 120 candidate regions associated with the differentiation of modern and historic pigs and analyzed genomic differences between the modern farms. The identified genes were associated with height, fitness, conformation, reproductive performance, and meat quality.

The landscape of gene-CDS-haplotype diversity in rice: Properties, population organization, footprints of domestication and breeding, and implications for genetic improvement.

Polymorphisms within gene coding regions represent the most important part of the overall genetic diversity of rice. We characterized the gene-coding sequence-haplotype (gcHap) diversity of 45 963 rice genes in 3010 rice accessions. With an average of 226 ± 390 gcHaps per gene in rice populations, rice genes could be classified into three main categories: 12 865 conserved genes, 10 254 subspecific differentiating genes, and 22 844 remaining genes. We found that 39 218 rice genes carry >255 179 major gcHaps of potential functional importance. Most (87.5%) of the detected gcHaps were specific to subspecies or populations. The inferred proto-ancestors of local landrace populations reconstructed from conserved predominant (ancient) gcHaps correlated strongly with wild rice accessions from the same geographic regions, supporting a multiorigin (domestication) model of Oryza sativa. Past breeding efforts generally increased the gcHap diversity of modern varieties and caused significant frequency shifts in predominant gcHaps of 14 266 genes due to independent selection in the two subspecies. Low frequencies of "favorable" gcHaps at most known genes related to rice yield in modern varieties suggest huge potential for rice improvement by mining and pyramiding of favorable gcHaps. The gcHap data were demonstrated to have greater power than SNPs for the detection of causal genes that affect complex traits. The rice gcHap diversity dataset generated in this study would facilitate rice basic research and improvement in the future.

Local ancestry prediction with PyLAE.

SUMMARY: We developed PyLAE, a new tool for determining local ancestry along a genome using whole-genome sequencing data or high-density genotyping experiments. PyLAE can process an arbitrarily large number of ancestral populations (with or without an informative prior). Since PyLAE does not involve estimating many parameters, it can process thousands of genomes within a day. PyLAE can run on phased or unphased genomic data. We have shown how PyLAE can be applied to the identification of differentially enriched pathways between populations. The local ancestry approach results in higher enrichment scores compared to whole-genome approaches. We benchmarked PyLAE using the 1000 Genomes dataset, comparing the aggregated predictions with the global admixture results and the current gold standard program RFMix. Computational efficiency, minimal requirements for data pre-processing, straightforward presentation of results, and ease of installation make PyLAE a valuable tool to study admixed populations. AVAILABILITY AND IMPLEMENTATION: The source code and installation manual are available at https://github.com/smetam/pylae.

Signatures of Selection From Candidate Gene Polymorphisms Panel in Five Zebu Breeds

ABSTRACT The signals of selection using candidate genes polymorphism were studied in five zebu breeds of Mexico. Three loci from GHRH and complementarily Steroyl Co Desaturase F762, Dopamine Beta Hydroxilase 17299, and LEP3272 were identified under selection. Findings depict Zebu selection pressure mainly on Brahman breed resulting in a divergent structure pattern.

Signatures of Selection for Environmental Adaptation and Zebu × Taurine Hybrid Fitness in East African Shorthorn Zebu.

The East African Shorthorn Zebu (EASZ) cattle are ancient hybrid between Asian zebu × African taurine cattle preferred by local farmers due to their adaptability to the African environment. The genetic controls of these adaptabilities are not clearly understood yet. Here, we genotyped 92 EASZ samples from Kenya (KEASZ) with more than 770,000 SNPs and sequenced the genome of a pool of 10 KEASZ. We observe an even admixed autosomal zebu × taurine genomic structure in the population. A total of 101 and 165 candidate regions of positive selection, based on genome-wide SNP analyses (meta-SS, Rsb, iHS, and ΔAF) and pooled heterozygosity (Hp) full genome sequence analysis, are identified, in which 35 regions are shared between them. A total of 142 functional variants, one novel, have been detected within these regions, in which 30 and 26 were classified as of zebu and African taurine origins, respectively. High density genome-wide SNP analysis of zebu × taurine admixed cattle populations from Uganda and Nigeria show that 25 of these regions are shared between KEASZ and Uganda cattle, and seven regions are shared across the KEASZ, Uganda, and Nigeria cattle. The identification of common candidate regions allows us to fine map 18 regions. These regions intersect with genes and QTL associated with reproduction and environmental stress (e.g., immunity and heat stress) suggesting that the genome of the zebu × taurine admixed cattle has been uniquely selected to maximize hybrid fitness both in terms of reproduction and survivability.

Towards a whole-genome sequence for rye (Secale cereale L.).

We report on a whole-genome draft sequence of rye (Secale cereale L.). Rye is a diploid Triticeae species closely related to wheat and barley, and an important crop for food and feed in Central and Eastern Europe. Through whole-genome shotgun sequencing of the 7.9-Gbp genome of the winter rye inbred line Lo7 we obtained a de novo assembly represented by 1.29 million scaffolds covering a total length of 2.8 Gbp. Our reference sequence represents nearly the entire low-copy portion of the rye genome. This genome assembly was used to predict 27 784 rye gene models based on homology to sequenced grass genomes. Through resequencing of 10 rye inbred lines and one accession of the wild relative S. vavilovii, we discovered more than 90 million single nucleotide variants and short insertions/deletions in the rye genome. From these variants, we developed the high-density Rye600k genotyping array with 600 843 markers, which enabled anchoring the sequence contigs along a high-density genetic map and establishing a synteny-based virtual gene order. Genotyping data were used to characterize the diversity of rye breeding pools and genetic resources, and to obtain a genome-wide map of selection signals differentiating the divergent gene pools. This rye whole-genome sequence closes a gap in Triticeae genome research, and will be highly valuable for comparative genomics, functional studies and genome-based breeding in rye.

Genetic diversity analyses reveal first insights into breed-specific selection signatures within Swiss goat breeds.

We used genotype data from the caprine 50k Illumina BeadChip for the assessment of genetic diversity within and between 10 local Swiss goat breeds. Three different cluster methods allowed the goat samples to be assigned to the respective breed groups, whilst the samples of Nera Verzasca and Tessin Grey goats could not be differentiated from each other. The results of the different genetic diversity measures show that Appenzell, Toggenburg, Valais and Booted goats should be prioritized in future conservation activities. Furthermore, we examined runs of homozygosity (ROH) and compared genomic inbreeding coefficients based on ROH (FROH ) with pedigree-based inbreeding coefficients (FPED ). The linear relationship between FROH and FPED was confirmed for goats by including samples from the three main breeds (Saanen, Chamois and Toggenburg goats). FROH appears to be a suitable measure for describing levels of inbreeding in goat breeds with missing pedigree information. Finally, we derived selection signatures between the breeds. We report a total of 384 putative selection signals. The 25 most significant windows contained genes known for traits such as: coat color variation (MITF, KIT, ASIP), growth (IGF2, IGF2R, HRAS, FGFR3) and milk composition (PITX2). Several other putative genes involved in the formation of populations, which might have been selected for adaptation to the alpine environment, are highlighted. The results provide a contemporary background for the management of genetic diversity in local Swiss goat breeds.

An interpretive review of selective sweep studies in Bos taurus cattle populations: identification of unique and shared selection signals across breeds.

This review compiles the results of 21 genomic studies of European Bos taurus breeds and thus provides a general picture of the selection signatures in taurine cattle identified by genome-wide selection-mapping scans. By performing a comprehensive summary of the results reported in the literature, we compiled a list of 1049 selection sweeps described across 37 cattle breeds (17 beef breeds, 14 dairy breeds, and 6 dual-purpose breeds), and four different beef-vs.-dairy comparisons, which we subsequently grouped into core selective sweep (CSS) regions, defined as consecutive signals within 1 Mb of each other. We defined a total of 409 CSSs across the 29 bovine autosomes, 232 (57%) of which were associated with a single-breed (Single-breed CSSs), 134 CSSs (33%) were associated with a limited number of breeds (Two-to-Four-breed CSSs) and 39 CSSs (9%) were associated with five or more breeds (Multi-breed CSSs). For each CSS, we performed a candidate gene survey that identified 291 genes within the CSS intervals (from the total list of 5183 BioMart-extracted genes) linked to dairy and meat production, stature, and coat color traits. A complementary functional enrichment analysis of the CSS positional candidates highlighted other genes related to pathways underlying behavior, immune response, and reproductive traits. The Single-breed CSSs revealed an over-representation of genes related to dairy and beef production, this was further supported by over-representation of production-related pathway terms in these regions based on a functional enrichment analysis. Overall, this review provides a comparative map of the selection sweeps reported in European cattle breeds and presents for the first time a characterization of the selection sweeps that are found in individual breeds. Based on their uniqueness, these breed-specific signals could be considered as "divergence signals," which may be useful in characterizing and protecting livestock genetic diversity.