A novel indel within goat casein alpha S1 gene is significantly associated with litter size

A novel 22-bp InDel within FGF7 gene is significantly associated with growth traits in goat

Fibroblast growth factor 7 (FGF7) is involved in lipid metabolism, which is considered as a candidate gene with close relation with muscle development by eGWAs and RNA-Seq analyses. To date, limited research has been conducted on the relationship between FGF7 gene and growth traits. The main objective of this work was to further investigate the association between novel InDel within FGF7 gene and growth traits in goat. Herein, FGF7 mRNA expression levels were investigated in various Fuqing goat tissues. We found that FGF7 gene was expressed in six adult goat tissues with the highest mRNA levels in adipose tissue. This result suggested that FGF7 gene might play a critical role in fat deposition. We also detected potential polymorphisms in Fuqing, Nubian and Jianyang Daer breeds. A 22-bp InDel polymorphism in FGF7 gene was detected in 396 goats and the three genotypes were designated as II, ID, and DD. Correlation analysis revealed that InDel polymorphism was significantly associated with growth traits (P < 0.05). Goats with genotypes ID and/or II had superior growth traits compared to those with genotype DD. In summary, our findings suggested that the 22-bp InDel within FGF7 gene could act as a molecular marker to improve the growth traits of goats in breeding programs.

Relationship between Indel Variants within the JAK2 Gene and Growth Traits in Goats

Janus kinase 2 (JAK2) plays a critical role in myoblast proliferation and fat deposition in animals. Our previous RNA-Seq analyses identified a close association between the JAK2 gene and muscle development. To date, research delving into the relationship between the JAK2 gene and growth traits has been sparse. In this study, we sought to investigate the relationship between novel mutations within the JAK2 gene and goat growth traits. Herein, two novel InDel (Insertion/Deletion) polymorphisms within the JAK2 gene were detected in 548 goats, and only two genotypes were designated as ID (Insertion/Deletion) and DD (Deletion/Deletion). The results indicate that the two InDels, the del19008 locus in intron 2 and del72416 InDel in intron 6, showed significant associations with growth traits (p < 0.05). Compared to Nubian and Jianzhou Daer goats, the del72416 locus displayed a more pronounced effect in the Fuqing breed group. In the Nubian breed (NB) group, both InDels showed a marked influence on body height (BH). There were strong linkages observed for these two InDels between the Fuqing (FQ) and Jianzhou (JZ) populations. The DD-ID diplotype was associated with inferior growth traits in chest width (ChW) and cannon circumference (CaC) in the FQ goats compared to the other diplotypes. In the NB population, the DD-DD diplotype exhibited a marked negative impact on BH and HuWI (hucklebone width index), in contrast to the other diplotypes. In summary, our findings suggest that the two InDel polymorphisms within the JAK2 gene could serve as valuable molecular markers for enhancing goat growth traits in breeding programs.

Relationships between the mutations of the goat GATA binding protein 4 gene and growth traits

Osteogenesis is a complex multilevel process regulated by multiple genes. The GATA binding protein 4 (GATA4) gene has been extensively studied for its pivotal role in bone genesis and bone differentiation. However, its relationship with the growth traits of Shaanbei white cashmere (SBWC) and Guizhou black (GB) goats remains unclear. This work aims to investigate the potential influence of genetic mutations in the GATA4 gene on the growth traits goats. Thus, two Insertion/deletion (InDel) polymorphisms (8-bp-InDel and 9-bp-InDel) were screened and detected in a total of 1161 goats (including 980 SBWC goats and 181 GB goats) using PCR and agarose gel electrophoresis. The analyses revealed that there were two genotypes (ID and DD) for these two loci. In SBWC goats, 8-bp-InDel and 9-bp-InDel loci were significantly associated with heart girth (HG) and hip width (HW). Notably, individuals with DD genotype of 8-bp-InDel locus were superior while those with DD genotype of 9-bp-InDel locus were inferior. Correlation analyses of the four combined genotypes revealed significant associations with cannon circumference (CC), body height (BH), HG and HW. This work provides a foundation for the application of molecular marker-assisted selection (MAS) in goat breeding programs. Furthermore, the findings highlight the potential of the GATA4 gene and its genetic variations as valuable indicators for selecting goats with desirable growth traits.

Goat CLSTN2 gene: tissue expression profile, genetic variation, and its associations with litter size

Calsyntenin-2 (CLSTN2) is involved in cell proliferation, differentiation, cell death, tumorigenesis, and follicular expression. Although CLSTN2 has been identified as a potential candidate gene for sheep prolificacy, no studies have been done on its effect on goat prolificacy. The purpose of this study was to identify mRNA expression and genetic variation within goat CLSTN2, and its association with prolificacy. Herein, we uncovered significant differences in mRNA levels of the CLSTN2 gene in different tissues in female goats (p < 0.01), including ovary tissue. Nine putative indels were designed to investigate their correlation to litter size, but only one 16-bp deletion was discovered in female Shaanbei white cashmere goats (n = 902). We discovered that a 16-bp deletion within the CLSTN2 gene was significantly correlated with first-born litter size (p = 0.0001). As shown by the chi-squared test, the genotypic II of single-lambs and multi-lambs was dramatically higher than with genotype ID (p = 0.005). Our findings suggest that indel within the CLSTN2 gene is a candidate gene affecting prolificacy in goats and may be used for Marker Assisted Selection (MAS) in goats.

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.

A genome-wide perspective on the diversity and selection signatures in indigenous goats using 53 K single nucleotide polymorphism array

Tibetan goats, Taihang goats, Jining grey goats, and Meigu goats are the representative indigenous goats in China, found in Qinghai-Tibet Plateau, Western pastoral area, Northern and Southern agricultural regions. Very few studies have conducted a comprehensive analysis of the genomic diversity and selection of these breeds. We genotyped 96 unrelated individuals, using goat 53 K Illumina BeadChip array, of the following goat breeds: Tibetan (TG), Taihang (THG), Jining grey (JGG), and Meigu (MGG). A total of 45 951 single nucleotide polymorphisms were filtered to estimate the genetic diversity and selection signatures. All breeds had a high proportion (over 95%) of polymorphic loci. The observed and excepted heterozygosity ranged from 0.338 (MGG) to 0.402 (JGG) and 0.339 (MGG) to 0.395 (JGG), respectively. Clustering analysis displayed a genetically distinct lineage for each breed, and their Fst were greater than 0.25, indicating that they had a higher genetic differentiation between groups. Furthermore, effective population size reduced in all four populations, indicating a loss of genetic diversity. In addition, runs of homozygosity were mainly distributed in 5-10 Mb. Lastly, we identified signature genes, which were closely related to high-altitude adaptation (ADIRF) and prolificity (CNTROB, SMC3, and PTEN). This study provides a valuable resource for future studies on genome-wide perspectives on the diversity and selection signatures of Chinese indigenous goats.

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.

Goat SMAD family member 1 (SMAD1): mRNA expression, genetic variants, and their associations with litter size

SMAD family member 1 (SMAD1) is phosphorylated and activated by the BMP receptors, which help regulate ovulation rate, cell growth, apoptosis, and development. Previously, the genome-wide association study revealed that it has been associated with fecundity in sheep. However, its effect on litter size has not been investigated in goats. Therefore, this study aimed to determine the level of SMAD1 mRNA expression in various tissues and to identify its polymorphisms and their association with litter size in Shaanbei white cashmere goat (SBWC). As a result, RT-qPCR analysis showed that SMAD1 was expressed in various tissues in female SBWC goats, including the ovary (P < 0.05). Importantly, the mRNA expression level in the ovaries of mothers of multi-lambs had a higher level than the mothers of single lambs (P < 0.05). Moreover, two InDels (18-bp and 7-bp) in intron 1 of SMAD1 were polymorphic among ten potential loci. Both 18-bp and 7-bp InDels were significantly correlated with litter size (P = 0.014) and (P = 0.0001), respectively. As shown by the chi-squared test, genotypic distributions of 18-bp and 7-bp were significantly distinct between single-lamb (P = 0.02) and multi-lamb mothers (P = 0.002). Our findings confirm that two InDels in SMAD1 were significantly associated with litter size and suggest that they could be used to improve fertility traits in goat breeding strategies.

Effect of Upregulation of Transcription Factor TFDP1 Binding Promoter Activity Due to RBP4 g.36491960G>C Mutation on the Proliferation of Goat Granulosa Cells

Retinol-binding protein 4 (RBP4), a member of the lipocalin family, is a specific carrier of retinol (vitamin A) in the blood. Numerous studies have shown that RBP4 plays an important role in mammalian embryonic development and that mutations in RBP4 can be used for the marker-assisted selection of animal reproductive traits. However, there are few studies on the regulation of reproduction and high-prolificacy traits by RBP4 in goats. In this study, the 5′ flanking sequence of RBP4 was amplified, and a G>C polymorphism in the promoter region -211 bp (g.36491960) was detected. An association analysis revealed that the respective first, second and third kidding number and mean kidding number of nanny goats with CC and GC genotypes (2.167 ± 0.085, 2.341 ± 0.104, 2.529 ± 0.107 and 2.189 ± 0.070 for CC and 2.052 ± 0.047, 2.206 ± 0.057, 2.341 ± 0.056 and 2.160 ± 0.039 for GC) were significantly higher (p < 0.05) than those with the GG genotype (1.893 ± 0.051, 2.027 ± 0.064, 2.107 ± 0.061 and 1.74 ± 0.05). The luciferase assay showed that luciferase activity was increased in C allele individuals compared with that in G allele individuals. A competitive electrophoretic mobility shift assay (EMSA) showed that individuals with the CC genotype had a stronger promoter region binding capacity than those with the GG genotype. In addition, transcription factor prediction software showed that the RBP4 g.36491960G>C mutation added a novel binding site for transcription factor DP-1 (TFDP1). RT−qPCR results showed that the expression of TFDP1 was significantly higher in the high-prolificacy group than in the low-prolificacy group, and the expression of RBP4 was higher in both the CC and GC genotypes than that in the GG genotype. TFDP1 overexpression significantly increased the expression of RBP4 mRNA (p < 0.05) and the expression of the cell proliferation factors cyclin-D1, cyclin-D2 and CDK4 (p < 0.05). The opposite trend was observed after interference with TFDP1. Both the EdU and CCK-8 results showed that TFDP1 expression could regulate the proliferation of goat ovarian granulosa cells. In summary, our results showed that RBP4 g.36491960G>C was significantly associated with fecundity traits in goats. The g.36491960G>C mutation enhanced the transcriptional activity of RBP4 and increased the expression of RBP4, thus improving the fertility of Yunshang black goats.

Genetic polymorphisms within the ETAA1 gene associated with growth traits in Chinese sheep breeds

Ewing tumor-associated antigen 1 (ETAA1) is an important candidate gene for fat deposition in sheep. This study aimed to investigate the indel variations in the ETAA1 gene and their associations with growth traits in the 1619 sheep. The polymorphism information content of this indel ranged from 0.308 to 0.375, with medium genetic diversity (0.25 ≤ polymorphism information content ≤ 0.5). The correlation analysis showed that an 8-bp insertion within the ETAA1 gene was significantly associated with growth traits in Luxi Blackhead sheep (LXBH), Lanzhou fat-tailed sheep, Hu sheep, Tong sheep, and Tan sheep (p < 0.05). Furthermore, the ETAA1 gene was significantly associated with several growth traits (p < 0.01), such as chest width and paunch girth of LXBH, Tong, and Tan rams. It was significantly related to the body morphometric traits of LXBH, Lanzhou fat-tailed sheep, Hu, Tong, and Tan ewes. In conclusion, the detected 8-bp indel within the ETAA1 gene was confirmed in sheep, significantly affecting the growth traits, and might be a potential DNA marker for the selection of high-quality individuals in marker-assisted selection for sheep breeding.

Genome-wide detection of genetic structure and runs of homozygosity analysis in Anhui indigenous and Western commercial pig breeds using PorcineSNP80k data

Background

Runs of homozygosity (ROH) are continuous homozygous regions typically located in the DNA sequence of diploid organisms. Identifications of ROH that lead to reduced performance can provide valuable insight into the genetic architecture of complex traits. Here, we systematically investigated the population genetic structure of five Anhui indigenous pig breeds (AHIPs), and compared them to those of five Western commercial pig breeds (WECPs). Furthermore, we examined the occurrence and distribution of ROHs in the five AHIPs and estimated the inbreeding coefficients based on the ROHs (F_ROH) and homozygosity (F_HOM). Finally, we identified genomic regions with high frequencies of ROHs and annotated candidate genes contained therein.

Results

The WECPs and AHIPs were clearly differentiated into two separate clades consistent with their geographical origins, as revealed by the population structure and principal component analysis. We identified 13,530 ROHs across all individuals, of which 4,555 and 8,975 ROHs were unique to AHIPs and WECPs, respectively. Most ROHs identified in our study were short (< 10 Mb) or medium (10–20 Mb) in length. WECPs had significantly higher numbers of short ROHs, and AHIPs generally had longer ROHs. F_ROH values were significantly lower in AHIPs than in WECPs, indicating that breed improvement and conservation programmes were successful in AHIPs. On average, F_ROH and F_HOM values were highly correlated (0.952–0.991) in AHIPs and WECPs. A total of 27 regions had a high frequency of ROHs and contained 17 key candidate genes associated with economically important traits in pigs. Among these, nine candidate genes (CCNT2, EGR2, MYL3, CDH13, PROX1, FLVCR1, SETD2, FGF18, and FGF20) found in WECPs were related to muscular and skeletal development, whereas eight candidate genes (CSN1S1, SULT1E1, TJP1, ZNF366, LIPC, MCEE, STAP1, and DUSP) found in AHIPs were associated with health, reproduction, and fatness traits.

Conclusion

Our findings provide a useful reference for the selection and assortative mating of pig breeds, laying the groundwork for future research on the population genetic structures of AHIPs, ultimately helping protect these local varieties.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08583-9.

A hypomorphic mutation in the mouse Csn1s1 gene generated by CRISPR/Cas9 pronuclear microinjection

Caseins are major milk proteins that have an evolutionarily conserved role in nutrition. Sequence variations in the casein genes affect milk composition in livestock species. Regulatory elements of the casein genes could be used to direct the expression of desired transgenes into the milk of transgenic animals. Dozens of casein alleles have been identified for goats, cows, sheep, camels and horses, and these sequence variants are associated with altered gene expression and milk protein content. Most of the known mutations affecting casein genes’ expression are located in the promoter and 3’-untranslated regions. We performed pronuclear microinjections with Cas9 mRNA and sgRNA against the first coding exon of the mouse Csn1s1 gene to introduce random mutations in the α-casein (Csn1s1) signal peptide sequence at the beginning of the mouse gene. Sanger sequencing of the founder mice identified 40 mutations. As expected, mutations clustered around the sgRNA cut site (3 bp from PAM). Most of the mutations represented small deletions (1–10 bp), but we detected several larger deletions as well (100–300 bp). Functionally most mutations led to gene knockout due to a frameshift or a start codon loss. Some of the mutations represented in-frame indels in the first coding exon. Of these, we describe a novel hypomorphic Csn1s1 (Csn1s1c.4-5insTCC) allele. We measured Csn1s1 protein levels and confirmed that the mutation has a negative effect on milk composition, which shows a 50 % reduction in gene expression and a 40–80 % decrease in Csn1s1 protein amount, compared to the wild-type allele. We assumed that mutation affected transcript stability or splicing by an unknown mechanism. This mutation can potentially serve as a genetic marker for low Csn1s1 expression.

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).

Population validation of reproductive gene mutation loci and association with the litter size in Nubian goat

Litter size is an important component trait of doe reproduction. By improving it, production efficiency and economic benefits can be significantly provided. Genetic marker-assisted selection (MAS) based on proven molecular indicators could enhance the efficacy of goat selection, as well as litter size trait. Many molecular markers have been identified that they can be used to improve litter size in different goat breeds. However, the presence and value of these markers vary among goat breeds. In the present study, we used the reported loci on other breeds of goat as candidate loci to detect whether these loci appear in this Nubian goat population; then we proceed to genotype and detect surrounding loci (50 bp) by multiplex PCR and sequencing technology. As a result, 69 mutations (59 SNPs and 10 indels) were screened out from 23 candidate genes in Nubian goat population, 12 loci were significantly associated with the litter size of first-parity individuals; 5 loci were significantly associated with the litter size of second-parity individuals; 3 loci were significantly associated with the litter size of third-parity individuals. In addition, five loci were significantly associated with the average litter size. The additive effect value of KITLG: g.18047318 G>A in first parity, KITLG: g.18152042G>A in third parity, KISS-1: g.1341674 C>G in first parity, and GHR: g.32134187G>A in second parity exceed more than 0.40, and the preponderant alleles are G, C, A and G, respectively. Further, linkage disequilibrium analysis of 21 mutation loci shows that 3 haplotype blocks are formed, and the litter size of combination type AACC in KISS-1 gene and AAGG in KITLG gene are significantly lower than that of other combinations genotype in first parity (P<0.05). These findings can provide effective candidate DNA markers for selecting superior individuals in Nubian goat breeding.

A deletion mutation within the goat AKAP13 gene is significantly associated with litter size

A-kinase anchoring protein 13 (AKAP13) is one of the AKAP protein family members, which is correlated with estrogen receptors (ERs) and progesterone receptor (PR) activity. Consequently, the AKAP13 gene is considered to be one of the candidate genes for regulating female fertility. Hence, the objectives of this study were to discover the potential insertion/deletion (indel) variants within the AKAP13 gene and evaluate their associations with litter size of Shaanbei white cashmere goats (SBWC) to screen candidate genes for the molecular marker-assisted selection (MAS). Ultimately, we found the 16-bp deletion of AKAP13 gene which displayed three genotypes (II, ID and DD). However, it was not confirmed to Hardy-Weinberg equilibrium (HWE) in the tested population. Statistical analysis demonstrated that this 16-bp indel locus was significantly associated with litter size in goats (p < 0.05), in which the ID genotype was a key genotype for increasing litter size in goats. Besides, independent χ² tests between different genotypes and litter size showed that high-prolific groups had higher frequency of the 'D' allele (p < 0.05). Briefly, AKAP13 gene is a candidate gene for improving fertility, and its 16-bp indel locus can be used as a valid DNA molecular marker for the MAS in goat breeding.

Detection of 15-bp Deletion Mutation within PLAG1 Gene and Its Effects on Growth Traits in Goats

Simple Summary

Goats have always served as an important domesticated livestock. PLAG1 is a major gene that affects the stature and growth of animals. Body size traits are very important for goats as it directly affects the economic characteristics of meat and cashmere production. This study showed that the 15-base pair (bp) InDel (rs637141549) can significantly affect growth traits such as body weight, height, height at hip cross, chest circumference, hip width and body index of goats through the detection of large samples (n = 1581) in four indigenous breeds. Accordingly, it is suggested that the deletion mutation can be used as a potential molecular marker that significantly affects goat growth traits. Moreover, the 15bp deletion mutation can be used as a potential molecular marker, which significantly affects the growth traits of goats and plays an important role in animal husbandry production.

Abstract

Stature and weight are important growth and development traits for animals, which also significantly affect the productivity of livestock. Polymorphic adenoma gene 1 (PLAG1) is located in the growth-related quantitative trait nucleotides (QTN), and its variation has been determined to significantly affect the body stature of bovines. This study found that novel 15-bp InDel could significantly influence important growth traits in goats. The frequencies of genotypes of the 15-bp mutation and relationship with core growth traits such as body weight, body height, height at hip cross, chest circumference, hip width and body index were explored in 1581 individuals among 4 Chinese native goat breeds. The most frequent genotypes of Shaanbei white Cashmere goat (SWCG), Inner Mongolia White Cashmere goat (IMCG) and Guanzhong Dairy goat (GZDG) were II genotypes (insertion/insertion), and the frequency of ID genotype (insertion/deletion) was found to be slightly higher than that of II genotype in Hainan Black goat (HNBG), showing that the frequency of the I allele was higher than that of the D allele. In adult goats, there were significant differences between 15-bp variation and body weight, chest circumference and body height traits in SWCG (p < 0.05). Furthermore, the locus was also found to be significantly correlated with the body index of HNBG (p = 0.044) and hip width in GZDG (p = 0.002). In regard to lambs, there were significant differences in height at the hip cross of SWCG (p = 0.036) and hip width in IMWC (p = 0.005). The corresponding results suggest that the 15-bp InDel mutation of PLAG1 is associated with the regulation of important growth characteristics of both adult and lamb of goats, which may serve as efficient molecular markers for goat breeding.

Shared Ancestry and Signatures of Recent Selection in Gotland Sheep

Gotland sheep, a breed native to Gotland, Sweden (an island in the Baltic Sea), split from the Gute sheep breed approximately 100 years ago, and since, has probably been crossed with other breeds. This breed has recently gained popularity, due to its pelt quality. This study estimates the shared ancestors and identifies recent selection signatures in Gotland sheep using 600 K single nucleotide polymorphism (SNP) genotype data. Admixture analysis shows that the Gotland sheep is a distinct breed, but also has shared ancestral genomic components with Gute (~50%), Karakul (~30%), Romanov (~20%), and Fjällnäs (~10%) sheep breeds. Two complementary methods were applied to detect selection signatures: A Bayesian population differentiation F_ST and an integrated haplotype homozygosity score (iHS). Our results find that seven significant SNPs (q-value < 0.05) using the F_ST analysis and 55 significant SNPs (p-value < 0.0001) using the iHS analysis. Of the candidate genes that contain significant markers, or are in proximity to them, we identify several belongings to the keratin genes, RXFP2, ADCY1, ENOX1, USF2, COX7A1, ARHGAP28, CRYBB2, CAPNS1, FMO3, and GREB1. These genes are involved in wool quality, polled and horned phenotypes, fertility, twining rate, meat quality, and growth traits. In summary, our results provide shared founders of Gotland sheep and insight into genomic regions maintained under selection after the breed was formed. These results contribute to the detection of candidate genes and QTLs underlying economic traits in sheep.

The mRNA expression profile of the goat prion protein testis-specific (PRNT) gene and its associations with litter size

The goat PRNT gene was initially identified as a testis-specific gene with a role in spermatogenesis. In this study, we used quantitative real-time PCR (qPCR) to first determine the mRNA expression profile of this gene in different goat tissues. Surprisingly, we found that PRNT was expressed not only in the testis but also in nine other tissues in goats. Moreover, PRNT was weakly expressed in the testis, while its expression was strongest in the ovary. These results, combined with those of other studies, led us to hypothesize that the goat PRNT gene has a role in both male and female reproduction. We further used direct DNA sequencing to detect potential SNPs within this gene in Shaanbei whit cashmere (SBWC) rams and ewes, and identified three SNPs within the PRNT gene, namely, c.-58C > T, c.71A > G (p.Alanine24Valine), and c.102C > T (synonymous). In rams, c.-58C > T and c.102C > T were strongly linked with each other (D' = 1.000, r² = 0.504), whereas no significant association (P > 0.05) was found between the three SNPs and semen quality, which was consistent with the low expression of the PRNT gene in the testis. Interestingly, in ewes (n = 502), c.-58C > T and c.71A > G were also strongly linked with each other (D' = 0.973, r² = 0.537). Additionally, the c.71A > G locus, especially the AA genotype, had a significant influence on litter size (P = 0.006), consistent with the high PRNT expression in the ovary. Combined, the results of the expression profiling and analysis of the association between the SNPs and reproductive traits showed that two strongly linked nucleotide sequence variants within PRNT were significantly associated with goat litter size. These findings provide potential DNA markers for use in the marker-assisted selection (MAS) of goats with high-fertility traits.

An insertion/deletion within the CREB1 gene identified using the RNA-sequencing is associated with sheep body morphometric traits

In a previous study, the cyclic AMP response element-binding protein 1 (CREB1) gene, which is likely involved in the regulation of fat metabolism in sheep adipose tissue, was identified using RNA sequencing. CREB1 is a transcription factor that participates in the regulation of cell proliferation, differentiation, and survival as well as energy metabolism. Therefore, based on preliminary studies, this study aimed to reveal the correlation between the insertion/deletion (indel) polymorphism of the CREB1 gene and sheep growth traits. One insertion variation of the ovine CREB1 gene, C3-ins-26 bp, was investigated in 1847 Chinese and Mongolian sheep breeds. The minor allele frequencies in the CREB1 gene varied from 0.021 to 0.938. Further, statistical analyses indicated that the C3-ins-26 bp indel in the CREB1 gene was significantly related to various body measurements (body length, height, and index; chest width, depth, and width index; cannon circumference index; and height at the hip cross) in a Tan sheep population (p < 0.05). Collectively, these findings may provide important insights into marker-assisted selection of sheep.

Whole-genome sequencing to identify candidate genes for litter size and to uncover the variant function in goats (Capra hircus)

To select candidate genes for goat prolificacy, we managed six multi- and six single-kid female goats at the same feeding level and in the same management mode over a 4-year period. These goats showed stable differences in litter size over five continuous parturition records. Whole-genome re-sequencing was used in all 12 to select candidate genes, namely, AURKA, ENDOG, SOX2, RORA, GJA10, RXFP2, CDC25C, and NANOS3, by the strength of their differentiation signals. Most of the selected genes were enriched in the coiled coil process and ovarian development, which suggests that the coiled coil process has a potential regulatory effect on fecundity. Detection of the distribution of variants and association analyses with litter size in 400 goats showed that NANOS3 exon mutations may lead to a transformation of the protein structure. The variation in CDC25C, ENDOG, and NANOS3 showed a significant association with litter size. These results can contribute to the improvement of reproduction traits in the artificial breeding of goats.

Relationships between novel nucleotide variants within the colony-stimulating factor 1 receptor (CSF1R) gene and mastitis indicators in sheep

Colony-stimulating factor 1 receptor (CSF1R) plays an important role in the process of innate immunity and inflammation, thus it was hypothesized that the CSF1R gene might affect the occurrence of mammalian mastitis. The purpose of this study was to investigate the association between nucleotide variations of CSF1R gene and mastitis in Australian white sheep (AUWs). Two indel variants (Intron5-27 bp and Intron5-22 bp) within the CSF1R gene have been found in AUWs. The Chi-square test for different mastitis symptoms demonstrated that individuals without symptoms of mastitis had higher 'I' allele frequencies and 'II' genotype frequencies (p < 0.01). We found strong correlation between mastitis and lactation score through Pearson correlation analysis. Therefore, we also analyzed the relationship between the two indel loci and lactation, we found that the lactation ability of individuals with type II was stronger than that of DD genotype at the Intron5-22 bp (p < 0.05). Additionally, we found that the combined genotype of the two loci was significantly associated with mastitis (p < 0.01). These findings indicated that CSF1R mutations were significantly associated with mastitis, and could affect lactation performance, suggesting that two deletion sites could be used as the effective molecular markers against mastitis in sheep breeding.

Determination of genetic variation within the DYRK2 gene and its associations with milk traits in cattle

To speed up the progress of marker-assisted selection (MAS) in cattle breeding, the dual-specificity tyrosine phosphorylation-regulated kinase 2 (DYRK2), cadherin 2 (CDH2), and kinesin family member 1A (KIF1A) genes were chosen based on our pervious genome-wide association study (GWAS) analysis results. DYRK2 is a kinase that may participate in cell growth and/or development; it shows phosphorylation activity toward serine, threonine, and tyrosine fragments of proteins, and it is different from other protein kinases. The CDH2 gene encodes a classic cadherin, which is a member of the cadherin superfamily. The protein encoded by KIF1A is a member of the kinesin family and plays a role in the transportation of membrane organelles along axon microtubules. We detected insertion/deletion (InDel) variation in these three candidate genes in 438 individual cattle (Xinjiang Brown cattle and Wagyu × Luxi crossbreed cattle). Only DYRK2-P3-11 bp was polymorphic and genotyped. The polymorphism information content of DYRK2-P3-11 bp was 0.336. Correlation analyses showed that InDel polymorphism was significantly associated with six different milk traits. These findings may aid future analyses of InDel genotypes in cattle breeds, and speed up the progress of MAS in cattle breeding.

Indel variants within the PRL and GHR genes associated with sheep litter size

Growth hormone and prolactin belong to the class of peptide hormones that have a wide range of regulatory functions. In this study, polymorphisms of growth hormone receptor (GHR) and prolactin (PRL) genes were analysed as candidate genes, which are responsible for the litter size in Australian White (AUW) sheep. According to the statistical analyses results, the polymorphism information content (PIC) values of the PRL-P1-ins-23 bp, GHR-P2-del-23 bp and GHR-P8-del-23 bp were 0.371, 0.366 and 0.375, respectively, which indicates the high genetic polymorphism in AUW sheep. Moreover, all indel loci are not conformed to the HWE (p < .05). Further, our findings revealed that the PRL-P1-ins-23 bp polymorphism in the ovine PRL gene was significantly related to the first parity litter size (p = .001) and the DD genotype displaying the highest genotypic mean. Meanwhile, the GHR-P2-del-23 bp and GHR-P8-23 bp indels in the ovine GHR gene were significantly correlated with first parity litter size (p < .05), and the individuals with the genotype II showed significantly higher litter size than others. Collectively, these results demonstrated that our findings could be useful for future sheep breeding strategies based on the molecular-assisted selection (MAS).

Two Novel Rare Strongly Linked Missense SNPs (P27R and A85G) Within the GDF9 Gene Were Significantly Associated With Litter Size in Shaanbei White Cashmere (SBWC) Goats

Growth differentiation factor 9 (GDF9) is a high-fertility candidate gene that plays a crucial role in early folliculogenesis in female mammals. In this study, direct sequencing was used to screen possible SNP loci in the goat GDF9 gene. Three SNP loci, p.proline27alanine (P27R), p.leucine61leucine (L61L), and p.alanine85glycine (A85G), were identified in Shaanbei white cashmere (SBWC) goats. Among the three SNPs, two rare missense SNP loci (P27R and A85G) were discovered to be strongly linked with each other (D′ value = 0.926, r² value = 0.703). Both P27R and A85G loci had two genotypes: wild type and heterozygous type. A85G exerted a significant effect on litter size (P = 0.029) in SBWC goats, and the heterozygous genotype was superior in comparison with the wild type. The heterozygous genotype was also superior in P27R but no significant association was found. However, the combination genotypes of P27R and A85G were identified to have superior effects on litter size (P = 3.8E−15). This information suggested that these two SNPs influenced litter size in goats synergistically. Combining this information with our previous studies, we propose that the GDF9 gene is the principal high-fertility candidate gene and that the A85G locus is a promising SNP that affects litter size in goats. These results may fill a research gap regarding rare mutations as well as provide crucial molecular markers that could be useful in marker-assisted selection (MAS) goat rearing when selecting superior individuals.

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.

Comprehensive analysis of mRNAs and miRNAs in the ovarian follicles of uniparous and multiple goats at estrus phase

Background

Fertility is an important economic trait in the production of meat goat, and follicular development plays an important role in fertility. Although many mRNAs and microRNAs (miRNAs) have been found to play critical roles in ovarian biological processes, the interaction between mRNAs and miRNAs in follicular development is not yet completely understood. In addition, less attention has been given to the study of single follicle (dominant or atretic follicle) in goats. This study aimed to identify mRNAs, miRNAs, and signaling pathways as well as their interaction networks in the ovarian follicles (large follicles and small follicles) of uniparous and multiple Chuanzhong black goats at estrus phase using RNA-sequencing (RNA-seq) technique.

Results

The results showed that there was a significant difference in the number of large follicles between uniparous and multiple goats (P < 0.05), but no difference in the number of small follicles was observed (P > 0.05). For the small follicles of uniparous and multiple goats at estrus phase, 289 differentially expressed mRNAs (DEmRNAs) and 16 DEmiRNAs were identified; and for the large follicles, 195 DEmRNAs and 7 DEmiRNAs were identified. The functional enrichment analysis showed that DE genes in small follicles were significantly enriched in ovarian steroidogenesis and steroid hormone biosynthesis, while in large follicles were significantly enriched in ABC transporters and steroid hormone biosynthesis. The results of quantitative real-time polymerase chain reaction were consistent with those of RNA-seq. Analysis of the mRNA-miRNA interaction network suggested that CD36 (miR-122, miR-200a, miR-141), TNFAIP6 (miR-141, miR-200a, miR-182), CYP11A1 (miR-122), SERPINA5 (miR-1, miR-206, miR-133a-3p, miR-133b), and PTGFR (miR-182, miR-122) might be related to fertility, but requires further research on follicular somatic cells.

Conclusions

This study was used for the first time to reveal the DEmRNAs and DEmiRNAs as well as their interaction in the follicles of uniparous and multiple goats at estrus phase using RNA-seq technology. Our findings provide new clues to uncover the molecular mechanisms and signaling networks of goat reproduction that could be potentially used to increase ovulation rate and kidding rate in goat.

A 5-bp mutation within MSTN/GDF8 gene was significantly associated with growth traits in Inner Mongolia White Cashmere goats

Myostatin (MSTN) gene, also known as growth and differentiation factor 8 (GDF8) gene, is a negative regulator of skeletal muscle growth and development, especially the number, size and type of muscle fibers. Its mutations contribute to the double-muscling (DBM) phenomenon which significantly increases the muscle mass. Hence, variations within MSTN/GDF8 gene receive so much attention in several kinds of species such as bovines, poultries, goats, sheep, horses. A 5-base pairs (bp) indel in the 5' untranslated region (5'UTR) of goat MSTN/GDF8 was verified to be significantly associated with growth traits except Inner Mongolia White Cashmere (IMWC) goats. Given that almost all sample sizes were below 150, we enlarged sample sizes to more than 500 to uncover the association between the 5-bp indel and growth traits in IMWC goats. Only two genotypes (deletion/deletion (DD) and insertion/deletion (ID)) were found, and DD genotypes were dominant genotypes. The detected locus displayed low genetic diversity (PIC = 0.090). Interestingly, the association analyses revealed that the 5-bp indel had a significant effect on the chest depth (p = 0.003), and DD genotypes were dominant genotypes. Hinted that the 5-bp indel could act as an effective marker in molecular marker-assisted selection (MAS) processes for selection of excellent goat individuals.

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.

Genetic effects of DSCAML1 identified in genome-wide association study revealing strong associations with litter size and semen quality in goat (Capra hircus)

The down syndrome cell adhesion molecule like 1 (DSCAML1), is associated with the development of the nervous system and neurologic diseases. Previous Genome-wide association studies have shown that it is associated with sperm morphology, suggesting it has a critical role in fecundity. In this study, expression profiles of goat DSCAML1 mRNA were analyzed. The results showed that its expression in the testis differ significantly between the mitotic stage and meiotic stage. Three insertion/deletion (indel) variants of goat DSCAML1 were determined in the Shaanbei White Cashmere Goat (SWCG, n = 2162). Based on the association analysis, two indels (P2-16bp, P14-15bp) were significantly related to sperm quality (sperm motility and sperm density) in male goat and three loci were markedly related to the first-birth litter size in female goat (P = 4.0 × 10^-6; P = 1.0 × 10^-6; P = 4.7 × 10^-2). In male goats, the different genotypes of P2-16bp and P14-15bp revealed a noticeable effect on the expression of DSCAML1. Moreover, the effects observed in the first-birth litter followed a similar trend, which may provide the basis for further research of DSCAML1 gene function and marker assisted selection (MAS) programs to improve reproductive traits.

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.

Correlation analysis of candidate gene SNP for high-yield in Liaoning cashmere goats with litter size and cashmere performance

This study was designed to identify the relationship of four genes (GDF9, BMPR-IB, FecB and ESR) polymorphisms in the 3'UTR region with litter size and cashmere performance of Liaoning cashmere goats (LCG, n = 1140). The ESR C463T and T575G loci of LCG were genotyped. The results of correlation analysis showed that five effective single nucleotide polymorphisms (SNPs) loci (C47T, C94T, C299T, C463T and T575G) were found in the four genes. The lambing number of CC and CT genotypic individuals at FecB C94T locus was significantly higher than that of TT genotypic individuals (45.7 and 46.8%, respectively); the lambing number of CC genotypic individuals at ESR C463T locus was significantly higher than that of CT, TT genotypic individuals (9 and 15%, respectively); There was a positive correlation between CC genotype at C463T locus and cashmere fineness. In this study, the relationship between FecB C94T and ESR C463T loci C alleles and lambing number in LCG was preliminarily revealed. These results further confirmed that FecB and ESR genes may be significantly correlated with high fecundity of LCG.

Recent advances in breeding and genetics for dairy goats

Goats (Capra hircus) were domesticated during the late Neolithic, approximately 10,500 years ago, and humans exerted minor selection pressure until fairly recently. Probably the largest genetic change occurring over the millennia happened via natural selection and random genetic drift, the latter causing genes to be fixed in small and isolated populations. Recent human-influenced genetic changes have occurred through biometrics and genomics. For the most part, biometrics has concentrated upon the refining of estimates of heritabilities and genetic correlations. Heritabilities are instrumental in the calculation of estimated breeding values and genetic correlations are necessary in the construction of selection indices that account for changes in multiple traits under selection at one time. Early genomic studies focused upon microsatellite markers, which are short tandem repeats of nucleic acids and which are detected using polymerase chain reaction primers flanking the microsatellite. Microsatellite markers have been very important in parentage verification, which can impact genetic progress. Additionally, microsatellite markers have been a useful tool in assessing genetic diversity between and among breeds, which is important in the conservation of minor breeds. Single nucleotide polymorphisms are a new genomic tool that have refined classical BLUP methodology (biometric) to provide more accurate genomic estimated breeding values, provided a large reference population is available.

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.

Relationship between SNPs of POU1F1 Gene and Litter Size and Growth Traits in Shaanbei White Cashmere Goats

POU (Pit-Oct-Unc) class 1 homeobox 1 (POU1F1, or Pit-1) is a transcription factor that directly regulates pituitary hormone-related genes, as well as affects the reproduction and growth in mammals. Thus, POU1F1 gene was investigated as a candidate gene for litter size and growth performance in goats. In the current study, using direct DNA sequencing, c.682G > T, c.723T > G and c.837T > C loci were genotyped in Shaanbei white cashmere (SBWC) goats (n = 609), but c.876 + 110T > C was monomorphic. Besides, the c.682G > T locus was first identified by HinfI (Haemophilus influenzae Rf) restriction endonuclease. Association analysis results showed that the c.682G > T, c.837T > C loci and diplotypes were significantly associated with goat litter size (p < 0.05). The positive genotypes were GT and TT for the two SNPs, respectively, and the optimal diplotype was H3H7 (GTTT-TTTT). On the other hand, the c.682G > T, c.723T > G and c.837T > C strongly affected growth traits and body measurement indexes in SBWC goats (p < 0.05). The positive genotypes or allele of these SNPs were GT, G and TT, respectively. Additionally, the goats with H3H7 diplotype also had a greater growth status than others (p < 0.05). Here, individuals with same genotype had both a better litter size and growth traits, showing a positive correlation between these economic traits. Meanwhile, the positive genotypes of four SNPs were combined to obtain the optimal diplotype, which was also H3H7. These SNPs, especially the diplotype, could be used for the genomic selection of excellent individuals with a greater litter size and better growth status in goat breeding.

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.

Two strongly linked single nucleotide polymorphisms (Q320P and V397I) in GDF9 gene are associated with litter size in cashmere goats

Growth differentiation factor 9 (GDF9) gene is an effective intra-ovarian regulator; it plays a crucial role in early folliculogenesis in female mammals. The non synonymous mutations: g.3905A > C (also known as p.Gln320Pro/Q320P) and g.4135G > A (also know as p.Val397Ile/V397I), are two well-known and controversial single nucleotide polymorphisms (SNPs) within GDF9 gene in goats with different prolificacy, and so far, there were no studies on linkage between Q320P and V397I. Therefore, the aim of this work was to study whether Q320P and V397I mutations have a significant effect on litter size, in Shaanbei white cashmere goats (SBWC, n = 1511), and to explore the specific relationship between these two SNPs. The results showed that both of Q320P and V397I mutations exhibited three genotypes; the minor allele frequencies (MAF) of the SNPs were 0.286 and 0.477, respectively; and these two SNPs were in strong linkage disequilibrium (D' = 0.976, r² = 0.348) in the studied goats. Moreover, association analyses revealed that Q320P was significantly associated with the first-born litter size in goats irrespective of the sample size (n = 1511; P = 0.008), while V397I significantly affected litter size until the sample size crossed 1300 (P = 0.015). Meanwhile, the diplotypes PP-II and QP-VI were observed to have a superior effect on litter size (P = 3.78 × 10^-5) to that of the haplotypes (P = 1.12 × 10^-7). Thus, the findings led us to assume that Q320P mutation was the major SNP affecting goat litter size. These findings can provide useful DNA markers for selecting superior individuals in marker-assisted selection (MAS) for breeding in relation to fecundity in goats.

The novel 22 bp insertion mutation in a promoter region of the <i>PITX2</i> gene is associated with litter size and growth traits in goats

Abstract. The paired-like homeodomain 2 (PITX2) gene plays a critical role in regulating development, reproduction, and growth traits in ruminants. Hence, the objective of this study was to explore the polymorphisms of this gene and to evaluate their associations with quantitative traits. Herein, a novel insertion in the promoter region of the PITX2 gene was reported in Shaanbei white cashmere (SBWC) goats (n=1012). The genotype distributions between mothers of single-kid and multi-kid groups within SBWC goats were significantly different (P<0.01), implying that this indel mutation might affect the litter size. Furthermore, association analysis found that this indel mutation was significantly associated with litter size (P=0.001). Individuals with genotype DD had a significantly smaller litter size than those with other genotypes (P<0.01). Besides, this indel was significantly associated with the body length (P=0.042) and the chest width (P=0.031). Especially, the individuals with genotype DD had a significantly lower body length than those with genotype II (P<0.05), which was consistent with the trend in litter size. These findings suggested that the new 22 bp indel mutation within the PITX2 gene is significantly associated with litter size and growth traits; this can be utilized as a functional molecular marker in goat breeding.

Detection of insertion/deletions (indels) of the <i>ATBF1</i> gene and their effects on growth-related traits in three indigenous goat breeds

Abstract. The AT motif-binding factor (ATBF1), also known as zinc finger homeobox 3 (ZFXH3), is necessary for activating the POU1F1 gene; thus, the ATBF1 gene greatly affects the grow traits in animals. The objective of this work was to explore novel indel (insertion/deletion) variations and their associations with growth traits in three native Chinese goat breeds. Two indels within the ATBF1 gene were found in the Shaanbei white cashmere goat (SWCG; n=581), the Guanzhong dairy goat (GZDG; n=334) and the Hainan black goat (HNBG; n=270) for the first time using 12 pairs of primers. Association analysis revealed that the P1-12-bp indel was consistently correlated with the body height of the three breeds, and individuals with ID (insertion/deletion) and DD (deletion/deletion) genotypes had a higher body weight than the II (insertion/insertion) genotype (P=0.036); the P11-6-bp indel was consistently correlated with chest circumference and hip width of the three breeds. Moreover, these two loci were associated with other several growth-related traits in different breeds. Hence, these findings indicated that the goat ATBF1 gene had marked effects on growth traits and the growth-trait-related loci, which would contribute to improving the growth-related traits of local breeds in the goat industry by implementing marker-assisted selection (MAS).