Whole-genome sequencing of a worldwide collection of sugarcane cultivars (Saccharum spp.) reveals the genetic basis of cultivar improvement.

Sugarcane is the main source of sugar worldwide, and 80% of the sucrose production comes from sugarcane. However, the genetic differentiation and basis of agronomic traits remain obscure. Here, we sequenced the whole-genome of 219 elite worldwide sugarcane cultivar accessions. A total of approximately 6 million high-quality genome-wide single nucleotide polymorphisms (SNPs) were detected. A genome-wide association study identified a total of 2198 SNPs that were significantly associated with sucrose content, stalk number, plant height, stalk diameter, cane yield, and sugar yield. We observed homozygous tendency of favor alleles of these loci, and over 80% of cultivar accessions carried the favor alleles of the SNPs or haplotypes associated with sucrose content. Gene introgression analysis showed that the number of chromosome segments from Saccharum spontaneum decreased with the breeding time of cultivars, while those from S. officinarum increased in recent cultivars. A series of selection signatures were identified in sugarcane improvement procession, of which 104 were simultaneously associated with agronomic traits and 45 of them were mainly associated with sucrose content. We further proposed that as per sugarcane transgenic experiments, ShN/AINV3.1 plays a positive role in increasing stalk number, plant height, and stalk diameter. These findings provide comprehensive resources for understanding the genetic basis of agronomic traits and will be beneficial to germplasm innovation, screening molecular markers, and future sugarcane cultivar improvement.

Genetic structural analysis of different breeds and geographical groups of Fenneropenaeus chinensis reveals population diversity.

Fenneropenaeus chinensis is a commercially important shrimp species cultured in China. This study investigated eight F. chinensis populations in China, including four geographical populations, three commercial breeds, and one wild population captured from the Yellow Sea. Population stratification analysis revealed that the Hebei geographical population and commercial breeding "Huanghai No. 4" were relatively independent and stable, reflecting a relatively closed breeding environment, whereas gene introgression was present between other populations. Selective signature analysis detected artificial selection for vision, growth, and disease resistance in the Hebei population. Neuronal development-related genes were detected to be under selection in the Changyi and Rizhao populations. Fertility of the Rizhao population was also investigated. Additionally, genes in the glycosaminoglycan biosynthesis-chondroitin sulfate/dermatan sulfate pathway were involved in the high pH tolerance of the "Huanghai No. 4" population. This study provided support for the genetic mechanism of parsing economic traits and the development of molecular breeding technologies.

Genome-wide comparative analyses highlight selection signatures underlying saline adaptation in Chilika buffalo.

Chilika, a native buffalo breed of the Eastern coast of India, is mainly distributed around the Chilika brackish water lake connected with the Bay of Bengal Sea. This breed possesses a unique ability to delve deep into the salty water of the lake and stay there to feed on local vegetation of saline nature. Adaptation to salinity is a genetic phenomenon; however, the genetic basis underlying salinity tolerance is still limited in animals, specifically in livestock. The present study explores the genetic evolution that unveils the Chilika buffalo's adaptation to the harsh saline habitat, including both water and food systems. For this study, whole genome resequencing data on 18 Chilika buffalo and for comparison 10 Murrah buffalo of normal habitat were generated. For identification of selection sweeps, intrapopulation and interpopulation statistics were used. A total of 709, 309, 468, and 354 genes were detected to possess selection sweeps in Chilika buffalo using the nucleotide diversity (θπ), Tajima's D, nucleotide diversity ratio (θπ-ratio), and FST methods, respectively. Further analysis revealed a total of 23 genes including EXOC6B, VPS8, LYPD1, VPS35, CAMKMT, NCKAP5, COMMD1, myosin light chain kinase 3 (MYLK3), and B3GNT2 were found to be common by all the methods. Furthermore, functional annotation study of identified genes provided pathways such as MAPK signaling, renin secretion, endocytosis, oxytocin signaling pathway, etc. Gene network analysis enlists that hub genes provide insights into their interactions with each other. In conclusion, this study has highlighted the genetic basis underlying the local adaptive function of Chilika buffalo under saline environment.NEW & NOTEWORTHY Indian Chilika buffaloes are being maintained on extensive grazing system and have a unique ability to convert local salty vegetation into valuable human food. However, adaptability to saline habitat of Chilika buffalo has not been explored to date. Here, we identified genes and biological pathways involved, such as MAPK signaling, renin secretion, endocytosis, and oxytocin signaling pathway, underlying adaptability of Chilika buffalo to saline environment. This investigation shed light on the mechanisms underlying the buffalo's resilience in its native surroundings.

Selection signatures of wool color in Gangba sheep revealed by genome-wide SNP discovery.

BACKGROUND: Gangba sheep as a famous breed of Tibetan sheep, its wool color is mainly white and black. Gangba wool is economically important as a high-quality raw material for Tibetan blankets and Tibetan serge. However, relatively few studies have been conducted on the wool color of Tibetan sheep. RESULTS: To fill this research gap, this study conducted an in-depth analysis of two populations of Gangba sheep (black and white wool color) using whole genome resequencing to identify genetic variation associated with wool color. Utilizing PCA, Genetic Admixture, and N-J Tree analyses, the present study revealed a consistent genetic relationship and structure between black and white wool colored Gangba sheep populations, which is consistent with their breed history. Analysis of selection signatures using multiple methods (FST, π ratio, Tajima's D), 370 candidate genes were screened in the black wool group (GBB vs GBW); among them, MC1R, MLPH, SPIRE2, RAB17, SMARCA4, IRF4, CAV1, USP7, TP53, MYO6, MITF, MC2R, TET2, NF1, JAK1, GABRR1 genes are mainly associated with melanin synthesis, melanin delivery, and distribution. The enrichment results of the candidate genes identified 35 GO entries and 19 KEGG pathways associated with the formation of the black phenotype. 311 candidate genes were screened in the white wool group (GBW vs GBB); among them, REST, POU2F1, ADCY10, CCNB1, EP300, BRD4, GLI3, and SDHA genes were mainly associated with interfering with the differentiation of neural crest cells into melanocytes, affecting the proliferation of melanocytes, and inhibiting melanin synthesis. 31 GO entries and 22 KEGG pathways were associated with the formation of the white phenotype. CONCLUSIONS: This study provides important information for understanding the genetic mechanism of wool color in Gangba, and provides genetic knowledge for improving and optimizing the wool color of Tibetan sheep. Genetic improvement and selective breeding to produce wool of specific colors can meet the demand for a diversity of wool products in the Tibetan wool textile market.

Identification of the genetic characteristics of copy number variations in experimental specific pathogen-free ducks using whole-genome resequencing.

BACKGROUND: Specific pathogen-free ducks are a valuable laboratory resource for waterfowl disease research and poultry vaccine development. High throughput sequencing allows the systematic identification of structural variants in genomes. Copy number variation (CNV) can explain the variation of important duck genetic traits. Herein, the genome-wide CNVs of the three experimental duck species in China (Jinding ducks (JD), Shaoxing ducks (SX), and Fujian Shanma ducks (SM)) were characterized using resequencing to determine their genetic characteristics and selection signatures. RESULTS: We obtained 4,810 CNV regions (CNVRs) by merging 73,012 CNVs, covering 4.2% of the duck genome. Functional analysis revealed that the shared CNVR-harbored genes were significantly enriched for 31 gene ontology terms and 16 Kyoto Encyclopedia of Genes and Genomes pathways (e.g., olfactory transduction and immune system). Based on the genome-wide fixation index for each CNVR, growth (SPAG17 and PTH1R), disease resistance (CATHL3 and DMBT1), and thermoregulation (TRPC4 and SLIT3) candidate genes were identified in strongly selected signatures specific to JD, SM, and SX, respectively. CONCLUSIONS: In conclusion, we investigated the genome-wide distribution of experimental duck CNVs, providing a reference to establish the genetic basis of different phenotypic traits, thus contributing to the management of experimental animal genetic resources.

Identification of genomic characteristics and selective signals in Guizhou black goat.

BACKGROUND: Guizhou black goat is one of the indigenous black goat breeds in the southwest region of Guizhou, China, which is an ordinary goat for mutton production. They are characterized by moderate body size, black coat, favorite meat quality with tender meat and lower odor, and tolerance for cold and crude feed. However, little is known about the genetic characteristics or variations underlying their important economic traits. RESULTS: Here, we resequenced the whole genome of Guizhou black goat from 30 unrelated individuals breeding in the five core farms. A total of 9,835,610 SNPs were detected, and 2,178,818 SNPs were identified specifically in this breed. The population structure analysis revealed that Guizhou black goat shared a common ancestry with Shaanbei white cashmere goat (0.146), Yunshang black goat (0.103), Iran indigenous goat (0.054), and Moroccan goat (0.002). However, Guizhou black goat showed relatively higher genetic diversity and a lower level of linkage disequilibrium than the other seven goat breeds by the analysis of the nucleotide diversity, linkage disequilibrium decay, and runs of homozygosity. Based on FST and θπ values, we identified 645, 813, and 804 selected regions between Guizhou black goat and Yunshang black goat, Iran indigenous goat, and cashmere goats. Combined with the results of XP-EHH, there were 286, 322, and 359 candidate genes, respectively. Functional annotation analysis revealed that these genes are potentially responsible for the immune response (e.g., CD28, CD274, IL1A, TLR2, and SLC25A31), humility-cold resistance (e.g., HBEGF, SOSTDC1, ARNT, COL4A1/2, and EP300), meat quality traits (e.g., CHUK, GAB2, PLAAT3, and EP300), growth (e.g., GAB2, DPYD, and CSF1), fertility (e.g., METTL15 and MEI1), and visual function (e.g., PANK2 and NMNAT2) in Guizhou black goat. CONCLUSION: Our results indicated that Guizhou black goat had a high level of genomic diversity and a low level of linkage disequilibrium in the whole genome. Selection signatures were detected in the genomic regions that were mainly related to growth and development, meat quality, reproduction, disease resistance, and humidity-cold resistance in Guizhou black goat. These results would provide a basis for further resource protection and breeding improvement of this very local breed.

Genomic patterns of selection in morphometric traits across diverse Indian cattle breeds.

This study seeks a comprehensive exploration of genome-wide selective processes impacting morphometric traits across diverse cattle breeds, utilizing an array of statistical methods. Morphometric traits, encompassing both qualitative and quantitative variables, play a pivotal role in characterizing and selecting livestock breeds based on their external appearance, size, and physical attributes. While qualitative traits, such as color, horn structure, and coat type, contribute to adaptive features and breed identification, quantitative traits like body weight and conformation measurements bear a closer correlation with production characteristics. This study employs advanced genotyping technologies, including the Illumina BovineSNP50 Bead Chip and next-generation sequencing methods like Reduced Representation sequencing, to identify genomic signatures associated with these traits. We applied four intra-population methods to find evidence of selection, such as Tajima's D, CLR, iHS, and ROH. We found a total of 40 genes under the selection signature, that were associated with morphometric traits in five cattle breeds (Kankrej, Tharparkar, Nelore, Sahiwal, and Gir). Crucial genes such as ADIPDQ, DPP6, INSIG1, SLC35D2 in Kankrej, LPL, ATP6V1B2, CDC14B in Tharparkar, HPSE2, PLAG1 in Nelore, PCSK1, PRKD1 in Sahiwal, and GNAQ, HPCAL1 in Gir were identified in our study. This approach provides valuable insights into the genetic basis of variations in body weight and conformation traits, facilitating informed selection processes and offering a deeper understanding of the evolutionary and domestication processes in diverse cattle breeds.

A comprehensive review of livestock development: insights into domestication, phylogenetics, diversity, and genomic advances.

Livestock plays an essential role in sustaining human livelihoods, offering a diverse range of species integral to food security, economic stability, and cultural traditions. The domestication of livestock, which began over 10,000 years ago, has driven significant genetic changes in species such as cattle, buffaloes, sheep, goats, and pigs. Recent advancements in genomic technologies, including next-generation sequencing (NGS), genome-wide association studies (GWAS), and genomic selection, have dramatically enhanced our understanding of these genetic developments. This review brings together key research on the domestication process, phylogenetics, genetic diversity, and selection signatures within major livestock species. It emphasizes the importance of admixture studies and evolutionary forces like natural selection, genetic drift, and gene flow in shaping livestock populations. Additionally, the integration of machine learning with genomic data offers new perspectives on the functional roles of genes in adaptation and evolution. By exploring these genomic advancements, this review provides insights into genetic variation and evolutionary processes that could inform future approaches to improving livestock management and adaptation to environmental challenges, including climate change.

Genome-wide analyses based on a novel donkey 40K liquid chip reveal the gene responsible for coat color diversity in Chinese Dezhou donkey.

Dezhou donkey is one of the representative local breeds in China, which is mainly divided into two strains: Sanfen and Wutou. There are obvious differences in coat color between the two strains. The former shows light points around the eyes, around the muzzle and under the belly, while the latter is completely solid black. In this study, genome-wide association analysis was performed for the differences in coat color traits between the Sanfen (n = 97) and Wutou (n = 108) strains using a novel donkey 40K liquid chip developed based on GenoBaits technology, to identify genomic regions and causal genes that could explain this variation. We also used FST and The cross-population composite likelihood ratio test (XPCLR) analyses to explore selected regions related to coat color differences. We identified one significant region on chromosome 15, with the most significant SNP located within the agouti signaling protein (ASIP) gene. At the same time, both FST and XPCLR methods detected the same selected region on chromosome 15, and ASIP was the gene with the strongest signal. ASIP and melanocortin 1 receptor (MC1R) control the ratio of eumelanin to pheomelanin through their protein activity. They are deeply involved in the process of melanosome organation and melanogenesis, thus affecting mammals' coat color variation. We used a range of genome-wide approach to identify the genetic basis of coat color variation in Dezhou donkeys. The results provide a supplement to the color variation study in Chinese donkeys at the genome-wide level, and preliminarily verified the reliability of the Molbreeding Donkey No. 1 40K liquid chip.

Genomic characterization of dryland adaptation in endangered Anxi cattle in China.

Local species exhibit distinctive indigenous characteristics while showing unique productive and phenotypic traits. However, the advent of commercialization has posed a substantial threat to the survival of indigenous species. Anxi cattle, an endangered native breed in China, have evolved unique growth and reproductive characteristics in extreme desert and semidesert ecosystems. In this study, we conducted a genomic comparison of 10 Anxi cattle genomes with those of five other global populations/breeds to assess genetic diversity and identify candidate genomic regions in Anxi cattle. Population structure and genetic diversity analyses revealed that Anxi cattle are part of the East Asian cattle clade, exhibiting higher genetic diversity than commercial breeds. Through selective sweep analysis, we identified specific genetic variations linked to the environmental adaptability of Anxi cattle. Notably, we identified several candidate genes, including CERS3 involved in regulating skin permeability and antimicrobial functions, RBFOX2 associated with cardiac development, SLC16A7 participated in the regulation of pancreatic endocrine function, and SPATA3 related to reproduction. Our findings revealed the distinctive genomic features of Anxi cattle in dryland environments, provided invaluable insights for further research and breed preservation, and had important significance for enriching the domestic cattle breeding gene bank.

Assessment of genetic diversity, inbreeding and collection completeness of Jersey bulls in the US National Animal Germplasm Program.

Genomic selection and extensive use of a few elite bulls through artificial insemination are leading to reduced genetic diversity in Jersey cattle. Conservation of genetic diversity through gene banks can protect a breed's genetic diversity and genetic gain, ensuring continued genetic advancement in the future. The availability of genomic information in the US National Animal Germplasm Program (NAGP) facilitates characterization of Jersey bulls in the germplasm collection. Therefore, in this study, we compared the genetic diversity and inbreeding between Jersey bulls in the NAGP and the national cooperator database (NCD). The NCD is maintained and curated by the Council on Dairy Cattle Breeding (CDCB). We found the genetic diversity to be marginally higher in NAGP (Ho = 0.34 ± 0.17) relative to the NCD population (Ho = 0.33 ± 0.16). The average pedigree and genomic inbreeding (FPED, FGRM, FROH > 2Mb) were similar between the groups, with estimates of 7.6% with FPED, 11.07% with FGRM and 20.13% with FROH > 2Mb. An increasing trend in inbreeding was detected, and a significantly higher level of inbreeding was estimated among the older bulls in the NAGP collection, suggesting an overrepresentation of the genetics from elite bulls. Results from principal component analyses (PCA) provided evidence that the NAGP collection is representative of the genetic variation found in the NCD population and a broad majority of the loci segregating (98.2%) in the NCD population were also segregating in the NAGP. Ward's clustering was used to assess collection completeness of Jerseys in the NAGP by comparison with top 1000 sires of bulls, top 1000 sires of cow, and bulls with high Lifetime Net Merit (NM$). All the clusters were represented in the NAGP suggesting that most of the genetic diversity in the US Jersey population is represented in the NAGP and confirmed the PCA results. The decade of birth was the major driver grouping bulls into clusters, suggesting the importance of selection over time. Selection signature analysis between the historic bulls in the NAGP with the newer bulls, born in the decade after implementation of genomic selection, identified selection for milk production, fat and protein yield, fertility, health, and reproductive traits. Cluster analysis revealed that the NAGP has captured allele frequency changes over time associated with selection, validating the strategy of repeated sampling and suggests that the continuation of a repeated sampling policy is essential for the germplasm collection to maintain its future utility. While NAGP should continue to collect bulls that have large influence on the population due to selection, care should be taken to include the entire breadth of bulls, including low merit bulls.

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 Scan Uncovers Candidate Genes Related to Milk Production Traits in Barka Cattle.

In this study, our primary aim was to explore the genomic landscape of Barka cattle, a breed recognized for high milk production in a semi-arid environment, by focusing on genes with known roles in milk production traits. We employed genome-wide analysis and three selective sweep detection methods (ZFST, θπ ratio, and ZHp) to identify candidate genes associated with milk production and composition traits. Notably, ACAA1, P4HTM, and SLC4A4 were consistently identified by all methods. Functional annotation highlighted their roles in crucial biological processes such as fatty acid metabolism, mammary gland development, and milk protein synthesis. These findings contribute to understanding the genetic basis of milk production in Barka cattle, presenting opportunities for enhancing dairy cattle production in tropical climates. Further validation through genome-wide association studies and transcriptomic analyses is essential to fully exploit these candidate genes for selective breeding and genetic improvement in tropical dairy cattle.

Integrating Multiple Database Resources to Elucidate the Gene Flow in Southeast Asian Pig Populations.

The domestic pig (Sus scrofa) and its subfamilies have experienced long-term and extensive gene flow, particularly in Southeast Asia. Here, we analyzed 236 pigs, focusing on Yunnan indigenous, European commercial, East Asian, and Southeast Asian breeds, using the Pig Genomics Reference Panel (PGRP v1) of Pig Genotype-Tissue Expression (PigGTEx) to investigate gene flow and associated complex traits by integrating multiple database resources. In this study, we discovered evidence of admixtures from European pigs into the genome of Yunnan indigenous pigs. Additionally, we hypothesized that a potential conceptual gene flow route that may have contributed to the genetic composition of the Diannan small-ear pig is a gene exchange from the Vietnamese pig. Based on the most stringent gene introgression scan using the fd statistic, we identified three specific loci on chromosome 8, ranging from 51.65 to 52.45 Mb, which exhibited strong signatures of selection and harbored the NAF1, NPY1R, and NPY5R genes. These genes are associated with complex traits, such as fat mass, immunity, and litter weight, in pigs, as supported by multiple bio-functionalization databases. We utilized multiple databases to explore the potential dynamics of genetic exchange in Southeast Asian pig populations and elucidated specific gene functionalities.

Genomic profiling of selective sweeps through haplotype differentiation unravelled genes associated with production and reproduction traits in Indian goat breeds.

A comprehensive genomic scan of selective sweeps was conducted in autochthonous Attappady Black and improved dual-purpose Malabari goat breeds in south India. High-throughput single nucleotide polymorphism (SNP) marker data, obtained through Illumina goat SNP50 BeadChip genotyping of 48 goats (24 each of Attappady Black and Malabari goats), were utilized for the analysis. Selection signature analysis, employing hapFLK analysis based on haplotype differentiation, identified seven significant sweep regions (p < 0.005). Notably, one of these regions encompassed the genomic area housing the casein cluster and quantitative trait loci associated with milk production on chromosome 6. Gene ontology enrichment analysis of 166 putative selective genes associated with these sweep regions revealed 13 significantly over-represented Panther pathways (p ≤ 0.05), including the TGF-beta signalling pathway and GNRHR pathway. The selective sweeps detected in this study contributed significantly to the phenotypic divergence observed between Attappady Black and Malabari goats in south India.

Deciphering climate resilience in Indian cattle breeds by selection signature analyses.

The signature of selection is a crucial concept in evolutionary biology that refers to the pattern of genetic variation which arises in a population due to natural selection. In the context of climate adaptation, the signature of selection can reveal the genetic basis of adaptive traits that enable organisms to survive and thrive in changing environmental conditions. Breeds living in diverse agroecological zones exhibit genetic "footprints" within their genomes that mirror the influence of climate-induced selective pressures, subsequently impacting phenotypic variance. It is assumed that the genomes of animals residing in these regions have been altered through selection for various climatic adaptations. These regions are known as signatures of selection and can be identified using various summary statistics. We examined genotypic data from eight different cattle breeds (Gir, Hariana, Kankrej, Nelore, Ongole, Red Sindhi, Sahiwal, and Tharparkar) that are adapted to diverse regional climates. To identify selection signature regions in this investigation, we used four intra-population statistics: Tajima's D, CLR, iHS, and ROH. In this study, we utilized Bovine 50 K chip data and four genome scan techniques to assess the genetic regions of positive selection for high-temperature adaptation. We have also performed a genome-wide investigation of genetic diversity, inbreeding, and effective population size in our target dataset. We identified potential regions for selection that are likely to be caused by adverse climatic conditions. We observed many adaptation genes in several potential selection signature areas. These include genes like HSPB2, HSPB3, HSP20, HSP90AB1, HSF4, HSPA1B, CLPB, GAP43, MITF, and MCHR1 which have been reported in the cattle populations that live in varied climatic regions. The findings demonstrated that genes involved in disease resistance and thermotolerance were subjected to intense selection. The findings have implications for marker-assisted breeding, understanding the genetic landscape of climate-induced adaptation, putting breeding and conservation programs into action.

Genome-wide landscape of runs of homozygosity and differentiation across Egyptian goat breeds.

Understanding the genomic features of livestock is essential for successful breeding programs and conservation. This information is scarce for local goat breeds in Egypt. In the current study, genomic regions with selection signatures were identified as well as runs of homozygosity (ROH), genomic inbreeding coefficients (FROH) and fixation index (FST) were detected in Egyptian Nubian, Damascus, Barki and Boer goat breeds. A total of 46,268 SNP markers and 337 animals were available for the genomic analyses. On average, 145.44, 42.02, 87.90 and 126.95 ROHs were detected per individual in the autosomal genome of the respective breeds. The mean accumulative ROH lengths ranged from 46.5 Mb in Damascus to 360 Mb in Egyptian Nubian. The short ROH segments (< 2 Mb) were most frequent in all breeds, while the longest ROH segments (> 16 Mb) were exclusively found in the Egyptian Nubian. The highest average FROH was observed in Egyptian Nubian (~ 0.12) followed by Boer (~ 0.11), while the lowest FROH was found in Damascus (~ 0.05) and Barki breed (~ 0.03). The estimated mean FST was 0.14 (Egyptian Nubian and Boer), 0.077 (Egyptian Nubian and Barki), 0.075 (Egyptian Nubian and Damascus), 0.071 (Barki and Boer), 0.064 (Damascus and Boer), and 0.015 (Damascus and Barki), for each pair of breeds. Interestingly, multiple SNPs that accounted for high FST values were observed on chromosome 6 in regions harboring ALPK1 and KCNIP4. Genomic regions overlapping both FST and ROH harbor genes related to immunity (IL4R, PHF23, GABARAP, GPS2, and CD68), reproduction (SPATA2L, TNFSF12, TMEM95, and RNF17), embryonic development (TCF25 and SOX15) and adaptation (MC1R, KDR, and KIT), suggesting potential genetic adaptations to local environmental conditions. Our results contribute to the understanding of the genetic architecture of different goat breeds and may provide valuable information for effective preservation and breeding programs of local goat breeds in Egypt.

Whole-genome sequencing revealed genetic diversity, structure and patterns of selection in Guizhou indigenous chickens.

BACKGROUND: The eight phenotypically distinguishable indigenous chicken breeds in Guizhou province of China are great resources for high-quality development of the poultry industry in China. However, their full value and potential have yet to be understood in depth. To illustrate the genetic diversity, the relationship and population structure, and the genetic variation patterns shaped by selection in Guizhou indigenous chickens, we performed a genome-wide analysis of 240 chickens from 8 phenotypically and geographically representative Guizhou chicken breeds and 60 chickens from 2 commercial chicken breeds (one broiler and one layer), together with 10 red jungle fowls (RJF) genomes available from previous studies. RESULTS: The results obtained in this present study showed that Guizhou chicken breed populations harbored higher genetic diversity as compared to commercial chicken breeds, however unequal polymorphisms were present within Guizhou indigenous chicken breeds. The results from the population structure analysis markedly reflected the breeding history and the geographical distribution of Guizhou indigenous chickens, whereas, some breeds with complex genetic structure were ungrouped into one cluster. In addition, we confirmed mutual introgression within Guizhou indigenous chicken breeds and from commercial chicken breeds. Furthermore, selective sweep analysis revealed candidate genes which were associated with specific and common phenotypic characteristics evolved rapidly after domestication of Guizhou local chicken breeds and economic traits such as egg production performance, growth performance, and body size. CONCLUSION: Taken together, the results obtained from the comprehensive analysis of the genetic diversity, genetic relationships and population structures in this study showed that Guizhou indigenous chicken breeds harbor great potential for commercial utilization, however effective conservation measures are currently needed. Additionally, the present study drew a genome-wide selection signature draft for eight Guizhou indigenous chicken breeds and two commercial breeds, as well as established a resource that can be exploited in chicken breeding programs to manipulate the genes associated with desired phenotypes. Therefore, this study will provide an essential genetic basis for further research, conservation, and breeding of Guizhou indigenous chickens.

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 FST 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.

Genomic regions underlying positive selection in local, Alpine cattle breeds.

We used genome-wide SNP data from 18 local cattle breeds from six countries of the Alpine region to characterize population structure and identify genomic regions underlying positive selection. The geographically close breeds Evolèner, Eringer, Valdostana Pezzata Nera, and Valdostana Castana were found to differ from all other Alpine breeds. In addition, three breeds, Simmental, and Original Braunvieh from Switzerland and Pinzgauer from Austria built three separate clusters. Of the 18 breeds studied, the intra-alpine Swiss breed Evolèner had the highest average inbreeding based on runs of homozygosity (FROH ) and the highest average genomic relationship within the breed. In contrast, Slovenian Cika cattle had the lowest average genomic inbreeding and the lowest average genomic relationship within the breed. We found selection signatures on chromosome 6 near known genes such as KIT and LCORL explaining variation in coat color and body size in cattle. The most prominent selection signatures were similar regardless of marker density and the breeds in the data set. In addition, using available high-density SNP data from 14 of the breeds we identified 47 genome regions as ROH islands. The proportion of homozygous animals was higher in all studied animals of local breeds than in Holstein and Brown Swiss cattle, the two most important commercial breeds in the Alpine region. We report ROH islands near genes related to thermoregulation, coat color, production, and stature. The results of this study serve as a basis for the search for causal variants underlying adaptation to the alpine environment and other specific characteristics selected during the evolution of local Alpine cattle breeds.