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

Potential Candidate Genes Associated with Litter Size in Goats: A Review

This review examines genetic markers associated with litter size in goats, a key reproductive trait impacting productivity in small ruminant farming. Goats play a vital socioeconomic role in both low- and high-income regions; however, their productivity remains limited due to low reproductive efficiency. Litter size, influenced by multiple genes and environmental factors, directly affects farm profitability and sustainability by increasing the output per breeding cycle. Recent advancements in genetic research have identified key genes and pathways associated with reproductive traits, including gonadotropin-releasing hormone (GnRH), inhibin (INHAA), Kit ligand (KITLG), protein phosphatase 3 catalytic subunit alpha (PPP3CA), prolactin receptor (PRLR), POU domain class 1 transcription factor 1 (POU1F1), anti-Müllerian hormone (AMH), bone morphogenetic proteins (BMP), growth differentiation factor 9 (GDF9), and KISS1 and suppressor of mothers against decapentaplegic (SMAD) family genes, among others. These genes regulate crucial physiological processes such as folliculogenesis, hormone synthesis, and ovulation. Genome-wide association studies (GWASs) and transcriptomic analyses have pinpointed specific genes linked to increased litter size, highlighting their potential in selective breeding programs. By incorporating genomic data, breeding strategies can achieve higher selection accuracy, accelerate genetic gains, and improve reproductive efficiency. This review emphasizes the importance of genetic markers in optimizing litter size and promoting sustainable productivity in goat farming.

Association analysis for SNPs of LIPE and ITGB4 genes with cashmere production performance, body measurement traits and milk production traits in Liaoning cashmere goats

Liaoning cashmere goat (LCG) is one of the excellent cashmere goat breeds in China. Because of its larger size, better cashmere, and better cashmere production performance, people pay special attention to it. This article mainly studied the relationship between SNP loci of LIPE gene and ITGB4 gene and milk production, cashmere production and body measurement traits of LCGs. We further identified potential SNP loci by PCR-Seq polymorphism detection and gene sequence comparison of LIPE and ITGB4 genes. Further, we use SPSS and SHEsis software to analyze their relationship to production performance. The consequence indicated that CC genotype of LIPE gene T16409C locus was dominant genotype in milk production and cashmere production, while CT genotype of LIPE gene T16409C locus was dominant in body size. The CT genotype of C168T locus of ITGB4 gene is the dominant genotype of body type and cashmere production, while the dominant genotype of milk production is TT genotype. Through joint analysis, in haploid combinations, H1H2:CCCT is the dominant haplotype combination in cashmere fineness. H3H4:TTCT is a dominant haplotype combination of milk production traits and body measurement traits. These dominant genotypes can provide a reliable basis for the study of production performance of LCG.

Association between the cashmere production performance, milk production performance, and body size traits and polymorphism of COL6A5 and LOC102181374 genes in Liaoning cashmere goats

The purpose of this study was to analyze the relationship between COL6A5 (collagen type VI alpha 5 chain) and LOC102181374 (alcohol dehydrogenase 1) genes and the production performance of Liaoning cashmere goats by single nucleotide polymorphism (SNP). We have searched for SNP loci of COL6A5 and LOC102181374 genes through sequence alignment and PCR experiments, and have used SPSS and SHEsis software to analyze production data. We obtained five SNP loci in total, including three SNP loci (G50985A, G51140T, G51175A) in COL6A5 gene and two SNP loci (A10067G, T10108C) in LOC102181374 gene. The genotypes G50985A (AG), G51140T (GT), G51175A (AA), A10067G (AA), and T10108C (CC) of these loci have certain advantages in improving the production performance of Liaoning cashmere goats. The haplotype combinations that can improve production performance in COL6A5 gene were H1H5:AGGGAG, H4H4:GGGGAA, and H4H4:GGGGAA. H3H3:GGCC and H2H4:AGTT were the dominant combinations in LOC102181374 gene. At G51175A and A10067G loci, we found that H1H2:AAAG and H1H3:AGAA have dominant effects. These results may provide some support for the molecular breeding of production traits in Liaoning cashmere goats.

Association analysis for SNPs of BAAT and COL1A1 genes with cashmere production performance and other production traits in Liaoning cashmere goats

This study aimed to investigate the relationship between the polymorphism of bile acid-CoA: amino acid N-acyltransferase (BAAT) and collagen type I alpha 1 chain (COL1A1) genes and the production performance of Liaoning Cashmere goat (LCG). The potential single nucleotide polymorphisms (SNPs) of LCG were detected by sequence comparison of BAAT and COL1A1 genes and PCR-Seq polymorphism, and the effect of SNPs on production performance was analyzed by SPSS software. The results showed that three SNPs loci were detected in BAAT gene: G7900A, T7967C, C7998T, and one SNP locus T6716C was detected in COL1AL gene. At G7900A locus, the dominant genotype for cashmere performance was GG, and the dominant genotype for body measurement traits and milk production traits was AG. At T7967C locus, the dominant genotype for cashmere performance was TT, and the dominant genotype for body measurement traits and milk production traits was CC. At C7998T locus, TT was the dominant genotype for cashmere performance, body measurement traits, and milk production traits. At the T6716C locus, TT was the dominant genotype for cashmere performance, body measurement traits, and milk production traits. H1H1: AACC is the dominant haplotype combination. Therefore, this study will provide a reliable reference for future research on cashmere production performance, body measurement traits, and milk production traits of LCG.

Association analysis for SNPs of MSTN and IGFBP-3 genes with body size and other production traits in Liaoning Cashmere Goats

Liaoning cashmere goat (LCG) have tall bones, high cashmere production and outstanding meat production performance. In recent years, good breeding progress has not been made in terms of body size, meat yield, milk yield and other properties in terms of production. The study focused on the correlation between the SNPs of MSTN and IGFBP-3 genes with the body size performance, cashmere production and milk performance. The MSTN and IGFBP-3 gene sequence alignment and PCR-Seq polymorphism were used to detect the potential SNPs, and the correlation with production performance was analyzed by SPSS and SHEsis software. The results showed that the TT genotype at the T1662G locus of the MSTN gene is dominant and has significant advantages in body measurements such as sacrum height, chest width, and waist height. The C allele at the C4021T locus of IGFBP-3 gene shows an advantage in the body measurement performance. Among the haplotype combinations, H2H2:TGTC is preponderant combination for body size performance, H2H2:TGTC and H1H2:TGCC are preponderant combinations for cashmere production performance, H1H3:GGCC is preponderant combination for milk production performance. It may be a molecular marker for future selection and breeding.

Combined analysis of mRNA-miRNA from testis tissue in Tibetan sheep with different FecB genotypes

Testis size is important for identifying breeding animals with adequate sperm production. The aim of this study was to survey the expression profile of mRNA and miRNA in testis tissue from rams carrying different FecB genotypes, including the wild-type and heterozygous genotypes in Tibetan sheep. Comparative transcriptome profiles for ovine testes were established for wild-type and heterozygote Tibetan sheep by next-generation sequencing. RNA-seq results identified 3,910 (2,034 up- and 1,876 downregulated) differentially expressed (DE) genes and 243 (158 up- and 85 downregulated) DE microRNAs (miRNAs) in wild-type vs heterozygote sheep, respectively. Combined analysis of mRNA-seq and miRNA-seq revealed that 20 miRNAs interacted with 48 true DE target genes in wild-type testes compared to heterozygous genotype testes. These results provide evidence for a functional series of genes operating in Tibetan sheep testis. In addition, quantitative real-time PCR analysis showed that the expression trends of randomly selected DE genes in testis tissues from different genotypes were consistent with high-throughput sequencing results.

Analysis of the Association of Two SNPs in the Promoter Regions of the PPP2R5C and SLC39A5 Genes with Litter Size in Yunshang Black Goats

Screening for candidate genes and genetic variants associated with litter size is important for goat breeding. The aim of this study was to analyze the relationship between single nucleotide polymorphisms (SNPs) in PPP2R5C and SLC39A5 and litter size in Yunshang black goats. KASP genotyping was used to detect the SNP genetic markers in the PPP2R5C and SLC39A5 in a population of 569 Yunshang black goats. The results show that there were two SNPs in the PPP2R5C and SLC39A5 promoter regions. Association analysis revealed that the polymorphisms PPP2R5C g.65977743C>T and SLC39A5 g.50676693T>C were significantly associated with the litter size of the third parity of Yunshang black goats (p < 0.05). To further explore the regulatory mechanism of the two genes, the expression of different genotypes of PPP2R5C and SLC39A5 was validated by RT-qPCR and Western blotting. The expression of PPP2R5C was significantly higher in individuals with the TT genotype than in those with the TC and CC genotypes (p < 0.05). The expression of SLC39A5 was also significantly higher in individuals with the TT genotype than in TC and CC genotypes (p < 0.05). Dual luciferase reporter analysis showed that the luciferase activity of PPP2R5C-C variant was significantly higher than that of PPP2R5C-T variant (p < 0.05). The luciferase activity of SLC39A5-T variant was significantly higher than that of SLC39A5-C variant (p < 0.05). Software was used to predict the binding of transcription factors to the polymorphic sites, and the results show that SOX18, ZNF418, and ZNF667 and NKX2-4 and TBX6 might bind to PPP2R5C g.65977743C>T and SLC39A5 g.50676693T>C, respectively. These results provide new insights into the identification of candidate genes for marker-assisted selection (MAS) in goats.

Are Copy Number Variations within the FecB Gene Significantly Associated with Morphometric Traits in Goats?

The Booroola fecundity (FecB) gene is a major fertility-related gene first identified in Booroola sheep. Numerous studies have investigated whether the FecB gene is a major fecundity gene in goats or whether there are other genes that play a critical role in goat fertility. Nevertheless, little attention has been paid to the role of the FecB gene in the body morphometric traits of goats, despite the positive relationship discerned between litter size and growth. We identified five copy number variations (CNVs) within the FecB gene in 641 goats, including 318 Shaanbei white cashmere (SBWC) goats, 203 Guizhou Heima (GZHM) goats, and 120 Nubian goats, which exhibited different distributions among these populations. Our results revealed that these five CNVs were significantly associated with goat morphometric traits (p < 0.05). The normal type of CNV3 was the dominant type and displayed superior phenotypes in both litter size and morphometric traits, making it an effective marker for goat breeding. Consequently, LD blocks in the region of 10 Mb upstream and downstream from FecB and potential transcription factors (TFs) that could bind with the CNVs were analyzed via bioinformatics. Although no significant LD block was detected, our results illustrated that these CNVs could bind to growth-related TFs and indirectly affect the growth development of the goats. We identified potential markers to promote litter size and growth, and we offer a theoretical foundation for further breeding work.

Association analysis of single-nucleotide polymorphism in prolactin and its receptor with productive and body conformation traits in Liaoning cashmere goats

The results of this study showed that the single-nucleotide polymorphism (SNP) sites of the PRL and PRLR genes have a certain association with the milk production performance, body size and cashmere performance of Liaoning cashmere goats (LCGs). Through our designed experiment, the potential SNPs of LCG were detected by sequence alignment, and two SNPs were found on two genes. The CC genotype of the PRL gene is the dominant genotype among the three genotypes. The GG genotype of the PRLR gene is the dominant genotype among the two genotypes. At the same time, the two genotypes also have good performance in cashmere production and body size. Through the screening of haplotype combination, the milk fat rate > 7.6 %, the milk protein rate > 5.6 %, the milk somatic cell number < 1500 × 103 mL-1, the cashmere fineness < 15.75 µm, the chest girth > 105 cm, the chest depth > 33 cm, and the waist height > 67.5 cm are considered as screening indexes for comprehensive production performance of Liaoning cashmere goats. It is concluded that the GCGC type is the dominant haplotype combination. According to our research data, we found that the biological indicators of Liaoning cashmere goat milk are higher than the national standards, so we think it is very significant to study the milk production performance of our experiment. Further research can be done on goat milk production and body conformation traits around PRL gene and PRLR gene.

Association analysis for SNPs of KRT26 and TCHH genes with cashmere production performance, body measurement traits and milk production traits in Liaoning cashmere goats

Cashmere fineness is getting thicker, which is one of the key problems in cashmere breeding, however, there have been no systematic studies on the molecular regulation of cashmere fineness. The aim of this study was to investigate the relationship between KRT26 and TCHH gene polymorphism and production performance in Liaoning cashmere goats (LCG). The potential single nucleotide polymorphisms (SNPs) of LCG were detected by sequence alignment and PCR-Seq polymorphism of KRT26 and TCHH genes and analyzed the effect of SNPs on production performance by SPSS software. Two SNPs sites (A559T and A6839G) of two genes were detected. The AA genotype of KRT26 A559T locus was the dominant genotype. AG and GG at TCHH A6839G locus were the dominant genotypes. AAAA was the dominant haplotype combination. The results showed that KRT26 and TCHH genes were associated with cashmere fineness of LCG, and A559T (AA) and A6839G (GG) genotypes were the preferred marker genotypes for cashmere fineness, which provided more theoretical basis for further research on cashmere fineness.

Detection of mRNA Expression and Copy Number Variations Within the Goat Fec B Gene Associated With Litter Size

The Booroola fecundity (Fec ^B ) gene, as the first major fecundity gene identified in Booroola sheep, has attracted careful attention. So far, previous research have uncovered the Fec^B mutation (Q249R) as the main mutation by virtue of which sheep exhibits multiple lambing phenomena. This mutation is now being intensively studied and widely used. However, such effect of the Fec^B mutation has not been applied to goats, and similar types of the Fec ^B gene in goats still need to be studied. Thus, the current study attempted to verify potential mutations in the goat Fec ^B gene as well as investigate their functions related to fecundity. First, Fec ^B expression was investigated in six different goat tissues, and we found that Fec ^B expression was highest in the mammary gland, followed by the ovary. Next, the influence of the Fec ^B gene was analyzed from a new perspective, where five potential copy number variations (CNVs) (CNV1-5) within the Fec ^B gene were identified for the first time, and then their effects on litter size were measured. Our results point out that CNV3 (P = 3.44E-4) and CNV5 (P = 0.034) could significantly influence the litter size of goats. Identically, the combination genotype of CNV3 and CNV5 which consisted of their dominant genotypes was also significantly associated with goat litter size (P = 7.80E-5). Hence, CNV3 and CNV5 could serve as potential DNA molecular markers applied to DNA editing and DNA microarray. Additionally, the abovementioned study has laid a theoretical foundation for the detection of potential fertility-related quantitative trait loci within the goat Fec ^B gene.

Polymorphism of fecundity genes (BMP15 and GDF9) and their association with litter size in Bangladeshi prolific Black Bengal goat

Goat farming in Bangladesh is primarily centred on indigenous Black Bengal goat, a highly prolific breed. Searching for genetic markers associated with prolificacy in this breed is vital for the country's goat breeding industry. However, there are no reports on polymorphisms associated with the fertility of Bangladeshi Black Bengal goats. This study investigated two major fecundity genes-bone morphogenetic protein 15 (BMP15) and growth differentiation factor 9 (GDF9) to detect any possible mutations in these two genes associated with litter size in Black Bengal goats. Blood samples were collected from 40 raised goats in Hathazari Government Goat Farm, Bangladesh. Genomic DNA was extracted; PCR amplification was performed; and sequencing of PCR products was performed to detect polymorphism loci in the target genes. Five SNPs viz. C735A, C743A, G754T, C781A and C808G were detected in exon 2 of BMP15 gene. A SNP (T1173A) was detected in GDF9 exon 2. Association results show that SNPs at the 735, 754 and 781 nucleotide positions of BMP15 exon 2 had a significant association with litter size in Black Bengal goat. The effect of parity was also highly significant (P < 0.001) on litter size. For the first time, this study explored SNP loci in fecundity genes in Bangladeshi prolific Black Bengal goats. Further studies with many genetically unrelated animals for assessing the association of these loci and others in the fecundity genes with litter size may be useful.

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.