Development of a touch-down multiplex PCR method for simultaneously rapidly detecting three novel insertion/deletions (indels) within one gene: an example for goat GHR gene

Detection of small sequence variations within the goat GHR gene and its effects on growth traits

Growth hormone receptor (GHR) gene is considered to be an important candidate gene in growth traits. Therefore, the purpose of this study was to detect whether there were potential indel variations in the GHR gene that were related to the growth traits of the Shaanbei white cashmere goats (SBWC). In this study, genomic DNA from 931 healthy SBWC individuals were used to verify the relationship between the indel of the GHR gene and growth traits. Two indel variants, P49-bp indel in intron 1 and P1410-bp indel in 3'-UTR, were confirmed. Association analyses demonstrated that these two indel polymorphism loci were associated with the chest circumference and chest width of SBWC. Additionally, for the P49-bp and P1410-bp indel loci, the ID and II genotypes were dominant genotypes, respectively. Moreover, the genotypic distributions of these two indel loci in SBWC were significantly different from those in three other Chinese indigenous goat breeds (HNBG, GZDG and IMWC) (p < 0.05). Taken together, two indel loci (P49-bp indel and P1410-bp indel) both significantly affected the growth traits of goats. This illustrated that these two indel loci might be the potential DNA marker for use in improving the selection and breeding of goats.

Expression profile of FASN gene and association of its polymorphisms with intramuscular fat content in Hu sheep

The content of intramuscular fat (IMF) is one of the most important factors that has a large impact on meat quality, and it is an effective way to improve IMF according to marker-assisted selection (MAS). Fatty-acid synthase (FASN) is a key gene in meat lipid deposition and fatty acid composition. Thus, this study was conducted to investigate the expression profile of FASN in mRNA and protein levels using real-time quantitative PCR (RT-qPCR) and western-blot methods. In addition, single nucleotide polymorphisms (SNPs) within FASN in 921 Hu rams with IMF content records were investigated using DNA-pooling sequencing and improved multiple ligase detection reaction (iMLDR) methods. Consequently, the highest mRNA expression level of FASN was observed in the perinephric fat, and the lowest in the liver among the 11 tissues analyzed, while no significant difference was found in mRNA and protein expression levels in longissimus dorsi among individuals with different IMF contents. A total of 10 putative SNPs were identified within FASN, and 9 of them can be genotyped by iMLDR method. Notably, two SNPs were significantly associated with IMF content, including NC_040262.1: g.5157 A > G in intron 5 (p = 0.046) and NC_040262.1: g.9413 T > C in intron 16 (p = 0.041), which supply molecular markers for improving meat quality in sheep breeding.

Insertion/deletion (Indel) variant of the goat RORA gene is associated with growth traits

RAR related orphan receptor A (RORA), which encodes the retinoid-acid-related orphan receptor alpha (RORα), is a clock gene found in skeletal muscle. Several studies have shown that RORα plays an important role in bone formation, suggesting that RORA gene may take part in the regulation of growth and development. The purpose of this research is to study the insertion/deletion (indel) variations of the RORA gene and investigate the relationship with the growth traits of Shaanbei white cashmere (SBWC) goats. Herein, the current study identified that the P4-11-bp and P11-28-bp deletion sites are polymorphic among 12 pairs of primers within the RORA gene in the SBWC goats (n = 641). Moreover, the P11-28-bp deletion locus was significantly related to the body height (p = 0.046), height at hip cross (p = 0.012), and body length (p = 0.003). Both of P4-11-bp and P11-28-bp indels showed the moderate genetic diversity (0.25 Collapse

Key Words

• RAR related orphan receptor A (RORA)
• clock gene
• goat
• growth traits
• insertion/deletion

Collapse

MESH Headings

• Pregnancy
• Female
• Animals
• Litter Size/genetics
• Goats/physiology
• INDEL Mutation/genetics
• Phenotype

Collapse

Grants Collapse

Affiliation(s)

Qian Zhou College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Huina Hu College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Yuta Yang College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Yuxin Kang College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Xianyong Lan College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Xianfeng Wu Institute of Animal Husbandry and Veterinary, Fujian Academy of Agricultural Sciences, Fuzhou, Fujian, China
Zhengang Guo College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China Animal Husbandry and Veterinary Science Institute of Bijie city, Bijie, Guizhou, China
Chuanying Pan College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China

Collapse

Investigation of the association between goat DNMT3B gene polymorphism and growth traits

The DNA methyltransferase 3 beta (DNMT3B) gene is key for DNA methylation and has been well recognized in regulating growth and development. A previous observation indicated that an 11-bp indel of DNMT3B affected the reproductive traits in goats, yet the effect of this polymorphism on body measurement traits in goats has not been reported. This study aims to investigate the associations between DNMT3B gene polymorphism and goat growth traits. We investigated this 11-bp indel in 2184 goats and three genotypes have been found in Shaanbei white cashmere goat (SBWC): insertion/insertion (II), deletion/deletion (DD) and insertion/deletion (ID). Only ID and DD genotypes were detected in Nubian goats and Guizhou heima goat (GZHM). The allele frequencies analyzed revealed that the 'D' allele frequencies were higher in all three goat breeds. Further association analysis demonstrated that this indel is markedly associated with the cannon circumference (CC) and cannon circumference index (CCI) of SBWC and cannon circumference (CC) of Nubian goats (p < .05). The CC and CCI are essential indicators to measure the growth status of goats. In summary, our study sheds some light on the potential impact of the 11-bp indel polymorphism of the DNMT3B gene on improving the growth traits in goats.

Detection distribution of CNVs of SNX29 in three goat breeds and their associations with growth traits

As a member of the SNX family, the goat sorting nexin 29 (SNX29) is initially identified as a myogenesis gene. Therefore, this study aimed to examine the polymorphism in the SNX29 gene and its association with growth traits. In this study, we used an online platform to predict the structures of the SNX29 protein and used quantitative real-time PCR to detect potential copy number variation (CNV) in Shaanbei white cashmere (SBWC) goats (n = 541), Guizhou black (GB) goats (n = 48), and Nubian (NB) goats (n = 39). The results showed that goat SNX29 protein belonged to non-secretory protein. Then, five CNVs were detected, and their association with growth traits was analyzed. In SBWC goats, CNV1, CNV3, CNV4, and CNV5 were associated with chest width and body length (P < 0.05). Among them, the CNV1 individuals with gain and loss genotypes were superior to those individuals with a median genotype, but CNV4 and CNV5 of individuals with the median genotype were superior to those with the loss and gain genotypes. In addition, individuals with the gain genotype had superior growth traits in CNV3. In brief, this study suggests that the CNV of SNX29 can be used as a molecular marker in goat breeding.

Goat Pleomorphic Adenoma Gene 1 (PLAG1): mRNA Expression, CNV Detection and Associations with Growth Traits

The pleomorphic adenoma gene 1 (PLAG1) gene, as the major gene responsible for growth, plays a vital role in myogenesis. Meanwhile, the relationship between copy number variation (CNV) of this gene and growth traits in goats remains unclear. Therefore, this study investigated four aspects: bioinformatics analysis, mRNA expression (n = 6), CNV detection (n = 224), and association analysis. The findings indicated that the gene had a large number of conserved motifs, and the gene expression level was higher in fetal goats than in adult goats. Three CNV loci were selected from the database, among which CNV1 was located in the bidirectional promoter region and was associated with goat growth traits. CNV analysis showed that CNV2 and CNV3 of the PLAG1 gene were associated with growth traits such as body weight, heart girth, height at hip cross, and hip width (p < 0.05), with CNV1 loss genotype being the superior genotype, and CNV2 and CNV3 median and gain genotypes of being superior genotypes. This finding further confirms that the PLAG1 gene is the dominant gene for growth traits, which will serve as theoretical guidance for goat breeding.

A 24-bp indel within the sheep AHR gene is associated with litter size

Aryl Hydrocarbon Receptor (AHR) is a member of the PER-ARNT-SIM (PAS) family, which could mediate various biological processes, for instance, the balance of the immune system, cell proliferation, differentiation, vascular tissue remodeling and reproduction ability regulation. A previous research showed that the AHR gene exerted important functions on the pig reproduction, implying that it could serve as a candidate gene related to animal reproductive traits. Here, the aim of this work was to identify potential insertion/deletion (indel) mutations of the AHR gene in three sheep breeds and analyze the associations between these mutations and reproductive traits. Results showed that a 24-bp indel was uncovered three genotypes (II, ID and DD) in the Australian White sheep (AuW) and Lanzhou fat-tail sheep (LZFT) population, while there were only two genotypes (ID and DD) in Luxi black-headed sheep (LXBH). Moreover, the Fisher's exact test showed that the 24-bp indel mutation was significantly associated with litter size and live litter size in AuW sheep (Fisher's p < 0.05). Therefore, the 24-bp indel of sheep AHR gene can contribute to sheep marker-assisted selection breeding and further improve the sheep reproductive performance.

The detection of mutation within goat cell division cycle 25 A and its effect on kidding number

Cell division cycle 25 A (CDC25A) accounts for an essential function on early folliculogenesis of female mammals, especially regulating the function of intra-ovarian, thus this gene is pinpointed as a candidate gene that influences the kidding number of goat. On this ground, the purpose of this study was to investigate whether the reported 20-nt nucleotide variants locus (rs639467625) of the CDC25A gene influences kidding number in Shaanbei white cashmere goat (SBWC). The χ²-test showed that there were more ID genotypes in mothers of multiple lambs than in mothers of single lambs. Interestingly, this indel locus was related to the first-born kidding number in the group of SBWC goats (p < 0.05). Similarly, the result of the t-test was consistent with the result of the χ²-test, showed the kidding number of ID genotype individuals was large than that of II individuals (p < 0.05). These findings proved that the different genotypes of CDC25A have impacts on goat kidding numbers. Thus, the results led us to speculate that the ID genotype of CDC25A was one of the main indel influencing goat kidding numbers. Simultaneously, this study was expected to provide useful DNA markers for superior individuals selection by marker-assisted selection (MAS) and make a contribution to goats breeding.

A 7-nt nucleotide sequence variant within the sheep KDM3B gene affects female reproduction traits

Lysine demethylase 3B (KDM3B) gene is a histone demethylase, demonstrating specific demethylation of the histone H3 lysine 9. It was detected as a sheep reproductive candidate gene by genome-wide scans, and related studies also showed its significance in female reproductive process. However, rare study researched its polymorphism. Herein, we hypothesized that the polymorphisms of KDM3B gene were associated with sheep reproduction traits. A 7-nt nucleotide sequence variant (rs1088697156) within KDM3B gene was identified in a total of 888 individuals, including the Australian White (AUW) sheep and Lanzhou Fat-tailed (LFT) sheep. II (insertion/insertion) and ID (insertion/deletion) genotypes of 7-nt variant were detected, which were at Hardy-Weinberg equilibrium (HWE) in detected breeds. Association analysis illustrated the 7-nt variant was significantly associated with the litter size, duration of pregnancy, live lamb number, live lamb rate, stillbirth number, stillbirth rate of average and different parity (P < 0.05) in AUW sheep. Moreover, 'ID' was the dominant genotype with excellent consistency in reproductive traits. It is instrumental to select individuals with ID genotype for improving the sheep reproduction traits. These findings suggest that the 7-nt variant within KDM3B gene can be used as a candidate marker of reproduction traits for sheep breeding improvement by marker-assisted selection.

Short variation of the sheep PDGFD gene is correlated with litter size

Platelet-derived growth factor (PDGF) family, exert plays a key role in embryonic development, cell proliferation, cell migration, angiogenesis and reproduction. Related studies about GWAS analyses have found that PDGFD significantly affected deposition of tail fat in sheep, but there are no studies on reproduction in animals. In this study, three breed of sheep were used to find insertion/deletion (indel) fragment polymorphism of PDGFD which including Australian white (AUW) sheep (Meat type, n = 932), Guiqian semi-fine wool (GSFW) sheep (wool type, n = 60) and East Friensian milk (EFM) sheep (dairy type, n = 60). Only a 18-bp variation was polymorphic in the study AUW sheep population and the genotypes of different sheep breed are also specific. Moreover, the association analysis indicated that this variant was associated with litter size of AUW sheep in the first parity (p < 0.05). The litter size of II genotype was significantly lower than other genotypes in the first parity (p < 0.05). We also revealed that the PDGFD gene was relatively conservative in eight species, PDGFD mRNA expression in 832 sheep samples implying this gene was related to reproduction traits. Hence, these finding demonstrated the one-cause multipotency of PDGFD gene. Collectively, these results suggest that this indel can be used as an effective marker for sheep breeding.

Goat SNX29: mRNA expression, InDel and CNV detection, and their associations with litter size

The sorting nexin 29 (SNX29) gene, a member of the SNX family, is associated with material transport and lipid metabolism. Previous studies have shown that lipid metabolism affects reproductive function in animals. Thus, we hypothesized there is a correlation between the SNX29 gene and reproductive trait. To date, studies on the relationship between the SNX29 gene and reproductive traits are limited. Therefore, the purpose of this study was to examine the polymorphism in the SNX29 gene and its correlation with litter size. Herein, the mRNA expression levels of SNX29 were assayed in various goat tissue. Surprisingly, we found that SNX29 was highly expressed in the corpus luteum, large and small follicles. This result led us to suggest that the SNX29 gene has a critical role in reproduction. We further detected potential polymorphisms in Shaanbei white cashmere (SBWC) goats, including insertion/deletion (InDel, n = 2,057) and copy number variation (CNV, n = 1,402), which were related to fertility. The 17 bp deletion (n = 1004) and the 20 bp deletion (n = 1,053) within the SNX29 gene were discovered to be significantly associated with litter size (P < 0.05), and individuals the ID genotype of P1-Del-17 bp and the DD genotype of P2-Del-20bp had larger litter size. Additionally, the four CNV loci had significant correlations with litter size (P < 0.01) in our detected population. In CNV5, individuals with the median genotype were superior compared to those with loss or gain genotype in term of litter size, and in other three CNVs showed better reproductive trait in the gain genotype. Briefly, these findings suggest that SNX29 could be used as a candidate gene for litter size in goat breeding through marker-assisted selection (MAS).

Rapid detection of InDel within the KDM3B gene in five sheep breeds using the mathematical expectation (ME) method

Lysine demethylase 3B (KDM3B), a candidate gene associated with bone formation and growth, and differentiation of osteoblast, might affect the animal growth traits. Herein, the insertion/deletion (InDel) of the KDM3B gene was quickly detected in 882 sheep from five breeds using the mathematical expectation (ME) method. The results showed that there were two genotypes of 7-bp variation in KDM3B, including II (insertion/insertion) and ID (insertion/deletion), and the frequency of two genotypes varied among the five sheep breeds. Association analysis results demonstrated that the 7-bp indel was significantly associated with chest depth of LFT sheep (P = 0.012), and body weight (P = 0.006), body height (P = 0.030), chest depth (P = 0.043), chest circumference (P = 0.016), abdominal width (P = 0.035) and height at hip cross (P = 0.022) in LXBH sheep. Moreover, II genotype was the predominant genotype with excellent consistency in sheep growth traits (P < 0.05). Collectively, the above results suggest that this locus can be used as an effective molecular marker to improve the sheep growth traits and provide a scientific basis for the development of sheep breeding.

Screen of small fragment mutations within the sheep thyroid stimulating hormone receptor gene associated with litter size

The thyroid stimulating hormone receptor (TSHR), a glycoprotein hormone receptor, plays an important role in metabolic regulation and photoperiod control in the time of reproduction in birds and mammals. Previous genome-wide association studies revealed that the TSHR gene was related to reproduction and its function was identified in female reproduction, but rare studies reported the polymorphism of TSHR gene. However, the molecular mutations of the TSHR gene in sheep have not been reported so far. Herein, we explored potential polymorphisms of the sheep TSHR gene, and a 29 bp nucleotide sequence variant (rs1089565492) was identified in the AUW sheep. There were three genotypes of the 29 bp variant locus detected which named 'II' 'DD' and 'ID' been identified. Association analysis results showed the 29 bp variant was significantly associated with the litter size of the AUW sheep (p < 0.05). This finding suggests that the 29 bp nucleotide sequence variant within TSHR gene could be a candidate marker of reproduction traits for sheep breeding improving through the marker-assisted selection (MAS).

Exploration of Genetic Variants within the Goat A-Kinase Anchoring Protein 12 (AKAP12) Gene and Their Effects on Growth Traits

Simple Summary

AKAP12, the family of A-kinase anchoring proteins (AKAPs), plays an important role in the regulation of growth and development. There have been no corresponding studies of the effect of the AKAP12 gene on growth traits in goats. In our previous study, 7 bp (intron 3) and 13 bp (3′UTR) indels within the AKAP12 gene significantly influenced AKAP12 gene expression. This study expected to identify the association between these two genetic variations and growth-related traits in 1405 Shaanbei white cashmere (SBWC) goats. The P1–7 bp indel locus was significantly correlated with height at hip cross (HHC; p < 0.05) and the P2–13 bp indel locus was associated with body weight, body length, chest depth, chest width, hip width, chest circumference and cannon (bone) circumference in SBWC goats (p < 0.05). These results prove that the AKAP12 gene plays an important role in the growth and development of goats.

Abstract

The A-kinase anchoring protein 12 gene (AKAP12) is a scaffold protein, which can target multiple signal transduction effectors, can promote mitosis and cytokinesis and plays an important role in the regulation of growth and development. In our previous study, P1–7 bp (intron 3) and P2–13 bp (3′UTR) indels within the AKAP12 gene significantly influenced AKAP12 gene expression. Therefore, this study aimed to identify the association between these two genetic variations and growth-related traits in Shaanbei white cashmere goats (SBWC) (n = 1405). Herein, we identified two non-linkage insertions/deletions (indels). Notably, we found that the P1–7 bp indel mutation was related to the height at hip cross (HHC; p < 0.05) and the P2–13 bp indel was associated with body weight, body length, chest depth, chest width, hip width, chest circumference and cannon (bone) circumference in SBWC goats (p < 0.05). Overall, the two indels’ mutations of AKAP12 affected growth traits in goats. Compared to the P1–7 bp indel, the P2–13 bp indel is more suitable for the breeding of goat growth traits.

The InDel variants of sheep IGF2BP1 gene are associated with growth traits

Insulin-like growth factor 2 mRNA binding protein 1 (IGF2BP1) plays positive roles in the growth, proliferation of cells and early embryos development by binding mRNA targets. Recently, it had been shown that some polymorphic loci within IGF2BP1 gene were associated with growth traits in animals, especially in goats. Therefore, it has been hypothesized that some variants within IGF2BP1 gene may be also involved in growth traits of sheep. Nine insertion/deletion (InDel) mutations within IGF2BP1 were identified and three loci were polymorphic. Meanwhile, the association analyses between three InDels and growth traits were carried out in 745 sheep. The results showed that all InDels included 5 bp InDel in downstream region, 9 bp InDel in intron 4 and 15 bp InDel in intron 2 within IGF2BP1 were significantly associated with growth traits (p<.05). Furthermore, at 5 and 9 bp InDel loci, the individuals of heterozygous genotype (ID) had superior growing performance especially at body weight (BW). In all, three InDels were crucial variants correlated with growth traits and could be applied in marker-assisted selection (MAS) in sheep.

Identification of three new microsatellites and their effects on body measurement traits in pigs using time of flight-mass spectrometry (TOF-MS)

The body status of livestock affects their physiological function and productive performances. Microsatellites, one of the most used DNA markers, have been found to be associated with pig productive traits. However, their identifications and effects on body measurement traits of the Chinese Qinghai Bamei pig still uncovered. According to our previous sequencing data, in this study, three novel microsatellites were found in this breed. Using time of flight-mass spectrometry (TOF-MS) method, these microsatellites were further identified in a large Bamei pig population. TOF-MS spectra showed that there are three microsatellites loci, named P1, P2 and P3. These microsatellites were linkage equilibrium based on the values of D' and r² tests. Association results demonstrated that P1 locus was associated with the body length, body height and chest width and the beneficial genotype was 150-/150-bp (p < 0.05); and P2 locus was associated with the body height (p < 0.05), and the 145-/145-bp, 145-/147-bp and 145-/149-bp were claimed as favorable genotypes and 145-bp allele was considered as the favorable allele. These findings suggested that P1 and P2 microsatellites might be considered as the candidate genetic markers to select pigs with superior body sizes, especially in local breed.

A 17-bp InDel (rs668420586) within goat CHCHD7 gene located in growth-related QTL affecting body measurement traits

The Coiled-Coil-Helix-Coiled-Coil-Helix Domain Containing 7 (CHCHD7) gene was located in a growth-related major QTL that participated in the process of bone cells metabolism in animals by regulating mitochondrial copper homeostasis and cytochrome C oxidase assembly. Therefore, we speculated that CHCHD7 gene might be involved in animal growth and body size. Herein, we discovered a 17-bp insertion/deletion (indel) within the goat CHCHD7 gene. Then, we detected this variation in Shaanbei White Cashmere (SBWC, n = 1055) goats and Inner Mongolia White Cashmere (IMWC, n = 743) goats (Alathai type) using the mathematical expectation (ME) method. We then analyzed the correlation between these genotypes and goat body measurement traits. The results showed that the minor allelic frequency (MAF) was 0.011 in SBWC, and 0.048 in IMWC. In SBWC and IMWC, the reaction time by ME method was reduced by 36.78% and 27.59%, respectively, compared to the traditional method of screening samples one by one. Moreover, in SBWC goats, the 17-bp indel was significantly associated with body measurement traits (e.g. body height, and body length) in adults. In IMWC goats, the 17-bp indel was correlated with body measurement traits (e.g. body height) in weaners. In SBWC and IMWC goat populations, the body measurement traits of the individuals homozygous for 17-bp indel were higher than those in heterozygous individuals, except for the case of cannon circumference in IMWC weaners. These findings showed that the 17-bp insertion mutation within the goat CHCHD7 gene significantly affected body morphometric traits, and could provide a basis for marker-assisted selection (MAS) breeding of cashmere goats.

A Deletion Downstream of the CHCHD7 Gene Is Associated with Growth Traits in Sheep

Simple Summary

The genes CHCHD7 and PLAG1 are located on the same growth-related major quantitative trait locus of sheep. PLAG1 affects sheep growth, but no corresponding studies have been conducted on CHCHD7. However, polymorphisms in the CHCHD7 gene are associated with carcass weight and muscle formation in cattle, body height in cattle and humans, and weaning weight in Duroc pigs. In this study, the mathematical expectation method was used to analyze an 8-bp deletion mutation located downstream of the CHCHD7 gene in 2350 individuals from seven sheep breeds. The associations between wild-type and deletion genotypes and growth traits in Tan sheep were also analyzed. The 8-bp deletion locus was significantly associated with body length (p = 0.032), chest depth (p = 0.015), and chest width (p = 0.047) of Tan sheep. Additionally, wild-type genotype carriers were more numerous than those heterozygous for the deletion genotype. Thus, the genotyped 8-bp deletion downstream of the CHCHD7 gene may be associated with growth and development traits in sheep.

Abstract

In sheep, the coiled-coil-helix-coiled-coil-helix domain containing 7 (CHCHD7) gene and the pleiomorphic adenoma gene 1 (PLAG1) are on the same growth-related major quantitative trait locus, positioned head-to-head approximately 420 bp apart on chromosome 9. PLAG1 affects sheep growth, but the effects of CHCHD7 have not been determined. In this study, an 8-bp deletion downstream of CHCHD7 was analyzed in 2350 sheep from seven breeds. The associations between the deletion and growth traits of Tan sheep were also determined. Both genotypes (homozygous wild-type and heterozygous) for the 8-bp deletion were found in Tan (TS), Luxi Blackhead (LXBH), Small-Tail Han (STHS), and Lanzhou Fat-Tail (LFTS) sheep. However, there were no polymorphic sites for the mutation in Hu (HS), Sartuul (SS), and Australian White (AUW) sheep. In TS, LXBH, STHS, and LFTS sheep, the deletion genotype was less frequent than the wild-type genotype, and the allele frequencies of the deletion variant were 0.007 (TS), 0.011 (LBXH), 0.008 (STHS), and 0.010 (LFTS). The 8-bp deletion was significantly associated with body length (p = 0.032), chest depth (p = 0.015), and chest width (p = 0.047) in Tan sheep. Thus, the 8-bp deletion downstream of the CHCHD7 gene might be associated with growth and development traits of sheep.

Detecting novel Indel variants within the GHR gene and their associations with growth traits in Luxi Blackhead sheep

The growth hormone is important in the regulation of metabolism and energy homeostasis and acts through a growth hormone receptor (GHR). In this work, genetic variations within the ovine GHR gene were identified and tested for associations with body morphometric traits in Chinese Luxi Blackhead (LXBH) sheep. Novel deletion loci in the LXBH GHR gene included P2-del-23 bp and P8-del-23 bp indel variants. The polymorphic information content (PIC) was 0.329 in P2-del-23 bp and 0.257 in P8-del-23 bp. Moreover, both indel polymorphisms were not at Hardy-Weinberg equilibrium (p < 0.05) in the LXBH population. Statistical analyses revealed that the P2-del-23 bp and P8-del-23 bp indels were significantly associated (p < 0.05) with several growth traits in rams and ewes, including body weight, body height, chest depth, chest width, chest circumference, cannon circumference, paunch girth and hip width. Among the tested sheep, the body traits of those with genotype DD were superior to those with II and ID genotypes, suggesting that the 'D' allele was responsible for the positive effects on growth traits. Thus, these results indicate that the P2-del-23 bp and P8-del-23 bp indel sites and the DD genotype can be useful in marker-assisted selection in sheep.

Effects of eight InDel variants in FHIT on milk traits in Xinjiang brown cattle

In our previous genome-wide association study (GWAS), we identified the fragile histidine triad diadenosine triphosphatase (FHIT) gene in Xinjiang brown cattle (XJBC) as a candidate gene associated with cattle productive traits, with potential application in mark-assisted selection (MAS) in cattle breeding. FHIT is a prototype of a class of tumor suppressor genes that contain genomic loci mapped to common fragile loci. Here, 388 healthy and unrelated XJBC were selected to identify insertion/deletion (InDel) variants in the bovine FHIT and assess their effects on milk traits. Eight of the thirteen InDel loci were found to be polymorphic in FHIT. The polymorphism information content of the eight loci ranged from 0.061 to 0.375. The correlation analysis showed that all the new InDel variants were significantly related to six different milk traits (p < 0.05). The following variants presented a significant relationship with productive traits: P2-23bp with the 305 milk yield (p = 0.005) in the sixth parity; P3-24bp with the milk fat yield (p = 0.009) in the third parity; P5-21bp with the somatic cell score (p = 0.001) in the first parity and with the milk protein percentage (p = 0.002) in the sixth parity; and P7-26bp with the somatic cell score (p = 0.003) in the sixth parity. These findings will help evaluate InDel genotypes, within and between cattle breeds and identify potential target loci to accelerate progress in MAS in cattle breeding.

A novel 23 bp indel mutation in PRL gene is associated with growth traits in Luxi Blackhead sheep

Prolactin is a highly versatile pituitary hormone with multiple biological functions. PRL expression is regulated by POU1F1 and the prophet of POU1F1 (PROP1). The aim of this study was to investigate the indel variations in ovine PRL and the directly related (PROP1 and POU1F1) genes, and their associations with growth traits in Luxi Blackhead (LXBH) sheep. A monomorphism in PROP1 and POU1F1 genes, and one novel 23-bp insertion mutation in the PRL gene, were identified in LXBH sheep. The 23 bp insertion mutation within PRL gene was significantly associated with several body measurements (e.g., body weight, body height) in sheep of different ages (p < 0.05). Ram lambs (p = 0.036) of genotype insertion/insertion (II) had significantly higher body weights. Weaners (p = 0.018) of genotypes insertion/insertion (II) and insertion/deletion (ID) also had significantly higher body weights compared with male sheep of deletion/deletion (DD) genotype. Moreover, among ewe lambs, individuals of genotype insertion/insertion (II) had a higher paunch girth compared to those with other genotypes (p = 0.044). These findings indicate that a 23 bp indel variant of the ovine PRL gene is correlated with body measurements in LXBH sheep. The findings have potential utility for sheep breeding programs based on marker-assisted selection.

Detection of a 4 bp Mutation in the 3'UTR Region of Goat Sox9 Gene and Its Effect on the Growth Traits

Simple Summary

The sex determining region Y (SRY)-type high mobility group (HMG) box 9 (Sox9) gene is critically important in the formation and development of cartilage and is considered the “main regulator” of chondrogenesis. Additionally, a large number of studies have shown that mutations in a single allele of human Sox9 can lead to campomelic dysplasia syndrome. Therefore, the mutations of Sox9 have been the subject of increasing interest among researchers. However, no studies to date have examined the association between Sox9 gene variants and growth traits in goats. Here, we detected a 4 bp indel in the 3′Untranslated Regions (3′UTR) region of Sox9 in Shaanbei white cashmere (SBWC) goats (n = 1109) and studied the association between this indel and growth traits. The 4 bp indel of Sox9 was significantly associated with body length, heart girth, hip width, and all body measurement indexes (p < 0.05) in SBWC goats. Thus, this deletion could be used as an effective molecular marker for maximizing the growth traits of goats in breeding programs.

Abstract

The SRY-type HMG box 9 (Sox9) gene plays an important role in chondrocyte development as well as changes in hypertrophic chondrocytes, indicating that Sox9 can regulate growth in animals. However, no studies to date have examined the correlation between variations in Sox9 and growth traits in goats. Here, we found a 4 bp indel in the 3′UTR of Sox9 and verified its association with growth traits in Shaanbei white cashmere goats (n = 1109). The frequencies of two genotypes (ID and II) were 0.397 and 0.603, respectively, and polymorphic information content (PIC) values showed that the indel had a medium PIC (PIC > 0.25). The 4 bp indel was significantly correlated with body length (p = 0.006), heart girth (p = 0.001), and hip width (p = 4.37 × 10 ⁻⁴). Notably, individuals with the ID genotype had significantly superior phenotypic traits compared with individuals bearing the II genotype. Hence, we speculated that the 4 bp indel is an important mutation affecting growth traits in goat, and may serve as an effective DNA molecular marker for marker-assisted selection in goat breeding programs.

Detection of rs665862918 (15-bp Indel) of the HIAT1 Gene and its Strong Genetic Effects on Growth Traits in Goats

Simple Summary

Growth traits are important in goats and can affect their body size and meat production. In this study, the hippocampus abundant transcript 1 (HIAT1) gene, which has been reported as a meat-associated trait in elite goat breeds and a dairy-associated trait in water buffalo, was chosen to detect its correlation with growth traits in goats. The results show that the rs665862918 polymorphism (a 15 bp insertion) in HIAT1 is associated with body length, chest width, chest depth, height at hip cross and heart girth in Shaanbei white cashmere goats (SBWC, n = 1013). Our results reveal that rs665862918 in HIAT1 is relevant to the growth traits of goats and could be used for marker-assisted selection (MAS) as a molecular marker in goat populations.

Abstract

The hippocampus abundant transcript 1 (HIAT1) gene, which was detected by the genome-wide identification of selective sweeps among elite goat breeds and water buffalo, is proposed to play an important role in meat characteristics. Four indels of the HIAT1 gene selected from the NCBI and Ensembl databases were detected via a pooling and sequencing strategy. A 15 bp insertion (rs665862918) in the first intron of HIAT1 was selected and classified on an electrophoresis platform in the Shaanbei white cashmere goat (SBWC) population. The correlation analysis revealed that rs665862918 is significantly highly associated with chest width (p = 1.57 × 10⁻⁵), chest depth (p = 8.85 × 10⁻⁵), heart girth (p = 1.05 × 10⁻⁷), body length (p = 0.022), and height at hip cross (p = 0.023) in the SBWC population (n = 1013). Further analysis revealed that individuals with a genotype insertion/insertion (II) of the rs665862918 locus exhibited better growth trait performance than individuals with an insertion/deletion (ID) or deletion/deletion (DD). These findings verify that HIAT1 affects the body size of goats and that rs665862918 could be a potential molecular marker for growth traits in goat breeding.

An 11-bp Indel Polymorphism within the CSN1S1 Gene Is Associated with Milk Performance and Body Measurement Traits in Chinese Goats

The casein alpha s1 (CSN1S1) gene encodes α-s1 casein, one of the proteins constituting milk, which affects milk performance, as well as improving the absorption of calcium and bone development in mammals. A previous study found that an 11-bp insertion/deletion (indel) of this gene strongly affected litter size in goats. However, to our knowledge, the relationships between this polymorphism and the milk performance and body measurement traits of goats have not been reported. In this paper, the previously identified indel has been recognized in three Chinese goat breeds, namely the Guanzhong dairy goat (GZDG; n = 235), Shaanbei white cashmere goat (SBWC; n = 1092), and Hainan black goat (HNBG; n = 278), and the following three genotypes have been studied for all of the breeds: insertion/insertion (II), deletion/deletion (DD), and insertion/deletion (ID). The allele frequencies analyzed signified that the frequencies of the "D" allele were higher (47.8%-65.5%), similar to the previous report, which indicates that this polymorphism is genetically stable in different goat breeds. Further analysis showed that this indel was markedly associated with milk fat content, total solids content, solids-not-fat content, freezing point depression, and acidity in GZDG (p < 0.05), and also affected different body measurement traits in all three breeds (p < 0.05). The goats with II genotypes had superior milk performance, compared with the others; however, goats with DD genotypes had better body measurement sizes. Hence, it may be necessary to select goats with an II or DD genotype, based on the desired traits, while breeding. Our study provides information on the potential impact of the 11-bp indel polymorphism of the CSN1S1 gene for improving the milk performance and body measurement traits in goats.

A novel duplicated insertion/deletion (InDel) of the CPT1a gene and its effects on growth traits in goat

Carnitine palmitoyltransferase 1a (CPT1a) is a rate-limited enzyme in the mitochondrial fatty acid β-oxidation pathway. It acts as a bridge between PPARα and the fatty acid oxidation pathways and is closely related to ruminant growth and development. In this study, one 12 bp InDel polymorphism of the CPT1a gene was identified in 700 goats, and we designated these three genotypes II, ID, and DD. Association analysis showed that the InDel polymorphism was closely associated with trunk index (p = 0.008) and body length index (p = 0.034) in Hainan black goats, and body length (p = 0.010), chest circumference (p = 0.004), chest depth (p = 0.029), and huckle bone width (p = 0.002) in Nubian goats, as well as the chest circumference (p = 0.016) in the Fuqing goat breed. In both kids and adult goats, qRT-PCR results showed that the CPT1a gene was expressed in all tissues, showing the highest mRNA levels in the liver, lung, spleen, and kidney, followed by the adipose tissue and brain. This indicates an association between the InDel of the CPT1a gene and growth traits in selected goat breeds, which may facilitate marker-assisted selection in goat genetics and breeding.

Polymorphisms within the Boule Gene Detected by Tetra-Primer Amplification Refractory Mutation System PCR (T-ARMS-PCR) are Significantly Associated with Goat Litter Size

As a gene contributing to spermatogenesis, the Boule gene (also called Boll), whose mutations result in azoospermia and sterility of flies and mice, was conserved in reductional maturation divisions. However, in goats, the polymorphisms of Boule, especially with regard to their fundamental roles in female reproduction traits, are still unknown. Therefore, the aims of this study were to detect a potential mutation (rs661484476: g.7254T>C) located in intron 2 of the Boule gene by tetra-primer amplification refractory mutation system PCR (T-ARMS-PCR) and to explore its potential association with the litter size of Shaanbei White-Cashmere goats (SBWGs). In this study, g.7254T>C was firstly detected. The TT genotype was the dominant genotype in the single-lamb group, and T was also the dominant allele in all tested groups. Additionally, the detected locus displayed moderate polymorphism with polymorphism information content (PIC) values among all studied goats ranging from 0.303 to 0.344. Notably, according to the χ² test, the distribution differences for the genotypic frequencies between the single- and multi-lamb groups was significant (p = 0.014). Furthermore, the polymorphisms of the goat Boule gene were significantly associated with the goat litter size in SBWGs (p < 0.05), which indicated that g.7254T>C could be a potential marker in the marker-assisted selection process for potential litter size in goats. These results also indicated that the Boule gene might exercise important functions in female goat reproduction, which provided new insight for female goat breeding and could accelerate the process of goat breeding.

A 21-bp indel within the LLGL1 gene is significantly associated with litter size in goat

The scribble cell polarity complex component (LLGL1) is part of the cytoskeletal network and is involved in maintaining cell polarity and epithelial integrity. Based on the whole-genome sequencing analysis in goat, LLGL1 gene is suggested as a putative important candidate gene affecting litter size in Shaanbei White Cashmere Goats (SBWC). Therefore, the objective of this study was to uncover the possible novel insertion/deletion (Indel) variant in goat LLGL1 gene and to evaluate its association with litter size of SBWC (n = 827). Using the PCR detection and DNA sequencing, the 21-bp indel in the upstream of LLGL1 was firstly founded and two genotypes were identified: II (insertion/insertion) and ID (insertion/deletion), respectively. Association analyses revealed that the 21-bp indel was significantly correlated with litter size (p = 0.017). Notably, the individuals with II genotype were significantly greater than that of the genotype ID, and the 'I' allele was dominant. Additionally, the remarkable influence of the indel on traits might be related to the change of DEAF-1-related (NUDR) binding site through bioinformatics analysis. Briefly, the 21-bp indel within the goat LLGL1 gene could be an effective DNA molecular marker and provide valuable theoretical basis for marker-assisted selection (MAS) in goat industry.

Novel lncRNA lncFAM200B: Molecular Characteristics and Effects of Genetic Variants on Promoter Activity and Cattle Body Measurement Traits

Skeletal muscle is one of the three major muscle types in an organism and has key roles in the motor system, metabolism, and homeostasis. RNA-Seq analysis showed that novel lncRNA, lncFAM200B, was differentially expressed in embryonic, neonatal, and adult cattle skeletal muscles. The main aim of this study was to investigate the molecular and expression characteristics of lncFAM200B along with its crucial genetic variations. Our results showed that bovine lncFAM200B was a 472 nucleotide (nt) non-coding RNA containing two exons. The transcription factor binding site prediction analysis found that lncFAM200B promoter region was enriched with SP1 transcription factor, which promotes the binding of myogenic regulatory factor MyoD and DNA sequence. The mRNA expression analysis showed that lncFAM200B was differentially expressed in embryonic, neonatal, adult bovine muscle tissues, and the lncFAM200B expression trend positively correlated with that of MyoG and Myf5 in myoblast proliferation and differential stages. To identify the promoter active region of lncFAM200B, we constructed promoter luciferase reporter gene vector pGL3-Basic plasmids containing lncFAM200B promoter sequences and transfected them into 293T, C2C12, and 3T3-L1 cells. Our results suggested that lncFAM200B promoter active region was from −403 to −139 (264 nt) of its transcription start site, covering 6 SP1 potential binding sites. Furthermore, we found a novel C-T variation, named as SNP2 (ERZ990081 in European Variation Archive) in the promoter active region, which was linked to the nearby SNP1 (rs456951291 in Ensembl database). The genotypes of SNP1 and combined genotypes of SNP1 and SNP2 were significantly associated with Jinnan cattle hip height. The luciferase activity analysis found that the SNP1-SNP2 haplotype CC had the highest luciferase activity, which was consistent with the association analysis result that the combined genotype CC-CC carriers had the highest hip height in Jinnan cattle. In conclusion, our data showed that lncFAM200B is a positive regulator of muscle development and that SNP1 and SNP2 could be used as genetic markers for marker-assisted selection (MAS) breeding of beef cattle.

Detection of Bovine TMEM95 p.Cys161X Mutation in 13 Chinese Indigenous Cattle Breeds

Simple Summary

Chinese indigenous cattle are economically important cattle breeds for animal husbandry development. The promotion and development of Chinese cattle breeds is essential. A previous study found that a nonsense mutation (rs378652941, c.483C > A, p.Cys161X) in bovine transmembrane protein 95 gene (TMEM95) seriously reduced reproductive performance in male Fleckvieh cattle; therefore, this locus was considered a candidate genetic marker in bovine marker-assisted selection (MAS) breeding. Until now, no study has identified this mutation in Chinese cattle breeds. Herein, we detected this c.483C > A mutation in 13 Chinese cattle breeds. Importantly, we found that this mutation did not exist at this locus in our analyzed breeds. Interestingly, we first identified a frameshift insertion/deletion (indel) mutation (NC_037346.1: g.27056998_27057000delCT) in the bovine TMEM95 gene in 11 cattle breeds. Together, the results of this study suggest that the mutation c.483C > A cannot serve as a genetic marker for molecular breeding among Chinese indigenous cattle breeds.

Abstract

Chinese indigenous cattle breeds have abundant genetic resources, which are valuable for the molecular breeding of cattle around the world. Thus, identifying important candidate genes and their genetic markers is very important for cattle molecular breeding. A previous study found that a nonsense mutation (rs378652941, c.483C > A, p.Cys161X) in the bovine transmembrane protein 95 gene (TMEM95) seriously reduced the reproductive performance in bulls, but few studies have detected this mutation in Chinese indigenous cattle breeds. Since the mutation c.483C > A may serve as a potential genetic marker for selecting higher fertility bulls, in the present study, using tetra-primer amplification refractory mutation system PCR (T-ARMS-PCR), forced PCR restriction fragment length polymorphism (forced PCR-RFLP), and DNA sequencing methods, the mutation c.483C > A was detected in 765 individuals from 13 Chinese cattle breeds. However, the results showed that this mutation did not exist at this locus in our analyzed breeds. Interestingly, we identified a newly frameshift insertion/deletion (indel) mutation (NC_037346.1: g.27056998_27057000delCT) in the bovine TMEM95 gene in 11 cattle breeds, which changed the location of the termination codon and changed the 16 amino acids in the C-terminal to 21 amino acids. Combined with previous studies, our study provides evidence that in Chinese cattle breeds the mutation c.483C > A cannot be used as a genetic marker in molecular breeding.

Two Insertion/Deletion Variants within SPAG17 Gene Are Associated with Goat Body Measurement Traits

Simple Summary

Sperm-associated antigen 17 (SPAG17) is a reproduction and skeletal development related gene. This study aimed to identify crucial insertion-deletion (indel) variations, which influence the body measurement traits of goats. Two intronic indels (14 bp and 17 bp indels) were identified by sequencing. In Shaanbei white cashmere goat (SBWC), the different genotypes of the 14 bp indel were markedly associated with goat body height, chest width, body length, and chest depth. The genotypes of the 17 bp indel were significantly associated with body height and chest width. The different combined genotypes were significantly associated with body height and chest width of SBWC and ten traits of Hainan black goat. These results suggested that the 14 and 17 bp indels within SPAG17 can be used in goat growth related traits marker-assisted selection breeding, especially body height.

Abstract

Sperm-associated antigen 17 (SPAG17) gene encodes a multifunctional cytoplasmic protein, which influences not only reproduction but also skeletal development related body measurement traits, especially body height. Thus, this study aimed to identify crucial insertion-deletion (indel) variations, which influence the body measurement traits of goats in large goat populations (n = 1725). As a result, two intronic indels (14 bp and 17 bp indel) were identified by sequencing. For the two indel loci, the distributions of genotypes and alleles were significantly different between the Shaanbei white cashmere goat (SBWC) and the Hainan black goat (HNBG). In SBWC goats, the different genotypes of the 14 bp indel were markedly associated with goat body height, chest width, body length and chest depth. The genotypes of the 17 bp indel were significantly related to body height and chest width. At the two loci, for all seven analyzed traits of SBWC goat, the growth data of DD homozygotes were the worst, which means that the 14 bp insertion and the 17 bp deletion were beneficial and detrimental variations, respectively. Moreover, the combined genotypes were significantly related to body height and chest width of SBWC goats and ten traits of HNBG. These results suggested that the 14 and 17 bp indels within SPAG17 can be used in goat growth related traits marker-assisted selection breeding, especially body height.

A 20-bp insertion/deletion (indel) polymorphism within the CDC25A gene and its associations with growth traits in goat

Cell division cycle 25A (CDC25A), a member of the CDC25 family of phosphatases, is required for progression from G1 to the S phase of the cell cycle. CDC25A provides an essential function during early embryonic development in mice, suggesting that it plays an important role in growth and development. In this study, we used mathematical expectation (ME) methods to identify a 20-bp insertion/deletion (indel) polymorphism of CDC25A gene in Shaanbei White Cashmere (SBWC) goats. We also investigated the association between this 20-bp indel and growth-related traits in SBWC goats. Association results showed that the indel was related to growth traits (height at hip cross, cannon circumference, and cannon circumference index) in SBWC goats. The height at hip cross of individuals with insertion/insertion (II) genotype was higher than those with insertion/deletion (ID) genotype (P=0.02); on the contrary, the cannon circumference and cannon circumference index of individuals with ID genotype were superior when compared with those with II genotype (P=0.017 and P=0.009). These findings suggest that the 20-bp indel in the CDC25A gene significantly affects growth-related traits, and could be utilized as a candidate marker for marker-assisted selection (MAS) in the cashmere goat industry.

Detection of polled intersex syndrome (PIS) and its effect on phenotypic traits in goats

Polled intersex syndrome (PIS), a physiological defect associates polledness and intersexuality, will directly affect the reproduction and other phenotypic traits of goats, thus, it is a hazard to the development of the goat industry. Recent studies have revealed that goat PIS was caused by an 11.7-kb deletion located in chromosome 1q43, but its DNA detection method was little reported. Herein, a total of 503 goats from three goat breeds were used to successfully establish a simple, rapid, accurate and effective method for the detection of goat PIS mutation. Furthermore, based on this assay, the associations between the PIS mutation and growth traits were analyzed in Guanzhong dairy goats. These findings would provide the potential practical application for molecular breeding to accelerate the development of the goat industry.

DNA methylation pattern of the goat PITX1 gene and its effects on milk performance

Paired-like homeodomain transcription factor 1 (PITX1) is a pivotal gene in the hypothalamic–pituitary–adrenal axis, which is a well-known pathway affecting lactation performance. The aim of this study was to analyze the DNA methylation profile of the PITX1 gene and its relevance to milk performance in Xinong Saanen dairy goats; thus, potential epigenetic markers of lactation performance were identified. A total of 267 goat blood samples were divided into “low” and “high” groups according to two milk traits: the average milk yield (AMY) and the average milk density (AMD). One CpG island in the 3′-flanking region of the PITX1 gene was identified as being related to milk performance. Fisher's exact test demonstrated that the methylation rates of the overall CpG island and the 3rd and 12th CpG-dinucleotide loci in the blood were significantly associated with the AMY, and the overall methylation rate of the high AMY group was relative hypomethylation compared with the low AMY group. The overall methylation rates of this CpG island in mammary gland tissue from dry and lactation periods again exhibited a significant difference: the lactation period showed relative hypomethylation compared with the dry period. Bioinformatic transcription factor binding site prediction identified some lactation performance related transcription factors in this CpG island, such as CTCF, STAT, SMAD, CDEF, SP1, and KLFS. Briefly, overall methylation changes of the CpG island in the PITX1 gene are relevant to lactation performance, which will be valuable for future studies and epigenetic marker-assisted selection (eMAS) in the breeding of goats with respect to lactation performance.

Effects of genetic variations within goat PITX2 gene on growth traits and mRNA expression

Paired-like homeodomain 2 (PITX2), a key gene in hypothalamic-pituitary-adrena axis, influences animal growth and development. The objective of this study was to identify the association of the functional genetic variations within goat PITX2 gene with growth traits and mRNA expression levels. According to the reported single nucleotide polymorphisms (SNPs) information in Guanzhong dairy goat (GZDG), we identified genotypes of the known SNPs in Hainan black goat (HNBG). Association analysis uncovered that the SNPs of AC_000163: g.18117T > C, g.18161C > G and g.18353T > C loci were significantly associated with several growth traits (e.g., body height and body length) in HNBG and GZDG breeds. According to the quantitative real-time PCR assay, β-Actin and ribosomal protein L19 (RPL19) were the most stable expressed housekeeping genes in heart and skeletal muscle, respectively; meanwhile, glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was the most stable expressed housekeeping gene in the other tissues. Based on the best housekeeping gene of varied tissues, we found the different genotypes of above loci were significantly associated with PITX2 mRNA expression in heart, muscle and small intestine. Briefly, the genetic variants of goat PITX2 associating with mRNA expression affected growth traits significantly, which would benefit for goat breeding in the future.

One 16 bp insertion/deletion (indel) within the KDM6A gene revealing strong associations with growth traits in goat

Lysine demethylase 6A (KDM6A), as a candidate gene associated with skeletal dysplasia and the differentiation of osteogenesis, might affect the growth traits in animals. In the previous study (Cui et al., 2018), one 16 bp intronic insertion/deletion (NW_017189516.1: g.138, 431_138,446delAATGTATAGCTTAAAA) within the KDM6A gene significantly influenced KDM6A gene expression. The objective of this work was to investigate the association between this 16 bp indel and growth-related traits in Shaanbei White Cashmere goat (SWCG). Association analyses showed the 16 bp indel was related to growth traits strongly (body height, chest depth, height across the hip, body length and chest circumference) in SWCG population (n = 1953, P < 0.05). Insertion/insertion (II) genotype individuals (n = 1502) had the best growth traits, by comparison of those of insertion/deletion (ID) (n = 410) and the deletion/deletion (DD) genotypes (n = 41). Particularly, the body weight of II genotype individuals were significantly higher than ID and DD genotypes (P < 0.01). Besides, the remarkable influence of this indel on traits might cause by the change of MEF2 binding site. These findings hinted that the 16 bp deletion mutation in KDM6A gene, which was significantly associated with growth-related traits, could be assigned to an effective molecular marker for growth traits in goat breeding.