Fine mapping of theFecLlocus influencing prolificacy in Lacaune sheep

Detection of Booroola Polymorphism of Bone Morphogenetic Protein Receptor 1b and Embrapa Polymorphism of Growth Differentiation Factor 9 in Sheep in Thailand

This study aimed to investigate the appearance and frequencies of the Booroola polymorphism of the bone morphogenetic protein receptor 1b (BMPR1B) gene (FecB) and the Embrapa polymorphism of the growth differentiation factor 9 (GDF9) gene (FecGE) in sheep in Thailand. A total of 454 crossbred sheep blood samples were collected from four provinces in Thailand during August 2022 to July 2023. The polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method was used to identify the FecB and FecGE genotypes. The history of ewe birth types was collected from the owners to analyze the association between fecundity (Fec) genotypes and the history of birth types. The genotypic frequencies of FecB for homozygous genotype (B/B), heterozygous genotype (+/B), and wildtype (+/+) were 0.22%, 1.54%, and 98.24%, respectively. Meanwhile, the genotypic frequencies of FecGE for homozygous genotype (E/E), heterozygous genotype (+/E), and wildtype (+/+) were 0.00%, 2.42%, and 97.58%, respectively. Furthermore, three ewes exhibited both FecB and FecGE genotypes. Fisher's exact test revealed that possession of the FecB genotype was associated with multiple births (p < 0.01). Both FecB and FecGE mutations were identified in crossbred sheep in Thailand. Sheep containing FecB allele could be alternative candidates to be selected to improve the prolificacy of crossbred sheep in Thailand.

An Identification of Functional Genetic Variants in B4GALNT2 Gene and Their Association with Growth Traits in Goats

β-1,4-N-acetylgalactosamine transferase 2 (B4GALNT2) is a vital candidate gene that affects the growth traits in sheep. However, whether it has the same function in goats remains to be investigated further. This study selected 348 Nanjiang Yellow goats, screened all exons, and conserved non-coding regions of the B4GALNT2 gene for single-nucleotide polymorphisms (SNPs). Our results revealed the presence of a synonymous mutation, rs672215506, within the exon of the B4GALNT2 gene in the Nanjiang Yellow goat population. The mutation resulted in a decrease in the mRNA stability of the B4GALNT2 gene. The results of SNP detection of the conserved non-coding region of the B4GALNT2 gene showed five potential regulatory SNPs in the Nanjiang Yellow goat population. Except for rs66095343, the ~500 bp fragments of the other four SNPs (rs649127714, rs649573228, rs652899012, and rs639183528) significantly increased the luciferase activity both in goat skeletal muscle satellite cells (MuSCs) and 293T cells. The genetic diversity indexes indicated low or intermediate levels for all six SNPs analyzed, and the genotype frequencies were in Hardy-Weinberg equilibrium. Association analysis showed that rs660965343, rs649127714, and rs649573228 significantly correlate with growth traits in the later stage of growth and development of Nanjiang Yellow goats. The haplotype combinations of H2H3 and H2H2 had higher body weight and greater body size. Moreover, H2H2 haplotype combinations significantly correlated with the litter size of the Nanjiang Yellow goats. The results of our study demonstrate the potential role of the B4GALNT2 gene as a functional genetic marker in the breeding programs of Nanjiang Yellow goats.

Bioinformatics analysis and the association of polymorphisms within the caprine GDF9 gene promoter with economically useful traits in Damani goats

The blood sample from 60 Damani does were collected and genomic DNA was extracted, and DNA integrity were investigated. A 447 bp promoter fragment of the GDF9 gene was amplified and Sanger sequenced for the identification of GDF9 gene polymorphism. Three novel SNPs were identified at positions g. 97(T > A), g. 142 (G > G) and g. 313(C > T) in the promoter region of the caprine GDF9 gene which significantly (P < 0.05) influenced litter size, body measurement, and milk production traits in Damani goats. The genotype CT of SNP1 significantly (P < 0.05) improved litter size, genotype GG of SNP2 significantly (P < 0.05) enhanced milk production, while the genotypes CC of SNP3 significant (P < 0.05) increased body measurement traits in Damani goats. Moreover, in SNP1 loss of 3 transcription factors (TF) binding sites occurred, SNP2 caused loss of two TFs binding sites, and SNP3 caused loss of a single TF binding site. Similarly, SNP1 and SNP2 caused the gain of three new potential TF binding sites, and SNP3 caused gain of two new TF binding sites. It is concluded that caprine GDF9 gene could be used as a candidate gene for litter size, milk production and body measurement traits in Damani goats through marker-assisted selection for future breeding program.

FecB mutation and litter size are associated with a 90-base pair deletion in BMPR1B in East Friesian and Hu crossbred sheep

Litter size is a critical economic trait in livestock, but only a few studies have focused on associated indel mutations in BMPR1B, a key regulator of ovulation and litter size in sheep. We evaluated the effects of BMPR1B mutations on the reproductive performance of sheep. We used Hu, East Friesian, and East Friesian/Hu crossbred sheep as experimental subjects and identified a novel 90 bp deletion in BMPR1B, which coincides with the c.746A > G (FecB mutation) genotype. The correlation between the two loci and litter size was then evaluated. We identified three genotypes for the Del-90bp locus, namely, II, ID, and DD, and three genotypes for the c.746A > G locus, namely ++, B+, and BB. Both Del-90bp and c.746A > G significantly affected the litter size of Hu and East Friesian/Hu crossbred sheep. Linkage disequilibrium analysis revealed a strong linkage disequilibrium between these loci in Hu sheep and the F1 population (r2 > 0.33), which suggests that detecting this 90 bp deletion might be a simple method to identify the likely carriers of c.746A > G. However, the function of this 90-bp deletion still needs further exploration. We provide genetic data that can be used as a reference for the breeding of improved prolific traits in sheep.

Effects of New Mutations in BMPRIB, GDF9, BMP15, LEPR, and B4GALNT2 Genes on Litter Size in Sheep

Prolificacy is a crucial characteristic of livestock, particularly for species such as sheep that have many births. The objectives of this study were as follows: (1) to investigate the genetic diversity of the 13 new and 7 known variants in the BMPRIB, GDF9, BMP15, LEPR, and B4GALNT2 genes in Ujimqin (UM), the F1 population of Dorper × Ujimqin crossbred (DPU), the F1 population of Suffolk × Ujimqin crossbred (SFKU), Sonid sheep (SN), Tan sheep (Tan), Hu sheep (Hu), and Small-tailed Han sheep (STH) sheep breeds/populations; (2) to perform an association analysis of the above 20 variants with litter size in 325 UM, 304 DPU, and 66 SFKU sheep populations; (3) to compare the frequencies of the litter-size-related alleles of these 20 variants among 8 sheep breeds/populations (the above seven sheep breeds + Mongolia sheep breed). With the use of the Sequenom MassARRAY®SNP assay technology, these 20 mutations were genotyped. The association analysis results showed that the c.746A>G (FecB) mutation in BMPR1B was significantly associated with the litter size of UM and DPU, the c.994A>G (FecGA) in GDF9 was significantly associated with the litter size of SFKU, and the c.31_33CTTinsdel (B1) in BMP15 was significantly associated with the litter size of UM. Our findings might provide valuable genetic markers for expanding sheep litter sizes.

Detection of polymorphism in growth differentiation factor 9 gene (GDF9) Exon1 and its association with litter size in local Iraqi goats

Litter size is one of the most important economic traits in goats. Growth Differential Factor 9 (GDF9) Gene is suggested as a functional candidate for fertility, fecundity and twining rate. This study was performed to identify SNPs in exon1 of the GDF9 gene and their association with fertility traits in local Iraqi goats. The exon 1 of GDF9 was sequenced, and the SNP (Single nucleotide polymorphism ) was determined in 36 local Iraqi goats. Then the association analyses between polymorphic locals of GDF9 and litter size were performed using the chi-square analysis procedure. The result shows that mutation (C 1902 G) was significantly associated with litter size in local Iraqi goats and the location of C1902 G mutation with three genotypes CC, CG and GG. A significant effect was found for birth weight, weaning weight and twin ratio, where the homozygous dominant genotype outperformed CC. A significant impact was seen for fertility and mortality, where CC topped CG and GG in fertility percentage and superiority of CG over CC and GG in mortality rate. The results preliminarily demonstrated that GDF9 was a critical gene affecting the fecundity of local Iraqi goats and that (C 1902 G) could be a potential genetic marker. Keywords: GDF9 gene, litter size, polymorphism, Iraqi goats

Identification of differentially expressed long noncoding RNAs in the ovarian tissue of ewes Shal and Sangsari using RNA-seq

Background

The ovary has an important role in reproductive function. Animal reproduction is dominated by numerous coding genes and noncoding elements. Although long noncoding RNAs (LncRNAs) are important in biological activity, little is known about their role in the ovary and fertility.

Methods

Three adult Shal ewes and three adult Sangsari ewes were used in this investigation. LncRNAs in ovarian tissue from two breeds were identified using bioinformatics analyses, and then target genes of LncRNAs were discovered. Target genes were annotated using the DAVID database, and their interactions were examined using the STRING database and Cytoscape software. The expression levels of seven LncRNAs with their target genes were assessed by real‐time PCR to confirm the RNA‐seq.

Results

Among all the identified LncRNAs, 124 LncRNAs were detected with different expression levels between the two breeds (FDR < 0.05). According to the DAVID database, target genes were discovered to be engaged in one biological process, one cellular component, and 21 KEGG pathways (FDR < 0.05). The PES1, RPS9, EF‐1, Plectin, SURF6, CYC1, PRKACA MAPK1, ITGB2 and BRD2 genes were some of the most crucial target genes (hub genes) in the ovary.

Conclusion

These results could pave the way for future efforts to address sheep prolificacy barriers.

Genome-wide detection of selective signals for fecundity traits in goats (Capra hircus)

Fecundity in livestock is an economically important complex quantitative trait that is influenced by both genetics and the environment. However, the underlying genetic mechanism of reproductive performance in goats has not been well investigated. To investigate the genomic basis of fecundity in goats, genomic sequencing data of the Jining grey goat (a high prolificacy breed in China) were collected, as well as data for other commonly available goat breeds, and a mass of genomic variants were generated after variation calling. We screened the Jining grey goat (20 individuals) using a selective sweep with the Asian wild goat population (5 individuals), and potential candidate genes were proposed, such as STIM1, ESR1, LRRC14B and SLC9A3. Among, STIM1 is a most promising one associated with high reproductive capacity. When compared to Chinese domestic goats with low fecundity (17 individuals), the genes including MLLT10, SPIRE2, TCF25, ZNF276 and FANCA were screened, and the SPIRE2 gene was thought to be associated with fecundity traits. Meanwhile, the functional enrichment of these candidate genes revealed that they were involved in biological processes of mammary gland morphogenesis, uterus development, gastrulation, mesoderm morphogenesis and formation, and blood vessel development, which might undergo natural or artificial selection during reproductive trait formation in goats. Thus, our findings could enrich the genetic basis of reproductive trait selection during goat domestication, which may serve to improve goat breeding practices.

Polymorphisms in BMPRIB gene affect litter size in Chinese indigenous sheep breed

The BMPRIB gene belongs to the TGF-β superfamily and is considered to be a regulator of sheep reproductive performance. Single nucleotide polymorphisms (SNPs) of BMPRIB gene in the Small Tail Han, Hu, Mongolian, Oula, Gansu Alpine Fine-wool, Dorper and Australian White sheep were detected by Sanger sequencing. Five SNPs (rs427897187 G > A, rs418841713 A > G, rs159952533 T > C, rs429416173 C > A and rs403555643 A > G) of BMPRIB gene were identified. For rs427897187 G > A, further analysis revealed that genotype GG and GA had 0.26 (p < 0.05) and 0.33 (p < 0.05) litter size less than those with genotype AA in Oula sheep. For rs403555643 A > G, further analysis revealed that genotype GG and AG had 0.65 (p < 0.05) and 0.38 (p < 0.05) litter size more than those with genotype AA in Oula sheep, and genotype GG had 0.56 (p < 0.05) litter size more than those with genotype AA in Mongolian sheep. The results showed that rs427897187 G > A and rs403555643 A > G are potential molecular markers wich could improve litter size of Chinese indigenous sheep and be used in Chinese indigenous sheep breeding.

Association between novel variants in BMPR1B gene and litter size in Mongolia and Ujimqin sheep breeds

Prolificacy is an important trait of animals, specifically for sheep. The Bone morphogenetic protein receptor 1B (BMPR1B) is a major gene affecting the litter size of many sheep breeds. The well-known FecB mutation (Q249R) was associated fully with the hyper prolific phenotype of Booroola Merino. However, the identification of variation in all exonic regions of BMPR1B was rare. In this study, we sequenced all exonic regions of BMPR1B gene of Mongolia sheep breed, and ten novel variants were detected by direct sequencing. Among them, the litter size of the Mongolia ewes with the CC genotype was significantly higher (0.34 additional lambs, p < .05) than those with the TT genotype of the g.29346567C>T single nucleotide polymorphism (SNP). The litter size of the Mongolia ewes with the TT genotype was significantly higher (0.19 additional lambs, p < .05 and .31 additional lambs, p < .01, respectively) than those with the GT and GG genotypes of the c.1470G>T SNP. The silent c.1470G>T mutation is predicted to increase the stability of the mRNA secondary structure through reducing minimum free energy and is predicted to change the mRNA secondary structure of BMPR1B. Our findings may give potentially useful genetic markers for increasing litter size in sheep.

Identification of sequence variation in the oocyte-derived bone morphogenetic protein 15 (BMP15) gene (BMP15) associated with litter size in New Zealand sheep (Ovis aries) breeds

Advances in the study of reproductive traits indicate that functional variation in fertility genes may be useful for improving sheep fertility. The aim of this study was to search for variation in the bone morphogenetic protein 15 gene (BMP15) and ascertain any association with litter size in purebred Finnish Landrace sheep (n = 148), Finnish Landrace × Texel-cross sheep (n = 45), and composite sheep (of varying breed background; n = 58) from New Zealand (NZ). A 482 bp and 312 bp fragment of exon 1 and 2, respectively, of BMP15 were analysed using polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP). The additive and dominance effect of BMP15 variation on litter size were estimated using animal and sire models. Two variants (A and B) were detected in exon 1; no sequence variation was detected in exon 2. Variant A had the nucleotide sequence CTT between positions c.31 and c.33, while variant B had a deletion (c.31_33del). The observed frequency for variant A in the Finnish Landrace sheep, Finnish Landrace × Texel-cross sheep and the composite sheep, was 0.77, 0.92, and 0.68, respectively while the frequency of variant B (c.31_33del) was 0.23, 0.08, and 0.32, respectively. An association between litter size and c.31_33del (P < 0.001) was observed in composite sheep. Analysis of more sheep will be required to confirm these results. Litter size did not differ significantly between sheep breeds regardless of the presence/absence of c.31_33del. Results suggested that c.31_33del might be a genetic marker for improving fecundity in some NZ sheep.

Genome-wide association study and pathway analysis identify NTRK2 as a novel candidate gene for litter size in sheep

Litter size is one of the most important economic traits in sheep. Identification of gene variants that are associated with the prolificacy rate is an important step in breeding program success and profitability of the farm. So, to identify genetic mechanisms underlying the variation in litter size in Iranian Baluchi sheep, a two-step genome-wide association study (GWAS) was performed. GWAS was conducted using genotype data from 91 Baluchi sheep. Estimated breeding values (EBVs) for litter size calculated for 3848 ewes and then used as the response variable. Besides, a pathway analysis using GO and KEGG databases were applied as a complementary approach. A total of three single nucleotide polymorphisms (SNPs) associated with litter size were identified, one each on OAR2, OAR10, and OAR25. The SNP on OAR2 is located within a novel putative candidate gene, Neurotrophic receptor tyrosine kinase 2. This gene product works as a receptor which is essential for follicular assembly, early follicular growth, and oocyte survival. The SNP on OAR25 is located within RAB4A which is involved in blood vessel formation and proliferation through angiogenesis. The SNP on OAR10 was not associated with any gene in the 1Mb span. Moreover, gene-set analysis using the KEGG database identified several pathways, such as Ovarian steroidogenesis, Steroid hormone biosynthesis, Calcium signaling pathway, and Chemokine signaling. Also, pathway analysis using the GO database revealed several functional terms, such as cellular carbohydrate metabolic, biological adhesion, cell adhesion, cell junction, and cell-cell adherens junction, among others. This is the first study that reports the NTRK2 gene affecting litter size in sheep and our study of this gene functions showed that this gene could be a good candidate for further analysis.

G1 point mutation in growth differentiation factor 9 gene affects litter size in Sudanese desert sheep

BACKGROUND AND AIM

Sudanese desert sheep encompass different sheep breeds named according to the different Sudanese tribes that rear them such as the Dubasi, Shugor, and Watish sheep. The objectives of this study were to screen for G1 point mutation in the polymorphic growth differentiation factor 9 (GDF9) gene, investigate its association with litter size, and construct the phylogeny of the different tribal breeds that belong to the Sudanese Desert sheep tribal types.

MATERIALS AND METHODS

Genomic DNA was extracted from whole blood of three tribal Desert sheep breeds (Dubasi, Watish, and Shugor) using the guanidine chloride method. Polymerase chain reaction-restriction fragment length polymorphism with HhaI restriction enzyme and sequencing techniques was used for genotyping the GDF9 locus for possible mutations associated with litter size in the three desert sheep tribal types.

RESULTS

G1 mutation in GDF9 caused the replacement of Arginine by Histidine at residue 87. The wild type allele (A) had the highest frequency, whereas the mutant type allele (a) had the lowest in all the sequenced subtypes. The genotype frequencies of the wild type ewes (AA) were higher than the heterozygous (Aa) and the mutant type (aa) frequencies in the three studied desert sheep types. No significant differences were found in the allele frequency between the three tribal types. Litter size was significantly influenced by the genotypes of GDF9 gene, parities, and subtypes (p≤0.01, 0.01, and 0.05, respectively). In the Watish sheep type, heterozygous sheep in their second parity recorded the highest litter size. Sequence alignment of GDF9 gene samples with the database entry indicated that all three tribal types were similar and identical to the reference sequence. The phylogenetic tree revealed that Shugor is the common ancestor of the studied types and Watish is more closely related to Shugor than Dubasi. This result mi ght partly explain the lower reproductive performance of Dubasi compared to Watish and Shugor.

CONCLUSION

The presence of one copy of GDF9 gene increased litter size in the studied Sudanese Desert sheep. This locus may be used as a biomarker for litter size improvement through genotypic selection and allele or gene introgression.

Novel Variants in GDF9 Gene Affect Promoter Activity and Litter Size in Mongolia Sheep

Litter size is an economically important trait in sheep breeding. The objectives of this study were as follows: (1) to ascertain if any of the 19 known variants in the BMPRIB, BMP15, and GDF9 genes are present and associated with the litter size of Mongolia sheep; (2) to identify novel variants in GDF9 and perform association analysis; and (3) to validate the effects of these GDF9 promoter variants on the activity of the gene. The results of the 19 known variants showed that the FecB^B affected the litter size of Mongolia sheep (p < 0.001). The association analysis results of novel variants showed that the g.46544883A>G (GenBank accession: NC_040256, the same below) in the 3' untranslated region (3' UTR), the c.1040T>C (Phe347Ser) in the exon 2, and the g.46547859C>T SNP in the promotor of GDF9 were significantly associated with litter size of Mongolia ewes (p < 0.01, p < 0.05, and p < 0.001, respectively). In addition, the GDF9 promoter activity analysis showed that the C allele at the -332 position (g.46547859C>T) could decrease luciferase activity compared with the T allele (p < 0.01). Our findings may facilitate effective marker-assisted selection to increase litter size in Mongolia sheep populations, as well as bring new insights into GDF9 expression.

Review on Genomic Regions and Candidate Genes Associated with Economically Important Production and Reproduction Traits in Sheep (Ovies aries)

Sheep (Ovis aries) is one of the most economically, culturally, and socially important domestic animals. They are reared primarily for meat, milk, wool, and fur production. Sheep were reared using natural selection for a long period of time to offer these traits. In fact, this production system has been slowing the productivity and production potential of the sheep. To improve production efficiency and productivity of this animal through genetic improvement technologies, understanding the genetic background of traits such as body growth, weight, carcass quality, fat percent, fertility, milk yield, wool quality, horn type, and coat color is essential. With the development and utilization of animal genotyping technologies and gene identification methods, many functional genes and genetic variants associated with economically important phenotypic traits have been identified and annotated. This is useful and presented an opportunity to increase the pace of animal genetic gain. Quantitative trait loci and genome wide association study have been playing an important role in identifying candidate genes and animal characterization. This review provides comprehensive information on the identified genomic regions and candidate genes associated with production and reproduction traits, and gene function in sheep.

Integrated Analysis of miRNA and mRNA Expression Profiles Reveals Functional miRNA-Targets in Development Testes of Small Tail Han Sheep

Small Tail Han Sheep is a highly valued local breed in China because of their precocity, perennial estrus, and high fecundity. The average annual lambing rate of ewes is as high as 180–270%, the semen of ram has characteristics of high yield, high density, and good motility. To reveal the key miRNAs and miRNA-targets underlying testis development and spermatogenesis in male Small Tail Han Sheep, integrated analysis of miRNA and mRNA expression profiles in 2-, 6-, and 12-month-old testes was performed by RNA-seq technology and bioinformatics methods. The results showed that total of 153 known sheep miRNAs and 2712 novel miRNAs were obtained in 2-,6 - and 12-month-old Small Tail Han Sheep testes; 5, 1, and 4 differentially expressed (DE) known sheep miRNAs, and 132, 105, and 24 DE novel miRNAs were identified in 2- vs. 6-, 6- vs. 12-, and 2- vs. 12-month-old testes, respectively. We combined miRNA results of this study and the mRNA results obtained in our previous study to predict the target mRNAs of DE known sheep miRNAs; 131, 10, and 15 target mRNAs of DE known sheep miRNAs and 76, 1, and 11 DE miRNA–targets were identified in the three groups, respectively. GO and KEGG analyses showed that: in 2- vs. 6-month-olds, the target genes of DE known sheep miRNAs were involved in 100 biological processes and 11 signaling pathways; in 6- vs. 12-month-olds, the target genes of DE known sheep miRNAs were involved in 4 biological processes; and in 2- vs. 12-month-olds, the target genes of DE known sheep miRNAs were involved in 17 biological processes and 4 signaling pathways. Three miR–target regulatory networks were constructed based on these DE miRNA–targets. The key miRNA-Targets involved in testis development and spermatogenesis were screened. 6 known sheep miRNAs and 6 novel miRNAs were selected to validate the accuracy of miRNA sequencing data by qRT-PCR. The binding sites of oar-miR-379-5p with WNT8A was validated by a dual luciferase reporter gene detection system.

The high prolificacy of D'man sheep is associated with the segregation of the FecLL mutation in the B4GALNT2 gene

Mutations in the FecL locus are associated with large variation in ovulation rate and litter size in the French Lacaune sheep breed. It has been shown that the B4GALNT2 gene within the FecL locus is most likely responsible for the high fecundity in the French breed. In this study, we have highlighted the segregation of the FecL^L mutation within the B4GALNT2 gene in North African sheep breeds and notably in the highly prolific D'man breed. Genotyping of a sample of 183 Tunisian D'man individuals revealed a high frequency (0.65) of the prolific allele FecL^L which was attributed to the adoption of a decades-old breeding strategy based on the selection of ewe lambs born from large litter size. Homozygous LL ewes showed a significantly increased litter size compared to heterozygous and non-carrier ewes (FecL^L /FecL^L  = 2.47 ± 0.09 vs. FecL^L /FecL⁺  = 2.23 ± 0.09, p < 0.05 and FecL⁺ /FecL⁺  = 1.93 ± 0.18, p < 0.01). The presence of the FecL^L polymorphism in both D'man and Lacaune breeds argues for an ancestral origin of this mutation and brings an answer to the old question of the genetic determinism of the extreme prolificacy of the D'man ewes. The results of this study can help to establish planned genotype-based mating allowing both higher profit for the breeders and an optimal management of the FecL^L mutation in D'man sheep populations.

Molecular Cloning of the B4GALNT2 Gene and Its Single Nucleotide Polymorphisms Association with Litter Size in Small Tail Han Sheep

Simple Summary

In French Lacaune sheep, the B4GALNT2 (beta-1, 4-N-acetyl-galactosaminyl transferase 2) gene was considered as the potential gene for a FecL (mutation), which regulates the ovine ovulation rate. Three specific mutation sites linked with the FecL mutation have not been previously found in 11 sheep breeds. However, two mutations of g.36946470C > T and g.36933082C > T in the exon of B4GALNT2 were found to have had a significant effect on the litter size in the first parity for Small Tail Han (STH) Sheep (p < 0.05). B4GALNT2, which is mainly expressed in ovine ovary, also plays an important role in sheep reproduction. Furthermore, we discovered two transcription start sites (TSS) of B4GALNT2 in its 5′-flanking region in ovine granule cells in vitro.

Abstract

A new fecundity gene named the FecL (mutation), which regulates the ovulation rate, was discovered in French Lacaune sheep. The B4GALNT2 (beta-1, 4-N-acetyl-galactosaminyl transferase 2) gene was considered as the potential FecL mutation gene. This study explores whether the effect of the FecL mutation exists in other sheep breeds, and the features of the B4GALNT2 gene in terms of the molecular structure and its expression profile. Using Sanger sequencing, we found that high and low fecundity breeds from among 11 measured sheep breeds all had no variation in the three specific mutation sites, which were linked with the FecL mutation. However, two mutations of g.36946470C > T and g.36933082C > T in the exon of B4GALNT2 had a significant effect on litter size in the first parity for Small Tail Han (STH) Sheep (p < 0.05). Two transcription start sites (TSS) of B4GALNT2 in its 5′-flanking region were discovered in ovine granule cells in vitro, through the RACE (Rapid amplification of cDNA ends) method. Except for in the kidney and oviduct, no significant difference in expression levels had been found between STH sheep and Tan sheep breeds. The B4GALNT2 gene, as a candidate for FecL, may have a relationship with the differences in litter size in STH sheep. B4GALNT2 is mainly expressed in the ovine ovary, which also suggests that B4GALNT2 plays an important role in sheep reproduction.

Detection of single nucleotide polymorphisms at major prolificacy genes in the Mehraban sheep and association with litter size

The present study aimed to investigate the presence of polymorphisms at four known genes controlling ovine prolificacy i.e. BMP15, GDF9, BMPR1B and B4GALNT2 in a sample of 115 Iranian Mehraban ewes and their association with litter size (LS) and lambs’ birth weight (BW) traits. Using Sanger sequencing of exons and polymorphism specific genotyping, ten SNPs (Single Nucleotide Polymorphisms) were observed in only two genes, GDF9 and BMPR1B. Seven SNPs were found in the GDF9 gene on the chromosome 5. Among them, six were already described in the coding sequence, and a new one (g.41840985C>T) was found in the 3’UTR. In the BMPR1B gene on the chromosome 6, three novel SNPs were detected in the exon 7 (g.29382184G>A; g.29382337G>A and g.29382340G>A). Allelic frequencies were established for six SNPs among the ten identified and they were in Hardy-Weinberg equilibrium. A significant association was found between the novel SNPs found in the exon 7 of BMPR1B and LS. Present results indicate the potential role of the BMPR1B locus in controlling prolificacy of Mehraban sheep and provide genetic markers for further exploitation in selection to improve reproductive efficiency.

Optimal mating strategies to manage a heterozygous advantage major gene in sheep

Some mutations (or 'major genes') have a desirable effect in heterozygous carriers but an undesirable effect in homozygous carriers. When these mutations affect a trait of significant economic importance, their eradication, depending on their effect and frequency, may be counterproductive. This is especially the case of major genes affecting the ovulation rate and thus the prolificacy in meat sheep populations. To manage such situations, a mating design based on the major genotypes of reproducers has to be optimized. Both the effect of the major gene and the cost of genotyping candidates at this locus influence the expected genetic progress and profitability of the breeding plan. The aim of this study was to determine the optimal combination of matings that maximizes profitability at the level of the whole population (nucleus + commercial flocks). A deterministic model was developed and, using sequential quadratic programming methodology, the optimal strategy (optimal combination of matings) that maximized the economic gain achieved by the population across a range of genotype effects and genotyping costs was determined. The optimal strategy was compared with simpler and more practical strategies based on a limited number of parental genotype mating types. Depending on the genotype effect and genotyping costs, the optimal strategy varied, such that either the heterozygous frequency and/or polygenic gain was maximized with a large number of animals genotyped, or when genotyping costs were higher, the optimization led to lower heterozygous frequency and/or polygenic gain with fewer animals genotyped. Comparisons showed that some simpler strategies were close to the optimal strategy. An overlapping model was then derived as an application of the real case of the French Lacaune meat sheep OVI-TEST breeding program. Results showed that a practical strategy based on mating non-carriers to heterozygous carriers was only slightly less effective than the optimal strategy, with a reduction in efficiency from 3% to 8%, depending on the genotyping costs. Based on only two different parental genotype mating types, this strategy would be easy to implement.

Genetic polymorphism of growth differentiation factor 9 (GDF9) gene related to fecundity in two Egyptian sheep breeds

PURPOSE

This study explores polymorphisms in the growth differentiation factor 9 (GDF9) gene (exon 1) with respect to fertility in Egyptian sheep.

METHODS

Blood samples were collected, and genomic DNA was extracted from 24 Saidi and 13 Ossimi ewes. A 710 bp portion of the GDF9 gene, was amplified using specific primers, and the sequence was analyzed to clarify the phylogenetic relationship of Egyptian breed sheep. In addition, the PCR-RFLP method using Pst1 or Msp1 restriction enzymes was used to mask polymorphisms of partial exon 1 of GDF9 gene to establish molecular markers for twinning.

RESULTS

The lambing rate percentage and litter size showed significant difference between ewes, which produce single and twin lamb for each breed individually, whereas the coefficient of variation of the Saidi breed is greater than that of the Ossimi breed. The results suggested that the GDF9 gene shared a similarity in sequence compared to six accession numbers of Ovis aries found in GenBank. Molecular phylogenetic analyses were performed based on nucleotide sequences in order to examine the position of the Egyptian breeds among many other sheep breeds. The results indicate that accession number AF078545 of O. aries is closely related with Saidi and Ossimi ewes that produce single or twin lamb using the unweighted pair group method with arithmetic mean (UPGMA) analysis. Results showed that Msp1 enzyme digestion revealed polymorphic restriction pattern consisting of one band with 710 bp for ewes producing single lamb and two bands with 710 and 600 bp for ewes producing twin lamb in Saidi sheep breed.

CONCLUSION

Sequence analysis and diversity of polymorphisms in the GDF9 gene (exon 1) have a novel base substitution (A-T) for detection of Fec^G mutations that serve as a molecular marker for twinning.

Association between single nucleotide polymorphisms and sexual precocity in Nellore heifers

The aim of this study was to determine the extent (r²) of linkage disequilibrium (LD) in the genome of Nellore cattle, and to examine associations between single nucleotide polymorphisms (SNP) and age at first calving (AFC) and early pregnancy (EP) using a panel of high-density SNPs and data from 1182 Nellore females. A total of 13 contemporary groups (CG) were used consisting of farm, season, and year of birth. For genome-wide association analysis, SNPs with a minor allele frequency (MAF)<0.05 and animals with a call rate<0.90 were excluded, totaling 431,885 SNPs. For statistical analysis, a linear model was used for AFC and a threshold model for EP. To estimate the significance of the associations for the two traits, the model included the categorical fixed effects of CG, SNPs, and sire. In addition, the polygenic effect was included in the analysis. The additive effects and dominance deviations of Bonferroni-adjusted significant SNPs for AFC and EP were estimated using orthogonal contrasts. The average estimate of r² for all autosomes was 0.18 at a distance of 4.8kb and the mean MAF was 0.25±0.13. The LD decreased as the distance between markers increased: 0.35 (1kb) to 0.12 (100kb). Eleven significant associations were detected in seven different chromosomes. Seven SNPs were associated with AFC and four were associated with EP. Three SNPs were significant for both traits. The identification of SNPs associated with AFC and EP may contribute for selecting sexually precocious animals.

Genome-Wide Association Study for Indicator Traits of Sexual Precocity in Nellore Cattle

The objective of this study was to perform a genome-wide association study (GWAS) to detect chromosome regions associated with indicator traits of sexual precocity in Nellore cattle. Data from Nellore animals belonging to farms which participate in the DeltaGen^® and Paint^® animal breeding programs, were used. The traits used in this study were the occurrence of early pregnancy (EP) and scrotal circumference (SC). Data from 72,675 females and 83,911 males with phenotypes were used; of these, 1,770 females and 1,680 males were genotyped. The SNP effects were estimated with a single-step procedure (WssGBLUP) and the observed phenotypes were used as dependent variables. All animals with available genotypes and phenotypes, in addition to those with only phenotypic information, were used. A single-trait animal model was applied to predict breeding values and the solutions of SNP effects were obtained from these breeding values. The results of GWAS are reported as the proportion of variance explained by windows with 150 adjacent SNPs. The 10 windows that explained the highest proportion of variance were identified. The results of this study indicate the polygenic nature of EP and SC, demonstrating that the indicator traits of sexual precocity studied here are probably controlled by many genes, including some of moderate effect. The 10 windows with large effects obtained for EP are located on chromosomes 5, 6, 7, 14, 18, 21 and 27, and together explained 7.91% of the total genetic variance. For SC, these windows are located on chromosomes 4, 8, 11, 13, 14, 19, 22 and 23, explaining 6.78% of total variance. GWAS permitted to identify chromosome regions associated with EP and SC. The identification of these regions contributes to a better understanding and evaluation of these traits, and permits to indicate candidate genes for future investigation of causal mutations.

Functional Characterization of a Single Nucleotide Polymorphism in the 3' Untranslated Region of Sheep DLX3 Gene

The Distal-less 3 (homeobox protein DLX-3), a transcription factor, is critical for the development of hair follicle and hair formation and regeneration. We previously identified and found that four SNPs (c. *118T>C, c. *228T>C, c. *688A>G and c. *1,038_1,039 insC) in 3' untranslated region (UTR) of sheep DLX3 were in high linkage disequilibrium with each other and significantly associated with wool crimp (P<0.05), however, the underlying mechanisms by which these SNPs affect the wool crimp remains unknown. In the present study, we performed association analysis between these four identified SNPs and DLX3 gene expression in sheep skin using quantitative real-time RT-PCR. The results showed that these SNPs were significantly associated with sheep skin DLX3 mRNA expression levels. Then, we constructed DLX3 3'UTR luciferase reporters and validated the association. The reporter assays showed that the three major haplotypes, derived from the four SNPs, had significantly different effects on luciferase reporter activity and the four SNPs also had significantly different allelic effects on the luciferase reporter activity (p < 0.05). Bioinformatics analysis showed that the SNP (c. *1,038_1,039 insC) was located within a potential miR-188 binding site of the 3'UTR of sheep DLX3 mRNA. This SNP may affect miR-188-mediated DLX3 gene expression and result in phenotypic variation. To test the hypothesis, we investigated the effects of miR-188 mimic and inhibitor on the activity of the DLX3 3'UTR luciferase reporter with different SNP alleles. The results showed that in both sheep fetal fibroblasts (SFFs) and human HaCaT cells, miR-188 mimic could significantly decrease the allele D (deletion) luciferase reporter activity (p < 0.05), but miR-188 inhibitor could increased the reporter activitiy. However, neither miR-188 mimc nor inhibitor could influence the allele I (insertion) reporter activity. In addition, transfection of miR-188 mimic dramatically decreased the endogenous expression of DLX3 in SFFs (p < 0.05). Taken together, we demonstrated that DLX3 is a target gene of miR-188 and the SNP (c. *1,038_1,039 insC) is a functional SNP, and affects miR-188-mediated gene regulation of sheep DLX3. Our finding may in part explain allelic difference in gene expression and wool crimp in our tested sheep population.

A Major Gene for Bovine Ovulation Rate

Half-sib daughters sired by a bull believed to be a carrier of a major gene for high ovulation rate were evaluated for ovulation rate and genotyped in an effort to both test the hypothesis of segregation of a major gene and to map the gene’s location. A total of 131 daughters were produced over four consecutive years at a University of Wisconsin-Madison research farm. All were evaluated for ovulation rate over an average of four estrous cycles using transrectal ultrasonography. The sire and all daughters were genotyped using a 3K SNP chip and the genotype and phenotype data were used in a linkage analysis. Subsequently, daughters recombinant within the QTL region and the sire were genotyped successively with 50K and 777K SNP chips to refine the location of the causative polymorphism. Positional candidate genes within the fine-mapped region were examined for polymorphism by Sanger sequencing of PCR amplicons encompassing coding and 5’ and 3’ flanking regions of the genes. Sire DNA was used as template in the PCR reactions. Strong evidence of a major gene for ovulation rate was observed (p<1x10^-28) with the gene localized to bovine chromosome 10. Fine-mapping subsequently reduced the location to a 1.2 Mb region between 13.6 and 14.8 Mb on chromosome 10. The location identified does not correspond to that for any previously identified major gene for ovulation rate. This region contains three candidate genes, SMAD3, SMAD6 and IQCH. While candidate gene screening failed to identify the causative polymorphism, three polymorphisms were identified that can be used as a haplotype to track inheritance of the high ovulation rate allele in descendants of the carrier sire.

Effects of the FecL major gene in the Lacaune meat sheep population

Background

The major prolificacy gene FecL was first described in the Lacaune sheep meat breed Ovi-Test in 1998. A few studies estimated the effect of this gene on prolificacy but little data is available. In 2010, the Ovi-Test cooperative started genotyping FecL in all of their replacement ewe lambs. Thanks to the large amount of genotyping data that is available now, gene effects on litter size and other relevant traits can be estimated more accurately.

Methods

Our study included 5775 ewes genotyped since 2010 and 1025 sires genotyped since 2002. Performances and pedigrees were extracted from the French national database for genetic evaluation and research. Analysis of the effect of the gene on different traits was based on linear or threshold genetic animal models using the ASReml software.

Results

The female population was composed of 71% homozygous wild type ewes (++), 27% heterozygous ewes for the FecL mutation (L+) and 2% homozygous mutant (LL) ewes. On average, L + ewes produced 0.5 more lambs per lambing than ++ ewes. The FecL gene not only affected the mean litter size but also its variability, which was lower for ++ than for L + ewes. Fertility after insemination was higher for L + ewes than for ++ ewes. Lambs from ++ dams were heavier (+300 g) than the lambs of L + dams and the mortality of twin lambs born from ++ dams was lower than those from L + dams. In addition, bias in estimated breeding values for prolificacy when ignoring the existence of this major gene was quantified.

Conclusions

The effect of the FecL gene on prolificacy was estimated more accurately and we show that this gene affects both the mean and the variability of litter size and other traits. This paper also shows that ignoring the existence of this major gene in genetic evaluation of prolificacy can lead to a large overestimation of polygenic breeding values.

Using sheep lines with mutations in single genes to better understand ovarian function

Livestock populations have been subjected to strong selection pressure to improve reproductive success, and this has led to the identification of lines of animals with increased fecundity. These animals provide a rich biological resource for discovery of genes and regulatory mechanisms that underpin improved reproductive success. To date, three genes, all related to the transforming growth factor β pathway, have been identified as having mutations that lead to alterations in ovulation in sheep. In addition, several other sheep lines have been identified with putative mutations in single genes with major effects on ovulation rate. This review is focused on the identification of the mutations affecting ovulation rate and how these discoveries have provided new insights into control of ovarian function.

The highly prolific phenotype of Lacaune sheep is associated with an ectopic expression of the B4GALNT2 gene within the ovary

Prolific sheep have proven to be a valuable model to identify genes and mutations implicated in female fertility. In the Lacaune sheep breed, large variation in litter size is genetically determined by the segregation of a fecundity major gene influencing ovulation rate, named FecL and its prolific allele FecL(L) . Our previous work localized FecL on sheep chromosome 11 within a locus of 1.1 Mb encompassing 20 genes. With the aim to identify the FecL gene, we developed a high throughput sequencing strategy of long-range PCR fragments spanning the locus of FecL(L) carrier and non-carrier ewes. Resulting informative markers defined a new 194.6 kb minimal interval. The reduced FecL locus contained only two genes, insulin-like growth factor 2 mRNA binding protein 1 (IGF2BP1) and beta-1,4-N-acetyl-galactosaminyl transferase 2 (B4GALNT2), and we identified two SNP in complete linkage disequilibrium with FecL(L) . B4GALNT2 appeared as the best positional and expressional candidate for FecL, since it showed an ectopic expression in the ovarian follicles of FecL(L) /FecL(L) ewes at mRNA and protein levels. In FecL(L) carrier ewes only, B4GALNT2 transferase activity was localized in granulosa cells and specifically glycosylated proteins were detected in granulosa cell extracts and follicular fluids. The identification of these glycoproteins by mass spectrometry revealed at least 10 proteins, including inhibin alpha and betaA subunits, as potential targets of B4GALNT2 activity. Specific ovarian protein glycosylation by B4GALNT2 is proposed as a new mechanism of ovulation rate regulation in sheep, and could contribute to open new fields of investigation to understand female infertility pathogenesis.

Genetic control of multiple births in low ovulating mammalian species

In mammals, litter size is a highly variable trait. Some species such as humans or cattle are monotocous, with one or sometimes two newborns per birth, whereas others, the polytocous species such as mice or pigs, are highly prolific and often produce a dozen newborns at each farrowing. In monotocous species, however, two or three newborns per birth may sometime be unwanted. In more polytocous species such as sheep or pigs, litter size is studied in order to increase livestock prolificacy. By contrast, twinning rates in humans or cattle may increase birth difficulties and health problems in the newborns. In this context, the aim of our review was to provide a clearer understanding of the genetic and physiological factors that control multiple births in low-ovulating mammalian species, with particular focus on three species: sheep, cattle, and humans, where knowledge of the ovulation rate in one may enlighten findings in the others. This article therefore reviews the phenotypic and genetic variability observed with respect to ovulation and twinning rates. It then presents the QTL and major genes that have been identified in each species. Finally, we draw a picture of the diversity of the physiological mechanisms underlying multiple ovulation. Although several major genes have been discovered in sheep, QTL detection methods in humans or cattle have suggested that the determinism of litter size is complex and probably involves several genes in order to explain variations in the number of ovulations.

Identification of a line of sheep carrying a putative autosomal gene increasing ovulation rate in sheep that does not appear to interact with mutations in the transforming growth factor beta superfamily

Sheep lines with mutations in single genes that have major effects on ovulation rate have been very useful in gaining a better understanding of pathways important in controlling follicular development and ovulation rate. To date however, all known mutations are in the transforming growth factor beta (TGFB) superfamily. Ovulation rates were measured in 720 progeny of 20 rams that were descendants of a single prolific ewe. Evaluation of ovulation rates of daughters of closely related sires suggests the presence of a segregating major gene Fecundity Davisdale (FECD) that increases ovulation rate between 0.4 and 0.8 in heterozygous daughters. Key features of mutations in genes of the TGFB superfamily pathway, such as synergistic interactions with other family members, infertility in homozygous carriers, and increased responsiveness to exogenous gonadotropins, were not observed in this line; thus, the mutation does not appear to be acting in the TGFB pathway. Hence, there is likely a novel mutation being carried in this line of sheep that alters ovulation rate. Future identification of the causative mutation may provide new insights into regulation of follicular development and ovulation rate.