Detection of CNV in the SH3RF2 gene and its effects on growth and carcass traits in chickens
BMC Genetics volume 21, Article number: 22 (2020)
The SH3RF2 gene is a protein-coding gene located in a quantitative trait locus associated with body weight, and its deletion has been shown to be positively associated with body weight in chickens.
In the present study, CNV in the SH3RF2 gene was detected in 4079 individuals from 17 populations, including the “Gushi ×Anka” F2 resource population and populations of Chinese native chickens, commercial layers, and commercial broilers. The F2 resource population was then used to investigate the genetic effects of the chicken SH3RF2 gene. The results showed that the local chickens and commercial layers were all homozygous for the wild-type allele. Deletion mutation individuals were detected in all of the commercial broiler breeds except Hubbard broiler. A total of, 798 individuals in the F2 resource group were used to analyze the effects of genotype (DD/ID/II) on chicken production traits. The results showed that CNV was associated with 2-, 6-, 10-, and 12-week body weight (P = 0.026, 0.042, 0.021 and 0.039 respectively) and significantly associated with 8-week breast bone length (P = 0.045). The mutation was significantly associated with 8-week body weight (P = 0.007) and 4-week breast bone length (P = 0.010). CNV was significantly associated with evisceration weight, leg muscle weight, carcass weight, breast muscle weight and gizzard weight (P = 0.032, 0.033, 0.045, 0.004 and 0.000, respectively).
CNV of the SH3RF2 gene contributed to variation in the growth and weight gain of chickens.
In recent years, with the development of the economy and the improvement of people’s living standards, consumer requirements regarding the quality of poultry products, especially flavor and taste, have increased. Chinese local chickens have excellent characteristics, such as tender meat, good taste and unique flavor, which are favored by consumers. However, as growth traits and carcass traits are the main economic traits of poultry, the slow growth rates and low feed utilization rates of local chicken breeds in China represent limits to production. Therefore, genetic improvement through activities such as cultivating new varieties to increase the growth rate and the rate of lean meat gain in chickens has been a focus of research . Thus, the use of modern molecular markers for marker-assisted selection and molecular breeding of chickens is important.
DNA molecular marker technology uses the gene library of the organism of interest without compromising the composition or expression of the genes . It is a kind of genetic marker technology that can be used to identify variation at the nucleic acid level that potentially reflects functional differences among individuals. Molecular marker technology is helpful in revealing differences in the composition or arrangement of the whole genome or variation at the nucleotide level within a gene, providing insight into DNA variability and polymorphism. It can also be used to identify individuals containing target genes by genotype analysis of closely linked genetic markers of target genes, which can help improve selection efficiency (e.g., by reducing blind search) and accelerate the breeding process . Compared with the use of traditional genetic markers, molecular marker-assisted selection (MAS) provides many marker loci, a large amount of genetic information, and strong repeatability of experiments, and it is not susceptible to environmental impacts and has no limitations regarding sex and age. Therefore, MAS allows early selection, shortens the generation interval, improves the selection intensity, and thus improves the efficiency and accuracy of selection. Due to these advantages, MAS has broad application prospects for animal genetic improvement . At present, the application of DNA molecular marker technology in poultry genetic breeding mainly comprises genetic diversity analysis, germplasm identification, genetic relationship research, genetic map construction, quantitative trait loci (QTL) mapping, genome wide association study (GWAS) and molecular marker-assisted breeding . A large number of genetic polymorphisms, including single nucleotide polymorphisms (SNPs), insertions/deletions (indels) and copy number variation (CNV), have been revealed in many species through whole-genome sequencing [6, 7]. CNV is an important source of genetic variation . CNV is the main form of genome structural variation, which refers to the insertion, deletion, duplication, translocation and derived chromosome structural variation of DNA fragments larger than 1 kb in the genome relative to the reference sequence of the genome . Because many CNVs contain entire genes, they are more difficult to identify and type than SNPs and indel copy number variants. As a result, they affect organisms to a greater extent than do these other types of polymorphisms. CNV is an important source of genetic variation complementary to SNP. CNV is associated with not only disease and abnormal development in livestock and poultry but also physical appearance and many economic traits [10,11,12,13]. For example, Wright  found that the first intron CNV of the SOX5 gene is related to crown type. The bean crown mutation in the chicken is a dominant mutation, which greatly reduces the size of the bean crown, thereby reducing heat loss and preventing frostbite, and it is an adaptive characteristic of the chicken in a cold environment . Dorshorst  found that the insertion of the endothelin 3 (EDN3) gene was the main cause of overstaining of the skin of the black chicken. Elferink  found that a segment comprising a chicken prolactin receptor in the Z-chromosome is associated with the growth of the chickenundefineds fast-and-slow feather. Gorla  provided that a chicken CNV map based on the 600 K SNP chip array jointly with a genome-wide gene copy number estimates in a native chicken population. Other reports have shown two copy number variable genes associated with Marek’s disease, namely frizzled family receptor 6 (FZD6) and LIM and senescent cell antigen-like domain 1 (LIMS1) [18, 19]. These findings indicate that CNV in a gene can have important influences on economic characteristics of poultry and enrich organismal diversity. To date, few studies have investigated the effects of CNV in the SH3 domain containing ring finger 2 (SH3RF2) gene on poultry growth and development. Rubin  resequenced the chicken genome and found a selective clearance region in many chicken genomes that included the SH3RF2 gene mutation site, which is located on chromosome 13, with all exons except the first exon being deleted. The total deletion length was 18,961 bp, and the gene was located in a body-related quantitative trait locus (QTL) range. The results of the study showed that the deletion mutation in the gene was positively associated with chicken body weight and fixed in the high-growth line of the broiler, appearing at lower frequencies in low-growth lines and commercial broilers, indicating that this deletion is closely related to broiler growth and body weight. Zhao  studied the frequency distribution of the indel mutation of the SH3RF2 gene in 15 local chicken species in China and found no local chickens with the deletion mutation. In this study, we tested the CNV genotypes of local chicken breeds, commercial broilers, commercial laying hens and the “Gushi × Anka” chicken F2 resource population to identify SH3RF2 gene mutations in the population and to analyze the association between CNV in this gene and chicken growth performance. The relationships between the CNV and growth and carcass traits of F2 resources in “Gushi × Anka” chickens showed that the SH3RF2 gene can be used as an effective marker for chicken production traits.
After amplification, the “ins” (I) allele yielded only one 305 bp fragment and was determined to be a nondeletion allele; individuals homozygous for this allele were assigned the II genotype. The “del” (D) allele yielded only one 416 bp fragment and was determined to be a deletion allele; homozygous individuals were assigned the DD genotype. Individuals exhibiting both 305 bp and 416 bp product bands were considered heterozygous and assigned the ID genotype (Fig. 1).
Analysis of the genetic parameters of SH3RF2 gene mutations in different populations
The distribution of SH3RF2 gene CNV genotypes and the gene frequencies in 4079 individuals from 17 populations were analyzed in this study (Additional file 1: Table S1). As shown in Table 1, there was a significant difference in genotypic frequency among the F2 resource population, AA broilers, R308 and CB broilers. ID genotypes and II genotypes were the predominant genotypes. The frequency of allele I was greater than that of allele D, and the frequency of II genotypes was higher than the frequencies of ID genotypes and DD genotypes in each population. The frequency of the DD genotype was highest in the AA broiler group. The expected heterozygosity (He) and number of allele (Ne) values were 0–0.35 and 1.00–1.55, respectively. According to the standard polymorphism information and polymorphism information content (PIC), PIC> 0.50 represents high polymorphism, PIC 0.25–0.50 represents moderate polymorphism, and PIC< 0.25 represents low polymorphism. The F2 resource group and AA broiler group showed moderate polymorphism, whereas the remaining populations showed low polymorphism.
Association analysis of SH3RF2 genotype and growth traits in the F2 resource population
As shown in Table 2, the growth traits of the F2 generation were significantly associated with SH3RF2 genotype. Genotype was significantly associated with BW at 2, 6, 8, 10 and 12 weeks (P < 0.05), BBL at 8 weeks (P < 0.05), 8-week BW (P < 0.01) and 4-week BBL (P < 0.01). The growth trait values of the homozygous genotype DD were higher than those of the other two genotypes each week. The 8-week BW and 4-week BBL of individuals with the DD genotype were 92.3 and 97.2% higher, respectively, than those of individuals with the II genotype, indicating that the DD genotype is an economically dominant genotype. The growth trait values of the ID genotype were higher than those of the II genotype each week, indicating that allele D was the dominant allele.
Association analysis of SH3RF2 genotype and carcass traits in the F2 resource population
The association analysis showed that SH3RF2 genotype was significantly related to SEW, EW, SEP, EP, LMW, LWR, GWR and CW (P < 0.05) and MW, GW and BMWR (P < 0.01). Table 3 shows that the slaughter index of the different genotypes followed the order DD genotype > ID genotype > II genotype, which is consistent with the growth trait results.
Analysis of SH3RF2 gene polymorphism and the genetic population structure of different breeds of chickens
The genetic resources of poultry species in China are diverse and contain great genetic variation and selection potential . The vigorous development of animal husbandry will promote breeding efforts to improve the performance of chicken production . Therefore, it is necessary to study the genetic variation of different chicken breeds. In this study, the frequencies of the DD genotype and D allele were lower than the other genotypes and I allele, respectively, in the studied populations, and few individuals with deletions were found among the populations of local and commercial layers. Whereas moderate polymorphism was found in the F2 resource population and the AA broiler, the remaining breeds showed low polymorphism. According to the polymorphism information content standard, this result indicates that the F2 resource population and AA broilers have greater selection potential in breeding and growth development efforts. In general, allele frequency reflects genetic diversity between populations. We found that allele D was the dominant allele in each population, suggesting that the gene did not undergo the same selection in the evolution of the different breeds of chicken. Moreover, artificial selection has a significant impact on the number of genes and distribution of genetic variation in different varieties . We found that the frequency of DD genotypes was higher in commercial broiler populations with fast, large and high yields than that in other breeds, indicating that targeted breeding of commercial broiler chickens promoted the fixation of DD genotypes, which may be related to the rapid growth of commercial broiler chickens. We found no mutant deletion in commercial layers, which may be related to the selective breeding of layers.
Associations between the SH3RF2 gene mutation and growth and carcass traits
Body weight and body size are important indicators of the body development of poultry and are closely related to important economic characteristics . To date, few studies on the SH3RF2 gene in poultry production have been conducted. Rubin  analyzed 400 chickens from F8 generations of high-growth and low-growth lines and found that the SH3RF2 gene deletion mutation had highly significant effects on chicken growth traits (P < 0.01). In the present study, the F2 resource population of Anka and Gushi chickens was used to analyze the relationships between the SH3RF2 gene mutation and growth and carcass traits of chickens. Consistent with the previous study, the results showed that the mutation had significant influences on the growth and carcass traits of chickens. As poultry muscles are mainly distributed in the chest and legs, the yield of chest and leg muscles is an important factor determining the slaughtering performance of poultry and is one of the traits considered in poultry breeding . In this study, SH3RF2 genotype was significantly associated with growth traits, such as BW and BBL, and slaughter indicators, such as GW, GW, LMW and CW (P < 0.05). In addition, individuals with the DD genotype had higher phenotype values than did individuals with II genotypes or ID genotype. The results show that the DD genotype has a significant dominant effect on traits of the F2 resource population.
This study found that CNV in the SH3RF2 gene was related to most of the studied growth traits and carcass traits in chickens. The results show that the mutation can be used for MAS breeding in chickens. Molecular breeding of chickens can reduce breeding costs, shorten the generation gap, and improve the efficiency of breeding to provide technical support for local chicken breeding in China.
To identify the distributions of the different genotypes in different populations, individuals of the Gushi-Anka chicken F2 generation resource population, 8 local varieties, 5 chicken varieties and 2 layer varieties were investigated, with DNA samples collected from a total of 4079 individuals (Table 4). All blood samples were collected through the wing vein, after which 1:300 multidimensional hormone was used to reduce the stress response. Whole genome DNA was extracted from whole blood using the phenol-chloroform method. All chickens used in the study were healthy animals raised in the same environment with ad libitum access to feed and water. To produce the F2 resource population, two hatchings were obtained from an F1 generation constructed via reciprocal crossing between Chinese native Gushi chickens (representing a slow-growing Chinese native chicken) and Anka broilers (representing a fast-growing broiler). Then according to a 1:9 ratio of males to females, an F2 generation was produced by mating the offspring with other family hens. The hens represented 7 families. The F2 generation was thus composed of 7 families, with Anka chicken as the male parent of 4 orthogonal lines and Gushi chicken as the male parent of the 3 reverse cross lines. We used an F2 resource family as previously described by Liang et al.  and Li et al. .
In henan provincial poultry germplasm resources innovation engineering research center, Gushi-Anka F2 resource population were housed in cages with 50 chickens per cage at 1 d, providing 392 cm2 per bird. At the age of 8-week, it was transferred to three birds per cage (448 cm2 / bird). A total of 798 individuals from the F2 resource population were euthanised at 84 d. In the laboratory of Henan Agricultural University, 5% Pentobarbital 2 mL (No. 57–33-0 of Chinese Academy of Sciences, Beijing Siyuan Technology Co., Ltd.) was injected intraperitoneally into the Gushi-Anka F2 resource population. After 2–3 min of no spontaneous respiration, carotid artery bleeding occurred. Intraperitoneal injection of pentobarbital sodium can lead to rapid and relatively stress-free death. From hatch to slaughter, several chicken growth traits, including body weight and body size indexes, were measured. Each chicken was weighed every 2 weeks. Shank length was measured at 0, 4, 8 and 12 weeks, and shank girth, chest depth, breast bone length, body slanting length and pelvis breadth were determined at 4, 8 and 12 weeks. Descriptions of the construction of the F2 population, feeding management and trait determination procedures can be found in a previous study [29, 30]. Information on 17 populations is presented in Table 4.
Primer and PCR amplification
The primers used in this study were all designed by Rubin  (Forward: 5′-TGCTTCGGGCTGAGCCTTCT-3′, Reverse1: 5′-CGCCCAAGCTGTGTCCT-3′, Reverse2: 5′-CTGTCGGGCACGTGAGTGAA-3′). Assays were performed by PCR in a total volume of 10 μL containing 5 μL of 2 × Taq Master Mix (Kangwei, Beijing, China), 0.5 μL of forward primer, 0.5 μL of reverse primer, 2.5 μL of ultrapure water and 1 μL of genomic DNA. The PCR amplification was performed as follows: 95 °C for 5 min followed by 30 cycles at 95 °C for 30 s, 64 °C for 30 s, and 72 °C for 10 min. The samples were then chilled at 4 °C. To determine genotype, aliquots from each reaction (7 μL) were subjected to electrophoresis in a 2% agarose gel.
Statistical analyses of the associations between genotype and the selected traits of the F2 chickens were performed using IBM SPSS (SPSS for Windows, Standard version 24; SPSS, USA). Genotype effects were analyzed by a multivariate linear model, and differences among genotypes were evaluated by Bonferroni’s multiple comparison method. The analyses followed a previous study . All data for each trait obtained by statistical analysis are presented as the mean ± standard error (mean ± SE). P < 0.05 was considered to indicate statistical significance .
Availability of data and materials
All the data supporting the conclusions of the study are included in the manuscript and Additional file 1.
Arbor Acres broiler
Breast muscle weight rate
Body slanting length
Breast muscle water loss rate
Copy number variation
F2 resource population
Frizzled family receptor 6
Genome wide association study
Gizzard weight rate
Hyline brown hen
Head weight percentage
Lohmann brown laying hen
Leg muscle weight
Liver weight rate
Breast muscle weight
Number of allele
Polymerase chain reaction
Polymorphism information content
Quantitative trait loci
Recessive white chicken
Semievisceration weight rate
SH3 domain containing ring finger 2
Single nucleotide polymorphisms
Tallentire CW, Leinonen I, Kyriazakis I. Breeding for efficiency in the broiler chicken: Areview. Agron Sustain Dev. 2016;36(4):66.
Vignal A, Milan D, Sancristobal M. André Eggen. A review on SNP and other types of molecular markers and their use in animal genetics. Genet Sel Evol. 2002;34(3):275–81.
Liu Z, Yang N, Yan YY, Li GQ, Liu AQ, Wu GQ, Sun CJ. Genome-wide association analysis of egg production performance in chickens across the whole laying period. BMC Genet. 2019;20(1):67.
Mahmoud KGM, Nawito MF. Molecular markers for fertility in farm animals. Iran J Appl Animal Sci. 2012;20(2):203–22.
Grover A, Sharma PC. Development and use of molecular markers: past and present. Crit Rev Biotechnol. 2014;36(2):1–13.
Mikkelsen TS, Hillier LW, Eichler EE, Zody MC, Jaffe DB, Yang SP, Enard W, Hellmann I, Lindblad-Toh K, Altheide TK, et al. Initial sequence of the chimpanzee genome and comparison with the human genome. Nature. 2005;437:69–87.
Begun DJ, Holloway AK, Stevens K, Hillier LW, Poh YP, Hahn MW, Nista PM, Jones CD, Kern AD, Dewey CN, et al. Whole-genome analysis of polymorphism and adivergence in Drosophila simulans. PLoS Biol. 2007;5:e3l0.
Conrad DF, Pinto D, Redon R, Feuk L, Gokcumen O, Zhang Y, Aerts J, Andrews TD, Barnes C, Campbell P, et al. Origins and functional impact of copy number variation in the human genome. Nature. 2010;464(7289):704–12.
Feuk L, Carson AR, Scherer SW. Structural variation in the human genome. Nat Rev Genet. 2006;7(2):85–97.
Gao YH, Jiang JP, Yang SH, Hou YL, Liu GE, Zhang SL, Zhang Q, Sun DX. CNV discovery for milk composition traits in dairy cattle using whole genome resequencing. BMC Genomics. 2017;18(1):265.
Moro C, Cornette R, Vieaud A, Vieaud A, Bruneau N, Gourichon D, Bed’hom B, Boichard M. Quantitative effect of a CNV on a morphological trait in chickens. PLoS One. 2015;10(3):e0118706.
Han RL, Yang PK, Tian YD, Wang DD, Zhang ZX, Wang LL, Li ZJ, Jiang RR, Kang XT. Identification and functional characterization of copy number variations in diverse chicken breeds. BMC Genomics. 2014;15(1):934.
Zheng L, Xu JW, Li JC, Wang DH, An QM, Xu LN, Ma YL, Wang J, Peng SJ, Lei CZ, et al. Distribution and association study in copy number variation of KCNJ12 gene across four Chinese cattle populations. Gene. 2019;689:90–6.
Wright D, Bije H, Meadows JRS, Bedhome B, Gourichon D, Vieaud A, Tixier-Boichard M, Rubin CJ, Imsland F, Hallböök F, et al. Copy number variation in intron 1 of SOX5 causes the pea-comb phenotype in chickens. PLoS Genet. 2009;5(6):e1000512.
Dorshorst B, Molin AM, Rubin CJ, Johansson AM, Strömstedt L, Pham MH, Chen CF, Hallböök F, Ashwell C, Andersson L. A complex genomic rearrangement involving the endothelin 3 locus causes dermal hyperpigmentation in the chicken. PLoS Genet. 2011;7(12):e1002412.
Elferink MG, Vallée AA, Jungerius AP, Crooijmans RP, Groenen MA. Partial duplication of the PRLR and SPEF2 genes at the late feathering locus in chicken. BMC Genomics. 2008;9(1):391.
Gorla E, Cozzi MC, Román-Ponce SI, Ruiz López FJ, Vega-Murillo VE, Cerolini S, Bagnato A, Strillacci MG. Genomic variability in Mexican chicken population using copy number variants. BMC Genet. 2017;18(1):61.
Luo J, Yu Y, Mitra A, Chang S, Zhang M, Liu G, Yang N. Genome-wide copy number variant analysis in inbred chickens lines with different susceptibility to Marek’s disease. G3 (Bethesda, Md.). 2013;3(2):217–23.
Yi G, Qu L, Liu J, Yan Y, Xu G, Yang N. Genome-wide patterns of copy number variation in the diversified chicken genomes using next-generation sequencing. BMC Genomics. 2014;15(1):962.
Rubin CJ, Zody MC, Eriksson J, Meadows JR, Sherwood E, Webster MT, Jiang L, Ingman M, Sharpe T, Ka S, et al. Whole-genome resequencing reveals loci under selection during chicken domestication. Nature. 2010;464(7288):587–93.
Zhao YJ, Zhang DX, Chen Y, Jia YX, Qu LJ. Study on the allele frequency distribution of an Indel of SH3RF2 gene in Chinese indigenous chicken breeds. Chinese Poultry. 2012;34(06):33–6.
Ma YL, Gu LT, Yang LB, Sun CH, Xie SS, Fang CC, Gong YZ, Li SJ. Identifying artificial selection signals in the chicken genome. PLoS One. 2018;13(4):e0196215.
Han RL, Li ZJ, Li MJ, Li JQ, Lan XY, Sun GR, Kang XT, Chen H. Novel 9-bp indel in visfatin gene and its associations with chicken growth. Br Poult Sci. 2011;52(1):52–67.
Han RL, Wang XN, Wang XL, Guo YP, Li DH, Li GX, Wang YB, Kang XT, Li ZJ. Chicken ZNF764L gene: mRNA expression profile, alternative splicing analysis and association analysis between first exon indel mutation and economic traits. Gene. 2019;695:92–8.
Allais S, Hennequet-Antier C, Berri C, Salles L, Demeure O, Le Bihan-Duval E. Mapping of QTL for chicken body weight, carcass composition, and meat quality traits in a slow-growing line. Poult Sci. 2019;98(5):1960–7.
Gaya LG, Ferraz JBS, Rezende FM. G B Mourão, Filho TM. Heritability and genetic correlation estimates for performance and carcass and body composition traits in a male broiler line. Poult Sci. 2006;85(5):837–43.
Liang K, Wang XN, Tian XX, Geng R, Li WY, Jing ZZ, Han RL, Tian YD, Liu XJ, Kang XT, et al. Molecular characterization and an 80-bp indel polymorphism within the prolactin receptor (PRLR) gene and its associations with chicken growth and carcass traits. 3. Biotech. 2019;9(8):296–306.
Li W, Liu D, Tang S, Li D, Han R, Tian Y, Li H, Li G, Li W, Liu X, Kang X, Li Z. A multiallelic indel in the promoter region of the Cyclin-dependent kinase inhibitor 3 gene is significantly associated with body weight and carcass traits in chickens. Poult Sci. 2019;98:556–65.
Wang XN, Li ZJ, Guo YP, Wang YB, Sun GR, Jiang RR, Kang XT, Han RL. Identification of a novel 43-bp indel in the heparan sulfate 6-O-sulfotransferase 3 (HS6ST3) gene and its associations with growth and carcass traits in chickens. Anim Biotechnol. 2019;30(3):252–9.
Wang XN, Wang XL, Chen BJ, Guo YP, Tang HH, Li DH, Liu DL, Wang YB, Li GX, Kang XT, et al. Association of a new 99-bp indel of the CEL gene promoter region with phenotypic traits in chickens. Sci Rep. 2020;10(1):1–10.
Ren TH, Li WY, Liu DL, Liang K, Wang XN, Li H, Jiang RR, Tian YD, Kang XT, Li ZJ. Two insertion/deletion variants in the promoter region of the QPCTL gene are significantly associated with body weight and carcass traits in chickens. Anim Genet. 2019;50(3):279–82.
Liu D, Han R, Wang X, Li W, Tang S, Li W, Wang Y, Jiang R, Yan F, Wang C, et al. A novel 86-bp indel of the motilin receptor gene is significantly associated with growth and carcass traits in Gushi-Anka F2 reciprocal cross chickens. Br Poult Sci. 2019;60(6):649–58.
We wish to thank our lab members for their helpful discussions.
This study was funded by resources from the Key Science and Technology Research Project of Henan Province (202102110085), National Natural Science Foundation of China-Henan joint grant (U1804107 and U1704233), National Natural Science Foundation of China (31872987), and Earmarked Fund for Modern Agro-Industry Technology Research Systems of China (No. CARS-40-K04). The funding bodies had no influence on the design, acquisition, analysis, interpretation, nor presentation of the data and results of this study.
Ethics approval and consent to participate
All animal experiments and animal care methods were approved by the Institutional Animal Care and Use Committee (IACUC) of Henan Agricultural University, Zhengzhou, P.R. China (Permit Number: 11–0085) and were performed in accordance with the protocols outlined in the “Guide for Care and Use of Laboratory Animals” (Henan Agricultural University).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Jing, Z., Wang, X., Cheng, Y. et al. Detection of CNV in the SH3RF2 gene and its effects on growth and carcass traits in chickens. BMC Genet 21, 22 (2020). https://doi.org/10.1186/s12863-020-0831-z
- SH3RF2 gene
- Growth traits
- Carcass traits
- Association analysis