Converting multi-trait breeding objectives into operative selection indexes to ensure genetic gains in low-input sheep and goat breeding programmes

Small ruminant breeding programmes in low-input production systems are best organised at the community level. Participant farmers have to agree on goal traits and their relative importance. When BLUP breeding values of goal traits are not available in time, appropriate selection indexes can be used to aid visual selection. Taking Ethiopian Abergelle goat and Bonga sheep community-based breeding programmes (CBBPs) as an example, breeding objective functions were defined and selection indexes were constructed and evaluated. Breeding goals for Abergelle goats included early sale weight, survival and milk production. Breeding goals for Bonga included the number of offspring born, sale weight and survival. Economic weights of objective traits can be used in several ways depending on measured traits and the reliability of their genetic parameters. Selection indexes included combinations of objective traits measured on candidates and their dams and situations when Abergelle communities prefer to restrict genetic changes in number of offspring born or adult weight and when Bonga communities prefer to restrict changes in adult weight. Genetic and economic gains were evaluated as well as sensitivity to feed cost assumptions and to repeated dam records. After independent culling on preponderant traits such as coat colour and horn/tail type, sires in Abergelle goat community breeding programmes should be selected on indexes including at least own early live weight and their dams average milk production records. Sires for Bonga sheep programmes should be selected on own early live weight and desirably also on their dam’s number of offspring born. Sensitivity to feed cost assumptions was negligible but repeated measurements of dam records improved index accuracies considerably. Restricting genetic changes in number of offspring born or adult weight is not recommended.     

Multiple genomic signatures of selection in goats and sheep indigenous to a hot arid environment

Goats and sheep are versatile domesticates that have been integrated into diverse environments and production systems. Natural and artificial selection have shaped the variation in the two species, but natural selection has played the major role among indigenous flocks. To investigate signals of natural selection, we analyzed genotype data generated using the caprine and ovine 50K SNP BeadChips from Barki goats and sheep that are indigenous to a hot arid environment in Egypt''s Coastal Zone of the Western Desert. We identify several candidate regions under selection that spanned 119 genes. A majority of the genes were involved in multiple signaling and signal transduction pathways in a wide variety of cellular and biochemical processes. In particular, selection signatures spanning several genes that directly or indirectly influenced traits for adaptation to hot arid environments, such as thermo-tolerance (melanogenesis) (FGF2, GNAI3, PLCB1), body size and development (BMP2, BMP4, GJA3, GJB2), energy and digestive metabolism (MYH, TRHDE, ALDH1A3), and nervous and autoimmune response (GRIA1, IL2, IL7, IL21, IL1R1) were identified. We also identified eight common candidate genes under selection in the two species and a shared selection signature that spanned a conserved syntenic segment to bovine chromosome 12 on caprine and ovine chromosomes 12 and 10, respectively, providing, most likely, the evidence for selection in a common environment in two different but closely related species. Our study highlights the importance of indigenous livestock as model organisms for investigating selection sweeps and genome-wide association mapping.     

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.     

Selection Practices of Bonga Sheep Reared in Southern Ethiopia

The study was conducted to assess the traits traditionally used to select the Bonga sheep reared in Southern Ethiopia. The breed was included for improvement under the community based breeding program (CBBP) launched in the year 2009. The results are based on the data collected from the project between 2009 till 2012. The findings are based on focus group discussions with the community elders who have been rearing the Bonga sheep for more than 20 years. The present findings are based on information obtained from 50% of the total respondents who were involved in CBBP. The Bonga sheep is a mutton type breed and the ewes are moderately prolific. In the past this breed of sheep were bartered to settle legal disputes and were traded against household commodities.. However, over the years (and also after the intervention of CBBP) the sheep have been selected for their body weight as well as for their distribution in many parts of Ethiopia The respondents selected the lambs at both pre weaning and post weaning stages. Traditionally the traits of rams are selected based on their body length, canon circumference, broad face, enlarged thyroid, while for the ewes, traits prolificacy, skin thickness and pelvic width are determining characters. Fat tail and brown coat color were preferred irrespective of the sexes in trait selection. The ram lambs selected for breeding purpose are locally known as “Dookoo”. These rams are selected based on some predetermined phenotypic traits and are initially selected at preweaning stage and further, again at the post weaning stage. These rams are preferentially cared and are provided with supplementary feed and comfortable housing.     

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.     

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.     

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.     

Community-based sheep breeding programs generated substantial genetic gains and socioeconomic benefits

Community-based breeding programs (CBBPs) for small ruminants have been suggested as alternatives to centralised, government-controlled breeding schemes which have been implemented in many developing countries. An innovative methodological framework on how to design, implement and sustain CBBPs was tested in three sites in Ethiopia: Bonga, Horro and Menz. In these CBBPs, the main selection trait identified through participatory approaches was 6-month weight in all three sites. In Horro and Bonga, where resources such as feed and water permitted larger litter sizes, twinning rate was included. Ten-year (2009 to 2018) performance data from the breeding programs were analysed using Average Information Restricted Maximum Likelihood method (AI-REML). Additionally, the socioeconomic impact of CBBPs was assessed. Results indicated that 6-month weight increased over the years in all breeds. In Bonga, the average increase was 0.21 ± 0.018 kg/year, followed by 0.18 ± 0.007 and 0.11 ± 0.003 kg/year in Horro and Menz, respectively. This was quite substantial in an on-farm situation. The birth weight of lambs did not improve over the years in Bonga and Horro sheep but significant increases occurred in Menz. Considering that there was no direct selection on birth weight in the community flock, the increased weights observed in Menz could be due to correlated responses, but this was not the case in Bonga and Horro. The genetic trend for prolificacy over the years in both Bonga and Horro flocks was positive and significant (P < 0.01). This increase in litter size, combined with the increased 6-month body weight, increased income by 20% and farm-level meat consumption from slaughter of one sheep per year to three. The results show that CBBPs are technically feasible, result in measurable genetic gains in performance traits and impact the livelihoods of farmers.     

Estimation of (co)variance components,genetic parameters,and genetic trends of growth traits in community-based breeding programs of Bonga sheep

Community-based sheep breeding programs (CBBPs) have been adopted strategically to improve Bonga sheep, the most popular sheep breed in Ethiopia. The present study was undertaken to estimate genetic parameters and genetic trends for growth traits and inbreeding levels in each Bonga sheep CBBP. Data pertaining to growth traits, spanning a period of seven years (2012–2017), were collected from 14 Bonga sheep CBBPs. Data were analyzed using the General Linear Model procedure of SAS to study the performance of the breed over the years. The genetic parameters were estimated by univariate and multivariate animal model using restricted maximum likelihood method of WOMBAT software. The genetic trends were estimated by the regression of the average breeding values of the animals on the year of birth. The overall least square means ± SE of BW (kg) were 3.10 ± 0.010, 16.1 ± 0.07, 24.7 ± 0.20, 30.4 ± 0.40 and 34.0 ± 0.84 for birth weight (BWT), weaning weight (WWT), six-month weight (SMWT), nine-month weight (NMWT) and yearling weight (YWT), respectively. Direct heritability estimates from selected models were 0.56 ± 0.030, 0.36 ± 0.030, 0.22 ± 0.040, 0.17 ± 0.070 and 0.13 ± 0.150 for BWT, WWT, SMWT, NMWT and YWT, respectively. Six-month weight was the selection trait and presented positive trends for 10 CBBPs, and negative trends for four CBBPs. Moderate to high heritability estimates and positive genetic trends indicated scope for further improvement of BW. Additionally, the positive and high correlation between BW traits indicated that selection for just one trait would also improve the other traits through correlated responses.     

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.     

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.     

Genotyping of Novel SNPs in BMPR1B,BMP15, and GDF9 Genes for Association with Prolificacy in Seven Indian Goat Breeds

Goats form the backbone of rural livelihood and financial security systems in India but their population is showing decreasing trend. Improvement of reproductive traits such as prolificacy offers a solution to stabilize the decreasing goat population and to meet the nutritional needs of growing human population. In the present study, six novel SNPs in three candidate genes for prolificacy (BMPR1B, BMP15, and GDF9) were genotyped in seven breeds of Indian goats to evaluate their association with litter size. Tetra primer ARMS-PCR and PCR-RFLP based protocols were developed for genotyping six novel SNPs, namely, T(-242)C in BMPR1B; G735A and C808G in BMP15; and C818T, A959C, and G1189A in GDF9 gene. The effect of breed was highly significant (p ≤ 0.01) on litter size but the effect of genotype was nonsignificant. The effect of parity on litter size was also significant in the prolific Black Bengal breed. The litter size differences observed between breeds are attributed to breed differences. Novel mutations observed at different loci in GDF9, BMP15, and BMPR1B genes do not contribute to the reproductive capability of the investigated breeds. Further studies with more number of breeds and animals exploring association of these novel SNPs with reproductive traits may be fruitful.     

Ovulation level and prolificacy in ewes depending on their age, birth type and percentage of prolific genotype

During December mating periods of 2000-2003, wool and meat breeds crossed with East Friesian milk sheep and rams of prolific breeds (Finn or Romanov) were examined laparoscopically. Ovine ovaries were examined between 4 and 9 days after mating. The litter size was recorded and compared with the number of ovulations. The fertility of the analyzed animals was 69.6%, with 192.6% prolificacy. The difference between the average number of corpora lutea found on ovaries and the mean number of lambs born was 0.35 and ranged from 0.32 to 0.41 according to age group, from 0.22 to 0.41 according to sheep birth type, and from 0 to 0.52 lamb per ewe according to percentage of prolific breeds.     

Freemartinism and FecX allele determination in replacement ewes of the Rasa Aragonesa sheep breed by duplex PCR

A new naturally occurring mutation in the fecundity gene BMP15 in the Rasa Aragonesa sheep breed (Ovis aries) has been found to affect prolificacy. This mutation (FecX^R allele) is a deletion of 17 base pairs that leads to an altered amino acid sequence, and this alteration increases prolificacy in heterozygous ewes but causes sterility in homozygous ewes. Selection of repository lambs with the FecX^R allele increases rates of twins and multiple lambing and thereby also increases the probability of lambing freemartins that will become sterile. In this sense, an accurate, reliable, and quick method was developed by duplex polymerase chain reaction (PCR) for sex, amplifying an ovine-specific Y chromosome repetitive fragment, and BMP15 genotype determination in replacement ewe lambs. The BMP15 fragment served as an internal control of the amplification and detected the FecX^R allele, avoiding a false negative and then a mistake in freemartin detection. This assay uncovered 6 freemartin females among 195 replacement ewes from 7 different commercial flocks and 1 experimental flock. Furthermore, 1554 rams from 64 commercial flocks were also analyzed to identify FecX^R rams. This analysis identified 103 rams hemizygous for the FecX^R allele and 1 heterozygous ram. Because this gene is located on the X chromosome, this heterozygous animal is a freemartin ram that is co-amplifying the DNA from XX and XY lymphocytes. These results confirm the usefulness of this multiplex PCR assay for detecting phenotypically sexed females, freemartins, and the BMP15 genotype to detect highly prolific ewes in commercial flocks and to assist breeders in selection of repository lambs.     

BMPR-IB和BMP15基因作为小尾寒羊多胎性能候选基因的研究

以控制BooroolaMerino羊多胎性能的BMPR IB基因 ,以及影响Invedale和Hanna羊排卵数的BMP15基因作为候选基因 ,从分子水平上对小尾寒羊的多胎机制进行研究 ,分析突变位点的特性 ,并通过大规模的群体检测统计推断其遗传效应。实验结果表明 :多胎品种小尾寒羊在BMPR IB基因的相应位置上发生了与BooroolaMerino羊相同的突变 (A74 6G) ,该基因的BB基因型在小尾寒羊群体内为优势基因型 ,且小尾寒羊初产和经产母羊的BB基因型比 ++基因型分别多产 0 97羔 (P <0 0 5 )和 1 5羔 (P <0 0 1) ,推测BMPR IB基因与控制小尾寒羊多胎性能的主效基因存在紧密的遗传连锁。而BMP15基因在小尾寒羊中不存在V31D或Q2 3Ter突变 ,说明小尾寒羊的多胎遗传机制与Romney羊不同 ,因此排除了BMP15突变影响小尾寒羊排卵数的可能性。     

贵州白山羊Bmp15和Gdf9基因多态性及其编码蛋白的进化分析(英文)

BMP15和GDF9是转化生长因子β（TGFβ）超家族的成员,对绵羊的繁殖性状有直接的调节作用,从中发现的多个高产突变位点直接提高了排卵数和产羔数。在之前的研究中,作者从贵州白山羊中找到了一个高产突变位点。为了进一步揭示Bmp15和Gdf9基因突变与繁殖性状之间的关系,对贵州白山羊Bmp15和Gdf9基因编码区进行了克隆,以人BMP7的晶体结构为模板构建了贵州白山羊BMP15和GDF9成熟肽的三维模型。贵州白山羊Bmp15和Gdf9基因分别编码394和453个氨基酸的蛋白前体。对BMP15和GDF9成熟肽序列进行分析发现,除了之前确认的BMP15中的FecX^B 突变（S99I）和GDF9中的V79I突变之外,还从贵州白山羊的BMP15和GDF9成熟肽分别发现7个和3个位点突变。其中,BMP15成熟肽的S32G、N66H、S99I/P99I和G107R突变可能影响二聚体与受体的结合;GDF9成熟肽的P78Q和V79I影响二聚体与I型受体的亲和力,将值得进一步深入研究。对Bmp15和Gdf9基因编码的蛋白前体序列进行聚类分析,结果显示在鱼类到哺乳类的进化过程中,BMP15出现长度逐渐增加的现象,以BMP15成熟肽N端长度增加为主。这种演变可能使BMP15对低排卵哺乳动物繁殖力的控制更为灵敏。该文的研究结果为贵州白山羊Bmp15和Gdf9基因变异与繁殖力的关系提出了合理的解释,并支持这两个因子是贵州白山羊高产性状重要调节因子的观点。     

A genome scan for selection signatures in Taihu pig breeds using next-generation sequencing

Taihu pig breeds are the most prolific breeds of swine in the world, and they also have superior economic traits, including high resistance to disease, superior meat quality, high resistance to crude feed and a docile temperament. The formation of these phenotypic characteristics is largely a result of long-term artificial or natural selection. Therefore, exploring selection signatures in the genomes of the Taihu pigs will help us to identify porcine genes related to productivity traits, disease and behaviour. In this study, we used both intra-population (Relative Extend Haplotype Homozygosity Test (REHH)) and inter-population (the Cross-Population Extend Haplotype Homozygosity Test (XPEHH); F-STATISTICS, F_ST) methods to detect genomic regions that might be under selection process in Taihu pig breeds. As a result, we found 282 (REHH) and 112 (XPEHH) selection signature candidate regions corresponding to 159.78 Mb (6.15%) and 62.29 Mb (2.40%) genomic regions, respectively. Further investigations of the selection candidate regions revealed that many genes under these genomic regions were related to reproductive traits (such as the TLR9 gene), coat colour (such as the KIT gene) and fat metabolism (such as the CPT1A and MAML3 genes). Furthermore, gene enrichment analyses showed that genes under the selection candidate regions were significantly over-represented in pathways related to diseases, such as autoimmune thyroid and asthma diseases. In conclusion, several candidate genes potentially under positive selection were involved in characteristics of Taihu pig. These results will further allow us to better understand the mechanisms of selection in pig breeding.     

Improving productivity of Malpura breed by crossbreeding with prolific Garole sheep in India

Garole, a prolific small sized sheep breed of India was introduced in the Mutton Project of the Institute during 1997 to explore the possibility of incorporating fecundity gene into Malpura (M) sheep of semi-arid region of Rajasthan. The data on reproduction and production traits of M and Garole × Malpura (GM) half-breds was used in the study to assess the improvement in evolved crossbred stock. The overall percent gain in prolificacy in GM was 52.38, which increased to 75.73 in third parity. The GM compared to M ewes also weaned more lambs per ewe lambed (1 versus 1.41). The lambing rate averaged 91.75 and 148.21% in M and GM ewes, respectively. The twin lambing percentage was only 4.71% in M flock while; it reached up to 51.10% in GM. Moreover, 5.72% of GM ewes produced triplets also. The genetic group non-significantly affect ewe efficiency (EE), but GM produced 3.3% more total lambs weight per ewe lambed at weaning and 6.15% more total lamb weight per ewe lambed at 6-month of age. The body weights at different ages of GM half-breds were lower (P < 0.01) compared to contemporary M lambs. The parity of ewes significantly affects body weight from birth weight to 12-month of age. The survivability of GM half-breds was almost at par with the local M sheep. Results revealed that the fecundity genes responsible for increasing prolificacy have been incorporated into the GM and it might prove a valuable germplasm towards evolving a new prolific strain of sheep. The crossing of Garole with M increased the ewe efficiency.