Genotyping and Whole-Genome Resequencing of Welsh Sheep Breeds Reveal Candidate Genes and Variants for Adaptation to Local Environment and Socioeconomic Traits

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.

Whole-genome sequencing reveals genomic diversity and selection signatures in Xia'nan cattle

BACKGROUND

The crossbreeding of specialized beef cattle breeds with Chinese indigenous cattle is a common method of genetic improvement. Xia'nan cattle, a crossbreed of Charolais and Nanyang cattle, is China's first specialized beef cattle breed with independent intellectual property rights. After more than two decades of selective breeding, Xia'nan cattle exhibit a robust physique, good environmental adaptability, good tolerance to coarse feed, and high meat production rates. This study analyzed the population genetic structure, genetic diversity, and genomic variations of Xia'nan cattle using whole-genome sequencing data from 30 Xia'nan cattle and 178 published cattle genomic data.

RESULT

The ancestry estimating composition analysis showed that the ancestry proportions for Xia'nan cattle were mainly Charolais with a small amount of Nanyang cattle. Through the genetic diversity studies (nucleotide diversity and linkage disequilibrium decay), we found that the genomic diversity of Xia'nan cattle is higher than that of specialized beef cattle breeds in Europe but lower than that of Chinese native cattle. Then, we used four methods to detect genome candidate regions influencing the excellent traits of Xia'nan cattle. Among the detected results, 42 genes (θπ and CLR) and 131 genes (FST and XP-EHH) were detected by two different detection strategies. In addition, we found a region in BTA8 with strong selection signals. Finally, we conducted functional annotation on the detected genes and found that these genes may influence body development (NR6A1), meat quality traits (MCCC1), growth traits (WSCD1, TMEM68, MFN1, NCKAP5), and immunity (IL11RA, CNTFR, CCL27, SLAMF1, SLAMF7, NAA35, and GOLM1).

CONCLUSION

We elucidated the genomic features and population structure of Xia'nan cattle and detected some selection signals in genomic regions potentially associated with crucial economic traits in Xia'nan cattle. This research provided a basis for further breeding improvements in Xia'nan cattle and served as a reference for genetic enhancements in other crossbreed cattle.

Whole-genome resequencing reveals the uniqueness of Subei yak

Subei yak is an essential local yak in the Gansu Province, which genetic resource has recently been discovered. It is a meat-milk dual-purpose variety with high fecundity and relatively stable population genetic structure. However, its population genetic structure and genetic diversity are yet to be reported. Therefore, this study aimed to identify molecular markers of Subei yak genome by whole-genome resequencing, and to analyze the population structure and genetic diversity of Subei yak. This study screened 12,079,496 single nucleotide polymorphism (SNP) molecular markers in the 20 Subei yaks genome using whole-genome resequencing technology. Of these SNPs, 32.09% were located in the intronic region of the genome. Principal component analysis, phylogenetic analysis, and population structure analysis revealed that the Subei yak belonged to an independent group in the domestic yak population. A selective clearance analysis was carried out on Subei yak and other domestic yaks, and the genes under positive selection were annotated. The functional enrichment analysis showed that Subei yak possessed prominent selection characteristics in terms of external environment perception, hypoxia adaptation, and muscle development. Furthermore, Subei yak showed excellent muscle fat deposition and meat quality traits. Thus, this study will serve as a reference for discovering population structure, genetic evolution, and other unique traits of Subei yak and for expanding the genetic variation catalog of yaks.

Assessing genetic diversity and defining signatures of positive selection on the genome of dromedary camels from the southeast of the Arabian Peninsula

Dromedary camels (Camelus dromedarius) are members of the Camelini tribe within the Camelidae family. They are distributed throughout North Africa, the Arabian Peninsula and Southeast Asia. This domestic species is characterized by its superior adaptability to the harsh desert environment. In this study, whole autosomal data of 29 dromedary samples from the Southeast Arabian Peninsula in Oman; 10 from Muscat, 14 from Al-Batinah, and 5 from Al-Sharqiya, were investigated to assess their genetic relationship and to define candidate signatures of positive selection. A minimal genetic distinction that separates Muscat dromedaries from the other two populations was observed, with a degree of genetic admixture between them. Using the de-correlated composite of multiple signals (DCMS) approach, a total of 47 candidate regions within the autosomes of these dromedary populations were defined with signatures of positive selection. These candidate regions harbor a total of 154 genes that are mainly associated with functional categories related to immune response, lipid metabolism and energy expenditure, optical and auditory functions, and long-term memory. Different functional genomic variants were called on the candidate regions and respective genes that warrant further investigation to find possible association with the different favorable phenotypes in dromedaries. The output of this study paves the way for further research efforts aimed at defining markers for use in genomic breeding programs, with the goal of conserving the genetic diversity of the species and enhancing its productivity.

Genetic diversity, population structure, and selective signature of sheep in the northeastern Tarim Basin

Local sheep in the northeastern Tarim Basin can adapt to dry and low-rainfall regional environments. In this study, three local sheep breeds in the northeastern Tarim Basin, LOP (LOP) sheep, Bayinbuluke (BYK) sheep, and Kunlun (KUN, also known as the Qiemo sheep) sheep, and three introduced sheep breeds, Suffolk (SUF) sheep, Dorset (APD) sheep, and Texel (TEX) sheep, were analyzed for genetic diversity, population structure, and selective signature using the Illumina OvineSNP50K BeadChip. We found that LOP, BYK, and KUN had lower observed heterozygosity and expected heterozygosity than TEX, SUF, and ADP, which were differentiated based on geographic distribution. We performed fixation index (FST) analysis on three local sheep breeds in the northeastern Tarim Basin (LOP, BYK, and KUN) and introduced sheep breeds (TEX, SUF, and ADP) to measure genetic differentiation. Nucleotide diversity (PI) analysis was performed on single-nucleotide polymorphism (SNP) data of LOP, BYK, and KUN. A total of 493 candidate genes were obtained by taking the intersection at a threshold of 5%. Among them, SMAD2, ESR2, and HAS2 were related to reproductive traits. PCDH15, TLE4, and TFAP2B were related to growth traits. SOD1, TSHR, and DNAJB5 were related to desert environmental adaptation. Analyzing the genetic patterns of local sheep in the northeastern Tarim Basin can protect the germplasm resources of local sheep and promote the development and utilization of sheep genetic resources.

Whole-genome resequencing reveals genetic diversity and selection signals in warm temperate and subtropical Sillago sinica populations

BACKGROUND

Genetic diversity and heterogeneous genomic signatures in marine fish populations may result from selection pressures driven by the strong effects of environmental change. Nearshore fishes are often exposed to complex environments and human activities, especially those with small ranges. However, studies on genetic diversity and population selection signals in these species have mostly been based on a relatively small number of genetic markers. As a newly recorded species of Sillaginidae, the population genetics and genomic selection signals of Sillago sinica are fragmented or even absent.

RESULTS

To address this theoretical gap, we performed whole-genome resequencing of 43 S. sinica individuals from Dongying (DY), Qingdao (QD) and Wenzhou (WZ) populations and obtained 4,878,771 high-quality SNPs. Population genetic analysis showed that the genetic diversity of S. sinica populations was low, but the genetic diversity of the WZ population was higher than that of the other two populations. Interestingly, the three populations were not strictly clustered within the group defined by their sampling location but showed an obvious geographic structure signal from the warm temperate to the subtropics. With further analysis, warm-temperate populations exhibited strong selection signals in genomic regions related to nervous system development, sensory function and immune function. However, subtropical populations showed more selective signalling for environmental tolerance and stress signal transduction.

CONCLUSIONS

Genome-wide SNPs provide high-quality data to support genetic studies and localization of selection signals in S. sinica populations. The reduction in genetic diversity may be related to the bottleneck effect. Considering that low genetic diversity leads to reduced environmental adaptability, conservation efforts and genetic diversity monitoring of this species should be increased in the future. Differences in genomic selection signals between warm temperate and subtropical populations may be related to human activities and changes in environmental complexity. This study deepened the understanding of population genetics and genomic selection signatures in nearshore fishes and provided a theoretical basis for exploring the potential mechanisms of genomic variation in marine fishes driven by environmental selection pressures.

Unravelling the effect of environment on the genome of Sarda breed ewes using Runs of Homozygosity

Natural adaptation and artificial selection have shaped the genome of modern livestock breeds. Among SNP-based metrics that are used to detect signatures of selection at genome-wide level, runs of homozygosity (ROH) are getting increasing popularity. In this paper, ROH distribution and features of a sample of 823 Sarda breed ewes farmed at different levels of altitude are analysed to investigate the effect of the environment on the patterns of homozygosity. A total of 46,829 (33,087 unique) ROH were detected. OAR2 exhibited the largest average number of ROH per animal. The most frequent ROH (OAR27, 38.9-44.2 Mb) was shared by 327. ROH length was statistically affected (p < 0.001) by both the altitude and temperature of the place where the flock was located. The highest probability of a SNP falling in a ROH was observed for hill ewes, whereas the smallest one for mountain. A total of 457 SNP exceeded the 99th percentile of the ROH count per SNP distribution and were considered significant. These markers mapped in eight chromosomes and they clustered into 17 ROH islands, where 80 candidate genes were mapped. Results of this study highlighted differences in the ROH distribution and features among sheep farmed in flocks located at different levels of altitude, confirming the role of environmental adaptability in shaping the genome of this breed.