Prolificacy genotypes at BMPR 1B,BMP15 and GDF9 genes in North African sheep breeds

The aim of this research was to investigate the genetic structure at BMPR 1B, BMP15 and GDF9 prolificacy genes in five sheep breeds reared in Tunisia: Barbarine, Queue Fine de L’Ouest, Noire de Thibar, Sicilo-Sarde and D’man. Genomic DNA of 204 sheep was investigated for the FecB^B (BMPR 1B), FecX^R, FecX^H, FecX^I, FecX^L, FecX^G, FecX^B (BMP15) and FecG^H (GDF9) mutations. The sequence variability of the different DNA fragments utilised for genotyping was further investigated by Single Stranded Conformation Polymorphism (SSCP) and sequencing. All the above-mentioned mutations were absent in the five sheep breeds examined. SSCP analysis and sequencing allowed the detection of two nucleotide variations. A non-functional mutation (T/C transition at nt 747 of BMP15 cDNA known as B3) was found at the BMP15 gene, in the Noire de Thibar breed; this mutation was first detected in the Belclare sheep. A new nucleotide change G/A at nt 1159 of BMP15 cDNA, causing the amino acid change A119T in the mature peptide, was detected in the Barbarine breed for the first time. The highly prolific D’man ewes were monomorphic for the absence of all the known prolificacy alleles.     

Screening of indigenous goats for prolificacy associated DNA markers of sheep

The present study was undertaken to explore the genetic basis of caprine prolificacy and to screen indigenous goats for prolificacy associated markers of sheep in BMPR1B, GDF9 and BMP15 genes. To detect the associated mutations and identify novel allelic variants in the candidate genes, representative samples were collected from the breeding tract of indigenous goat breeds varying in prolificacy and geographic distribution. DNA was extracted and PCR amplification was done using primers designed or available in literature for the coding DNA sequence of candidate genes. Direct sequencing was done to identify the genetic variations. Mutations in the candidate genes associated with fecundity in sheep were not detected in Indian goats. Three non-synonymous SNPs (C818T, A959C and G1189A) were identified in exon 2 of GDF9 gene out of which mutation A959C has been associated with prolificacy in exotic goats. Two novel SNPs (G735A and C808G) were observed in exon 2 of BMP15 gene.     

Polymorphism of fecundity genes (BMPR1B, BMP15 and GDF9) in the Indian prolific Black Bengal goat

The Black Bengal is a prolific goat breed in India. Natural mutations in prolific sheep breeds have shown that the transforming growth factor beta (TGF-β) super family ligands such as growth differentiation factor 9 (GDF9), bone morphogenetic protein 15 (BMP15) and their type I receptor (bone morphogenetic protein receptor, BMPR1B) are crucial for ovulation and as well as for increasing litter size. Mutations in any of these genes increased prolificacy in sheep. Based on the known mutation information in sheep PCR primers were designed to screen known polymorphism in 88 random Black Bengal goats. Only the BMPR1B gene was polymorphic. Three genotypes of animals were detected in tested animals with mutant (FecB^B) and wild type (FecB⁺) alleles were 0.57 and 0.43, respectively. Non-carrier, heterozygous carrier and homozygous carrier Black Bengal does had 2.7, 3.04 and 3.11 kids, respectively. All known point mutations of BMP15 and GDF9 genes were monomorphic in the animals tested. These results preliminarily showed that the BMPR1B gene might be a major gene that influences prolificacy of Black Bengal goats.     

Genome‐wide association study to identify genomic regions affecting prolificacy in Lori‐Bakhtiari sheep

Several causative mutations in candidate genes affecting prolificacy have been detected in various sheep breeds. A genome‐wide association study was performed on estimated breeding values for litter size in Lori‐Bakhtiari sheep. Prolific ewes with twinning records and others with only singleton records were genotyped using the medium‐density Illumina Ovine SNP50 array. Four single nucleotide polymorphisms (SNPs) associated with litter size were identified on chromosomes 3, 6 and 22. The region on sheep chromosome 3 between 75 739 167 and 75 745 152 bp included two significant SNPs (s52383.1 and OAR3_80038014_X.1) in high linkage disequilibrium with each other. The region that surrounds these SNPs contains a novel putative candidate gene: luteinizing hormone/choriogonadotropin receptor (LHCGR), known to be involved in ovarian steroidogenesis and organism‐specific biosystem pathways in sheep. Known prolificacy genes BMPR1B, BMP15 and GDF9 were not associated with litter size in Lori‐Bakhtiari sheep, suggesting that other biological mechanisms could be responsible for the trait's variation in this breed.     

Polymorphism of fecundity genes (FecB, FecX, and FecG) in the Indian Bonpala sheep

The present study was designed for screening polymorphism of known fecundity genes in prolific Indian Bonpala sheep. Employing tetra-primer amplification refractory mutation system PCR, 11-point mutations of BMP1B, BMP15, and GDF9 genes of 97 Bonpala ewes were genotyped. The FecB locus of the BMPR1B gene and two loci (G1 and G4) of GDF9 gene were found to be polymorphic. In FecB locus, three genotypes, namely, wild type (Fec++, 0.02), heterozygous (FecB+, 0.23), and mutant (FecBB, 0.75) were detected. At G1 locus of GDF9 gene, three genotypes, namely, wild type (GG, 0.89), heterozygous (GA, 0.10), and mutant (AA, 0.01) were detected. At G4 locus of GDF9 gene, three genotypes, namely, wild type (AA, 0.01), heterozygous (AG, 0.14), and mutant (GG, 0.85) were detected. Statistically no significant correlation of polymorphism of FecB, G1, and G4 loci and litter size was found in this breed. All five loci of BMP15 and three loci of GDF 9 genes were monomorphic. This study reports Bonpala sheep as the first sheep breed where concurrent polymorphism at three important loci (FecB, G1, and G4) of two different fecundity genes (BMPR1B and GDF9) has been found.     

GDF9 as a candidate gene for prolificacy of Small Tail Han sheep

Growth differentiation factor 9 (GDF9) which controls the fecundity of Belclare, Cambridge, Santa Ines, Moghani, Ghezel and Thoka ewes was studied as a candidate gene for the prolificacy of Small Tail Han sheep. According to the sequence of ovine GDF9 gene, six pairs of primers were designed to detect single nucleotide polymorphisms of two exons of GDF9 gene in both high fecundity breed (Small Tail Han sheep) and low fecundity breed (Dorset sheep) by polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP). Only the products amplified by primers 2-1 and 2-2 displayed polymorphisms. For primer 2-1, three genotypes (AA, AB and BB) were detected in both sheep breeds. Sequencing revealed one silent mutation (G477A) in exon 2 of GDF9 gene in the BB genotype in comparison with the AA, which was known as G3 mutation of GDF9 gene in Belclare and Cambridge ewes. The relationship of least squares means for litter size was AA?>?AB?>?BB in Small Tail Han sheep (P?>?0.05). For primer 2-2, two genotypes (CC and CD) were detected in both sheep breeds. Sequencing revealed one novel single nucleotide mutation (G729T) in exon 2 of GDF9 gene in the CD genotype in comparison with the CC, which resulted in an amino acid change (Gln243His). The ewes with mutation heterozygous genotype CD had 0.77 (P?相似文献   

Patterns of expression of messenger RNAs encoding GDF9, BMP15, TGFBR1, BMPR1B, and BMPR2 during follicular development and characterization of ovarian follicular populations in ewes carrying the Woodlands FecX2W mutation

Woodlands sheep have a putative genetic mutation (FecX2(W)) that increases ovulation rate. At present, the identity of FecX2(W) is unknown. The trait does not appear to be due to the previously described mutations in bone morphogenetic protein 15 (BMP15), growth differentiation factor 9 (GDF9), or bone morphogenetic protein receptor type 1B (BMPR1B) that affect ovulation rate in sheep. Potentially, FecX2(W) could be an unidentified genetic mutation in BMP15 or in the closely related GDF9, which interacts with BMP15 to control ovarian function. Alternatively, FecX2(W) may affect ovulation rate by changing the expression patterns in the molecular pathways activated by genes known to regulate ovulation rate. The objectives of these experiments were to sequence the complete coding region of the BMP15 and GDF9 genes, determine the patterns of expression of mRNAs encoding GDF9, BMP15, TGFBR1, BMPR1B, and BMPR2 during follicular development, and characterize the follicular populations in ewes heterozygous for the Woodlands mutation and their wild-type contemporaries. No differences in the coding sequences of BMP15 or GDF9 genes were identified that were associated with enhanced ovulation rate. The expression patterns of GDF9 and BMPR2 mRNAs were not different between genotypes. However, expression of BMP15 mRNA was less in oocytes of FecX2(W) ewes in large preantral and antral follicles. Expression of ALK5 mRNA was significantly higher in the oocytes of FecX2(W) ewes, whereas expression of BMPR1B was decreased in both oocytes and granulosa cells of FecX2(W) ewes. FecX2(W) ewes also had increased numbers of antral follicles <1 mm in diameter. These follicles were smaller in average diameter, with the oocytes also being of a smaller mean diameter. Given that a mutation in BMP15 or BMPR1B results in increased ovulation rates in sheep, the differences in expression levels of BMP15 and BMPR1B may play a role in the increase in ovulation rate observed in Woodlands ewes with the FecX2(W) mutation.     

Mutations in the genes for oocyte-derived growth factors GDF9 and BMP15 are associated with both increased ovulation rate and sterility in Cambridge and Belclare sheep (Ovis aries)

Belclare and Cambridge are prolific sheep breeds, the origins of which involved selecting ewes with exceptionally high litter size records from commercial flocks. The variation in ovulation rate in both breeds is consistent with segregation of a gene (or genes) with a large effect on this trait. Sterile ewes, due to a failure of normal ovarian follicle development, occur in both breeds. New naturally occurring mutations in genes for the oocyte-derived growth factors growth differentiation factor 9 (GDF9) and bone morphogenetic protein 15 (BMP15) are described. These mutations are associated with increased ovulation rate in heterozygous carriers and sterility in homozygous carriers in both breeds. This is the first time that a mutation in the gene for GDF9 has been found that causes increased ovulation rate and infertility in a manner similar to inactivating mutations in BMP15, and shows that GDF9 is essential for normal folliculogenesis in sheep. Furthermore, it is shown, for the first time in any species, that individuals with mutations in both GDF9 and BMP15 have a greater ovulation rate than sheep with either of the mutations separately.     

The proregion of mouse BMP15 regulates the cooperative interactions of BMP15 and GDF9

Bone morphogenetic protein 15 (BMP15) and growth and differentiation factor 9 (GDF9) are secreted by the mammalian oocyte and are essential for ovarian follicular development, ovulation, and fertility. However, the secreted forms of the BMP15 and GDF9 proteins and the nature of cooperative molecular interactions between BMP15 and GDF9 previously reported have not been fully characterized. In this study, we found that recombinant mouse BMP15 and GDF9 are secreted as cleaved mature and proregion proteins, with BMP15 also secreted as uncleaved promature protein. Noncovalent interactions were identified between the mature and proregion proteins of each growth factor. Moreover, GDF9 mature protein was found to coimmunoprecipitate with the BMP15 proregion, suggestive of a heteromeric association between BMP15 and GDF9. Mouse GDF9 was found to exist mostly as a dimer of mature protein, in both the presence and absence of BMP15. In contrast, BMP15 formed mostly multimers of proregion and mature protein when combined with GDF9, providing further evidence for heteromeric interaction. Mouse BMP15 was found to act cooperatively with GDF9 in a rat granulosa cell thymidine incorporation bioassay and to signal through the BMPR2 and ACVR1B/TGFBR1/ACVR1C receptor-mediated pathways. Immunoneutralization experiments using GDF9 mature protein antibody indicated that these cooperative interactions are species specific. Additionally, immunoneutralization with proregion antibodies highlighted the involvement of the BMP15 proregion in BMP15/GDF9 cooperative interactions. Taken together, these findings support a novel hypothesis where the extracellular cooperative interactions of recombinant mouse BMP15 and GDF9 are multimeric, involving the proregion of BMP15, and may well be species specific.     

A 17 bp deletion in the Bone Morphogenetic Protein 15 (BMP15) gene is associated to increased prolificacy in the Rasa Aragonesa sheep breed

Different mutations in the Bone Morphogenetic Protein 15 (BMP15) and the Growth Differentiation Factor 9 (GDF9) genes cause increased ovulation rate and infertility in a dosage-sensitive manner in sheep. They cause increased ovulation rate and twin and triplet births in heterozygotes, and complete primary ovarian failure in homozygotes resulting in total infertility. We are here presenting a novel mutation in the second exon of the ovine BMP15 gene, found in the Spanish breed Rasa Aragonesa. It consists of a 17 bp deletion resulting in displacement of the open reading frame and premature stop codons. As a consequence, nearly 85% of the sequence of the wild type aminoacidic chain in the second exon of the BMP15 pro-protein is modified or suppressed as only the first 45 amino acids are conserved of the 245 original. The mature peptide is lost. The ewes heterozygous for this deletion present very high prolificacy (2.66 lambs/birth) when compared to a mean flock prolificacy of 1.36 lambs. The deletion causes a complete lack of functionality of the second exon of BMP15, comparable to the effect of premature stop codons in other mutations. Therefore, homozygous females for the deletion are expected to present primary ovarian failure. DNA sequence analysis of the GDF9 coding regions detected only a synonymous Single Nucleotide Polymorphism (SNP), apparently not linked to changes in prolificacy.     

甘肃肉羊新品种选育群GDF9基因外显子2多态性及其与产羔性状关联分析

根据GenBank发布的绵羊GDF9基因外显子2的序列设计4对引物,采用PCR-SSCP技术分析GDF9基因外显子2在甘肃内羊新品种选育群羊中的单核苷酸多态性,并与产羔性状进行关联分析.结果表明,GDF9基因的扩增片段在所检测的新品种群羊中存在PCR-SSCP多态性,检测到3种基因型(AA、AB和BB),而在32只无角陶赛特母羊群中只检测到AA和AB基因型.测序结果显示,GDF9基因编码区第978位碱基发生A→G突变,但没有导致氨基酸的改变;第994位碱基发生G→A突变,导致Ⅴ变成Ⅰ(缬氨酸→异亮氨酸).新品种选育群羊产羔数的最小二乘均值关系为AB＞ AA＞ BB,统计分析结果初步表明3种基因型之间差异不显著(P＞0.05).故该区域可能不是影响新品种群羊繁殖力的功能结构区城.     

Cloning and expression levels of genes relating to the ovulation rate of the Yunling black goat

This study investigated the nucleotide sequences and tissue expression levels of genes relating to the ovulation rate in Yunling black goats, a famous Yunnan province, China, local breed with low fecundity. Five genes, FSHB, FSHR, BMP15, BMPR1B, and ESR2, were investigated; the complete cDNA sequences of these genes were 390, 2088, 1185, 1509, and 1585 bp, respectively, and compared with Boer goats (a more fecund breed), the sequence identities were 99%, 99%, 99%, 100%, and 99%, respectively. There were two base differences in FSHB and BMP15, four in FSHR, and three in ESR2. There were fewer follicles and oocytes in Yunling black goats than in Boer goats. Expression levels of FSHB, FSHR, and BMP15 genes in Yunling black goats were lower, and expression levels of BMPR1B and ESR2 genes were higher. Serum FSH content was lower in Yunling black goats, but serum estrogen content was higher. Protein expression levels of FSHR, BMP15, BMPR1B, and ESR2 in ovary and pituitary correlated positively with gene mRNA expression levels. In Yunling black goats, the mRNA expression levels of FSHB, FSHR, and BMP15 positively correlated with litter size, but those of BMPR1B and ESR2 correlated negatively. Together, base changes and variations in mRNA and protein expression levels of genes relating to the ovulation rate result in low fecundity in the Yunling black goat. Reduced BMP15 and FSHR levels may be related to the observed fewer oocytes and, consequently, fewer follicles.     

Polymorphism and nucleotide sequencing of BMPR1B gene in prolific Assam hill goat

Assam hill goat (Capra hircus) is a prolific local goat in India. bone morphogenetic protein receptor (BMPR1B) gene was studied as a candidate gene for the prolificacy of goats. The objective of the present study was to detect the incidence of mutation in the exonic region of BMPR1B gene of Assam hill goat. Total 90 blood samples were collected randomly from different parts of Assam and genomic DNA were extracted using phenol–chloroform method. The quantity and quality of extracted DNA was examined by spectrophotometry and gel electrophoresis, respectively. PCR amplicon showed a product of 140 bp fragment of BMPR1B gene. The purified product upon digestion with AvaII showed monomorphic banding pattern and revealed wild type alleles with AA genotype. Nucleotide sequencing showed one new mutation 773 (G→C) which is found to be unique in Assam hill goat. Construction of tree at nucleotide level generates from the present experiment lies in common cluster which differs from the other breeds of goat. The analysis of polymorphism for BMPR1B in Assam hill goat indicates that the genetic factor responsible for prolificacy or multiple kidding rates is not related to the reported mutated alleles of BMPR1B gene. Therefore, attempts to be made to detect other SNPs for BMPR1B gene or otherwise effort should be made towards other fecundity gene which might be responsible for the prolificacy of Assam hill goat.     

Association between PCR-SSCP of bone morphogenetic protein 15 gene and prolificacy in Jining Grey goats

The Jining Grey is a prolific local goat breed in P.R. China. Bone morphogenetic protein 15 (BMP15) gene that controls high fecundity of Inverdale, Hanna, Lacaune, Belclare, Cambridge, and Small Tailed Han ewes was studied as a candidate gene for the prolificacy of Jining Grey goats. According to the sequence of ovine BMP15 gene, six pairs of primers were designed to detect single nucleotide polymorphisms in exon 1 and exon 2 of the BMP15 gene in both high fecundity breed (Jining Grey goats) and low fecundity breeds (Boer, Liaoning Cashmere, and Inner Mongolia Cashmere goats) by single strand conformation polymorphism (SSCP). Two pairs of primers (F1/R1 and F2/R2) were used to amplify the exon 1. Four pairs of primers (F3/R3, F4/R4, F5/R5, and F6/R6) were used to amplify the exon 2. Only the products amplified by primer F5/R5 displayed polymorphism. Results indicated that two genotypes (AA and AB) were detected in prolific Jining Grey goats and only one genotype (AA) was detected in low fecundity goat breeds. In Jining Grey goats frequencies of genotypes AA and AB were 0.10 and 0.90, respectively. Sequencing revealed two point mutations (G963A and G1050C) of BMP15 gene in the AB genotype in comparison to the AA genotype. In Jining Grey goats the heterozygous AB does had 1.13 (p < 0.01) kids more than the homozygous AA does. These results preliminarily showed that the BMP15 gene is either a major gene that influences the prolificacy of Jining Grey goats or a molecular genetic marker in close linkage with such a gene.     

Polymorphism identification in goat GNRH1 and GDF9 genes and their association analysis with litter size

This study investigated the polymorphisms of GNRH1 and GDF9 genes in 641 goats of three breeds: Xinong Saanen, Guanzhong and Boer. Two allelic variants were identified in the GNRH1 gene (JN645280:g.3548A>G and JN645281:g.3699G>A) and one allelic variant was found in the GDF9 gene (JN655693:g.4093G>A). Furthermore, g.4093G>A was a missense mutation (p.Val397Ile of GDF9). Results of an association analysis indicated that SNPs g.3548A>G and g.4093G>A had significant effects on litter size (P < 0.05). The combination genotypes of SNPs g.3548A>G, g.3699G>A and g.4093G>A also affected litter size with the C₅ genotype having the highest litter size in the first, third, fourth and average parity. Hence, the biochemical and physiological functions, together with the results obtained in our investigation, suggest that the GNRH1 and GDF9 genes could serve as genetic markers for litter size in goat breeding.     

Association between PCR-SSCP of growth differentiation factor 9 gene and high prolificacy in Small Tail Han sheep

Small Tail Han sheep that has significant characteristics of high prolificacy and nonseasonal ovulatory activity is an excellent local sheep breed in P.R. China. The lambing percentage averaged 260% in Small Tail Han sheep. Growth differentiation factor 9 (GDF9) gene, which was essential for growth and differentiation of early ovarian follicles, was considered as a possible candidate gene for litter size in Small Tail Han sheep. The genetic polymorphism of a part of the GDF9 gene was detected in 130 ewes of Small Tail Han sheep by PCR-SSCP. The results indicated that there were two genotypes (AA and AB) detected by two primer pairs. In both exon 1 and exon 2 of the GDF9 gene in Small Tail Han sheep, frequencies of AA genotype were 0.846 and 0.908, frequencies of AB genotype were 0.154 and 0.092, frequencies of A allele were 0.923 and 0.954, and frequencies of B allele were 0.077 and 0.046, respectively. The results of chi2 fitness test indicated that both exon 1 and exon 2 of the GDF9 gene were in Hardy-Weinberg equilibrium (p > 0.05) in Small Tail Han sheep. Least squares means of litter size in the first and the second parity for genotype AA were 0.30 (p <0.05) and 0.77 (p <0.0001) more than those for genotype AB detected in exon 1 of the GDF9 gene in Small Tail Han sheep, respectively. Fragments detected in exon 2 of the GDF9 gene had no significant effect (p > 0.05) on litter size in both the first and the second parity in Small Tail Han sheep. Litter size in sheep is lowly heritable, expressed only in females, and manifested relatively late in life. Access to genetic markers would thus be advantageous in selection programs.     

Maximum‐likelihood approaches reveal signatures of positive selection in BMP15 and GDF9 genes modulating ovarian function in mammalian female fertility

Bone morphogenetic proteins (BMPs) and the growth factors (GDFs) play an important role in ovarian folliculogenesis and essential regulator of processes of numerous granulosa cells. BMP15 gene variations linked to various ovarian phenotypic consequences subject to the species, from infertility to improved prolificacy in sheep, primary ovarian insufficiency in women or associated with minor subfertility in mouse. To study the evolving role of BMP15 and GDF9, a phylogenetic analysis was performed. To find out the candidate gene associated with prolificacy in mammals, the nucleotide sequence of BMP15 and GDF9 genes was recognized under positive selection in various mammalian species. Maximum‐likelihood approaches used on BMP15 and GDF9 genes exhibited a robust divergence and a prompted evolution as compared to other TGFβ family members. Furthermore, among 32 mammalian species, we identified positive selection signals in the hominidae clade resulting to 132D, 147E, 163Y, 191W, and 236P codon sites of BMP15 and 162F, 188K, 206R, 240A, 244L, 246H, 248S, 251D, 253L, 254F and other codon sites of GDF9. The positively selected amino acid sites such as Alanine, Lucien, Arginine, and lysine are important for signaling. In conclusion, this study evidences that GDF9 and BMP15 genes have rapid evolution than other TGFß family members and was subjected to positive selection in the mammalian clade. Selected sites under the positive selection are of remarkable significance for the particular functioning of the protein and consequently for female fertility.     

Polymorphisms of caprine <Emphasis Type="Italic">GDF</Emphasis>9 gene and their association with litter size in Jining Grey goats

The exons 1, 2 and flanking region of growth differentiation factor 9 (GDF9) gene in five randomly selected does of Jining Grey, Boer and Liaoning Cashmere goats were amplified and analyzed. Thirteen nucleotide differences were identified in GDF9 gene between sheep (AF078545) and goats. Four SNPs (G3288A in intron 1, G423A, A959C [Gln320Pro] and G1189A [Val397Ile] in exon 2) were detected in four goat breeds with different prolificacy, in which G3288A was a new SNP in goats. The results showed that loci 3288, 423 and 1189 in Boer goats, loci 3288 and 423 in Guizhou White goats, loci 423 and 1189 in Liaoning Cashmere goats were all in complete linkage disequilibrium (D′ = 1, r ² = 1), respectively. In moderate (Boer goat) and low prolificacy (Liaoning Cashmere goat) breeds, linkage analysis indicated that there were more fervent linkage disequilibrium among loci 3288, 423 and 1189 than high prolificacy (Jining Grey and Guizhou White goats) breeds. For the 959 locus, the genotype distribution showed obvious difference between high prolificacy breeds and moderate or low prolificacy breeds (P < 0.05 or P < 0.01). The Jining Grey goat does with genotype CC or AC had 0.81 (P < 0.01) or 0.63 (P < 0.01) kids more than those with genotype AA, respectively. The present study preliminarily showed an association between allele C at 959 locus of GDF9 gene and high litter size in Jining Grey goats. These results provide further evidence that the GDF9 gene may be significantly correlated with high prolificacy in goats.     

Whole-genome sequencing identifies potential candidate genes for reproductive traits in pigs

Reproductive performance is a complex quantitative trait, that is determined by multiple genes, regulatory pathways and environmental factors. A list of major genes with large effect have been detected, although multiple QTLs are identified. To identify candidate genes for pig prolificacy, whole genome variants from five high- and five low-prolificacy Yorkshire sows were collected using whole-genome resequencing. A total of 13,955,609 SNPs and 2,666,366 indels were detected across the genome. Common differential SNPs and indels were identified between the two groups of sows. Genes encoding components of the TGF-beta signaling pathway were enriched with the variations, including BMP5, BMP6, BMP7, ACVR1, INHBA, ZFYVE9, TGFBR2, DCN, ID4, BAMBI, and ACVR2A. Several differential variants within these genes related to reproductive traits were identified to be associated with litter size. A comparison of selective regions and published QTL data suggests that NEDD9, SLC39A11, SNCA, and UNC5D are candidate genes for reproduction traits.