Two mutations in the 5'-flanking region of the KITLG gene are associated with litter size of dairy goats

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.

Single-nucleotide polymorphisms in <i>FLT3</i>, <i>NLRP5</i>, and <i>TGIF1</i> are associated with litter size in Small-tailed Han sheep

Previous studies have indicated that FLT3, NLRP5, and TGIF1 play a pivotal role in sheep fecundity. Nevertheless, little is known about the association of the polymorphisms of these genes with litter size (LS). In this study, the selected single-nucleotide polymorphisms (SNPs) were genotyped using a Sequenom MassARRAY® platform, and the distribution of different genotypes of the SNPs in the seven sheep breeds (Small-tailed Han, Hu, Cele Black, Suffolk, Tan, Prairie Tibetan, and Sunite sheep) were analyzed. The reliability of the estimated allele frequency for all seven SNPs was at least 0.9545. Given the association of the TGIF1 g.37866222C > T polymorphism with LS in Small-tailed Han sheep (p<0.05), fecundity differences might be caused by the change in amino acid from proline (Pro) to serine (Ser), which has an impact on secondary, tertiary protein structures with concomitant TGIF1 functionality changes. The FLT3 rs421947730 locus has a great effect on the LS (p<0.05), indicating that the locus of FLT3 in synergy with KILTG is likely to facilitate ovarian follicle maturation and ovulation. Moreover, NLRP5 rs426897754 is associated with the LS of the second and third parities (p<0.05). We speculate that a synonymous variant of NLRP5 may be involved in folliculogenesis accompanied by BMP15, FSHR, BMPR1B, AMH, and GDF9, resulting in the different fecundity of Small-tailed Han sheep. Our studies provide valuable genetic markers for sheep breeding.

Whole-genome sequence analysis unveils different origins of European and Asiatic mouflon and domestication-related genes in sheep

The domestication and subsequent development of sheep are crucial events in the history of human civilization and the agricultural revolution. However, the impact of interspecific introgression on the genomic regions under domestication and subsequent selection remains unclear. Here, we analyze the whole genomes of domestic sheep and their wild relative species. We found introgression from wild sheep such as the snow sheep and its American relatives (bighorn and thinhorn sheep) into urial, Asiatic and European mouflons. We observed independent events of adaptive introgression from wild sheep into the Asiatic and European mouflons, as well as shared introgressed regions from both snow sheep and argali into Asiatic mouflon before or during the domestication process. We revealed European mouflons might arise through hybridization events between a now extinct sheep in Europe and feral domesticated sheep around 6000-5000 years BP. We also unveiled later introgressions from wild sheep to their sympatric domestic sheep after domestication. Several of the introgression events contain loci with candidate domestication genes (e.g., PAPPA2, NR6A1, SH3GL3, RFX3 and CAMK4), associated with morphological, immune, reproduction or production traits (wool/meat/milk). We also detected introgression events that introduced genes related to nervous response (NEURL1), neurogenesis (PRUNE2), hearing ability (USH2A), and placental viability (PAG11 and PAG3) into domestic sheep and their ancestral wild species from other wild species.

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.

Combined approaches identify known and novel genes associated with sheep litter size and non-seasonal breeding

Improvement of ewe reproduction is considerable by appropriately increasing litter size and sustaining non-seasonal breeding. However, their genetic makeups have not been entirely elucidated. Genome-wide analyses of 821 individuals were performed by combining three genomic approaches (genome-wide association study, XP-nSL, and runs of homozygosity). Consequently, 35 candidate genes including three domestication genes (TSHR, GTF2A1, and KITLG) were identified. Other than the FecB mutation at BMPR1B, we described a significant association of a missense mutation rs406686139 at seasonal lambing-associated TSHR gene with litter size. Some promising novel genes may be relevant for sheep reproduction by multitude biological processes, such as FETUB functioning in fertilization, HNRNPA1 in oogenesis, DCUN1D1 in spermatogenesis, and HRG in fertility outcome. The present study suggests that improvement of ewe reproduction is attributed to selective breeding, and casts light on the genetic basis and improvement of sheep reproduction.

Identification of candidate genomic regions for chicken egg number traits based on genome-wide association study

Background

Since the domestication of chicken, various breeds have been developed for food production, entertainment, and so on. Compared to indigenous chicken breeds which generally do not show elite production performance, commercial breeds or lines are selected intensely for meat or egg production. In the present study, in order to understand the molecular mechanisms underlying the dramatic differences of egg number between commercial egg-type chickens and indigenous chickens, we performed a genome-wide association study (GWAS) in a mixed linear model.

Results

We obtained 148 single nucleotide polymorphisms (SNPs) associated with egg number traits (57 significantly, 91 suggestively). Among them, 4 SNPs overlapped with previously reported quantitative trait loci (QTL), including 2 for egg production and 2 for reproductive traits. Furthermore, we identified 32 candidate genes based on the function of the screened genes. These genes were found to be mainly involved in regulating hormones, playing a role in the formation, growth, and development of follicles, and in the development of the reproductive system. Some genes such as NELL2 (neural EGFL like 2), KITLG (KIT ligand), GHRHR (Growth hormone releasing hormone receptor), NCOA1 (Nuclear receptor coactivator 1), ITPR1 (inositol 1, 4, 5-trisphosphate receptor type 1), GAMT (guanidinoacetate N-methyltransferase), and CAMK4 (calcium/calmodulin-dependent protein kinase IV) deserve our attention and further study since they have been reported to be closely related to egg production, egg number and reproductive traits. In addition, the most significant genomic region obtained in this study was located at 48.61–48.84 Mb on GGA5. In this region, we have repeatedly identified four genes, in which YY1 (YY1 transcription factor) and WDR25 (WD repeat domain 25) have been shown to be related to oocytes and reproductive tissues, respectively, which implies that this region may be a candidate region underlying egg number traits.

Conclusion

Our study utilized the genomic information from various chicken breeds or populations differed in the average annual egg number to understand the molecular genetic mechanisms involved in egg number traits. We identified a series of SNPs, candidate genes, or genomic regions that associated with egg number, which could help us in developing the egg production trait in chickens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07755-3.

Genomic Signatures of Selection Associated With Litter Size Trait in Jining Gray Goat

Litter size (LS), an important economic trait in livestock, is so complicate that involves many aspects of reproduction, the underlying mechanism of which particularly in goat has always been scanty. To uncover the genetic basis of LS, the genomic sequence of Jining Gray goat groups (one famous breed for high prolificacy in China) with LS 1, 2, and 3 for firstborn was analyzed, obtaining 563.67 Gb sequence data and a total of 31,864,651 high-quality single nucleotide polymorphisms loci were identified. Particularly, the increased heterozygosity in higher LS groups, and large continuous homozygous segments associated with lower LS group had been uncovered. Through an integrated analysis of three popular methods for detecting selective sweeps (Fst, nucleotide diversity, and Tajima’s D statistic), 111 selected regions and 42 genes associated with LS were scanned genome wide. The candidate genes with highest selective signatures included KIT, KCNH7, and KMT2E in LS2 and PAK1, PRKAA1, and SMAD9 in LS3 group, respectively. Meanwhile, functional terms of programmed cell death involved in cell development and regulation of insulin receptor signaling pathway were mostly enriched with 42 candidate genes, which also included reproduction related terms of steroid metabolic process and cellular response to hormone stimulus. In conclusion, our study identified novel candidate genes involving in regulation of LS in goat, which expand our understanding of genetic fundament of reproductive ability, and the novel insights regarding to LS would be potentially applied to improve reproductive performance.

Goat CMTM2: mRNA expression profiles of different alternative spliced variants and associations analyses with growth traits

CKLF like MARVEL transmembrane domain containing 2 (CMTM2) plays crucial roles in spermiogenesis, skeletogenous, growth, and development through PI3K/Akt and other pathways. The purpose of this study was to explore the expression profile and variation of different spliced CMTM2 gene in Shaanbei white cashmere goats, as well as to find the relationships between a CMTM2 promoter region 14 bp genetic variant and growth traits in 1366 Shaanbei white cashmere goats. In this study, we identified alternative CMTM2 splicing and detected the effects of the spliced variants on mRNA expression levels in tissues. Meanwhile, an unreported spliced variant of CMTM2 in goat was identified using in CDS cloning and RT-PCR, namely, CMTM2-AS2. Compared with the normal transcript (CMTM2-AS1), the novel variant had the higher expression level in muscle and liver tissues, indicating that it plays an effective role in growth traits. Furthermore, a 14 bp deletion was detected within CMTM2 promoter region, and the different genotypes were significantly associated with growth traits (e.g., body length, circumference of cannon bone) in the large group of 1366 individuals in Shaanbei white cashmere goats. We found that the body length of the individuals with II (n = 571) genotype had better phenotypes than those with DD (n = 118) and ID (n = 650) genotypes. These results have direct guiding significance for goat breeding in the future and provide a new idea for studying the characteristics and functions of CMTM2 gene in goats.

Association of Polymorphisms in Candidate Genes with the Litter Size in Two Sheep Breeds

Hu sheep and Small-tailed Han sheep are the most widely raised and most famous maternal sheep breeds in China, which are known for precocious puberty, perennial oestrus and high fecundity (1-6 lambs each parity). Therefore, it is crucial to increase litter size of these two breeds for intensive sheep industry. The objective of this study was to identify potential genetic markers linked with sheep litter size located at ten genes. This study collected blood sample of 537 Hu sheep and 420 Small-tailed Han sheep with litter size of first parity. The average litter sizes in Hu sheep and Small-tailed Han sheep were 2.21 and 1.93. DNA-pooling sequencing method was used for detecting the potential single nucleotide polymorphisms (SNPs) in ten genes related to follicle development and female reproduction. SNPscan® was used for individually genotyping. As a result, a total of 78 putative SNPs in nine out of ten candidate genes (except NOG) were identified. In total, 50 SNPs were successfully genotyped in Hu sheep and Small-tailed Han sheep. After quality control, a total of 42 SNPs in Hu sheep and 44 SNPs in Small-tailed Han sheep were finally used for further analysis. Association analysis revealed that nine SNPs within six genes (KIT: g.70199073A>G, KITLG: g.124520653G>C, ADAMTS1: g.127753565T>C, ADAMTS1: g.127754640G>T, NCOA1: g.31928165C>T, NCOA1: g.32140565G>A, LIFR: g.35862868C>T, LIFR: g.35862947G>T and NGF: g.91795933T>C) were significantly associated with litter size in Hu sheep or Small-tailed Han sheep. A combined haplotypes analysis of the two loci (LIFR: g.35862868C>T and LIFR: g.35862947G>T) revealed that H2H3 (CTTT) combined haplotypes had the largest litter size than the rest combined haplotypes and more than those with either mutation alone in Small-tailed Han sheep. Taken together, our study suggests that nine significant SNPs in six genes can be served as useful genetic markers for MAS in sheep.

Convergent genomic signatures of domestication in sheep and goats

The evolutionary basis of domestication has been a longstanding question and its genetic architecture is becoming more tractable as more domestic species become genome-enabled. Before becoming established worldwide, sheep and goats were domesticated in the fertile crescent 10,500 years before present (YBP) where their wild relatives remain. Here we sequence the genomes of wild Asiatic mouflon and Bezoar ibex in the sheep and goat domestication center and compare their genomes with that of domestics from local, traditional, and improved breeds. Among the genomic regions carrying selective sweeps differentiating domestic breeds from wild populations, which are associated among others to genes involved in nervous system, immunity and productivity traits, 20 are common to Capra and Ovis. The patterns of selection vary between species, suggesting that while common targets of selection related to domestication and improvement exist, different solutions have arisen to achieve similar phenotypic end-points within these closely related livestock species.

The sheep and goat were domesticated ~10,500 years ago in the same region of the Middle-East. Here, Alberto et al compare the genomes of wild Asiatic mouflon and Bezoar ibex with that of domestics from local, traditional and improved breeds and find common targets of selection related to domestication and improvement in sheep and goats.

Analysis of the SNP loci around transcription start sites related to goat fecundity trait base on whole genome resequencing

Genome flanking regions surrounding transcription start sites (TSSs) are critical for the regulation of gene expression, containing many translational regulatory elements. To investigate whether critical single nucleotide polymorphisms (SNPs) exist around TSSs in the dairy goat genome, we performed high throughput DNA sequencing to compare two dairy goat groups with discrepant litter sizes. After genome mapping, SNP calling, and annotation, we screened the SNPs within 2kb scales surrounding annotated TSSs in high fecundity (HF) and low fecundity (LF) groups, respectively. We attempted to identify distinct SNPs and motifs near the TSSs in both groups. The SNPs near the TSSs most were consistent; 318 new SNPs were uncovered in the HF group, of which 305 were heterozygote SNPs, 13 were homozygote SNPs, and majority of which were distributed on chromosome 2 and 29. After validation by Sanger sequencing we found that a SNP in CHI16: 27612330 C>A in the PSEN2 gene presented an A/A genotype in the HF group and an A/A or A/C genotype in the LF group. In conclusion, our study provides insightful information into the dairy goat genomic variations surrounding TSSs, which may contribute to enhanced litter size. Based on comparison studies of SNPs exist around transcription start sites between high fecundity group and low fecundity group. Our finding provides insights concerning the goat litter size phenotypic and will promote future goat breeding.

Association of polymorphisms at the microRNA binding site of the caprine KITLG 3'-UTR with litter size

This study identified three novel single nucleotide polymorphisms (SNPs) (c.1389C > T, c.1457A > C and c.1520G > A) in the caprine KITLG 3'-UTR through DNA sequencing. The three SNP loci were closely linked in Guanzhong dairy (GD) goats. Two alleles of the c.1457A > C SNP introduced two miRNA sites (chi-miR-204-5p and chi-miR-211). Individuals with combined genotype TT-CC-AA had a higher litter size compared with those with combined genotypes CC-AA-GG, TC-CC-GA and CC-AC-GG (P < 0.05). Luciferase assays showed that chi-miR-204-5p and chi-miR-211 suppressed luciferase expression in the presence of allele 1457A compared with negative control (NC) and allele 1457C (P < 0.05). Western blot revealed that KITLG significantly decreased in the granulosa cells (GCs) of genotype AA compared with that in the GCs of genotype CC and NC (P < 0.05). The KITLG mRNA levels of the CC-AA-GG carriers significantly decreased compared with those of the TT-CC-AA, TC-CC-GA and CC-AC-GG carriers. In addition, cell proliferation was reduced in haplotype C-A-G GCs compared with that in haplotype T-C-A GCs. These results suggest that SNPs c.1389C > T, c.1457A > C and c.1520G > A account for differences in the litter size of GD goats because chi-miR-204-5p and chi-miR-211 could change the expression levels of the KITLG gene and reduce GC proliferation.