Next Article in Journal
Association between Cardiac Auscultation and Echocardiographic Findings in Warmblood Horses
Previous Article in Journal
Alternatives in Education—Evaluation of Rat Simulators in Laboratory Animal Training Courses from Participants’ Perspective
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Polymorphisms of SORBS1 Gene and Their Correlation with Milk Fat Traits of Cattleyak

1
Key Laboratory of Qinghai-Tibetan Plateau Animal Genetic Reservation and Exploitation of Ministry of Education, Southwest Minzu University, Chengdu 610041, China
2
Key Laboratory of Animal Science of National Ethnic Affairs Commission, Southwest Minzu University, Chengdu 610041, China
*
Author to whom correspondence should be addressed.
Animals 2021, 11(12), 3461; https://doi.org/10.3390/ani11123461
Submission received: 1 November 2021 / Revised: 29 November 2021 / Accepted: 3 December 2021 / Published: 5 December 2021
(This article belongs to the Section Animal Reproduction)

Abstract

:

Simple Summary

Increasing milk fat rate has a good effect on the milk quality of cattleyak. SNPs can help us find potential molecular markers for the milk fat traits of cattleyak, and they can be screened according to molecular markers when they are young. It provides a reference for cultivating high milk fat cattle population in the future. The results of this study suggest that the SORBS1 gene polymorphism is closely related to the milk fat traits of cattleyak, which could be used as a candidate genetic marker for milk fat trait selection in cattleyak. This study provides a new molecular marker and theoretical basis for screening the milk fat traits of cattleyak. It has a certain reference value for the research and improvement of milk quality.

Abstract

This study aimed to find the SNPs in the SORBS1 gene of cattleyak, analyze the relationship between its polymorphisms and the milk fat traits, and find potential molecular markers for the milk fat traits of cattleyak. The polymorphism of the SORBS1 gene in 350 cattleyak from Hongyuan County (Sichuan, China) were detected by PCR and DNA sequencing, and the correlation between these SNPs and the milk production traits of cattleyak was analyzed. The results showed that there were nine SNPs in the CDS and their adjacent non-coding regions of the SORBS1 gene, and all SNPs have three genotypes. The correlation analysis found that the genotypes with superior milk fat traits in the other eight alleles were homozygous genotypes with a high genotype frequency except the g.96284 G > A (c.3090 G > A) (p < 0.05). However, at locus g.96284 G > A, the milk fat percentage, monounsaturated fatty acids (MUFAs), polyunsaturated fatty acids (PUFAs) and saturated fatty acids (SFAs) of the GA genotype were significantly higher than that of GG and AA genotypes (p < 0.05). Among these SNPs, three SNPs (g.6256 C > T (c.298 C > T), g.24791 A > G (c.706 A > G) and g.29121 A > G (c.979 A > G)) caused the amino acids change. The genotypes of the three SNPs consist of three haplotypes and four diplotypes. The amino acid mutation degree of diplotype H1–H1 (CCAAAA) was the highest, and its milk fat percentage, MUFAs, PUFAs and SFAs were also the highest (p < 0.05). Taken together, we found nine SNPs in the SORBS1 gene that are closely related to the milk fat traits of cattleyak. Moreover, the mutation of amino acids caused by SNPs had positive effects on the milk fat traits of cattleyak. H1-H1 is the dominant diplotype which significantly related to the milk fat traits of cattleyak. This study provides a new molecular marker and theoretical basis for screening the milk fat traits of cattleyak.

1. Introduction

Yak (Bos grunniens), a wild bovine species that lives 3500 m above sea level, is an iconic symbol of the Qinghai–Tibetan Plateau region and the nearby areas. Moreover, yak raising has become the main livelihood in the local area. However, the reproductive efficiency and production performance of yak are low, which cannot satisfy the local herdsman needs and industrial development [1]. As the hybrid offspring of yak (♀) and cattle (Bos taurus) (♂), the female cattleyak has obvious heterosis compared to yak, especially in its lactation performance [2,3]. Dairy products made from fermented and processed milk are popular with local herdsmen and tourists. As the main energy substance in milk, milk fat only accounts for 3~5% of the milk content, but it determines the nutritional value of milk [4,5]. Therefore, it is greatly meaningful to increase the percentage of milk fat and improve the composition of fatty acids in cattleyak.
The relevance between polymorphism and milk quality has been widely applied to the breeding of lactation traits, and the relationship between milk composition and genetic variation in milk production has a far-reaching impact on milk production [6,7,8,9,10,11,12]. In genetics, milk fatty acids have been found to be heritable, and the heritability was estimated at 0.22 to 0.71 [13]. In the past few years, with the rapid development of biotechnology, the mapping of candidate genes and quantitative trait locus (QTL) has been gradually applied to analyze lactation traits and breeding. Therefore, some promising SNPs have been shown to be significantly associated with lactation and its milk components [14,15,16].
The sorbin and SH3 domain-containing 1 (SORBS1) gene is a member of the SORBS family, and it encodes a Cbl-associated protein (CAP) that plays key roles in the signaling and stimulation of insulin [17,18,19]. Previous studies found that after a postprandial lipid (PPL) challenge, fatty acids and sterols related to cholesterol absorption were increased, while sterols related to cholesterol synthesis were decreased, and they detected that two SNPs (rs.12247017 and rs.12240292) in the SORBS1 gene were correlated with b-Sitosterol after correcting for multiple testing [20]. In addition, SORBS1 blocks the insulin stimulation of glucose transport in 3T3-L1 adipocytes [21]. The T228A polymorphism of SORBS1 is correlated with obesity and diabetes [22]. Recently, it was found that the SORBS1 gene has an effect on cell proliferation and migration, which can inhibit tumor metastasis and improve the sensitivity of tumors to chemotherapeutic drugs [23]. Above all, SORBS1 has multiple functions, but there are still many potential mechanisms that need to be studied.
A genome-wide association study revealed that the SORBS1 gene was one of 20 new candidate genes that were closely related to milk fat traits in Chinese Holstein cows [24]. However, there are few studies on the breeding of the milk fat traits of cattleyak, and there is no report on the association between cattleyak and the SORBS1 gene. Therefore, whether the SORBS1 gene is related to the creamy traits of cattleyak remains to be determined. The cattleyak is the subject of this experiment, the detected SNPs of the SORBS1 gene and its association with cattleyak milk fat traits was analyzed, and we looked for the molecular markers for the assisted selection of the milk fat traits of cattleyak.

2. Materials and Methods

2.1. Animal and Preparation

The study samples and data were collected in Ngawa Tibetan and Qiang Autonomous Prefecture Hongyuan County, Sichuan Province, China (longitudinal 32°78′ N to 32°82′ N, latitudinal 102°52′ E to 102°55′ E). In total, 350 female cattleyak (Jersey cattle ♂ × yak ♀) were selected in this experiment. The female parent of our cattleyak comes from a group of relatively closed and natural mating yak, while the male parent comes from the frozen semen of the same batch of Jersey cattle. All animals were collected from August, and they have the same parity and lactation, a similar growth environment and equivalent nutritional conditions. The 25 mL milk samples of each cattleyak were collected and sent to the Hongya Yangping branch of New Hope Co., Ltd. The physical and chemical indexes (DHI) of each sample were analyzed by a milk composition analyzer. The 10 mL blood samples were collected from the jugular vein of each cattleyak by using a sterile blood collection device, and then put into the 4 °C ice box. All blood samples were taken back to the laboratory and placed in a freezer at −20 °C.

2.2. DNA Extraction from Blood

The DNA was extracted from whole-blood samples using the Blood Genomic DNA Extraction Kit (Solarbio, Beijing, China). The purity and concentration of the DNA were detected by a nucleic acid concentration detector, and the integrity of the DNA was then detected by agarose gel electrophoresis.

2.3. Primer Design and Sequencing

The standard sequence is from NCBI, which has 31 exons and 31 introns (Accession Number: NW_005394038.1). Based on the DNA sequence, 12 pairs of primers were designed for the PCR amplification of the SORBS1 gene CDS and their adjacent regions (Table 1). The PCR reaction was 25 µL: 12.5 µL 2× Rapid Taq Master Mix, l µL upstream and 1 µL downstream primers, 1 µL DNA template and add ddH2O to 25 µL. PCR reaction conditions were as follows: denaturation for 3 min at 95 °C, 35 cycles for 15 s at 95 °C, annealing temperature for 15 s and 72 °C for 15 s, and extension at 72 °C for 5 min. PCR products were detected by 2% agarose gel electrophoresis, and sent to Tsingke Biotechnology Co., Ltd. for sequencing.

2.4. Genotyping and Statistical Analysis

Based on the sequencing results, Chromas 2 software was used to analyze and find SNP sites. Using Excel 2019 to calculate the genotype frequency and allele frequency, the polymorphic information content (PIC) was calculated with different genotype frequency and allele frequency. The formula is as follows:
P I C = 1 H 0 i = 1 n 1 j i + 1 n 2 P i 2 P j 2
Pi and Pj are the frequencies of the i and j alleles, and n is the number of alleles. The Standard Amino acid sequence was derived from NCBI and compared with the pre- and post-mutation amino acid sequence in BLAST to test if the SNPs changed the amino acid sequence. Haploview 4.2 software was used to analyze the haplotype. The effect of different genotypes and diplotypes on the milk fat traits of cattleyak was analyzed by the Analyze–General linear model univariate in SPSS 19.0 software. The formula is as follows:
Yi = µ + Gi + e
Among them, Yi is the phenotypic value of milk fat traits, and µ is the population mean of milk fat traits, Gi is the total effect and e is the random residual effect. Differences were considered significant at p < 0.05.

3. Results

3.1. PCR Amplification and SNP Screening

The concentration and purity of genomic DNA in cattleyak blood were detected by a nucleic acid concentration detector, and the samples were all qualified. Using 2% agarose gel electrophoresis to detect four amplified products, clear and single-purpose bands were obtained. The specificity is good and can be used for subsequent sequencing analysis (Figure 1).
A total of nine SNPs in the CDS and their adjacent regions of the SORBS1 gene were detected in the cattleyak population by PCR and sequencing (Figure 2). g.6256 C > T (c.298 C > T) and g.24791 G > A (c.706 A > G) were detected by primer 1 and primer 4 PCR product, g.29029 C > T, g.29050 A > G, g.29121 A > G (c.979 A > G), g.29245 T > C, g.29305 T > C and g.29347 T > C were detected by primer 5 PCR product, g.96284 G > A (c.3090 G > A) was detected by primer 10 PCR product. Additionally, four SNPs (g.6256 C > T, g.24791 G > A, g.29121 A > G and g.96284 G > A) are located in the CDS region, while the other five SNPs are located in the intron region adjacent of the SORBS1 gene.

3.2. SNP Genotyping of the SORBS1 Gene

The genotype frequency, allele frequency and polymorphic information content of the nine SNP loci were shown in Table 2. The results showed that all the nine SNP loci had three genotypes in the cattleyak population. The PIC of g.6256 C > T, g.24791 G > A, g.29029 C > T and g.29050 A > G were 0.361, g.29121 A > G was 0.346, g.29245 T > C, g.29305 T > C, g.29347 T > C and g.96284 G > A were 0.358. The PIC of all SNPs was between 0.25~0.5, and all of them were moderately polymorphic. All the other SNPs conformed to the Hardy–Weinberg balance except g.29121 A > G site (p > 0.05).

3.3. Haplotype and Diplotype Analysis of Each Amino Acid Mutation Site of SORBS1

During the comparison of the amino acid residue sequence alignment of the SORBS1 gene before and after mutation, we found that three SNPs (g.6256 C > T, g.24791 A > G and g.29121 A > G) caused amino acid changes among the four SNPs located in the CDS region. Moreover, no amino acid changes were observed at g.96284 G > A. The three mutations were Ser100Pro, Arg229His and Ala327Thr (Table 2). Using the three SNPs, a haplotype analysis was carried out using Haploview software (Table 3). The haplotype analysis showed that the three haplotypes were in fact found in the population of cattleyak, of which H1 (C-A-A) and H2 (T-G-G) were the major haplotypes, accounting for 61.9% and 33.3%, respectively. The frequency of H3 (T-G-A) was 4.9%. Based on these three haplotypes, we further obtained four combinations, in the form of H1–H1 (CCAAAA), H1–H2 (CTAGAG), H1–H3 (CTAGAA) and H2–H2 (TTGGGG) (Figure 3), of which H1–H1 (38.3%) and H2–H2 (37.4%) belonged to dominant diplotypes.

3.4. Association Analysis between SNP Genotypes and Milk Fat Traits

Based on the data of the SNP genotyping, the correlation analysis of a single locus with the milk fat percentage and fatty acids of cattleyak was carried out. As shown in Table 4, in terms of milk fat percentage and fatty acid content, most of the homozygous (CC genotype at g.6256 C > T, AA genotype at g.24791 A > G, TT genotype at g.29029 T > C, AA genotype at g.29050 A > G, AA genotype at g.29121 A > G, CC genotype at g.29245 C > T, CC genotype at g.29305 C > T, TT genotype at g.29347 T > C and GA genotype at g.96284 G > A) were significantly higher than the other genotypes (p < 0.05). At locus g.96284 G > A, the milk fat percentage, monounsaturated fatty acids (MUFAs), polyunsaturated fatty acids (PUFAs) and saturated fatty acids (SFAs) of the GA genotype were significantly higher than that in GG and AA genotypes (p < 0.05).

3.5. Association Analysis between Genotype Combinations and the Milk Fat Traits of Cattleyak

The genotype combination of H2–H2 did not cause the change of amino acid. There is a possibility of amino acid variation in the part of the cattleyak whose genotype combinations are H1–H2 and H1–H3. The base mutation in the combination of H1–H1 resulted in the mutation of all three amino acids. The results of the ANOVA analysis showed that H1–H1 had the highest milk fat percentage, MUFAs, PUFAs and SFAs, while H2–H2 had the lowest milk fat percentage and fatty acid content (Table 5). The differences were statistically significant (p < 0.05).

4. Discussion

Recently, the study of the polymorphism has been used to optimize the performance and composition of dairy cows [25,26,27,28]. The GWAS data in Li and his colleagues’ study suggested that a SNP in the SORBS1 gene can significantly influence the composition of fatty acids in milk, especially MUFAs and PUFAs [24]. The fatty acids in milk have positive effects on human health. For example, previous study showed that MUFAs had a positive effect on human health due to their cholesterol-lowering properties [29]. PUFAs are essential nutrients which play the important role in plasma lipids and endothelial function. Additionally, PUFAs can be also used in the prevention and treatment of coronary heart disease [30]. Hereby, we detected the potential SNP in the SORBS1 gene of cattleyak, and we revealed the effect of the SNPs in the SORBS1 CDS and their adjacent non-coding regions on the milk fat traits of cattleyak.
In this study, a total of nine SNPs located in the CDS and their adjacent non-coding regions of the SORBS1 gene in cattleyak were detected. After genotyping, all nine SNPs were found that have three genotypes in the population of cattleyak, and the dominant genotypes were mainly heterozygotes. In addition, our results showed that the genetic richness was high and had the potential of genetic variation, which lead to a more abundant selection effect. The natural selection, artificial selection intervention, genetic drift and other factors could lead to changes in some animal nucleotides so that the biological genetic traits become more colorful [31]. The size of the population genetic variation represents the level of the population genetic richness. Generally speaking, the higher the population genetic variation is, the higher is the population genetic richness. Thus, the PIC is often used as an indicator [32,33]. With the exception of g.29121 A > G, the other eight SNPs have reached the Hardy–Weinberg balance, which suggested that these SNPs of cattleyak were not strongly selected and might not be influenced by mutation, selection and genetic drift. It is also possible that the alleles and genotypes have a genetic advantage in their adaptability that rebalanced after a long period of artificial selection and breeding [34]. However, g.29121 A > G was out of balance, which means that traces of the artificial selection can be found. Thus, the other eight alleles are more likely to be rebalanced after a long period of artificial selection and reproduction.
In a previous study, it was generally assumed that only SNPs in exons were associated with biological traits, and these SNPs in introns did not affect biological traits [35,36,37]. In this study, five SNPs (g.29029 T > C, g.29050 A > G, g.29245 C > T, g.29305 C > T and g.29347 T > C), located in the introns of the SORBS1 gene were found to be significantly associated with the cattleyak milk fat traits. A recent study found that introns were involved in the gene splicing in eukaryotes, and that the mutation of introns might change the splicing efficiency, which resulted in the change of amino acid coding and eventually affected the biological traits [38,39]. Even if the intron mutation has no effect on the amino acid sequence, it may also result in the inactivation of the splicing site, or it may contain enhancer sequences to promote the specific transcription of the gene, affecting the animal phenotype [40,41]. The correlation between the SNPs on the introns and the milk fat traits of cattleyak may have a causal mutation, or it may be that these five SNPs are closely related to the SNPs on the exons. The exact cause needs further functional verification in a future study.
The effects of a specific gene on biological traits might be related to a combination of SNPs in that gene [42]. This study focuses on the genotype combination of SNPs causing amino acid changes. It was found that three SNPs (g.6256 C > T, g.24791 A > G, g.29121 A > G) caused amino acids changes and we obtained three haplotypes (Figure S1). Based on these three haplotypes, four combinations were obtained. The amino acids of the three SNPs in H1-H1 were changed compared with the standard sequence. Moreover, the milk fat content and fatty acid content of the H1-H1 population were significantly higher than those of other diplotypes. The milk fat content and fatty acid content of the H2-H2 population without amino acid mutation was the lowest. Therefore, it is clear that the mutation of amino acids has a positive effect on the percentage of milk fat and the content of fatty acids.

5. Conclusions

In conclusion, nine SNPs were detected in the CDS and their adjacent non-coding regions of the SORBS1 gene of the cattleyak, and all of these SNPs had a significant correlation with their milk fat traits. Additionally, H1-H1 (CCAAAA) is the dominant diplotype significantly related to the milk fat traits of cattleyak. These sites provide new molecular markers and a theoretical basis for the screening of the milk fat traits of cattleyak.

Supplementary Materials

The following are available online at www.mdpi.com/article/10.3390/ani11123461/s1. Figure S1. The mutation of SNPs caused the change of amino acid sequence (red maker).

Author Contributions

Experiment design, L.Y. and X.X.; collection of blood samples, L.Y., X.M. and Y.H.; software, X.M. and H.Y.; formal analysis, L.Y. and M.Y.; writing—original draft preparation, L.Y.; writing—review and editing, X.X., Y.Z. and J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Science and Technology Department of Sichuan province Key R and D Program (grant No. 2021YFN0121), the Fund of China Scholarship Council (No. 202008510085) and the Fundamental Research Funds for the Central Universities of Southwest Minzu University (grant No. 2021PTJS20).

Institutional Review Board Statement

The experimental animals and operating procedures in this study obeyed the policies of the national and local animal welfare agencies, and they were approved by the Animal Care and Use Committee of Southwest Minzu University and performed in accordance with animal welfare and ethics.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

No conflict of interest exists in the submission of this manuscript, and the manuscript is approved by all authors for publication.

References

  1. Zi, X.D. Reproduction in female yaks (Bos grunniens) and opportunities for improvement. Theriogenology 2003, 59, 1303–1312. [Google Scholar] [CrossRef]
  2. Barsila, S.R.; Devkota, N.R.; Kreuzer, M.; Marquardt, S. Effects of different stocking densities on performance and activity of cattle × yak hybrids along a transhumance route in the Eastern Himalaya. Springerplus 2015, 4, 398. [Google Scholar] [CrossRef] [Green Version]
  3. Shah, M.A.; Xu, C.; Wu, S.; Zhao, W.; Luo, H.; Yi, C.; Liu, W.; Cai, X. Isolation and characterization of spermatogenic cells from cattle, yak and cattleyak. Anim. Reprod. Sci. 2018, 193, 182–190. [Google Scholar] [CrossRef]
  4. Filipczak-Fiutak, M.; Pluta-Kubica, A.; Domagała, J.; Duda, I.; Migdał, W. Nutritional value and organoleptic assessment of traditionally smoked cheeses made from goat, sheep and cow’s milk. PLoS ONE 2021, 16, e0254431. [Google Scholar] [CrossRef]
  5. Park, C.H.; Ranaraja, U.; Dang, C.G.; Kim, J.J.; Do, C.H. Genetic parameters for milk fatty acid composition of Holstein in Korea. Asian-Australas. J. Anim. Sci. 2020, 33, 1573–1578. [Google Scholar] [CrossRef]
  6. Liang, Y.; Gao, Q.; Zhang, Q.; Arbab, A.; Li, M.; Yang, Z.; Karrow, N.A.; Mao, Y. Polymorphisms of the ACSL1 Gene Influence Milk Production Traits and Somatic Cell Score in Chinese Holstein Cows. Animals 2020, 10, 2282. [Google Scholar] [CrossRef]
  7. Sanchez, M.P.; Ferrand, M.; Gelé, M.; Pourchet, D.; Miranda, G.; Martin, P.; Brochard, M.; Boichard, D. Short communication: Genetic parameters for milk protein composition predicted using mid-infrared spectroscopy in the French Montbéliarde, Normande, and Holstein dairy cattle breeds. J. Dairy Sci. 2017, 100, 6371–6375. [Google Scholar] [CrossRef]
  8. Prakapenka, D.; Liang, Z.; Jiang, J.; Ma, L.; Da, Y. A Large-Scale Genome-Wide Association Study of Epistasis Effects of Production Traits and Daughter Pregnancy Rate in U.S. Holstein Cattle. Genes 2021, 12, 1089. [Google Scholar] [CrossRef]
  9. Sun, X.; Liang, Y.; Gao, Q.; Guo, J.; Tang, C.; Shi, K.; Yang, Z.; Mao, Y. AGPAT3 Gene polymorphisms are associated with milk production traits in Chinese Holstein cows. J. Dairy Res. 2021, 15, 1–6. [Google Scholar] [CrossRef]
  10. Brzáková, M.; Rychtářová, J.; Čítek, J.; Sztankóová, Z. A Candidate Gene Association Study for Economically Important Traits in Czech Dairy Goat Breeds. Animals 2021, 11, 1796. [Google Scholar] [CrossRef]
  11. Yang, F.; Zhang, M.; Rong, Y.; Liu, Z.; Yang, S.; Zhang, W.; Li, J.; Cai, Y. A Novel SNPs in Alpha-Lactalbumin Gene Effects on Lactation Traits in Chinese Holstein Dairy Cows. Animals 2019, 10, 60. [Google Scholar] [CrossRef] [Green Version]
  12. Pegolo, S.; Mota, L.; Bisutti, V.; Martinez-Castillero, M.; Giannuzzi, D.; Gallo, L.; Schiavon, S.; Tagliapietra, F.; Revello Chion, A.; Trevisi, E.; et al. Genetic parameters of differential somatic cell count, milk composition, and cheese-making traits measured and predicted using spectral data in Holstein cows. J. Dairy Sci. 2021, 104, 10934–10949. [Google Scholar] [CrossRef]
  13. Stoop, W.M.; van Arendonk, J.A.; Heck, J.M.; van Valenberg, H.J.; Bovenhuis, H. Genetic parameters for major milk fatty acids and milk production traits of Dutch Holstein-Friesians. J. Dairy Sci. 2008, 91, 385–394. [Google Scholar] [CrossRef]
  14. Naserkheil, M.; Miraie-Ashtiani, S.R.; Sadeghi, M.; Nejati-Javaremi, A.; Park, C.W.; Min, K.S.; Lee, D. Exploring novel single nucleotide polymorphisms and haplotypes of the diacylglycerol O-acyltransferase 1 (DGAT1) gene and their effects on protein structure in Iranian buffalo. Genes Genom. 2019, 41, 1265–1271. [Google Scholar] [CrossRef]
  15. Schennink, A.; Heck, J.M.; Bovenhuis, H.; Visker, M.H.; van Valenberg, H.J.; van Arendonk, J.A. Milk fatty acid unsaturation: Genetic parameters and effects of stearoyl-CoA desaturase (SCD1) and acyl CoA: Diacylglycerol acyltransferase 1 (DGAT1). J. Dairy Sci. 2008, 91, 2135–2143. [Google Scholar] [CrossRef] [Green Version]
  16. Schennink, A.; Stoop, W.M.; Visker, M.H.; van der Poel, J.J.; Bovenhuis, H.; van Arendonk, J.A. Short communication: Genome-wide scan for bovine milk-fat composition. II. Quantitative trait loci for long-chain fatty acids. J. Dairy Sci. 2009, 92, 4676–4682. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Chang, T.J.; Wang, W.C.; Hsiung, C.A.; He, C.T.; Lin, M.W.; Sheu, W.H.; Chang, Y.C.; Quertermous, T.; Chen, I.; Rotter, J.; et al. Genetic variation in the human SORBS1 gene is associated with blood pressure regulation and age at onset of hypertension: A SAPPHIRe Cohort Study. Medicine 2016, 95, e2970. [Google Scholar] [CrossRef]
  18. Chang, T.J.; Wang, W.C.; Hsiung, C.A.; He, C.T.; Lin, M.W.; Sheu, W.H.; Chang, Y.C.; Quertermous, T.; Chen, Y.I.; Rotter, J.I.; et al. Genetic variation of SORBS1 gene is associated with glucose homeostasis and age at onset of diabetes: A SAPPHIRe Cohort Study. Sci. Rep. 2018, 8, 10574. [Google Scholar] [CrossRef]
  19. Germain, M.; Pezzolesi, M.G.; Sandholm, N.; McKnight, A.J.; Susztak, K.; Lajer, M.; Forsblom, C.; Marre, M.; Parving, H.H.; Rossing, P.; et al. SORBS1 gene, a new candidate for diabetic nephropathy: Results from a multi-stage genome-wide association study in patients with type 1 diabetes. Diabetologia 2015, 58, 543–548. [Google Scholar] [CrossRef] [Green Version]
  20. Irvin, M.R.; Zhi, D.; Aslibekyan, S.; Claas, S.A.; Absher, D.M.; Ordovas, J.M.; Tiwari, H.K.; Watkins, S.; Arnett, D.K. Genomics of post-prandial lipidomic phenotypes in the Genetics of Lipid lowering Drugs and Diet Network (GOLDN) study. PLoS ONE. 2014, 9, e99509. [Google Scholar] [CrossRef]
  21. Baumann, C.A.; Ribon, V.; Kanzaki, M.; Thurmond, D.C.; Mora, S.; Shigematsu, S.; Bickel, P.E.; Pessin, J.E.; Saltiel, A.R. CAP defines a second signalling pathway required for insulin-stimulated glucose transport. Nature 2000, 407, 202–207. [Google Scholar] [CrossRef]
  22. Lin, W.H.; Chiu, K.C.; Chang, H.M.; Lee, K.C.; Tai, T.Y.; Chuang, L.M. Molecular scanning of the human sorbin and SH3-domain-containing-1 (SORBS1) gene: Positive association of the T228A polymorphism with obesity and type 2 diabetes. Hum. Mol. Genet. 2001, 10, 1753–1760. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Cho, W.C.; Jang, J.E.; Kim, K.H.; Yoo, B.C.; Ku, J.L. SORBS1 serves a metastatic role via suppression of AHNAK in colorectal cancer cell lines. Int. J. Oncol. 2020, 56, 1140–1151. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Li, C.; Sun, D.; Zhang, S.; Wang, S.; Wu, X.; Zhang, Q.; Liu, L.; Li, Y.; Qiao, L. Genome wide association study identifies 20 novel promising genes associated with milk fatty acid traits in Chinese Holstein. PLoS ONE 2014, 9, e96186. [Google Scholar]
  25. Cao, M.; Shi, L.; Peng, P.; Han, B.; Liu, L.; Lv, X.; Ma, Z.; Zhang, S.; Sun, D. Determination of genetic effects and functional SNPs of bovine HTR1B gene on milk fatty acid traits. BMC Genom. 2021, 22, 575. [Google Scholar] [CrossRef]
  26. Singh, A.; Kumar, A.; Gondro, C.; da Silva Romero, A.R.; Karthikeyan, A.; Mehrotra, A.; Pandey, A.K.; Dutt, T.; Mishra, B.P. Identification of genes affecting milk fat and fatty acid composition in Vrindavani crossbred cattle using 50 K SNP-Chip. Trop. Anim. Health Prod. 2021, 53, 347. [Google Scholar] [CrossRef]
  27. Zwierzchowski, L.; Ostrowska, M.; Żelazowska, B. Bagnicka, E. Single nucleotide polymorphisms in the bovine SLC2A12 and SLC5A1 glucose transporter genes—The effect on gene expression and milk traits of Holstein Friesian cows. Anim. Biotechnol. 2021, 6, 1–11. [Google Scholar] [CrossRef]
  28. Bordonaro, S.; Tumino, S.; Marletta, D.; De Angelis, A.; Di Paola, F.; Avondo, M.; Valenti, B. Effect of GH p.L127V Polymorphism and Feeding Systems on Milk Production Traits and Fatty Acid Composition in Modicana Cows. Animals 2020, 10, 1651. [Google Scholar] [CrossRef]
  29. Labonté, M.È.; Jenkins, D.J.; Lewis, G.F.; Chiavaroli, L.; Wong, J.M.; Kendall, C.W.; Hogue, J.C.; Couture, P.; Lamarche, B. Adding MUFA to a dietary portfolio of cholesterol-lowering foods reduces apoAI fractional catabolic rate in subjects with dyslipidaemia. Br. J. Nutr. 2013, 110, 426–436. [Google Scholar] [CrossRef] [Green Version]
  30. Golanski, J.; Szymanska, P.; Rozalski, M. Effects of Omega-3 polyunsaturated fatty acids and their metabolites on haemostasis-current perspectives in cardiovascular disease. Int. J. Mol. Sci. 2021, 22, 2394. [Google Scholar] [CrossRef]
  31. Mäki-Tanila, A.; Webster, L. Heritability, SNP, inbreeding, dairy cattle, genomic selection-and other keywords. J. Anim. Breed. Genet. 2019, 136, 1–2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Azimu, W.; Manatbay, B.; Li, Y.; Kaimaerdan, D.; Wang, H.E.; Reheman, A.; Muhatai, G. Genetic diversity and population structure analysis of eight local chicken breeds of Southern Xinjiang. Br. Poult. Sci. 2018, 59, 629–635. [Google Scholar] [CrossRef]
  33. Santos, C.P.; Aguiar, A.F.; Giometti, I.C.; Mariano, T.B.; de Freitas, C.; Nai, G.A.; de Freitas, S.Z.; Pai-Silva, M.D.; Pacagnelli, F.L. High final energy of gallium arsenide laser increases MyoD gene expression during the intermediate phase of muscle regeneration after cryoinjury in rats. Lasers Med. Sci. 2018, 33, 843–850. [Google Scholar] [CrossRef] [Green Version]
  34. Huang, Y.Z.; He, H.; Sun, J.J.; Wang, J.; Li, Z.J.; Lan, X.Y.; Lei, C.Z.; Zhang, C.L.; Zhang, E.P.; Wang, J.Q.; et al. Haplotype combination of SREBP-1c gene sequence variants is associated with growth traits in cattle. Genome 2011, 54, 507–516. [Google Scholar] [CrossRef] [PubMed]
  35. He, X.; Zhang, Z.; Chu, M. The effect of SNP rs400827589 in exon 2 of the MTNR1B gene on reproductive seasonality and litter size in sheep. Vet. Med. Sci. 2020, 6, 804–812. [Google Scholar] [CrossRef]
  36. Jones, M.; Sergeant, C.; Richardson, M.; Groth, D.; Brooks, S.; Munyard, K. A non-synonymous SNP in exon 3 of the KIT gene is responsible for the classic grey phenotype in alpacas (Vicugna pacos). Anim. Genet. 2019, 50, 493–500. [Google Scholar] [CrossRef]
  37. Mesrian Tanha, H.; Rahgozar, S.; Mojtabavi Naeini, M. ABCC4 functional SNP in the 3′ splice acceptor site of exon 8 (G912T) is associated with unfavorable clinical outcome in children with acute lymphoblastic leukemia. Cancer Chemother. Pharmacol. 2017, 80, 109–117. [Google Scholar] [CrossRef]
  38. Xiang, G.; Ren, J.; Hai, T.; Fu, R.; Yu, D.; Wang, J.; Li, W.; Wang, H.; Zhou, Q. Editing porcine IGF2 regulatory element improved meat production in Chinese Bama pigs. Cell. Mol. Life Sci. 2018, 75, 4619–4628. [Google Scholar] [CrossRef]
  39. Yadav, T.; Magotra, A.; Bangar, Y.C.; Kumar, R.; Yadav, A.S.; Garg, A.R.; Bahurupi, P.; Kumar, P. Effect of BsaA I genotyped intronic SNP of leptin gene on production and reproduction traits in Indian dairy cattle. Anim. Biotechnol. 2021, 9, 1–7. [Google Scholar] [CrossRef] [PubMed]
  40. Mi, T.; Liu, K.; Guo, T.; Li, L.; Wang, Y.; Li, C.; Cui, Y.; Dai, J.; Zhang, Y.; Hu, S. Analysis of the eighth intron polymorphism of NR6A1 gene in sheep and its correlation with lumbar spine number. Anim. Biotechnol. 2021, 4, 1–7. [Google Scholar] [CrossRef]
  41. Nair, V.; Sankaranarayanan, R.; Vasavada, A.R. Deciphering the association of intronic single nucleotide polymorphisms of crystallin gene family with congenital cataract. Indian J. Ophthalmol. 2021, 69, 2064–2070. [Google Scholar] [PubMed]
  42. Magotra, A.; Gupta, I.D.; Ahmad, T.; Alex, R. Polymorphism in DNA repair gene BRCA1 associated with clinical mastitis and production traits in indigenous dairy cattle. Res. Vet. Sci. 2020, 133, 194–201. [Google Scholar] [CrossRef] [PubMed]
Figure 1. PCR amplification products (12 pairs of primers) for CDS and their adjacent regions of the SORBS1 gene.
Figure 1. PCR amplification products (12 pairs of primers) for CDS and their adjacent regions of the SORBS1 gene.
Animals 11 03461 g001
Figure 2. Sequencing diagram of heterozygous genotypes in the SNP loci of the SORBS1 gene (red box marker). Four SNPs (A,B,E,I) are located in the CDS region of the SORBS1 gene, c. system naming indicates the position of the site on the cDNA sequence. Five SNPs (C,D,F,G,H) are located in the intron region adjacent to the CDS region of the SORBS1 gene.
Figure 2. Sequencing diagram of heterozygous genotypes in the SNP loci of the SORBS1 gene (red box marker). Four SNPs (A,B,E,I) are located in the CDS region of the SORBS1 gene, c. system naming indicates the position of the site on the cDNA sequence. Five SNPs (C,D,F,G,H) are located in the intron region adjacent to the CDS region of the SORBS1 gene.
Animals 11 03461 g002
Figure 3. Diplotypes and frequencies of the three mutations in the SORBS1 gene.
Figure 3. Diplotypes and frequencies of the three mutations in the SORBS1 gene.
Animals 11 03461 g003
Table 1. Conventional PCR amplification primers for CDS and their adjacent regions of the SORBS1 gene.
Table 1. Conventional PCR amplification primers for CDS and their adjacent regions of the SORBS1 gene.
Primer NameSequence(5′-3′)PCR Product Size (bp)Region (bp)Annealing Temperature (°C)
Primer 1F: CACTTGCTCTCCCCTTCCTG
R: CAACGTTCAGCCTCTGGACT
6755791–675062
Primer 2F: ATGCCCTGTGCTGTCAACTT
R: TACAGTGGTCGCTGCCATAC
6427741–874060
Primer 3F: GGACAGGAGAGTTCTGTGGC
R: AAGGACAGAGCTGCTGGAAC
93918,181–19,18063
Primer 4F: AGAGTGCCTCACTGCATGTC
R: ACAGACTGGTGAACAGCCAC
68424,361–25,36059
Primer 5F: ACCGGATTGAGCCACAGTTT
R: GGCACCAAGATTTTCCCAGC
78928,681–29,68060
Primer 6F: ACTGAGGTCTCTCAGCCAGT
R: TACAGTGGTCGCTGCCATAC
95043,111–44,11059
Primer 7F: CTGTCTGACCCTGCTCTGTG
R: GCCGGTGAGAAACTCAGGAA
68079,641–80,64061
Primer 8F: TGCCATCTCCTCCCTACACA
R: GTCCACACCATGGCCACTAA
72185,681–86,68061
Primer 9F: CCAAGATGAGCACGGAAGGT
R: GGGATTGTGGTGGTACCCAG
64994,081–95,08061
Primer 10F: TCTCCAGACATCCCGTGTGA
R: GGTCTTGTGGGCATCCACTT
95296,061–97,06061
Primer 11F: GTTGAACGGATCTCCCCCAA
R: GCAACTGGAAACTGCCCTTC
530112,141–113,14063
Primer 12F: AAGCCCCTAACCTTGGTGTG
R: AGAGCACGTGCAGGCTAAAT
790114,001–115,00061
Abbreviations: F = Forward primer; R = Reverse primer.
Table 2. SNP analysis of the population genetic parameters of SORBS1.
Table 2. SNP analysis of the population genetic parameters of SORBS1.
LociGenomic RegionAmino AcidGenotypeGenotype Frequency (%)AlleleAllele Frequency (%)HoHeNePICp Value
g.6256 C > T
(c.298 C > T)
Exon 3Ser100ProCC38.3C61.90.5280.4721.8940.3610.054
CT47.1T38.1
TT14.6
g.24791 A > G
(c.706 A > G)
Exon 7Arg229HisAA38.3A61.90.5280.4721.8940.3610.054
AG47.1G38.1
GG14.6
g.29029 T > CIntron 7-TT38.3T61.90.5280.4721.8940.3610.054
TC47.1C38.1
CC14.6
g.29050 A > GIntron 7-AA38.3A61.90.5280.4721.8940.3610.054
AG47.1G38.1
GG14.6
g.29121 A > G
(c.979 A > G)
Exon 8Ala327ThrAA48.0A66.70.5560.4441.7990.3460.017
AG37.4G33.3
GG14.6
g.29245 C > TIntron 8-CC40.6C63.00.5340.4661.8730.3580.073
CT44.8T37.0
TT14.6
g.29305 C > TIntron 8-CC40.6C63.00.5340.4661.8730.3580.073
CT44.8T37.0
TT14.6
g.29347 T > CIntron 8-TT40.6T63.00.5340.4661.8730.3580.073
TC44.8C37.0
CC14.6
g.96284 G > A
(c.3090 G > A)
Exon 28Ser1030SerGG39.7G62.70.5320.4681.8790.3580.057
GA46.0A37.3
AA14.3
Abbreviations: Ho = Homozygosity; He = Heterozygosity; Ne = Effective number of alleles; PIC = Polymorphic information content. p > 0.05 suggested that the population gene is in Hardy–Weinberg balance and the sample comes from the same mendel population.
Table 3. Haplotype and haplotype frequencies of the three mutations in the SORBS1 gene.
Table 3. Haplotype and haplotype frequencies of the three mutations in the SORBS1 gene.
HaplotypePloymorphism Sites of SORBS1 GeneFrequency (%)
g.6256 C > T
(c.298 C > T)
g.24791 A > G
(c.706 A > G)
g.29121 A > G
(c.979 A > G)
H1CAA61.8
H2TGG33.3
H3TGA4.9
Table 4. Association analysis between genotypes and the milk fat of cattleyak.
Table 4. Association analysis between genotypes and the milk fat of cattleyak.
LociGenotypeMilk Fat Percentage (%)MUFAsPUFAsSFAs
g.6256 C > TCC5.649 ± 0.499 a1.605 ± 0.212 a0.180 ± 0.018 a3.771 ± 0.388 a
CT4.783 ± 0.732 b1.305 ± 0.177 b0.162 ± 0.023 b3.268 ± 0.479 b
TT3.280 ± 0.677 c1.033 ± 0.231 c0.130 ± 0.011 c1.834 ± 0.529 c
g.24791 A > GAA5.649 ± 0.499 a1.605 ± 0.212 a0.180 ± 0.018 a3.771 ± 0.388 a
AG4.783 ± 0.732 b1.305 ± 0.177 b0.162 ± 0.023 b3.268 ± 0.479 b
GG3.280 ± 0.677 c1.033 ± 0.231 c0.130 ± 0.011 c1.834 ± 0.529 c
g.29029 T > CTT5.649 ± 0.499 a1.605 ± 0.212 a0.180 ± 0.018 a3.771 ± 0.388 a
TC4.783 ± 0.732 b1.305 ± 0.177 b0.162 ± 0.023 b3.268 ± 0.479 b
CC3.280 ± 0.677 c1.033 ± 0.231 c0.130 ± 0.011 c1.834 ± 0.529 c
g.29050 A > GAA5.649 ± 0.499 a1.605 ± 0.212 a0.180 ± 0.018 a3.771 ± 0.388 a
AG4.783 ± 0.732 b1.305 ± 0.177 b0.162 ± 0.023 b3.268 ± 0.479 b
GG3.280 ± 0.677 c1.033 ± 0.231 c0.130 ± 0.011 c1.834 ± 0.529 c
g.29121 A > GAA5.516 ± 0.621 a1.568 ± 0.225 a0.178 ± 0.021 a3.662 ± 0.516 a
AG4.729 ± 0.715 b1.274 ± 0.152 b0.160 ± 0.020 b3.278 ± 0.401 b
GG3.280 ± 0.677 c1.033 ± 0.231 c0.130 ± 0.011 c1.834 ± 0.529 c
g.29245 C > TCC5.611 ± 0.588 a1.587 ± 0.231 a0.178 ± 0.021 a3.746 ± 0.465 a
CT4.773 ± 0.696 b1.305 ± 0.167 b0.163 ± 0.022 a3.625 ± 0.428 b
TT3.280 ± 0.677 c1.033 ± 0.231 c0.130 ± 0.011 c1.834 ± 0.529 c
g.29305 C > TCC5.611 ± 0.588 a1.587 ± 0.231 a0.178 ± 0.021 a3.746 ± 0.465 a
CT4.773 ± 0.696 b1.305 ± 0.167 b0.163 ± 0.022 a3.625 ± 0.428 b
TT3.280 ± 0.677 c1.033 ± 0.231 c0.130 ± 0.011 c1.834 ± 0.529 c
g.29347 T > CTT5.611 ± 0.588 a1.587 ± 0.231 a0.178 ± 0.021 a3.746 ± 0.465 a
TC4.773 ± 0.696 b1.305 ± 0.167 b0.163 ± 0.022 a3.625 ± 0.428 b
CC3.280 ± 0.677 c1.033 ± 0.231 c0.130 ± 0.011 c1.834 ± 0.529 c
g.96284 G > AGG4.335 ± 0.772 c1.297 ± 0.377 b0.152 ± 0.029 b2.829 ± 0.679 b
GA5.334 ± 0.612 a1.483 ± 0.167 a0.174 ± 0.021 a3.565 ± 0.467 a
AA5.040 ± 0.283b1.278 ± 0.088b0.168 ± 0.011 a3.417 ± 0.024 a
Abbreviations: MUFAs = Monounsaturated fatty acids; PUFAs = Polyunsaturated fatty acids; SFAs = Saturated fatty acids. a–c Values within a row with different superscripts differ significantly at p < 0.05.
Table 5. Association analysis of different genotype combinations of amino acid mutation sites with milk fat traits.
Table 5. Association analysis of different genotype combinations of amino acid mutation sites with milk fat traits.
DiplotypeCombinatorial GenotypeNumber of Amino Acid MutationsMilk Fat Percentage (%)MUFAsPUFAsSFAs
H1–H1CCAAAA35.649 ± 0.499 a1.605 ± 0.161 a0.180 ± 0.018 a3.771 ± 0.388 a
H1–H2CTAGAG0–34.783 ± 0.595 b1.305 ± 0.143 c0.162 ± 0.020 c3.268 ± 0.396 b
H1–H3CTAGAA1–34.992 ± 0.465 b1.421 ± 0.122 b0.171 ± 0.027 b3.321 ± 0.516 b
H2–H2TTGGGG03.280 ± 0.652 c1.033 ± 0.153 d0.130 ± 0.011 d1.834 ± 0.262 c
Abbreviations: MUFAs = Monounsaturated fatty acids; PUFAs = Polyunsaturated fatty acids; SFAs = Saturated fatty acids. a–d Values within a row with different superscripts differ significantly at p < 0.05.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, L.; Min, X.; Zhu, Y.; Hu, Y.; Yang, M.; Yu, H.; Li, J.; Xiong, X. Polymorphisms of SORBS1 Gene and Their Correlation with Milk Fat Traits of Cattleyak. Animals 2021, 11, 3461. https://doi.org/10.3390/ani11123461

AMA Style

Yang L, Min X, Zhu Y, Hu Y, Yang M, Yu H, Li J, Xiong X. Polymorphisms of SORBS1 Gene and Their Correlation with Milk Fat Traits of Cattleyak. Animals. 2021; 11(12):3461. https://doi.org/10.3390/ani11123461

Chicago/Turabian Style

Yang, Luyu, Xingyu Min, Yanjin Zhu, Yulei Hu, Manzhen Yang, Hailing Yu, Jian Li, and Xianrong Xiong. 2021. "Polymorphisms of SORBS1 Gene and Their Correlation with Milk Fat Traits of Cattleyak" Animals 11, no. 12: 3461. https://doi.org/10.3390/ani11123461

APA Style

Yang, L., Min, X., Zhu, Y., Hu, Y., Yang, M., Yu, H., Li, J., & Xiong, X. (2021). Polymorphisms of SORBS1 Gene and Their Correlation with Milk Fat Traits of Cattleyak. Animals, 11(12), 3461. https://doi.org/10.3390/ani11123461

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