Next Article in Journal
Identification and Characterization of SOG1 (Suppressor of Gamma Response 1) Homologues in Plants Using Data Mining Resources and Gene Expression Profiling
Next Article in Special Issue
The Effect of EGR1 on the Proliferation of Dermal Papilla Cells
Previous Article in Journal
MicroRNA-132-3p, Downregulated in Myeloid Angiogenic Cells from Hereditary Hemorrhagic Telangiectasia Patients, Is Enriched in the TGFβ and PI3K/AKT Signalling Pathways
Previous Article in Special Issue
MicroRNA-181a Regulates the Proliferation and Differentiation of Hu Sheep Skeletal Muscle Satellite Cells and Targets the YAP1 Gene
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 13
Issue 4
10.3390/genes13040666
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Genetic Polymorphisms of IGF1 and IGF1R Genes and Their Effects on Growth Traits in Hulun Buir Sheep

by
Ning Ding
1,2,3,4,5,†,
Dehong Tian
1,3,†,
Xue Li
1,2,3,
Zhichao Zhang
1,2,3,
Fei Tian
1,3,
Sijia Liu
1,3,
Buying Han
1,2,3,
Dehui Liu
1,2,3 and
Kai Zhao
1,3,*
1
Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining 810001, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Qinghai Provincial Key Laboratory of Animal Ecological Genomics, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining 810008, China
4
Hulun Buir State Farm, Hulun Buir 021000, China
5
Hulun Buir Ecological Industry Academy of Technology, Hulun Buir 021000, China
*
Author to whom correspondence should be addressed.
†
These authors contribute equally to this work.
Genes 2022, 13(4), 666; https://doi.org/10.3390/genes13040666
Submission received: 2 March 2022 / Revised: 25 March 2022 / Accepted: 7 April 2022 / Published: 9 April 2022
(This article belongs to the Special Issue Molecular Mechanism Analysis of Important Traits and Breeding of New Breeds of Sheep)
Browse Figures
Review Reports Versions Notes

Abstract

:
The identification of candidate genes and genetic variations associated with growth traits is important for sheep breeding. Insulin like growth factor 1 (IGF1) and insulin like growth factor 1 receptor (IGF1R) are well-accepted candidate genes that affect animal growth and development. The current study attempted to assess the association between IGF1 and IGF1R genetic polymorphisms and growth traits in Hulun Buir sheep. To achieve this goal, we first identified three and ten single nucleotide polymorphisms (SNPs) in exons of IGF1 and IGF1R in Hulun Buir sheep and then constructed six haplotypes of IGF1R based on linkage disequilibrium, respectively. Association studies were performed between SNPs and haplotypes of IGF1 and IGF1R with twelve growth traits in a population encompassing 229 Hulun Buir sheep using a general linear model. Our result indicated three SNPs in IGF1 were significantly associated with four growth traits (p < 0.05). In IGF1R, three SNPs and two haplotype blocks were significantly associated with twelve growth traits (p < 0.05). The combined haplotype H5H5 and H5H6 in IGF1R showed the strong association with 12 superior growth traits in Hulun Buir sheep (p < 0.05). In conclusion, we identified SNPs and haplotype combinations associated with the growth traits, which provided genetic resources for marker-assisted selection (MAS) in Hulun Buir sheep breeding.

1. Introduction

Growth traits are among the most important economic attributes in sheep breeding and are of great concern to breeding experts. Growth traits, including body weight, average daily gain and body size greatly influence meat productivity, which influences production and profitability in the mutton sheep industry [1]. Studies have revealed that many candidate genes are related to growth traits, among which IGF1 and IGF1R genes are well-accepted candidate genes that affect growth and production performance in livestock [2,3]. Insulin-like growth factor 1 (IGF1) is an endocrine growth factor involved in normal growth and development [4,5,6], fetal development and metabolism [7,8]. Insulin-like growth factor 1 receptor (IGF1R) is encoded by the IGF1R gene and is a receptor tyrosine kinase that mediates the actions of IGF1 [9,10].
Significant associations were identified between single nucleotide polymorphisms (SNPs) of the two genes and growth performance in diverse farm animals, including cattle [11,12,13], buffaloes [14], pigs [15,16] and goats [17,18,19,20]. In sheep, it has been reported that SNPs of IGF1 and IGF1R are related to meat production and growth [21,22,23]. Using the PCR restriction fragment length polymorphism (PCR-RFLP) method, Grochowska et al. found a highly significant effect of SNPs in the 5' untranslated (5’ UTR) region of IGF1 on carcass traits and meat compositions in local sheep breeds in Poland Merino sheep [24]. Negahdary et al. found a significant effect of the 5’ UTR region of the IGFI gene on birth weight, weaning weight, 6-month weight, average daily gain from birth to weaning and average daily gain from 6 to 9 months in Makooei sheep [25]. A mutation in intron 12 of the IGF1R gene was significantly associated with body weight and growth rate in Pomeranian Coarse wool ewes [26]. Later, associations were found between SNP in exon 3 of IGF1R and daily gain in the early developmental stage of Colored Polish Merino sheep [23]. The discovery of associations between genetic polymorphisms and growth traits provides useful information for the genetic improvement in sheep breeding. SNPs in the exon of genes are important because they may cause potentially functional variations, which lead to phenotypic changes in livestock. Since most identified SNPs in IGF1 and IGF1R were located in the 5’ flanking regions, we paid particular attention to genetic variations in the exons of the two genes.
Hulun Buir sheep are one of the representative indigenous sheep breeds in northern China, characterized by their high-grade meat quality and outstanding resistance to stress, such as cold and roughage. A lack of advanced breeding methods leads to poor growth performance compared to commercial breeds. To conduct the genetic improvement and breeding in Hulun Buir sheep, several works have been performed to identify genetic variations that were associated with economic traits in Hulun Buir sheep. It has been shown that the somatostatin receptor 1 (SSTR1) gene harbors two SNPs that were remarkably associated with growth traits of Hulun Buir sheep [27]. Based on Genome-wide association studies (GWAS), six SNP loci from 526,225 autosomal markers were greatly associated with carcass traits and chest girth [28]. Candidate genes and SNPs have been reported to be associated with fat deposition and fat metabolism [29,30]. However, no systematic investigations have been reported on the association between genetic polymorphism and the early growth traits in Hulun Buir sheep.
To improve the growth performance of Hulun Buir sheep, we investigated the genetic polymorphisms of IGF1 and IGF1R and their associations with twelve growth traits. By scanning exons of IGF1 and IGF1R, we identified thirteen SNPs in the IGF1 and IGF1R genes and two haplotype blocks involving six haplotypes in 229 Hulun Buir sheep. Among these SNPs and haplotype blocks, six SNPs and two haplotype blocks were remarkably associated with growth traits in Hulun Buir sheep. Notably, we identified two combined haplotypes that demonstrated a strong association with twelve greater phenotypic traits. Conclusively, our study provided useful information and laid the foundation for the genetic breeding of Hulun Buir sheep.

2. Materials and Methods

2.1. Animals and Data Collection

In total, 229 Hulun Buir lambs (male = 106, female = 123), which were born in March 2019 on Hulun Buir sheep farms (Hulun Buir city, Inner Mongolia, China), were investigated. The animals were grazed in identical conditions. The birth weight (BW), weaning weight adjusted at 4-month-old (WW) and body weight at 9-month-old (NBW) were recorded. Meanwhile, average daily gains (ADG) during birth to weaning, weaning to 9-month-old and birth to 9-month-old periods were calculated. Body height (BH), body length (BL) and chest girth (CG) were measured at weaning and at 9 months of age, respectively. Approximately 1 cm3 marginal ear tissues were collected and preserved in 95% ethanol. Genomic DNA of Hulun Buir sheep was extracted using a Tiangen DNA extraction kit (Tiangen Biotech Co., Ltd, Beijing, China) and stored at −20 °C for PCR amplification.

2.2. SNP Identification and Genotyping

Four and twenty-one pairs of primers were designed for all exons of IGF1 and IGF1R genes based on the published mRNA sequences (Gene ID: 443318, GenBank No. NM_001009774 (IGF1), Gene ID: 443515, GenBank No. XM_027957015 (IGF1R)), using Primer Premier 5.0 (Premier Biosoft, Palo Alto, Santa Clara, CA, USA), respectively. The primer information is listed in Supplementary materials Table S1. The PCR contained 100 ng template DNA, 10 pM of each primer, 3.5 μL 10 × PCR buffer, 2.5 mM dNTP, 1 U of Taq DNA polymerase (Takara Biotechnology Co., Ltd., Beijing, China) and double-distilled water (ddH2O), to make up a volume of 35 μL. PCR was performed in Thermocycler System (ABI 9700, Applied Biosystems, Waltham, MA, USA) with the following reaction procedure: predenaturation at 94 °C for 5 min, followed by 35 cycles at 94 °C for 30 s, 53–60 °C for 30 s and 72 °C for 40 s, with a final extension at 72 °C for 10 min. PCR products were separated by gel electrophoresis (1.5% agarose), purified using magnetic beads (Agencourt AMPure XP, Beckman Coulter, Krefeld, Germany) and sequenced in an Agilent 3730 sequencer (Agilent Technologies, Santa Clara, CA, USA). The sequencing results were aligned to published sheep IGF1 and IGF1R genes using Chromas 2.0 and SeqMan (DNASTAR software, version 7.1) to identify potential SNPs.

2.3. Population Genetics of IGF1 and IGF1R Genes

Genotypic and allelic frequencies were estimated with the direct counting method. Hardy–Weinberg equilibrium (HWE), observed heterozygosity (Ho), expected heterozygosity (He) and effective allele numbers (Ne) were analyzed according to the genotype frequencies of SNPs [31]. Cervus (version 3.0) was used to calculate the polymorphic information content (PIC) of each mutation site [32].

2.4. Linkage Disequilibrium Analysis and Haplotype Construction

The extent of linkage disequilibrium (LD) between each pair of SNPs in IGF1 and IGF1R was analyzed according to the value of r2 using Haploview software (version 4.2) [33]. Haplotype blocks with strong LD of SNPs (r2 > 0.33) were defined based on the confidence intervals methods [34].

2.5. Statistical Analyses

SAS software (version 13.0, SAS Institute) was applied for statistical analyses, and the results were expressed as the mean ± SE (standard error). The associations were carried out between the genotypes and individual growth traits using general linear model (GLM):
Yij = μ + Gi + Sj + Gi × Sj + ε
where Yij is a growth trait measured on an individual animal (BW, WW, NBW, ADG, BH, BL and CG); µ is the mean value of overall; Gi is the fixed effect of genotypes of the population (i = 3 levels, except rs600896367 of IGF1 gene and c.244C>T, rs162159917, rs601806812 and rs193644211 of IGF1R gene where i = 2 levels); Sj is the fixed effect of sex (j = 2 levels); Gi × Sj is the interaction effect between sex and genotypes; if the difference of interaction effect between the sex and genotypes is not significant, the general linear model should be reduced as:
Yi = μ + Gi + ε
The association analysis between haplotype combinations and individual growth traits was analyzed by the following GLM:
Yij = μ + Hi + Sj + Hi × Sj + ε
where Yij is a growth trait measured on an individual animal (BW, WW, NBW, ADG, BH, BL and CG); µ is the mean value of overall; Hi is the fixed effect of haplotype combinations of the population (i = 5 levels); Sj is the fixed effect of sex (j = 2 levels); Hi × Sj is the interaction effect between sex and haplotype combinations; if the difference of interaction effects between the sex and haplotype combinations is not significant, the general linear model should be reduced as:
Yi = μ + Hi + ε
ε is the random error in the above models. Tukey’s test and Bonferroni corrections were performed for multiple pairwise comparisons between genotypes or haplotype combinations based on SNPs. The p value of 0.05 was defined as statistical significance.

3. Results

3.1. SNP Detection of IGF1 and IGF1R Genes in Hulun Buir Sheep

We detected three SNPs in the IGF1 gene and ten SNPs in the IGF1R gene in 229 Hulun Buir sheep (Figure 1, Figure 2, Table 1). All of the detected SNPs were transition mutations except for SNP13 (transversion mutation) in exon 19 of IGF1R. SNP4 in IGF1R was a nonsynonymous mutation resulting in a substitution of Cys for Arg in amino acid sequence, and the rest were synonymous mutations. By searching in the dbSNP database of NCBI, we found that SNP4 and SNP8 in IGF1R were two novel single-nucleotide mutations in sheep and will be uploaded to the SNP data bank (Table 1).

3.2. Population Genetic Analyses

3.2.1. Genotyping, Genotypic and Allelic Frequencies

Among all the SNPs, the wild types were dominant alleles compared with the mutants (Table 2). Genotyping results showed that SNP1 in IGF1 as well as SNP4, SNP5, SNP9 and SNP12 in IGF1R displayed two genotypes: wild-type homozygotes and mutant heterozygotes, and the remaining eight SNPs showed three different genotypes: wild-type homozygotes, mutant heterozygotes and mutant homozygotes (Table 2). In SNP6–8, heterozygotes showed the highest genotype frequencies compared with wild-type and mutant homozygous. In the remaining 10 SNPs, the wild-type homozygotes had the highest genotype frequencies compared with mutant heterozygotes and homozygotes (Table 2).

3.2.2. Genetic Diversity and Hardy–Weinberg Equilibrium

The allelic frequencies of all 13 SNPs obey the HWE law (p > 0.05). The Ne values of SNP2 in IGF1 and SNP6–SNP8 in IGF1R were close to 2. The PIC value showed that the five SNP loci (SNP2, SNP3, SNP6–SNP8) exhibited low polymorphism (PIC < 0.25), while the remaining eight SNPs showed moderate polymorphism in the Hulun Buir sheep population (0.25 < PIC < 0.5) (Table 2).

3.3. Effects of Genotypes on Growth Traits

Association analysis was performed between genotypes of the SNPs and growth traits on 229 Hulun Buir sheep. The statistical results were listed in Supplementary Tables S2–S5.

3.3.1. Effects of SNP Genotypes in IGF1 on Growth Traits

The GA genotype of SNP1 had significantly greater WCG and NBL than the GG genotype (p < 0.05). At the SNP2 locus, the higher NCG was observed in TC genotype than that in the CC genotype but not in the TT genotype (p < 0.05). The GG and GA genotypes of SNP3 were significantly associated with greater 4–9 ADG than the AA genotype (p < 0.05, Figure 3).

3.3.2. Effects of SNP Genotypes in IGF1R on Growth Traits

The mutant homozygotes (CC) of SNP6 had significantly longer NBL than those individuals with the TC genotype (p < 0.05, Figure 4A). Significant differences (p < 0.05) and extremely significant differences (p < 0.01) were found between genotypes of the SNP8 locus with the 11 growth traits out of 4–9 ADG (Figure 4B,C). The genotypes containing the wild–type allele had better phenotypic values than mutant homozygotes. At the SNP13 locus, the individuals with the CC genotype had greater NBW, 0–9 ADG, WCG, NBH and NCG than those with the GG genotype (p < 0.05); the CC and CG genotypes were associated with significantly longer NBL than the GG genotype (p < 0.05, Figure 4D,E). No significant effects were detected among the remaining seven SNP loci and early growth traits of Hulun Buir sheep (p > 0.05).

3.4. Linkage Disequilibrium and Haplotype Analysis

A strong linkage disequilibrium (r2 > 0.33) was observed among SNP5, SNP9 and SNP11, and between SNP6 and SNP7, as well as SNP8 and SNP9 loci in the IGF1R gene (Figure 5). In particular, SNP6 to SNP9 loci formed two haplotype blocks. The first haplotype block was composed of SNPs 6 and 7, including three common haplotypes. The haplotypes H1 (TC), H2 (CT) and H3 (CC) occurred at frequencies of 0.537, 0.389 and 0.074, respectively, and five haplotype combinations were generated (Table 3). The second haplotype block was composed of SNP8 and SNP9, including three common haplotypes. The haplotypes H4 (CG), H5 (TG) and H6 (CA) occurred at frequencies of 0.321, 0.581 and 0.098, respectively, and generated five haplotype combinations (Table 4). We did not detect the linkage disequilibrium among three SNP loci (r2 < 0.33) in the IGF1 gene (Figure 6).

3.5. Effects of Haplotype Combinations on Growth Traits

Association analysis was performed between haplotypes in the IGF1R gene and growth traits of 229 Hulun Buir sheep populations. The statistical results were shown in Supplementary Tables S6–S9. The haplotype block 1 was only significantly associated with NBL, in which H1H3 (TCCC) haplotype combination had significantly longer NBL than those individuals with the H2H3 (CTCC) haplotype combination (p < 0.05) (Figure 7A). For haplotype block 2, the sheep with H5H6 (TGCA) haplotype combination was significantly heavier than that of the H4H4 (CGCG) haplotype combination of BW (p < 0.05). The individuals with the H5H5 (TGTG) and H5H6 (TGCA) haplotype combinations had significantly greater WW, NBW, 0–4 ADG, 4–9 ADG, 0–9 ADG, WBL, WCG, NBH, NBL and NCG than those with the H4H6 (CGCA) haplotype combination (p < 0.05). H5H5 (TGTG) and H5H6 (TGCA) with the wild-type allele T were the predominant haplotype combinations in the experimental population (Figure 7B,C). Therefore, haplotype combinations H1H3 (TCCC), H5H5 (TGTG) and H5H6 (TGCA) can be used as candidate markers for better growth traits of Hulun Buir sheep.

4. Discussion

The growth of the animal was subject to growth hormone (GH)-IGF1 somatrotropic axis, in which GH acts as a major regulator for development, growth and anabolic processes. IGF1 modulates the biological actions of GH by binding to its receptor (IGF1R) [35]. IGF system includes IGF ligands and their receptors, which influences glycogenesis, glucogenesis and protein synthesis through the regulation of downstream gene expression and signaling pathways [36]. Among IGF ligands and receptors, IGF1 and IGF1R proteins are crucial regulators of cell growth and metabolism [37,38]. Genetic variation may have an impact on the phenotypic characteristics of animals by influencing the expression and function of the genes [39,40]. Therefore, we inferred that the genetic variation in IGF1 and IGF1R may also influence the growth traits of sheep.
In the present study, we discovered genetic polymorphisms of the IGF1 and IGF1R genes and evaluated their effects on growth traits in Hulun Buir sheep. Our results indicated that IGF1 and IGF1R exhibited low to medium genetic diversity, and some of the genetic variations exhibited a significant association with the growth performance in Hulun Buir sheep. This observation provided SNP marker information, which has potential feasibility for MAS in Hulun Buir sheep breeding schemes.
The Hardy–Weinberg equilibrium of all 13 SNPs indicated the absence of artificial selection of Hulun Buir sheep [41]. In the current study, two novel single-nucleotide polymorphisms were identified, including a nonsynonymous mutation of SNP4. A growing body of evidence has shown that the synonymous mutations could influence phenotypic performance by influencing gene expression through the regulation of mRNA stability and protein expression [42,43,44,45]. Maria et al. reported that the synonymous mutation rs159876393 SNP1 of IGF1 was associated with milk protein and casein contents in Sarda sheep [46]. A synonymous mutation SNP2 (rs159876393) in exon 2 of IGF1 was associated with variations in carcass traits of New Zealand Romney Sheep, including carcass weight, backfat thickness and the lean meat percentage [47]. Consistent with previous reports on other sheep breeds, we also identified a strong association of SNP1 and SNP2 with growth traits in Hulun Buir sheep, which indicated that SNP1 and SNP2 of the IGF1 gene might be related to multiple traits in sheep. A remarkable association was found between SNP3 (rs400398060) of the IGF1 gene and average daily gain from 4–9 months of age (4–9 ADG) in the present study. This mutant site was also detected in Egyptian Barki sheep and was not correlated with growth traits, indicating that its association might be dependent on the genetic backgrounds of sheep breeds [48]. Few studies reported the association between genetic polymorphisms of the IGF1R gene and growth traits in sheep. A significant correlation was detected between average daily gain and an SNP of the IGF1R gene in local sheep breed in Poland Merino sheep [23]. The present study reported 10 SNPs in the IGF1R gene, and SNP6, SNP8 and SNP13 were significantly associated with growth traits in Hulun Buir sheep. In addition, the sheep with homozygous wild genotype TT of SNP8 and CC of SNP13 had superior growth traits than those with homozygous mutant genotypes CC and GG, suggesting that they could serve as the predominant genotypes.
Generally, linked SNP loci are of much concern because of the existence of substantial LD between causal SNPs [49]. Haplotype combinations involving multiple linked SNP loci may provide more precise information than single SNP markers for association analysis [50,51,52]. In this study, the strong LD suggested that these alleles were tightly linked; thus, we carried out an association analysis between the haplotypes and growth traits. The association and multiple comparison analyses demonstrated that the H5H5 (TGTG) haplotype combination with wild-type alleles was the dominant haplotype. This was consistent with the result that the wild-type allele T of SNP8 was related to better growth traits. Additionally, SNP6-SNP9 formed two haplotype blocks, which displayed a remarkably significant effect on growth traits. Based on the results above, we inferred that the four SNPs did not act independently [53], and SNP6 and SNP8 of IGF1R may be causal mutations that affect phenotypic traits [54].

5. Conclusions

Conclusively, our analysis showed that SNP1, SNP2 and SNP3 of the IGF1 gene, SNP6, SNP8 and SNP13, as well as haplotype block 1 and haplotype block 2 of the IGF1R gene can be used as candidate markers for early growth traits in MAS of Hulun Buir sheep. Further studies will be conducted to investigate the effects of these SNPs on other economic traits in Hulun Buir sheep. The wild-type alleles of SNP8, haplotype combinations H5H5 (TGTG) and H5H6 (TGCA) in the IGF1R gene showed superior growth traits during the early stage. Overall, our study provided important genetic variations, which could serve as potential markers for growth trait selection in Hulun Buir sheep.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes13040666/s1, Table S1: Primer information of IGF1 and IGF1R in Hulun Buir sheep; Table S2: Associations for the SNPs of IGF1 gene with body weight traits and ADG traits in Hulun Buir sheep; Table S3: Associations for the SNPs of IGF1 gene with body size traits in Hulun Buir sheep; Table S4: Associations for the SNPs of IGF1R gene with body weight traits and ADG traits in Hulun Buir sheep; Table S5: Associations for the SNPs of IGF1R gene with body size traits in Hulun Buir sheep; Table S6: Associations for the haplotype combinations (block 1) of IGF1R gene with body weight traits and ADG traits in Hulun Buir sheep; Table S7: Associations for the haplotype combinations (block 1) of IGF1R gene with body size traits in Hulun Buir sheep; Table S8: Associations for the haplotype combinations (block 2) of IGF1R gene with body weight traits and ADG traits in Hulun Buir sheep; Table S9: Associations for the haplotype combinations (block 2) of IGF1R gene with body size traits in Hulun Buir sheep.

Author Contributions

N.D. and D.T. contributed equally to this paper. Conceptualization, K.Z. and F.T.; methodology, N.D., D.T. and F.T.; software, S.L. and Z.Z.; formal analysis, N.D., D.T. and Z.Z.; investigation, N.D., X.L., B.H. and D.L.; resources, N.D. and D.T.; data curation, N.D., D.T. and K.Z.; writing—original draft preparation, N.D. and D.T.; writing—review and editing, F.T.; visualization, N.D., X.L., B.H. and D.L.; supervision, K.Z. and D.T.; project administration, D.T.; funding acquisition, K.Z. and D.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Strategic Priority Research Program of Chinese Academy of Sciences, Grant No. XDA (26040303-01); “Wang Kuancheng Leading Talents Program for Industrial and Research Talents support Project” of Chinese Academy of Sciences; Project of Ecological Grassland Animal Husbandry Engineering Laboratory, Chinese Academy of Sciences (KFJ-PTXM-007).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of Northwest Institute of Plateau Biology, Chinese Academy of Sciences (NWIPB2021311).

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank Dejun Fu and Xiaoliang Song for their kind help in data collection of growth traits in Hulun Buir sheep.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Luo, W.W.; Zhou, Y.; Wang, J.R.; Yu, X.M.; Tong, J.G. Identifying Candidate Genes Involved in the Regulation of Early Growth Using Full-Length Transcriptome and RNA-Seq Analyses of Frontal and Parietal Bones and Vertebral Bones in Bighead Carp (Hypophthalmichthys nobilis). Front. Genet. 2021, 11, 603454. [Google Scholar] [CrossRef] [PubMed]
  2. XU, S.S.; Li, M.H. Recent advances in understanding genetic variants associated with economically important traits in sheep (Ovis aries) revealed by high-throughput screening technologies. Front. Agric. Sci. Eng. 2017, 4, 25–34. [Google Scholar] [CrossRef] [Green Version]
  3. Valencia, C.P.L.; Franco, L.; Herrera, D.H. Association of single nucleotide polymorphisms in the CAPN, CAST, LEP, GH, and IGF-1 genes with growth parameters and ultrasound characteristics of the Longissimus dorsi muscle in Colombian hair sheep. Trop. Anim. Health Prod. 2022, 54, 82. [Google Scholar] [CrossRef] [PubMed]
  4. Kindler, J.M.; Pollock, N.K.; Laing, E.M.; Jenkins, N.T.; Oshri, A.; Isales, C.; Hamrick, M.; Lewis, R.D. Insulin Resistance Negatively Influences the Muscle-Dependent IGF-1-Bone Mass Relationship in Premenarcheal Girls. J. Clin. Endocrinol. Metab. 2016, 101, 199–205. [Google Scholar] [CrossRef] [PubMed]
  5. Cannata, D.; Lann, D.; Wu, Y.; Elis, S.; Sun, H.; Yakar, S.; Lazzarino, D.A.; Wood, T.L.; Leroith, D. Elevated circulating IGF-I promotes mammary gland development and proliferation. Endocrinology 2010, 151, 5751–5761. [Google Scholar] [CrossRef]
  6. Estany, J.; Tor, M.; Villalba, D.; Bosch, L.; Gallardo, D.; Jiménez, N.; Altet, L.; Noguera, J.L.; Reixach, J.; Amills, M.; et al. Association of CA repeat polymorphism at intron 1 of insulin-like growth factor (IGF-I) gene with circulating IGF-I concentration, growth, and fatness in swine. Physiol. Genom. 2007, 31, 236–243. [Google Scholar] [CrossRef]
  7. Hu, F.; Liu, F. Targeting tissue-specific metabolic signaling pathways in aging: The promise and limitations. Protein Cell 2014, 5, 21–35. [Google Scholar] [CrossRef]
  8. Shen, W.; Wisniowski, P.; Ahmed, L.; Boyle, D.W.; Denne, S.C.; Liechty, E.A. Protein anabolic effects of insulin and IGF-I in the ovine fetus. Am. J. Physiol. Endocrinol. Metab. 2003, 284, E748–E756. [Google Scholar] [CrossRef] [Green Version]
  9. Alves, A.P.N.R.; Fernandes, J.C.; Fenerich, B.A.; Coelho-Silva, J.L.; Scheucher, P.S.; Simoes, B.P.; Rego, E.M.; Ridley, A.J.; Machado-Neto, J.A.; Traina, F. IGF1R/IRS1 targeting has cytotoxic activity and inhibits PI3K/AKT/mTOR and MAPK signaling in acute lymphoblastic leukemia cells. Cancer Lett. 2019, 456, 59–68. [Google Scholar] [CrossRef]
  10. Zhang, J.Y.; Liu, M.Q.; Huang, M.H.; Chen, M.F.; Zhang, D.; Luo, L.P.; Ye, G.N.; Deng, L.J.; Peng, Y.H.; Wu, X.; et al. Ginsenoside F1 promotes angiogenesis by activating the IGF-1/IGF1R pathway. Pharm. Res. 2019, 144, 292–305. [Google Scholar] [CrossRef]
  11. Reyna, X.F.; Montoya, H.M.; Castrellón, V.V.; Rincón, A.M.; Bracamonte, M.P.; Vera, W.A. Polymorphisms in the IGF1 gene and their effect on growth traits in Mexican beef cattle. Genet. Mol. Res. 2010, 9, 875–883. [Google Scholar] [CrossRef] [PubMed]
  12. Szewczuk, M.; Zyh, S.; Wójcik, J.; Czerniawska-Piątkowska, E. Association of two SNPs in the coding region of the insulin-like growth factor 1 receptor (IGF1R) gene with growth-related traits in Angus cattle. J. Appl. Genet. 2013, 54, 305–308. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Putra, D.E. Polymorphism of Insulin-like Growth Factor 1 Gene (IGF1/TasI, IGF1/SnaBI, IGF1/RsaI) and the Association with Daily Gain of Pesisir Cattle Local Breed from West Sumatera, Indonesia. Pak. J. Biol. Sci. 2017, 20, 210–216. [Google Scholar] [CrossRef]
  14. El-Magd, M.A.; Abbas, H.E.; El-kattawy, A.M.; Mokhbatly, A. Novel polymorphisms of the IGF1R gene and their association with average daily gain in Egyptian buffalo (Bubalus bubalis). Domest. Anim. Endocrinol. 2013, 45, 105–110. [Google Scholar] [CrossRef]
  15. Wang, B.; Li, P.; Zhou, W.; Gao, C.; Liu, H.; Li, H.; Niu, P.; Zhang, Z.; Li, Q.; Zhou, J.; et al. Association of Twelve Candidate Gene Polymorphisms with the Intramuscular Fat Content and Average Backfat Thickness of Chinese Suhuai Pigs. Animals 2019, 9, 858. [Google Scholar] [CrossRef] [Green Version]
  16. Yue, M.; Tian, Y.G.; Wang, Y.J.; Gu, Y.; Bayaer, N.; Hu, Q.; Gu, W.W. Associated analysis of single nucleotide polymorphisms found on exon 3 of the IGF-1 gene with Tibetan miniature pig growth traits. Genet. Mol. Res. 2014, 13, 1263–1269. [Google Scholar] [CrossRef]
  17. Luo, J.; Qin, F.; Deng, C.; Li, F.; Li, W.; Yue, X. Polymorphisms of IGF-IR gene and their association with economic traits in two indigenous Chinese dairy goat breeds. Gene 2019, 695, 51–56. [Google Scholar] [CrossRef]
  18. Wang, W.J.; Hui, K.; Su, X.F.; Xu, M.S.; Chen, X.; Guan, S. Polymorphism of Insulin-like Growth Factor 1 Receptor Gene in 12 Pig Breeds and Its Relationship with Pig Performance Traits. Asian-Australas. J Anim. Sci. 2006, 19, 1541–1545. [Google Scholar] [CrossRef]
  19. Naicy, T.; Venkatachalapathy, R.T.; Aravindakshan, T.V.; Kurian, E. Association of a Cac8I polymorphism in the IGF1 gene with growth traits in Indian goats. J. Genet. Eng. Biotechnol. 2017, 15, 7–11. [Google Scholar] [CrossRef]
  20. Lestari, D.A.; Oikawa, T.; Sutopo, S.; Purbowati, E.; Setiaji, A.; Kurnianto, E. Effect of insulin-like growth factor 1 gene on growth traits of Kejobong goat and its growth analysis. Vet. World 2020, 13, 127–133. [Google Scholar] [CrossRef] [Green Version]
  21. Meira, A.N.; Montenegro, H.; Coutinho, L.L.; Mourão, G.B.; Azevedo, H.C.; Muniz, E.N.; Machado, A.L.; Sousa, L.P., Jr.; Pedrosa, V.B.; Pinto, L.F.B. Single nucleotide polymorphisms in the growth hormone and IGF type-1 (IGF1) genes associated with carcass traits in Santa Ines sheep. Animal 2019, 13, 460–468. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. He, J.N.; Zhang, B.Y.; Chu, M.X.; Wang, P.Q.; Feng, T.; Cao, G.L.; Di, R.; Fang, L.; Huang, D.W.; Tang, Q.Q.; et al. Polymorphism of insulin-like growth factor 1 gene and its association with litter size in Small Tail Han sheep. Mol. Biol. Rep. 2012, 39, 9801–9807. [Google Scholar] [CrossRef] [PubMed]
  23. Grochowska, E.; Lisiak, D.; Akram, M.Z.; Adeniyi, O.O.; Lühken, G.; Borys, B. Association of a polymorphism in exon 3 of the IGF1R gene with growth, body size, slaughter and meat quality traits in Colored Polish Merino sheep. Meat. Sci. 2021, 172, 108314. [Google Scholar] [CrossRef] [PubMed]
  24. Grochowska, E.; Borys, B.; Janiszewski, P.; Knapik, J.; Mroczkowski, S.J.A.A.B. Effect of the IGF-I gene polymorphism on growth, body size, carcass and meat quality traits in Coloured Polish Merino sheep. Arch. Anim. Breed. 2017, 60, 161–173. [Google Scholar] [CrossRef] [Green Version]
  25. Negahdary, M.; Hajihosseinlo, A.; Ajdary, M. PCR-SSCP Variation of IGF1 and PIT1 Genes and Their Association with Estimated Breeding Values of Growth Traits in Makooei Sheep. Genet. Res. Int. 2013, 2013, 272346. [Google Scholar] [CrossRef]
  26. Proskura, W.S.; Szewczuk, M.J.P.V.J. The Polymorphism in the IGF1R Gene is Associated with Body Weight and Average Daily Weight Gain in Pomeranian Coarsewool Ewes. Pak. Vet. J. 2014, 34, 514–517. [Google Scholar] [CrossRef]
  27. Li, X.; Ding, N.; Zhang, Z.Z.; Tian, D.H.; Han, B.Y.; Liu, S.J.; Liu, D.H.; Tian, F.; Zhao, K. Identification of Somatostatin Receptor Subtype 1 (SSTR1) Gene Polymorphism and Their Association with Growth Traits in Hulun Buir Sheep. Genes 2021, 13, 77. [Google Scholar] [CrossRef]
  28. Zhang, T.Y.; Gao, H.D.; Sahana, G.; Zan, Y.Z.; Fan, H.Y.; Liu, J.X.; Shi, L.Y.; Wang, H.W.; Du, L.X.; Wang, L.X.; et al. Genome-wide association studies revealed candidate genes for tail fat deposition and body size in the Hulun Buir sheep. J. Anim. Breed. Genet. 2019, 136, 362–370. [Google Scholar] [CrossRef]
  29. Zhi, D.F.; Da, L.; Liu, M.N.; Cheng, C.; Zhang, Y.K.; Wang, X.; Li, X.N.; Tian, Z.P.; Yang, Y.Y.; He, T.Y.; et al. Whole Genome Sequencing of Hulunbuir Short-Tailed Sheep for Identifying Candidate Genes Related to the Short-Tail Phenotype. G3 Genes Genomes Genet. 2018, 8, 377–383. [Google Scholar] [CrossRef] [Green Version]
  30. Fan, H.Y.; Hou, Y.L.; Sahana, G.; Gao, H.D.; Zhu, C.Y.; Du, L.X.; Zhao, F.P.; Wang, L.X. A Transcriptomic Study of the Tail Fat Deposition in Two Types of Hulun Buir Sheep According to Tail Size and Sex. Animals 2019, 9, 655. [Google Scholar] [CrossRef] [Green Version]
  31. Nei, M.; Roychoudhury, A.K. Sampling Variances of Heterozygosity and Genetic Distance. Genetics 1974, 76, 379–390. [Google Scholar] [CrossRef] [PubMed]
  32. Kalinowski, S.T.; Taper, M.L.; Marshall, T.C. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol. Ecol. 2007, 16, 1099–1106. [Google Scholar] [CrossRef] [PubMed]
  33. Barrett, J.C.; Fry, B.; Maller, J.; Daly, M.J. Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics 2005, 21, 263–265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Gabriel, S.B.; Schaffner, S.F.; Nguyen, H.; Moore, J.M.; Roy, J.; Blumenstiel, B.; Higgins, J.; DeFelice, M.; Lochner, A.; Faggart, M.; et al. The structure of haplotype blocks in the human genome. Science. 2002, 296, 2225–2229. [Google Scholar] [CrossRef] [Green Version]
  35. Renaville, R.; Hammadi, M.; Portetelle, D. Role of the somatotropic axis in the mammalian metabolism. Domest. Anim. Endocrinol. 2002, 23, 351–360. [Google Scholar] [CrossRef]
  36. Wrigley, S.; Arafa, D.; Tropea, D. Insulin-like growth factor 1: At the crossroads of brain development and aging. Front. Cell Neurosci. 2017, 11, 14. [Google Scholar] [CrossRef] [Green Version]
  37. Byun, S.O.; Zhou, H.; Hickford, J.G. Polymorphism of the ovine insulin-like growth factor I receptor (IGFIR) gene. Mol. Cell. Probes 2008, 22, 131–132. [Google Scholar] [CrossRef]
  38. Kurokawa, M.; Sato, F.; Aramaki, S.; Soh, T.; Yamauchi, N.; Hattori, M.A. Monitor of the myostatin autocrine action during differentiation of embryonic chicken myoblasts into myotubes: Effect of IGF-I. Mol. Cell. Biochem. 2009, 331, 193–199. [Google Scholar] [CrossRef]
  39. Wyszyńska-Koko, J.; Pierzchała, M.; Flisikowski, K.; Kamyczek, M.; Rózycki, M.; Kurył, J. Polymorphisms in coding and regulatory regions of the porcine MYF6 and MYOG genes and expression of the MYF6 gene in m. longissimus dorsi versus productive traits in pigs. J. Appl. Genet. 2006, 47, 131–138. [Google Scholar] [CrossRef]
  40. Bonafè, M.; Barbieri, M.; Marchegiani, F.; Olivieri, F.; Ragno, E.; Giampieri, C.; Mugianesi, E.; Centurelli, M.; Franceschi, C.; Paolisso, G. Polymorphic variants of insulin-like growth factor I (IGF-I) receptor and phosphoinositide 3-kinase genes affect IGF-I plasma levels and human longevity: Cues for an evolutionarily conserved mechanism of life span control. J. Clin. Endocrinol. Metab. 2003, 88, 3299–3304. [Google Scholar] [CrossRef]
  41. Liu, X.Y.; Lu, R.; Xia, Y.L.; Sun, J. Global analysis of the eukaryotic pathways and networks regulated by Salmonella typhimurium in mouse intestinal infection in vivo. BMC Genom. 2010, 11, 722. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Supek, F.; Minana, B.; Valcarcel, J.; Gabaldon, T.; Lehner, B. Synonymous Mutations Frequently Act as Driver Mutations in Human Cancers. Cell 2014, 156, 1324–1335. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Diederichs, S.; Bartsch, L.; Berkmann, J.C.; Frose, K.; Heitmann, J.; Hoppe, C.; Iggena, D.; Jazmati, D.; Karschnia, P.; Linsenmeier, M.; et al. The dark matter of the cancer genome: Aberrations in regulatory elements, untranslated regions, splice sites, non-coding RNA and synonymous mutations. EMBO Mol. Med. 2016, 8, 442–457. [Google Scholar] [CrossRef] [PubMed]
  44. Keller, T.E.; Mis, S.D.; Jia, K.E.; Wilke, C.O. Reduced mRNA Secondary-Structure Stability Near the Start Codon Indicates Functional Genes in Prokaryotes. Genome Biol. Evol. 2012, 4, 80–88. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Chu, D.; Wei, L. Nonsynonymous, synonymous and nonsense mutations in human cancer-related genes undergo stronger purifying selections than expectation. BMC Cancer 2019, 19, 359. [Google Scholar] [CrossRef]
  46. Dettori, M.L.; Pazzola, M.; Paschino, P.; Amills, M.; Vacca, G.M. Association between the GHR, GHRHR, and IGF1 gene polymorphisms and milk yield and quality traits in Sarda sheep. J. Dairy Sci. 2018, 101, 9978–9986. [Google Scholar] [CrossRef] [Green Version]
  47. Li, S.B.; Zhou, H.T.; Zhao, F.F.; Fang, Q.; Wang, J.Q.; Liu, X.; Luo, Y.Z.; Hickford, J.G.H. Nucleotide Sequence Variation in the Insulin-Like Growth Factor 1 Gene Affects Growth and Carcass Traits in New Zealand Romney Sheep. DNA Cell Biol. 2020, 40, 265–271. [Google Scholar] [CrossRef]
  48. Abousoliman, I.; Reyer, H.; Oster, M.; Murani, E.; Mourad, M.; Rashed, M.A.S.; Mohamed, I.; Wimmers, K. Analysis of Candidate Genes for Growth and Milk Performance Traits in the Egyptian Barki Sheep. Animals 2020, 10, 197. [Google Scholar] [CrossRef] [Green Version]
  49. Routtu, J.; Hall, M.D.; Albere, B.; Beisel, C.; Bergeron, R.D.; Chaturvedi, A.; Choi, H.; Colbourne, J.; De Meester, L.; Stephens, M.T.; et al. An SNP-based second-generation genetic map of Daphnia magna and its application to QTL analysis of phenotypic traits. BMC Genom. 2014, 15, 1033. [Google Scholar] [CrossRef] [Green Version]
  50. Akey, J.; Jin, L.; Xiong, M.M. Haplotypes vs single marker linkage disequilibrium tests: What do we gain? Eur. J. Hum. Genet. 2001, 9, 291–300. [Google Scholar] [CrossRef] [Green Version]
  51. Martin, E.R.; Lai, E.H.; Gilbert, J.R.; Rogala, A.R.; Afshari, A.J.; Riley, L.; Finch, K.L.; Stevens, F.; Livak, K.J.; Slotterbeck, B.D.; et al. SNPing away at complex diseases: Analysis of single-nucleotide polymorphisms around APOE in Alzheimer disease. Am. J. Hum. Genet. 2000, 67, 383–394. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Stephens, J.C.; Schneider, J.A.; Tanguay, D.A.; Choi, J.; Acharya, T.; Stanley, S.E.; Jiang, R.H.; Messer, C.J.; Chew, A.; Han, J.H.; et al. Haplotype variation and linkage disequilibrium in 313 human genes. Science 2001, 293, 489–493. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Ryckman, K.K.; Simhan, H.N.; Krohn, M.A.; Williams, S.M. Molecular Human Reproduction: Predicting risk of bacterial vaginosis: The role of race, smoking and corticotropin-releasing hormone-related genes. Mol. Hum. Reprod. 2009, 15, 131–137. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Wang, T.J.; Zhang, F.; Richards, J.B.; Kestenbaum, B.; Meurs, J.B.V.; Berry, D.; Kiel, D.P.; Streeten, E.A.; Ohlsson, C.; Koller, D.L.; et al. Common Genetic Determinants of Vitamin D Insufficiency: A Genome-Wide Association Study. Lancet 2011, 376, 180–188. [Google Scholar] [CrossRef] [Green Version]
Genes 13 00666 g001 550
Figure 1. The sequencing peaks for three SNP loci of the IGF1 gene in Hulun Buir sheep. SNP1: c.144G>A (rs600896367); the arrow indicates the G–A mutation site. SNP2: c.150T>C (rs159876393); the arrow indicates the T–C mutation site. SNP3: c.495G>A (rs400398060); the arrow denotes to the G–A mutation site.
Figure 1. The sequencing peaks for three SNP loci of the IGF1 gene in Hulun Buir sheep. SNP1: c.144G>A (rs600896367); the arrow indicates the G–A mutation site. SNP2: c.150T>C (rs159876393); the arrow indicates the T–C mutation site. SNP3: c.495G>A (rs400398060); the arrow denotes to the G–A mutation site.
Genes 13 00666 g001
Genes 13 00666 g002 550
Figure 2. The sequencing peak maps for the ten detected SNP loci of the IGF1R gene in Hulun Buir sheep. SNP4: c.244C>T; the arrow denotes to the C–T mutation site. SNP5: c.714G>A (rs162159917); the arrow demonstrates the G–A mutation site. SNP6: c.924T>C (rs161166969); the arrow indicates the T–C mutation site, reverse sequenced as an A–G change. SNP7: c.939C>T (rs162159919); the arrow pinpoints the C–T mutation site, reverse sequenced as a G–A change. SNP8: c.1305T>C; the arrow points to the T–C mutation site. SNP9: c.1320G>A (rs601806812); the arrow indicates the G–A mutation site. SNP10: c.1401A>G (rs161166977); the arrow indicates the A–G mutation site. SNP11: c.1722T>C (rs161166984); the arrow demonstrates the T–C mutation site, reverse sequenced as an A–G change. SNP12: c.2253C>T (rs193644211); the arrow indicates the C–T mutation site. SNP13: c.2634C>G (rs161167008); the arrow pinpoints the C–G mutation site.
Figure 2. The sequencing peak maps for the ten detected SNP loci of the IGF1R gene in Hulun Buir sheep. SNP4: c.244C>T; the arrow denotes to the C–T mutation site. SNP5: c.714G>A (rs162159917); the arrow demonstrates the G–A mutation site. SNP6: c.924T>C (rs161166969); the arrow indicates the T–C mutation site, reverse sequenced as an A–G change. SNP7: c.939C>T (rs162159919); the arrow pinpoints the C–T mutation site, reverse sequenced as a G–A change. SNP8: c.1305T>C; the arrow points to the T–C mutation site. SNP9: c.1320G>A (rs601806812); the arrow indicates the G–A mutation site. SNP10: c.1401A>G (rs161166977); the arrow indicates the A–G mutation site. SNP11: c.1722T>C (rs161166984); the arrow demonstrates the T–C mutation site, reverse sequenced as an A–G change. SNP12: c.2253C>T (rs193644211); the arrow indicates the C–T mutation site. SNP13: c.2634C>G (rs161167008); the arrow pinpoints the C–G mutation site.
Genes 13 00666 g002
Genes 13 00666 g003 550
Figure 3. Associations for the SNPs of IGF1 gene with growth traits in Hulun Buir sheep. (A) The comparison of growth traits in SNP1 genotypes of IGF1 gene; WCG = chest girth at weaning (4-month-old); NBL = body length at 9-months-old. (B) The comparison of growth traits in SNP2 genotypes of IGF1 gene; NCG = chest girth at 9-months-old. (C) The comparison of growth traits in SNP3 genotypes of IGF1 gene; 4–9 ADG = average daily gain from 4 to 9-months-old. Different letters (small letters: p < 0.05) above the column indicate significant differences among the different genotypes.
Figure 3. Associations for the SNPs of IGF1 gene with growth traits in Hulun Buir sheep. (A) The comparison of growth traits in SNP1 genotypes of IGF1 gene; WCG = chest girth at weaning (4-month-old); NBL = body length at 9-months-old. (B) The comparison of growth traits in SNP2 genotypes of IGF1 gene; NCG = chest girth at 9-months-old. (C) The comparison of growth traits in SNP3 genotypes of IGF1 gene; 4–9 ADG = average daily gain from 4 to 9-months-old. Different letters (small letters: p < 0.05) above the column indicate significant differences among the different genotypes.
Genes 13 00666 g003
Genes 13 00666 g004 550
Figure 4. Associations for the SNPs of IGF1R gene with growth traits in Hulun Buir sheep. (A) Association analysis for different genotypes of SNP6 in the IGF1R gene with growth traits; NBL = body length at 9-months-old. (B) The comparison of body weight traits in SNP8 genotypes of IGF1R gene; BW = birth weight; WW = weaning weight (4-month-old); NBW = body weight at 9-months-old; 0–4 ADG = average daily gain from birth to 4-months-old; 0–9 ADG = average daily gain from birth to 9-months-old. (C) Association analyses for different genotypes of SNP8 in IGF1R with body size traits; WBH = body height at 4-months-old; WBL = body length at 4-months-old; WCG = chest girth at weaning (4-months-old); NBH = body height at 9-months-old; NBL = body length at 9-months-old; NCG = chest girth at 9-months-old. (D) The comparison of body weight traits in SNP13 genotypes of IGF1R gene; NBW = body weight at 9-months-old; 0–9 ADG = average daily gain from birth to 9-months-old. (E) Association analyses for different genotypes of SNP13 in the IGF1R with body size traits; WCG = chest girth at weaning (4-month-old); NBH = body height at 9-months-old; NBL = body length at 9-months-old; NCG = chest girth at 9-months-old. Different letters (small letters: p < 0.05; capital letters: p < 0.01) above the column indicate significant differences among the different genotypes.
Figure 4. Associations for the SNPs of IGF1R gene with growth traits in Hulun Buir sheep. (A) Association analysis for different genotypes of SNP6 in the IGF1R gene with growth traits; NBL = body length at 9-months-old. (B) The comparison of body weight traits in SNP8 genotypes of IGF1R gene; BW = birth weight; WW = weaning weight (4-month-old); NBW = body weight at 9-months-old; 0–4 ADG = average daily gain from birth to 4-months-old; 0–9 ADG = average daily gain from birth to 9-months-old. (C) Association analyses for different genotypes of SNP8 in IGF1R with body size traits; WBH = body height at 4-months-old; WBL = body length at 4-months-old; WCG = chest girth at weaning (4-months-old); NBH = body height at 9-months-old; NBL = body length at 9-months-old; NCG = chest girth at 9-months-old. (D) The comparison of body weight traits in SNP13 genotypes of IGF1R gene; NBW = body weight at 9-months-old; 0–9 ADG = average daily gain from birth to 9-months-old. (E) Association analyses for different genotypes of SNP13 in the IGF1R with body size traits; WCG = chest girth at weaning (4-month-old); NBH = body height at 9-months-old; NBL = body length at 9-months-old; NCG = chest girth at 9-months-old. Different letters (small letters: p < 0.05; capital letters: p < 0.01) above the column indicate significant differences among the different genotypes.
Genes 13 00666 g004
Genes 13 00666 g005 550
Figure 5. Linkage disequilibrium plot (r2) and haplotype blocks for SNPs of the IGF1R gene in Hulun Buir sheep. The values in boxes are pairwise SNP correlations (r2); dark red boxes indicate strong LD (r2 > 0.33) and light red boxes without numbers represent very weak LD (r2 < 0.001).
Figure 5. Linkage disequilibrium plot (r2) and haplotype blocks for SNPs of the IGF1R gene in Hulun Buir sheep. The values in boxes are pairwise SNP correlations (r2); dark red boxes indicate strong LD (r2 > 0.33) and light red boxes without numbers represent very weak LD (r2 < 0.001).
Genes 13 00666 g005
Genes 13 00666 g006 550
Figure 6. Linkage disequilibrium plot (r2) and haplotype blocks for SNPs of the IGF1 gene in Hulun Buir sheep. The values within boxes are pairwise SNP correlations (r2) and light red boxes represent very weak LD (r2 < 0.001).
Figure 6. Linkage disequilibrium plot (r2) and haplotype blocks for SNPs of the IGF1 gene in Hulun Buir sheep. The values within boxes are pairwise SNP correlations (r2) and light red boxes represent very weak LD (r2 < 0.001).
Genes 13 00666 g006
Genes 13 00666 g007 550
Figure 7. Associations for the haplotype combinations of SNPs in the IGF1R gene with growth traits in Hulun Buir sheep. (A) Association analysis for the haplotype combinations (block 1) of the IGF1R gene with growth traits; NBL = body length at 9–months–old. (B) The comparison of body weight traits for the haplotype combinations (block 2) of IGF1R gene in Hulun Buir sheep; BW = birth weight; WW = weaning weight (4–month–old); NBW = body weight at 9–months–old; 0–4 ADG = average daily gain from birth to 4–months–old; 4–9 ADG = average daily gain from 4 to 9–months–old; 0–9 ADG = average daily gain from birth to 9–months–old. (C) Association analyses for the haplotype combinations (block 2) of IGF1R gene with body size traits in Hulun Buir sheep; WBH = body height at 4 months of age; WBL = body length at 4–months–old; WCG = chest girth at weaning (4–month–old); NBH = body height at 9–months–old; NBL = body length at 9–months–old; NCG = chest girth at 9–months–old. Different letters (small letters: p < 0.05; capital letters: p < 0.01) above the column indicate significant differences among the different haplotype combinations.
Figure 7. Associations for the haplotype combinations of SNPs in the IGF1R gene with growth traits in Hulun Buir sheep. (A) Association analysis for the haplotype combinations (block 1) of the IGF1R gene with growth traits; NBL = body length at 9–months–old. (B) The comparison of body weight traits for the haplotype combinations (block 2) of IGF1R gene in Hulun Buir sheep; BW = birth weight; WW = weaning weight (4–month–old); NBW = body weight at 9–months–old; 0–4 ADG = average daily gain from birth to 4–months–old; 4–9 ADG = average daily gain from 4 to 9–months–old; 0–9 ADG = average daily gain from birth to 9–months–old. (C) Association analyses for the haplotype combinations (block 2) of IGF1R gene with body size traits in Hulun Buir sheep; WBH = body height at 4 months of age; WBL = body length at 4–months–old; WCG = chest girth at weaning (4–month–old); NBH = body height at 9–months–old; NBL = body length at 9–months–old; NCG = chest girth at 9–months–old. Different letters (small letters: p < 0.05; capital letters: p < 0.01) above the column indicate significant differences among the different haplotype combinations.
Genes 13 00666 g007
Table
Table 1. The information of SNP in IGF1 and IGF1R in Hulun Buir Sheep.
Table 1. The information of SNP in IGF1 and IGF1R in Hulun Buir Sheep.
GeneMutant LociSNPsRefSNPRegionAlleleAmino Acid VariationMutation Type
AB
IGF1c.144G>ASNP1rs600896367exon2GAAlasynonymous
c.150T>CSNP2rs159876393exon2TCProsynonymous
c.495G>ASNP3rs400398060exon5GAThrsynonymous
IGF1Rc.244C>TSNP4-exon3CTp.Arg81Cysnonsynonymous
c.714G>ASNP5rs162159917exon6GALyssynonymous
c.924T>CSNP6rs161166969exon8TCAspsynonymous
c.939C>TSNP7rs162159919 exon8CTCyssynonymous
c.1305T>CSNP8-exon11TCAspsynonymous
c.1320G>ASNP9rs601806812exon11GAThrsynonymous
c.1401A>GSNP10 rs161166977exon11AGAlasynonymous
c.1722T>CSNP11rs161166984exon12TCSersynonymous
c.2253C>TSNP12rs193644211exon17CTAlasynonymous
c.2634C>GSNP13rs161167008exon19CGGlysynonymous
Table
Table 2. Genetic diversity of the SNP loci within IGF1 and IGF1R genes in Hulun Buir sheep population.
Table 2. Genetic diversity of the SNP loci within IGF1 and IGF1R genes in Hulun Buir sheep population.
GeneSNPsGenotype FrequencyAllele FrequencyNeHoHePICP (HWE)
Wild TypeHybrid SubtypeMutant TypeWild TypeMutant Type
AAABBBAB
IGF1SNP10.9840.01600.9920.0081.0170.0160.0160.0160.057
SNP20.4900.4360.0740.7080.2921.7050.4360.4140.3280.604
SNP30.6460.3290.0250.8110.1891.4430.3290.3070.2600.485
IGF1RSNP40.9480.05200.9740.0261.0530.0520.0510.0490.085
SNP50.8100.19000.9050.0951.2080.1900.1720.1570.685
SNP60.3070.4510.2420.5320.4681.9920.4500.4980.3740.841
SNP70.3680.4930.1390.6150.3851.9000.4940.4740.3610.818
SNP80.3200.5110.1690.5760.4241.9550.5110.4890.3690.562
SNP90.7970.20300.8980.1021.2240.2030.1830.1660.085
SNP100.6930.2810.0260.8330.1671.3850.2810.2780.2390.685
SNP110.7230.2470.0300.8460.1541.3520.2470.2600.2260.841
SNP120.9310.06900.9650.0351.0720.0690.0670.0650.818
SNP130.4930.4160.0910.7010.2991.7210.4160.4190.3310.562
P (HWE) = P value of Hardy-Weinberg equilibrium, PIC < 0.25 demonstrates low polymorphism, 0.25 < PIC < 0.5 demonstrates medium polymorphism, PIC > 0.5 demonstrates high polymorphism.
Table
Table 3. Haplotype and haplotype combination analyses of SNPs (block1) in IGF1R gene.
Table 3. Haplotype and haplotype combination analyses of SNPs (block1) in IGF1R gene.
HaplotypeSNP6SNP7FrequencyHaplotype CombinationFrequency
H1 (TC)TC0.537H1H10.310
H2 (CT)CT0.389H1H20.402
H3 (CC)CC0.074H1H30.096
H2H20.052
H2H30.140
Table
Table 4. Haplotype and haplotype combination analyses of SNPs (block2) in IGF1R gene.
Table 4. Haplotype and haplotype combination analyses of SNPs (block2) in IGF1R gene.
HaplotypeSNP8SNP9FrequencyHaplotype CombinationFrequency
H4 (CG)CG0.321H1H10.114
H5 (TG)TG0.581H1H20.367
H6 (CA)CA0.098H1H30.048
H2H20.323
H2H30.148
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ding, N.; Tian, D.; Li, X.; Zhang, Z.; Tian, F.; Liu, S.; Han, B.; Liu, D.; Zhao, K. Genetic Polymorphisms of IGF1 and IGF1R Genes and Their Effects on Growth Traits in Hulun Buir Sheep. Genes 2022, 13, 666. https://doi.org/10.3390/genes13040666

AMA Style

Ding N, Tian D, Li X, Zhang Z, Tian F, Liu S, Han B, Liu D, Zhao K. Genetic Polymorphisms of IGF1 and IGF1R Genes and Their Effects on Growth Traits in Hulun Buir Sheep. Genes. 2022; 13(4):666. https://doi.org/10.3390/genes13040666

Chicago/Turabian Style

Ding, Ning, Dehong Tian, Xue Li, Zhichao Zhang, Fei Tian, Sijia Liu, Buying Han, Dehui Liu, and Kai Zhao. 2022. "Genetic Polymorphisms of IGF1 and IGF1R Genes and Their Effects on Growth Traits in Hulun Buir Sheep" Genes 13, no. 4: 666. https://doi.org/10.3390/genes13040666

APA Style

Ding, N., Tian, D., Li, X., Zhang, Z., Tian, F., Liu, S., Han, B., Liu, D., & Zhao, K. (2022). Genetic Polymorphisms of IGF1 and IGF1R Genes and Their Effects on Growth Traits in Hulun Buir Sheep. Genes, 13(4), 666. https://doi.org/10.3390/genes13040666

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop