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.

1 Introduction

The varying altitude gradient on Earth leads to complex and diverse climates, and China’s topography exhibits a terraced structure due to these differences in altitude (Wu et al., 2023). Atushi City in Kezhou, Xinjiang, located in the southwest of the Tarim Basin, lies at the southern foot of the Tianshan Mountains with altitudes ranging from about 1,500 to 2,000 m. This area is characterized by undulating hills and peaks, with the highest point reaching 7,719 m (Qi et al., 2021). Kirghiz sheep, native to high-altitude regions, exhibit strong environmental adaptability and stable growth traits. Predominantly found in Wuqia County, Kezilesu Autonomous Prefecture, in the southern Tianshan Mountains of Xinjiang, these sheep are known for their tolerance to rough feeding, disease resistance, strong stress resistance, and rapid weight gain. They have adapted well to the cold high-altitude environment of the Pamirs in southern Xinjiang. Animals in high-altitude regions undergo complex physiological and behavioral adjustments to adapt to their harsh environment (Guo et al., 2021). These animals need to cope with extreme climatic conditions such as cold, hypoxia, and intense ultraviolet radiation to ensure their survival and reproduction (Castiglione et al., 2017).Adaptations to high-altitude environments were analyzed by genome-wide scanning of natural selection traits in each species. Kirghiz sheep as one of Xinjiang’s meat sheep,In the 1970s,Kirghiz sheep were endangered until 1974 when China restored the Kirghiz Autonomous Prefecture sheep farm, allowing for the protection of this breed.Currently, there are approximately 640,000 Kirghiz sheep in the Kirghiz Autonomous Prefecture. Although breeding efforts for Kirghiz sheep have been ongoing, there have been relatively few studies on their genetic structure and genetic diversity, resulting in incomplete population genetic structures and significant degradation of breed characteristics. Through in-depth research and identification of the gene types adapted to the environment (Porcelli et al., 2015), new breeds that can adapt to extreme environment and have excellent production performance can be effectively screened and developed (Williams et al., 2007), and the genetic potential of individual animals can be more accurately assessed (Wearing and Scott, 2022).

Gene chip technology plays an important role in gene selection and mapping, genome breeding, breed identification,and parentage testing (Li Y. et al., 2014).The 50 K SNP chip for sheep can evenly cover the entire sheep genome and includes important trait loci.Wu (Han and Xiong, 2001) discovered candidate genes related to growth traits and high-altitude adaptability in Cashmere goats from two different regions using the Fst and PI methods.Wang (Guo et al., 2019)compared and analyzed Tibetan sheep, Altay sheep, Duolang sheep,Hu sheep,and Mongolian sheep using sheep genome re-sequencing and found genes related to hemoglobin levels and red blood cell counts in certain selection regions, such as the CYP17 gene on chromosome 22 and the DNAJB5 gene on chromosome 2 of Duolang sheep.Caiye Zhu (Wang et al., 2019)studied seven sheep breeds from different regions using the 50 K SNP chip technology and identified candidate genes associated with high-altitude adaptability through Fst and XP-EHH analysis, such as EPAS1,CRYAA, LONP1,NF1,DPP4,SOD1,PPARG, SOCS2 and detected 31 significant SNPs associated with tail type traits in indigenous Chinese sheep and identified BMP2 and PDGFD as candidate genes for tail type traits.Based on the identified candidate genes for tail type traits.Baazaouil (Caiye et al., 2023) conducted whole-genome sequencing of semi-arid sheep using this chip and identified candidate genes PROKR1 and BMP2 related to tail fat. Caiye Zhu (Baazaoui et al., 2021) use a total of 31 significant SNPs related to tail type traits were detected in Chinese native sheep, and BMP2 and PDGFD were identified as tail type trait candidate genes. Based on the above identified tail type trait candidate genes, BMP2 and PDGFD genes were selected to study the relationship between tail SNPs of Altay sheep and Tibetan sheep. (Wei et al., 2016) conducted genome selection on 140 individuals of indigenous Chinese sheep breeds and determined that PDGFD may affect fat deposition in fat-type sheep.

Population genetic diversity and selection signature are among the fundamental methods for studying sheep genomes and can reveal how genes respond to environmental selection. Therefore,the aim of this study was to use genome selection and haplotype statistics to study the population structure and genetic diversity of Kirghiz sheep, identify genetic markers associated with production traits,and provide theoretical basis for the conservation of sheep germplasm resources in southern Xinjiang,to promote animal husbandry sustainable development.

2 Materials and methods

2.1 Animal care

The study followed the guidelines of the Ethics Committee of Tarim University of Science and Technology (SYXK 2020-009) for animal experiments.

2.2 Experimental animals

In this study, 189 blood samples from adult Qira black sheep (QR) and 42 venous blood samples from adult Kirghiz sheep (KE) were randomly collected from Jinken Aoqun Agriculture and Animal Husbandry Technology Limited in Qira County, Hotan District and Bozheng Sheep Industry Technology Limited in Wuqia County, Kezhou. In addition, data from 160 Chinese merino sheep (CM) and 24 Tien-Shan in Kyrgyzstan (TNSH) were used.The data are from:https://doi.org/10.5061/dryad.37pvmcvff (Deniskova et al., 2019a).

2.3 DNA extraction and identification

The remaining samples of genomes DNA (gDNA) underwent agarose gel electrophoresis and nanodrop ND-2000 (Thermo Scientific) concentration analysis to ensure their quality and concentration.After adjusting the gDNA concentration to 50 ng/L,whole-genome amplification was performed. Subsequently,the gDNA was fragmented and precipitated before being re-suspended in hybridization buffer. The re-suspended DNA fragments were then applied to the chip for hybridization.Following hybridization,non-specifically bound DNA was removed, leaving behind specifically bound sites. These specific binding sites were then subjected to single-base extension and staining, followed by scanning using the Illumina iScan Reader. After completing the scan,the extracted DNA underwent quality control to obtain detailed quality control reports. Finally,the raw data scanned by the iScan system were imported into the Illumina official data analysis pipeline to generate PLINK files.

2.4 Genetic diversity and population structure

We utilized liftover (https://genome.ucsc.edu/cgi-bin/hgLiftOver) to convert the genome coordinates of 37 Tien-Shan sheep, then merged the files and performed quality control on the SNP data using PLINK software (version 1.90). We used the following criteria: inclusion detection rate less than 90% and hardy-Weinberg balance test p-value less than 10^–5 (Zhang et al., 2022).

2.5 Genetic diversity and selection signature

The PLINK (V1.90) software was used to perform principal component analysis based on the variance-standardized relationship matrix of the quality-controlled data. The NJ matrix was calculated using VCF2Dis (https://github.com/BGI-shenzhen/VCF2Dis) to draw the evolutionary tree, and the sparse non-negative matrix factorization algorithm sNMF (http://membres-timc.imag) was employed to generate estimates of ancestral proportions.

2.6 Fixation index

Fst is a statistical test used to measure the degree of differentiation between populations, mainly used to study the extent of genetic variation between different populations,as well as population structure and genetic diversity.The calculation formula is as follows:

$Fst=\frac{MSP-MSG}{MSP+\left(n_c-1\right)MSG}$

In this formula,MSG rep resents the mean square error of the detected intra-group sites, while MSP stands for the mean square variance of the inter-group sites, indicating the corrected average sample size between groups (Zhang W. et al., 2023). In this study,SNP sites retained were calculated using the unbiased estimation Fst method, yielding comparisons between Kirghiz sheep and Qira black sheep, as well as Kirghiz sheep and Chinese Merino sheep.

2.7 Haplotype analysis

We conducted haplotype statistics and inference on the PDGFD gene region genotype data of Kirghiz sheep and Tien-Shan sheep using the genehap R package in R scripts (Zhang R. et al., 2023).With this software, we combined continuous SNP data into individual haplotypes and generated corresponding haplotype information for each individual. We employed the genehap R software for phenotype-associated haplotype statistics in the vicinity of the PDGFD gene region.

2.8 Nucleotide diversity (PI)

Nucleotide diversity (PI) refers to the average number of different nucleotide at the same position of a random sequence taken from the DNA of multiple samples in a population, representing the degree of nucleotide polymorphisms within the population (Brouillette et al., 2000). The calculation formula is as follows:

$PI={\displaystyle \sum _{j=i}^S}{\mathrm{h}}_j$

In the above formula, S represents the number of segregating sites, and ${\mathrm{h}}_j$ represents the heterozygosity of the j segregating site. We utilized the vcftools software to compute population nucleotide diversity.

2.9 Enrichment analysis of candidate genes

We conducted statistical analysis of the p-values for Fst and PI calculated for the Kirghiz sheep (KE), Chinese Merino sheep (CM), and Qira black sheep (QR) populations,and plotted Manhattan distribution graphs. The top 5% of selected sites were referenced against two databases: the sheep genome Ovis Oar_v4.0 for annotation and the NCBI database (http://www.ncbi.nlm.nih.gov/gene) for enrichment analysis using cluster profile (Yu et al., 2012). The KEGG and GO enrichment analyses covered three aspects:Biological Process, Cellular Component, and Molecular Function.

3 Results

3.1 Population structure

Cluster analysis and principal component analysis (PCA) were conducted on SNP data from 415 sheep. The four populations were clustered into four groups, and significant differentiation among them was observed using PC1, PC2, and, PC3. TNSH showed slight differentiation compared to other populations, while QR exhibited closer proximity to TNSH in terms of principal components, consistent with the results of the evolutionary tree. The four sheep breeds QR, KE, CM, and TNSH formed separate branches, with KE being distantly separated from the other three sheep breeds.QR and TNSH populations showed closer structural proximity, with QR extending outward, indicating the separation of KE from the other three sheep breeds (Figure 1B). Considering the different evolutionary processes among the four sheep breeds, we further inferred the ancestral proportions of the four distinct sheep populations using the sNMF software (Figure 1C). When K = 5, each of the four sheep breeds exhibited distinct ancestral components, consistent with the separation observed in the PCA results for TNSH.The evolutionary tree results revealed a close genetic relationship between TNSH and QR, while KE and CM were situated on the same evolutionary branch (Figure 1A).

Figure 1

Figure 1. (A) NJ evolutionary tree. (B) Principal Component Analysis of Four Sheep Breeds (X-axis represents PC1, Y-axis represents PC2, Z-axis represents PC3. (C) sNMF Ancestral Proportion Estimation (minimized CV_error when K = 5).

3.2 Selective sweeps

3.2.1 Fixation index

Initially, we ranked the Fst values of TNSH in descending order and selected the top 5% as selection regions. After Fst filtering, we identified a total of 3048 genes. When examining the degree of differentiation between KE and TNSH populations using Fst, we discovered the gene PDGFD, which influences tail fat deposition. Fst tests for QR, CM, and KE revealed 2723,3264 candidate genes, respectively. (Supplementary Data 1).

3.2.2 Haplotype statistics

Subsequently,we generated tail trait phenotype-associated haplotypes in the PDGFD gene region,and the results showed that the H001 and H003 haplotypes were positively correlated with tail fat deposition in the PDGFD gene (Figure 2, Supplementary Data 2).

Figure 2

Figure 2. (A) Fst results of KE and TNSH, the red line in the figure represents the 0.05 threshold (B) Haplotype analysis of the region near PDGFD (chr15:3800010-4200000) based on the genehap R package. (C, D) Phenotypic correlation analyses of the fat tail of Kirghiz sheep with the thin tail of TNSH).

3.2.3 Nucleotide diversity

We sorted the PI values of KE in ascending order and selected the top 5% as the selection region, revealing a total of 3484 genes (Supplementary Data 3). Following the cross-analysis of genes filtered by both Fst and PI, we identified 639 intersecting candidate genes (Figure 3B, Supplementary Data 3).

Figure 3

Figure 3. (A) Plot of Fst and PI results for KE vs. QR and CM,with different colored dots representing whether the locus is under selection, using a significance threshold of 0.05. (B) Intersection of genes annotated from nucleotide diversity in KE and the top 5% Fst loci from KE. QR and CM comparisons).

Through Fst and PI analysis of KE, CM, and QR, we identified several genes associated with Kirghiz sheep. These genes are related to the high-altitude adaptation of organisms, such as ASIP, FMO1, FMO2, FMO4, as well as genes related to animal growth and reproduction, such as PTK6, PAG11, and CDC16.

3.2.4 Candidate gene enrichment analysis

The GO enrichment pathway analysis was performed on 590 candidate genes selected by QR and KE. There were 12 significant p-value<0.05 pathways, among which the following were related to high altitude adaptation:GO:0050732∼negative regulation of peptidyl-tyrosine phosphorylation,GO:0006805∼xenobiotic metabolic process, KEGG enrichment pathways showed 18 significant p-value<0.05 pathways, among which OAS00430-Taurine and hypotaurine metabolism were associated with high altitude adaptation.oas04610∼Complement and coagulation cascades; oas04672∼Intestinal immune network for IgA production (Supplementary Data 4, Figure 4).

Figure 4

Figure 4. significance threshold of the top twenty key genes for the interaction of enrichment pathway.

4 Discussion

4.1 Genetic diversity and population structure

Principal component analysis of Kirghiz sheep, Tien-shan sheep, Qira black and Merino sheep showed no significant aggregation in these four populations. However, Qira black and Tien-shan sheep exhibited closer proximity, suggesting a potential shared ancestry, possibly influenced by historical breeding practices.Notably, despite inhabiting the Naryn region of Kyrgyzstan (Deniskova et al., 2019b), Tien-shan sheep demonstrated unstable population structure, likely influenced by human-mediated selection over recent decades (Eydivandi et al., 2021).Conversely, Qira black sheep from Hotan cele County displayed greater ancestral resemblance to Tien-shan sheep, possibly attributable to trade dynamics in livestock (Valerio et al., 2020).

The phylogenetic tree reconstructed analysis indicated no significant kinship among these populations, corroborating the PCA findings and suggesting independent evolutionary trajectories in relatively isolated environments.

However, ancestral proportion analysis revealed a close genetic affinity between Kirghiz and Tien-shan sheep, as well as Qira black sheep. This likely stems from shared ancestral components and genetic distances, potentially shaped by prolonged cohabitation in similar environments and the introduction of foreign breeds, contributing to genetic homogeneity among these populations.

4.2 Adaptation mechanism in plateau environment

ASIP (agouti signaling protein) located on chromosome 13 of sheep, from physical distance 63063228 to 63063230 bp, ASIP, which is an autosomal gene, is directly related to a pathway that regulates melanin production (Goutte et al., 2022).ASIP plays a crucial role in regulating changes in the color of the animal’s back skin, It not only affects the muscle tone of the back, but also directly controls the production of melanin in the ventral skin area. It can regulate the process of protein formation in animal skin pigment, by precisely regulating the distribution and activity of melanocytes. Thus, it can be said that ASIP has a decisive influence in shaping the color pattern of the animal’s back and the absence or overexpression of melanin in the ventral skin (Li MH. et al., 2014). ASIP leads to the synthesis of eumelanin in hair follicle melanocytes, rather than the black or brown pigment true melanin, Gene switching events at ASIP sites in sheep may play an important role in the evolution of pigmentation in sheep (Norris and Whan, 2008). The pleiotropic effects expressed in animal models include obesity (Kempf et al., 2022),increased susceptibility to tumors, and premature infertility (Lima et al., 2015). Therefore,we can infer that the formation of ASIP is closely related to environmental factors. At higher altitudes, Ultraviolet radiation levels increase significantly, and sheep in these areas may adapt to this change by developing specific physiological mechanisms to reduce Ultraviolet radiation damage to them, which can help sheep better withstand Ultraviolet radiation stress in high-altitude environments. FMO1 (flavin containing dimethylaniline monoxygenase 1) is an enzyme containing flavin, located on chromosome 12 of sheep, with a physical position between 36894203 and 36894334 bp. Together with P450 enzymes, FMO1 participates in the reduction of TNO (N-oxidized dimethylaniline) to TAM (dimethylaniline), and the oxidation of TAM to TNO (Parte and Kupfer, 2005). FMO is the primary enzyme oxidant in this process,it binds to NADPH and exerts its catalytic activity (Eswaramoorthy et al., 2006). FMO can regulate cellular stress resistance through various cellular energy metabolism activities such as mitochondrial respiration pathway and glycolysis (Huang et al., 2021). When the levels of mean erythrocyte hemoglobin and mean erythrocyte hemoglobin concentration in the blood are elevated, it may indicate that sheep exhibit enhanced oxygen carrying capacity. FMO1 as a regulator of energy homeostasis (Veeravalli et al., 2014), is highly expressed in the liver and kidney of rats (Lattard et al., 2002) but can cause dysregulation of lipid metabolism (Zou et al., 2023), suggesting improved adaptation to high altitude hypoxia (Zhang Y. et al., 2021). This helps the Kirghiz adapt to the low oxygen and cold conditions at high altitudes. FMO2 (flavin containing dimethylaniline monoxygenase 2) gene, located on chromosome 12 of sheep, has a physical position between 36833449 and 36833580 bp.The human FMO2 gene has been confirmed to regulate oxidative stress levels (Henderson et al., 2008), thus playing a role in innate immunity against microbial infections, including tuberculosis (Mekonnen and Bekele, 2017). the FMO2 gene as an immune regulatory factor, may play a role in their resistance to tuberculosis. In addition,the transcription of this gene is more than 50% expressed in the lung (Siddens et al., 2008), but not in the liver or kidney (Yueh et al., 1997). The lung pressure, blood oxygen-carrying capacity and lipid metabolism of Kirghiz sheep living in high altitude have important effects on their adaptation to high altitude environment. These results suggest that FMO1 and FMO2 genes may play an important role in the cold anoxic environment and immune response mechanism under the plateau, and further studies may reveal their specific functions and mechanisms in Kirghiz sheep population.

4.3 Tail fat-related genes

Tail fat plays an important role in the adaptability of the Kirghiz sheep, not only being economically valuable (Zhao et al., 2022), but also storing energy and acting as a buffer to help the sheep survive extremely cold and arid environments (Wang et al., 2021). The deposition of fat in the tail can store energy and act as a buffer, thereby protecting the organism from the effects of extreme environments (Moradi et al., 2022). Additionally, fat is involved in regulating key physiological and biochemical responses in the organism (Zhang W. et al., 2021), crucial for the adaptation of Kirghiz sheep to high-altitude environments.Through selection signal analysis, we discovered candidate genes related to sheep tail fat, located in the CDS region of the PDGFD gene area, from chr15:3800010 to 4200000 bp. Three haplotypes were identified: H001 (chr15:3972342) with a frequency of 14, H002 with a frequency of 2, and H003 (chr15:4131705) with a frequency of 1 (Figure 2). The expression of the PDGFD region is related to the maturation of adipocytes (Luo et al., 2021). Members of the PDGF (prostaglandin-related factors) family have been shown to promote the proliferation of certain types of preadipocytes and effectively inhibit differentiation into mature adipocytes (Ma et al., 2018) (Pan et al., 2019), playing a key role in the process of body fat adipogenesis in humans and mice, and controlling the shape of the fat-tailed sheep to a certain extent (Dong et al., 2020). This gene is important in regulating fat deposition in the thin-tailed sheep (Li et al., 2020). We identified functional genes and selective gene regions associated with fat deposition in Kirghiz sheep breeds, specifically in the PDGFD region.

5 Conclusion

In conclusion,we employed the Illumina Ovine SNP50K Bead Chip for targeted signal analysis of Kirghiz sheep thriving in high-altitude environments.Our investigation unveiled candidate genes such as ASIP, FMO1, FMO2, implicated in the growth and development of Kirghiz sheep and their adaptation to elevated altitudes, along with PDGFD,a pivotal gene governing fat deposition in sheep tails. Our research outcomes are instrumental for pinpointing candidate genes associated with crucial traits across diverse sheep breeds, addressing the challenges posed by global climate change, cultivating novel sheep varieties endowed with robust cold tolerance,and offering insights into the development of new sheep breeds boasting desirable tail characteristics. This study establishes a solid groundwork for breeding superior new sheep germplasm.

Data availability statement

The data presented in the study are deposited in the Figshare repository, accession link: https://doi.org/10.6084/m9.figshare.25744491.v1.

Ethics statement

The animal study was approved by the Tarim University of Science and Technology (SYXK 2020-009) for animal experiments. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

XL: Conceptualization, Writing–original draft. LiZ: Conceptualization, Writing–review and editing. C-LZ: Writing–review and editing, Investigation. XW: Formal Analysis, Writing–review and editing. YL: Formal Analysis, Writing–review and editing. WZ: Formal Analysis, Writing–review and editing. ZH: Methodology, Writing–review and editing. RY: Formal Analysis, Writing–review and editing. YP: Formal Analysis, Writing–review and editing. YH: Formal Analysis, Writing–review and editing. LuZ: Formal Analysis, Writing–original draft, Writing–review and editing. SL: Conceptualization, Writing–review and editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Natural Science Foundation of China (NO: 32060743), Bingtuan Science and Technology Program (NO: 2022CB001-09), Autonomous Region Pastoral Area Indigenous Meat Sheep Breeding and Promotion Technology System (NO: MQTX-11), and 2022 Autonomous Region’s Scientific and Technological Action Plan for Rural Revitalization and Industrial Development (NO: 2022NC110).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2024.1432105/full#supplementary-material

References

Baazaoui, I., Bedhiaf-Romdhani, S., Mastrangelo, S., and Ciani, E. (2021). Genome-wide analyses reveal population structure and identify candidate genes associated with tail fatness in local sheep from a semi-arid area. Animal 15 (4), 100193. doi:10.1016/j.animal.2021.100193

PubMed Abstract | CrossRef Full Text | Google Scholar

Brouillette, J. A., Andrew, J. R., and Venta, P. J. (2000). Estimate of nucleotide diversity in dogs with a pool-and-sequence method. Mamm. Genome 11 (12), 1079–1086. doi:10.1007/s003350010220

PubMed Abstract | CrossRef Full Text | Google Scholar

Caiye, Z., Song, S., Li, M., Huang, X., Luo, Y., and Fang, S. (2023). Genome-wide DNA methylation analysis reveals different methylation patterns in Chinese indigenous sheep with different type of tail. Front. Vet. Sci. 10, 1125262. doi:10.3389/fvets.2023.1125262

PubMed Abstract | CrossRef Full Text | Google Scholar

Castiglione, G. M., Hauser, F. E., Liao, B. S., Lujan, N. K., Van Nynatten, A., Morrow, J. M., et al. (2017). Evolution of nonspectral rhodopsin function at high altitudes. Proc. Natl. Acad. Sci. U. S. A. 114 (28), 7385–7390. doi:10.1073/pnas.1705765114

PubMed Abstract | CrossRef Full Text | Google Scholar

Deniskova, T., Dotsev, A., Lushihina, E., Shakhin, A., Kunz, E., Medugorac, I., et al. (2019a). Population structure and genetic diversity of sheep breeds in the Kyrgyzstan. Front. Genet. 10, 1311. doi:10.3389/fgene.2019.01311

PubMed Abstract | CrossRef Full Text | Google Scholar

Deniskova, T., Dotsev, A., Lushihina, E., Shakhin, A., Kunz, E., Medugorac, I., et al. (2019b). Population structure and genetic diversity of sheep breeds in the Kyrgyzstan. Front. Genet. 10, 1311. doi:10.3389/fgene.2019.01311

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, K., Yang, M., Han, J., Ma, Q., Han, J., Song, Z., et al. (2020). Genomic analysis of worldwide sheep breeds reveals PDGFD as a major target of fat-tail selection in sheep. BMC Genomics 21 (1), 800. doi:10.1186/s12864-020-07210-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Eswaramoorthy, S., Bonanno, J. B., Burley, S. K., and Swaminathan, S. (2006). Mechanism of action of a flavin-containing monooxygenase. Proc. Natl. Acad. Sci. U. S. A. 103 (26), 9832–9837. doi:10.1073/pnas.0602398103

PubMed Abstract | CrossRef Full Text | Google Scholar

Eydivandi, S., Roudbar, M. A., Ardestani, S. S., Momen, M., and Sahana, G. (2021). A selection signatures study among Middle Eastern and European sheep breeds. J. Anim. Breed. Genet. 138 (5), 574–588. doi:10.1111/jbg.12536

PubMed Abstract | CrossRef Full Text | Google Scholar

Goutte, S., Hariyani, I., Utzinger, K. D., Bourgeois, Y., and Boissinot, S. (2022). Genomic analyses reveal association of ASIP with a Recurrently evolving adaptive color pattern in frogs. Mol. Biol. Evol. 39 (11), msac235. doi:10.1093/molbev/msac235

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, J., Zhong, J., Li, L., Zhong, T., Wang, L., Song, T., et al. (2019). Comparative genome analyses reveal the unique genetic composition and selection signals underlying the phenotypic characteristics of three Chinese domestic goat breeds. Genet. Sel. Evol. 51 (1), 70. doi:10.1186/s12711-019-0512-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, T., Zhao, H., Yuan, C., Huang, S., Zhou, S., Lu, Z., et al. (2021). Selective sweeps uncovering the genetic basis of horn and adaptability traits on fine-wool sheep in China. Front. Genet. 12, 604235. doi:10.3389/fgene.2021.604235

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, A. J., and Xiong, M. (2001). Gene chip technology and its advances in medical science. Zhongguo Yi Xue Ke Xue Yuan Xue Bao 23 (5), 528–531.

PubMed Abstract | Google Scholar

Henderson, M. C., Siddens, L. K., Morre, J. T., Krueger, S. K., and Williams, D. E. (2008). Metabolism of the anti-tuberculosis drug ethionamide by mouse and human FMO1, FMO2 and FMO3 and mouse and human lung microsomes. Toxicol. Appl. Pharmacol. 233 (3), 420–427. doi:10.1016/j.taap.2008.09.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, S., Howington, M. B., Dobry, C. J., Evans, C. R., and Leiser, S. F. (2021). Flavin-containing monooxygenases are conserved regulators of stress resistance and metabolism. Front. Cell Dev. Biol. 9, 630188. doi:10.3389/fcell.2021.630188

PubMed Abstract | CrossRef Full Text | Google Scholar

Kempf, E., Landgraf, K., Stein, R., Hanschkow, M., Hilbert, A., Abou Jamra, R., et al. (2022). Aberrant expression of agouti signaling protein (ASIP) as a cause of monogenic severe childhood obesity. Nat. Metab. 4 (12), 1697–1712. doi:10.1038/s42255-022-00703-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Lattard, V., Lachuer, J., Buronfosse, T., Garnier, F., and Benoit, E. (2002). Physiological factors affecting the expression of FMO1 and FMO3 in the rat liver and kidney. Biochem. Pharmacol. 63 (8), 1453–1464. doi:10.1016/s0006-2952(02)00886-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M. H., Tiirikka, T., and Kantanen, J. (2014b). A genome-wide scan study identifies a single nucleotide substitution in ASIP associated with white versus non-white coat-colour variation in sheep (Ovis aries). Hered. (Edinb). 112 (2), 122–131. doi:10.1038/hdy.2013.83

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Q., Lu, Z., Jin, M., Fei, X., Quan, K., Liu, Y., et al. (2020). Verification and analysis of sheep tail type-associated PDGF-D gene polymorphisms. Anim. (Basel) 10 (1), 89. doi:10.3390/ani10010089

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Wu, D. D., Boyko, A. R., Wang, G. D., Wu, S. F., Irwin, D. M., et al. (2014a). Population variation revealed high-altitude adaptation of Tibetan mastiffs. Mol. Biol. Evol. 31 (5), 1200–1205. doi:10.1093/molbev/msu070

PubMed Abstract | CrossRef Full Text | Google Scholar

Lima, F. A., Gonçalves, F. T., and Fridman, C. (2015). SLC24A5 and ASIP as phenotypic predictors in Brazilian population for forensic purposes. Leg. Med. (Tokyo) 17 (4), 261–266. doi:10.1016/j.legalmed.2015.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, R., Zhang, X., Wang, L., Zhang, L., and Zheng, Z. (2021). GLIS1,a potential candidate gene affect fat deposition in sheep tail. Mol. Biol. Rep. 48 (5), 4925–4931. doi:10.1007/s11033-021-06468-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, L., Li, Z., Cai, Y., Xu, H., Yang, R., and Lan, X. (2018). Genetic variants in fat-and short-tailed sheep from high-throughput RNA-sequencing data. Anim. Genet. 49 (5), 483–487. doi:10.1111/age.12699

PubMed Abstract | CrossRef Full Text | Google Scholar

Mekonnen, E., and Bekele, E. (2017). An ancestral human genetic variant linked to an ancient disease: a novel association of FMO2 polymorphisms with tuberculosis (TB) in Ethiopian populations provides new insight into the differential ethno-geographic distribution of FMO2*1. PLoS One 12 (10), e0184931. doi:10.1371/journal.pone.0184931

PubMed Abstract | CrossRef Full Text | Google Scholar

Moradi, M. H., Nejati-Javaremi, A., Moradi-Shahrbabak, M., Dodds, K. G., Brauning, R., and McEwan, J. C. (2022). Hitchhiking mapping of candidate regions associated with fat deposition in Iranian thin and fat tail sheep breeds suggests new insights into molecular aspects of fat tail selection. Anim. (Basel) 12 (11), 1423. doi:10.3390/ani12111423

CrossRef Full Text | Google Scholar

Norris, B. J., and Whan, V. A. (2008). A gene duplication affecting expression of the ovine ASIP gene is responsible for white and black sheep. Genome Res. 18 (8), 1282–1293. doi:10.1101/gr.072090.107

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, Z., Li, S., Liu, Q., Wang, Z., Zhou, Z., Di, R., et al. (2019). Rapid evolution of a retro-transposable hotspot of ovine genome underlies the alteration of BMP2 expression and development of fat tails. BMC Genomics 20 (1), 261. doi:10.1186/s12864-019-5620-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Parte, P., and Kupfer, D. (2005). Oxidation of tamoxifen by human flavin-containing monooxygenase (FMO) 1 and FMO3 to tamoxifen-N-oxide and its novel reduction back to tamoxifen by human cytochromes P450 and hemoglobin. Drug Metab. Dispos. 33 (10), 1446–1452. doi:10.1124/dmd.104.000802

PubMed Abstract | CrossRef Full Text | Google Scholar

Porcelli, D., Butlin, R. K., Gaston, K. J., Joly, D., and Snook, R. R. (2015). The environmental genomics of metazoan thermal adaptation. Hered. (Edinb). 114 (5), 502–514. doi:10.1038/hdy.2014.119

PubMed Abstract | CrossRef Full Text | Google Scholar

Qi, Y., Zhao, W., Li, Y., and Zhao, Y. (2021). Environmental and geological changes in the Tarim Basin promoted the phylogeographic formation of Phrynocephalus forsythii (Squamata:Agamidae). Gene 768, 145264. doi:10.1016/j.gene.2020.145264

PubMed Abstract | CrossRef Full Text | Google Scholar

Siddens, L. K., Henderson, M. C., Vandyke, J. E., Williams, D. E., and Krueger, S. K. (2008). Characterization of mouse flavin-containing monooxygenase transcript levels in lung and liver,and activity of expressed isoforms. Biochem. Pharmacol. 75 (2), 570–579. doi:10.1016/j.bcp.2007.09.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Valerio, V. C., Walther, O. J., Eilitta, M., Cissé, B., Muneepeerakul, R., and Kiker, G. A. (2020). Network analysis of regional livestock trade in West Africa. PLoS One 15 (5), e0232681. doi:10.1371/journal.pone.0232681

PubMed Abstract | CrossRef Full Text | Google Scholar

Veeravalli, S., Omar, B. A., Houseman, L., Hancock, M., Gonzalez Malagon, S. G., Scott, F., et al. (2014). The phenotype of a flavin-containing monooyxgenase knockout mouse implicates the drug-metabolizing enzyme FMO1 as a novel regulator of energy balance. Biochem. Pharmacol. 90 (1), 88–95. doi:10.1016/j.bcp.2014.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W., Zhang, X., Zhou, X., Zhang, Y., La, Y., Zhang, Y., et al. (2019). Deep genome resequencing reveals artificial and natural selection for visual Deterioration,Plateau adaptability and high prolificacy in Chinese domestic sheep. Front. Genet. 10, 300. doi:10.3389/fgene.2019.00300

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Fang, C., He, H., Cao, H., Liu, L., Jiang, L., et al. (2021). Identification of key genes in sheep fat tail evolution Based on RNA-seq. Gene 781, 145492. doi:10.1016/j.gene.2021.145492

PubMed Abstract | CrossRef Full Text | Google Scholar

Wearing, O. H., and Scott, G. R. (2022). Evolved reductions in body temperature and the metabolic costs of thermoregulation in deer mice native to high altitude. Proc. Biol. Sci. 289, 20221553. doi:10.1098/rspb.2022.1553

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, C., Wang, H., Liu, G., Zhao, F., Kijas, J. W., Ma, Y., et al. (2016). Genome-wide analysis reveals adaptation to high altitudes in Tibetan sheep. Sci. Rep. 6, 26770. doi:10.1038/srep26770

PubMed Abstract | CrossRef Full Text | Google Scholar

Williams, R. B., Chan, E. K., Cowley, M. J., and Little, P. F. R. (2007). The influence of genetic variation on gene expression. Genome Res. 17 (12), 1707–1716. doi:10.1101/gr.6981507

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, C., Ma, S., Zhao, B., Qin, C., Wu, Y., Di, J., et al. (2023). Drivers of plateau adaptability in cashmere goats revealed by genomic and transcriptomic analyses. BMC Genomics 24 (1), 428. doi:10.1186/s12864-023-09333-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, G., Wang, L. G., Han, Y., and He, Q. Y. (2012). clusterProfiler:an R package for comparing biological themes among gene clusters. OMICS 16 (5), 284–287. doi:10.1089/omi.2011.0118

PubMed Abstract | CrossRef Full Text | Google Scholar

Yueh, M. F., Krueger, S. K., and Williams, D. E. (1997). Pulmonary flavin-containing monooxygenase (FMO) in rhesus macaque:expression of FMO2 protein,mRNA and analysis of the cDNA. Biochim. Biophys. Acta 1350 (3), 267–271. doi:10.1016/s0167-4781(97)00004-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C. L., Liu, C., Zhang, J., Zheng, L., Chang, Q., Cui, Z., et al. (2022). Analysis on the desert adaptability of indigenous sheep in the southern edge of Taklimakan Desert. Sci. Rep. 12 (1), 12264. doi:10.1038/s41598-022-15986-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, R., Jia, G., and Diao, X. (2023b). geneHapR:an R package for gene haplotypic statistics and visualization. BMC Bioinforma. 24, 199. doi:10.1186/s12859-023-05318-9

CrossRef Full Text | Google Scholar

Zhang, W., Jin, M., Lu, Z., Li, T., Wang, H., Yuan, Z., et al. (2023a). Whole genome resequencing reveals selection signals related to wool color in sheep. Anim. 13 (20), 3265. doi:10.3390/ani13203265

CrossRef Full Text | Google Scholar

Zhang, W., Xu, M., Wang, J., Wang, S., Wang, X., Yang, J., et al. (2021b). Comparative transcriptome analysis of key genes and pathways activated in response to fat deposition in two sheep breeds with distinct tail phenotype. Front. Genet. 12, 639030. doi:10.3389/fgene.2021.639030

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Xue, X., Liu, Y., Abied, A., Ding, Y., Zhao, S., et al. (2021a). Genome-wide comparative analyses reveal selection signatures underlying adaptation and production in Tibetan and Poll Dorset sheep. Sci. Rep. 11 (1), 2466. doi:10.1038/s41598-021-81932-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Y., Zhang, D., Zhang, X., Li, F., Xu, D., Zhao, L., et al. (2022). Expression features of the ovine FTO gene and association between FTO polymorphism and tail fat deposition related-traits in Hu sheep. Gene 826, 146451. doi:10.1016/j.gene.2022.146451

PubMed Abstract | CrossRef Full Text | Google Scholar

Zou, L., Shi, Q., Li, Y., Yuan, Z., Peng, L., Lu, J., et al. (2023). FMO1 promotes nonalcoholic fatty liver disease progression by regulating PPARα activation and inducing ferroptosis. Discov. Med. 35 (177), 612–622. doi:10.24976/Discov.Med.202335177.60

PubMed Abstract | CrossRef Full Text | Google Scholar