Direct responses to six generations of selection for ovulation rate or prenatal survival in Large White pigs

Mechanisms of Oogenesis-Related Long Non-coding RNAs in Porcine Ovaries Treated With Recombinant Pig Follicle-Stimulating Hormone

Reproductive efficiency is of significant importance in pork production for it has a great impact on economic success. Ovulation rate is an early component of reproduction efficiency of pigs, and it contributes to the upper limit of litter size. In this study, we used the newly developed recombinant pig follicle stimulating hormone (rpFSH) instead of traditional PMSG to increase ovulation rate of pigs in order to achieve higher litter size, for it was better at stimulating ovulation, and showed more cheaper and greener. However, relatively little is known about the underlying genetic bases and molecular mechanisms. Consequently, an experiment was carried out in ovaries of replacement gilts to screen the key genes and lncRNAs that affect the fecundity of pigs by RNA-seq technology. Twenty gilts were divided into two groups, including 10 rpFSH treatment pigs and 10 control animals. After slaughtering and collecting the phenotypic data, ovaries of five pigs in each group were selected for RNA-seq. Total RNA was extracted to construct the library and then sequence on an Illumina Hiseq 4000 system. A comprehensive analysis of mRNAs and long non-coding RNAs (lncRNAs) from 10 samples was performed with bioinformatics. The phenotypic data showed that rpFSH treatment groups had the higher (P < 0.01) ovarian weight and more mature follicles. The RNA-seq results showed that a total of 43,499 mRNAs and 21,703 lncRNAs were identified, including 21,300 novel lncRNAs and 403 known lncRNAs, of which 585 mRNAs and 398 lncRNAs (P < 0.05) were significantly differentially expressed (DE) between the two groups of rpFSH treatment group and controlled group. GO and KEGG annotation analysis indicated that the target genes of DE lncRNAs and DE mRNAs were related to prolactin receptor activity, mitophagy by induced vacuole formation, and meiotic spindle. Moreover, we found that NR5A2 (nuclear receptor subfamily 5, group A, member 2), a target gene of lncRNA MSTRG.3902.1, was involved in regulating follicular development, ovulation, and estrogen production. Our study provided a catalog of lncRNAs and mRNAs associated with ovulation of rpFSH treatment, and they deserve further study to deepen the understanding of biological processes in the regulation of ovaries of rpFSH treatment pigs.

Litter size components traits in two Algerian rabbit lines

<p>The aim of this study was to estimate the limiting litter size components in rabbit females from a Synthetic line (n=32) and a Local population (n=34). Ovulation rate, number of implanted and live embryos were counted by laparoscopy at 12 d after mating. Prolificacy (total newborn, number born alive and mortality) and embryonic, foetal and prenatal survival at day of birth of the 3^rd gestation were measured. The analysed traits were body weight of the female at mating, ovulation rate, implanted, live and resorbed embryos, embryonic, foetal and prenatal survival, as well as total newborn, number born alive and mortality at birth. Synthetic line females had a higher ovulation rate compared to the Local population (11.03±0.23 vs. 8.41±0.23 corpora lutea; <em>P</em>&lt;0.0001). Synthetic line displayed a higher number of implanted embryos (10.00±0.25 vs. 7.85±0.25 embryos; <em>P</em>&lt;0.0001). No difference was found between groups for number of resorbed embryos. Similar embryonic, foetal and prenatal survival rates were reported between the Synthetic line and the Local population. Additionally, total newborn was higher in the Synthetic line than in the Local population (+1.46 kits; <em>P</em>&lt;0.05). A principal components analysis was performed. The first four principal components (PC) explained more than 90% of the total variation in both lines. Total newborn, number born alive and live embryos were the main variables defining the 1^st PC. Resorbed embryos and foetal survival were located in the 2^nd PC. Ovulation rate and embryonic survival were the predominant variables defining the 3^rd PC. The body weight of females was located in the 4^th PC. The phenotypic correlation between total newborn and its components were high and positive in both lines, except for ovulation rate and total newborn, where it was moderate in Synthetic line. In conclusion, the females from Synthetic line have a higher total newborn than those from Local population, as a consequence of a higher number of released oocytes and embryos that successfully reach implantation. However, a higher uterine crowding in Synthetic line seems to limit survival of foetuses that reach term of gestation, while ovulation rate is the principal limiting factor of total newborn in Local population.</p>

Ovarian Transcriptomic Analysis Reveals Differential Expression Genes Associated with Cell Death Process after Selection for Ovulation Rate in Rabbits

Simple Summary

Transcriptomic analysis showed nineteen potential biomarkers in ovarian tissue from females belonged to a rabbit line selected for ovulation rate for 10 generations and the control line. These females differed not only in ovulation rate but also in prenatal survival since similar litter size were observed.

Abstract

Litter size is an essential trait in rabbit meat production but with low heritability. A selection experiment for ovulation rate has been performed for 10 generations to improve litter size in rabbits. The selected line increased two ova more than the control line but nevertheless a negative correlation was observed with prenatal survival. A transcriptomic study was performed, using microarrays, in ovarian tissue from females belonging to the selected line and the control line. Our results showed 1357 differential expressed genes and nineteen potential biomarkers associated with prenatal mortality, which could explain differences between litter size in rabbits. Cell death was the most relevant process.

Systematic Analysis of Long Non-Coding RNAs and mRNAs in the Ovaries of Duroc Pigs During Different Follicular Stages Using RNA Sequencing

The dynamic process involving the selection and maturation of follicles is regulated and controlled by a highly synchronized and exquisitely timed cascade of gene expression. Studies have shown that long non-coding RNA (lncRNA) is essential for the normal maintenance of animal reproductive function and has an important regulatory function in ovarian development and hormone secretion. In this study, a total of 2076 lncRNAs (1362 known lncRNAs and 714 new lncRNAs) and 25,491 mRNAs were identified in libraries constructed from Duroc ovaries on days 0, 2 and 4 of follicle development. lncRNAs were shorter, had fewer exons, exhibited a shorter ORF (Open Reading Frame) length and lower expression levels, and were less conserved than mRNAs. Furthermore, 1694 transcripts (140 lncRNAs and 1554 mRNAs) were found to be differentially expressed in pairwise comparisons. A total of 6945 co-localized mRNAs were detected in cis in 2076 lncRNAs. The most enriched GO (Gene Ontology) terms were related to developmental processes. KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis revealed that the differentially expressed lncRNAs targeted mRNAs, and the differentially expressed mRNAs were related to the TGF-β signaling pathway, the PI3K-Akt signaling pathway, the Retinol metabolic pathway and the Wnt signaling pathway. This study deepened our understanding of the genetic basis and molecular mechanisms of follicular development in pigs.

Relationship Between the Estimated Breeding Values for Litter Traits at Birth and Ovarian and Embryonic Traits and Their Additive Genetic Variance in Gilts at 35 Days of Pregnancy

We investigated (1) the relationship between the estimated breeding values (EBVs) for litter traits at birth and ovulation rate (OR), average corpora luteal weight, uterine length and embryonic survival and development traits in gilts at 35 days of pregnancy by linear regression, (2) the genetic variance of OR, average corpora lutea (CL) weight, uterine length and embryonic survival and development traits at 35 days of pregnancy, and (3) the genetic correlations between these traits. Landrace (n = 86) and Yorkshire × Landrace (n = 304) gilts were inseminated and slaughtered at 35 days of pregnancy. OR was assessed by dissection of the CL on both ovaries. Individual CL was weighed and the average CL weight calculated. The number of embryos (total and vital) were counted and the vital embryos were individually weighed for calculation of within litter average and standard deviation (SD) of the embryo weight. Length of the uterine implantation site of the vital embryos was measured and the average per gilt calculated. Results suggests that increasing the EBV for total number of piglets born would proportionally increase OR and number of embryos, while decreasing the average CL weight. On the contrary, increasing the EBV for average piglet birth weight and for within litter birth weight standard deviation would increase the average CL weight. There was no relationship between the EBVs for BW and for BWSD and vital embryonic weight at 35 days of pregnancy. OR, average CL weight, number of embryos, average weight and implantation length of the vital embryos had all moderate to high heritabilities, ranging from 0.36 (±0.18) to 0.70 (±0.17). Thus, results indicate that there is ample genetic variation in OR, average CL weight and embryonic development traits. This knowledge could be used to optimize the balance between selection for litter size, average piglets birth weight and within litter birth weight uniformity.

Correlated response in litter size components in rabbits selected for litter size variability

A divergent selection experiment for the environmental variability of litter size (Ve) over seven generations was carried out in rabbits at the University Miguel Hernández of Elche. The Ve was estimated as the phenotypic variance within the female, after correcting for year-season and parity-lactation status. The aim of this study was to analyse the correlated responses to selection in litter size components. The ovulation rate (OR) and number of implanted embryos (IE) in females were measured by laparoscopy at 12 day of the second gestation. At the end of the second gestation, the total number of kits born was measured (TB). Embryonic (ES), foetal (FS) and prenatal (PS) survival were computed as IE/OR, TB/IE and TB/OR, respectively. A total of 405 laparoscopies were performed. Data were analysed using Bayesian methodology. The correlated response to selection for litter size environmental variability in terms of the litter size components was estimated as either genetic trends, estimated by computing the average estimated breeding values for each generation and each line, or the phenotypic differences between lines. The OR was similar in both lines. However, after seven generations of selection, the homogenous line showed more IE (1.09 embryos for genetic means and 1.23 embryos for phenotypic means) and higher ES than the heterogeneous one (0.07 for genetic means and 0.08 for phenotypic means). The probability of the phenotypic differences between lines being higher than zero (p) was 1.00 and .99, respectively. A higher uterine overcrowding of embryos in the homogeneous line did not penalize FS; as a result, this line continued to show a greater TB (1.01 kits for genetic means and 1.30 kits for phenotypic means, p = .99, in the seventh generation). In conclusion, a decrease in litter size variability showed a favourable effect on ES and led to a higher litter size at birth.

Estimation of the effects of selection on French Large White reproductive performance using frozen semen

Genetic trends for male and female sexual development and components of litter size and weight at birth (LB) as well as traits related to sow lifetime productivity were estimated in a French Large White (LW) pig population using frozen semen. Two experimental groups (EXP = L77 and L98) were produced by inseminating LW sows with either stored frozen semen from 17 LW boars born in 1977 (EXP = L77) or with semen from 23 LW boars born in 1998 (EXP = L98). In each group, about 100 males were measured for testicular development and bulbo-urethral gland development, and 90 gilts were checked for puberty, among which about 50 gilts were measured for their first ovulation rate at puberty (ORP). The females were then mated to 15 randomly chosen males for each group for 6 successive parities and measured for total number born (TNB), number born alive (NBA), litter weight at birth (LWB), and interval from weaning to first estrus as well as ovulation rate at fertilization (ORF = ovulation rate of gestating sows) on a subset of litters. The data were analyzed using mixed linear models including the fixed effects of EXP, of finishing (sexual development traits) or farrowing group, and parity (LB traits); the random effects of birth litter (sexual development traits) or sow (LB traits); and, when relevant, age or BW as covariates. The homogeneity of residual variances across EXP was also tested. For each trait, the genetic trend was estimated as twice the difference between L98 and L77. Limited differences were observed for sexual development traits, except for ORF, which increased by 0.17 ± 0.06/yr. The residual SD of NBA (2.7 vs. 2.0; = 0.004) and LWB (3.1 vs. 2.5; = 0.03) were larger in L98 sows than in L77 sows. Positive trends were observed for all litter traits (0.22 ± 0.07 ova, 0.12 ± 0.04 piglet, 0.09 ± 0.04 piglet, and 0.23 ± 0.06 kg/yr for ORF, TNB, NBA, and LWB, respectively) but at the expense of a degradation of prenatal and farrowing survival. A negative trend was also observed for length of productive life (-8.9 ± 3.9 d/yr). These results suggest that gains in productivity were associated with a loss of robustness.

Effect of increased ovulation rate on embryo and foetal survival as a model for selection by ovulation rate in rabbits

Selection for ovulation rate in prolific species has not improved litter size, due to an increase in prenatal mortality, with most mortality observed in the foetal period. The aim of this study was to investigate the magnitude and timing of embryo and early foetal survival in females with high ovulation rate using hormonal treatment as a model for selection by ovulation rate. Two groups of females (treated and untreated) were used. Treated females were injected with 50 IU equine chorionic gonadotropin 48 h before mating. Females were slaughtered at 18 d of gestation. Ovulation rate (OR), number of implanted embryos (IE), number of live foetuses at 12 and 18 d (LF₁₂ and LF₁₈, respectively) were recorded. In addition, embryo survival (ES=IE/OR), foetal survival at 18 d of gestation (FS_LF18=LF₁₈/IE), foetal survival between 12 and 18 d of gestation (FS_LF18/LF12=LF₁₈/LF₁₂) and prenatal survival (PS_LF18=LF₁₈/OR) were estimated. For each female, the mean and variability of the weight for live foetuses (LFWm and LFWv, respectively) and their placentas (LFPWm and LFPWv, respectively) were calculated. Treated females had a higher ovulation rate (+3.02 ova) than untreated females, with a probability of 0.99. An increase in the differences (D) between treated and untreated females was observed from implantation to 18 d of gestation (D=–0.33, –0.70 and –1.28 for IE, LF₁₂ and LF₁₈, respectively). These differences had a low accuracy and the probability that treated females would have a lower number of foetuses also increased throughout gestation (0.60, 0.70 and 0.86 for IE, LF₁₂ and LF₁₈, respectively). According to the previous results for OR and LF₁₈, treated females showed a lower survival rate from ovulation to 18 d of gestation (D=–0.12, P=0.98 for PS_LF18). Treated females also had lower embryo and foetal survival (D=–0.10 and P=0.94 for ES and D=–0.08 and P=0.93 for FS_LF18). Main differences in foetal survival appeared from 12 to 18 d of gestation (D=–0.09 and P=0.98 for FS_LF18/LF12). Unexpectedly, treated females showed similar foetus weight and higher foetal placenta weight than untreated females (D=0.25 g, P=0.98) and lower variability for these traits (D=–0.02 g, P=0.72 for LFWv and D=–0.05 g, P=0.83 for LFPWv). These results are not related to a lower number of IE or LF₁₈. Thus, the effect of increasing by three ova in rabbits leads to a lower embryo and early foetal survival. There seems to be no relationship between foetal mortality and foetus weight.

Analysis of 28 generations of selection for reproduction, growth, and carcass traits in swine

Selection (28 generations, G) in a Large White-Landrace composite population for traits aimed at increasing live pigs born per litter (BA), with additional selection for increased 180-d weight (WT180) and longissimus muscle area (LMA) and decreased back fat (BF10) in the last 8 generations, was practiced. Objectives herein were to estimate genetic and phenotypic responses and genetic parameters (n = 1,883 to 54,174) and to investigate whether a plateau in response for BA occurred. Line 2 (L2) was selected for an index of ovulation rate and embryo survival (G0 to G11), fully formed pigs (FF) per litter (G12 to 14), and BA and pig birth weight (PBW, G15 to G19), and its control line (LC1) was selected randomly (G0 to G21). Line 4 (L4), derived from L2, and line 5 (L5), derived from LC1, at G8 were selected in 2 stages for ovulation rate and FF (G9 to G16) and BA and PBW (G17 to G19), and their control (LC6) was selected randomly. At G20, L4 and L5 were crossed to form L45, and L4 and L2 were crossed to continue L2; L2 and L45 were subsequently selected for BA, WT180, LMA, and BF10 (G21 to G28). At G21, LC1 and LC6 were reciprocally crossed to form LC16, control for L2, and LC61, control for L45. Selection in L2 and L45 was first for BA and then for other traits among pigs selected for BA. Line sizes were 40 to 60 litters by 15 to 20 sires/G. Cumulative selection differentials (CSD) were calculated. MTDFREML was used to estimate variance components, EBV, and responses. Genetic changes at G28 in L2 were 4.63 FF and 3.66 BA, with 72% (FF) and 86% (BA) of the change occurring after G11. Two-stage selection produced similar responses (P < 0.01) in FF in L4 and L5 (0.27 and 0.29 pigs/G) but a greater response in BA in L5 (0.19 vs. 0.28 pigs/G). Genetic change in L45 from G20 to G28 was 0.17 pigs/G for both FF and BA (P < 0.01). Genetic changes at G28 in L45 were 4.16 FF and 3.68 BA. Genetic correlations of reproductive and growth traits were near zero, ranging from -0.43 (stillborn pigs/litter with BF10) to 0.21 (mummies/litter with LMA). Selection for growth traits along with litter size selection during G19 to G28 resulted in responses consistent with the selection applied and the heritability of the trait. No evidence for a selection plateau existed; selection differentials and variances of FF and BA in selection lines during G20 to G28 were similar to those in earlier generations. Over all generations, heritability of BA was 0.20 ± 0.03 and remained at approximately 0.17 in selection lines in later generations.

Incorporation of genetic technologies associated with applied reproductive technologies to enhance world food production

Animal breeding and reproductive physiology have been closely related throughout the history of animal production science. Artificial insemination provides the best method of increasing the influence of sires with superior genetics to improve production traits. Multiple ovulation embryo transfer (MOET) provides some ability to increase the genetic influence of the maternal line as well. The addition of genetic technologies to this paradigm allows for improved methods of selecting sires and dams carrying the best genes for production and yield of edible products and resistance to diseases and parasites. However, decreasing the number of influential parents within a population also increases the risk of propagating a recessive gene that could negatively impact the species (Reprod Domest Anim 44:792-796, 2009; BMC Genomics 11:337, 2010). Furthermore, antagonistic genotypic relationships between production traits and fertility (Anim Prod Sci 49:399-412, 2009; Anim Genet 43:442-446, 2012) suggest that care must be taken to ensure that increasing the frequency of genes with a positive influence on production does not negatively impact the fertility of the replacement females entering the herd.

Genetic and non-genetic parameters related to embryo production in superovulated Large White (LW) gilts

The aim of this study was to identify genetic and non genetic factors which might affect results of embryo production of Large White (LW) cyclic gilts from data collected in one herd during 6 years. Donors (n=1060) were synchronized with a progestogen treatment and luteolysis was induced 13-15 days later by 2 injections of cloprostenol. To stimulate follicular development 800IU eCG was then injected 24h later, followed by 500IU hCG 48h later. Donors were inseminated twice; depending on the onset of oestrus, the interval between hCG treatment and first insemination (hCGAI1) was either 24 or 41 h. Embryos were collected at 5-6 days after the 1st AI by flushing uterine horns. Traits of interest were the number of corpora lutea (CL), the number of flushed embryos (FE), the number of transferable embryos (TE) and the number of unfertilized embryos (UE). The average number of TE was 18.8 ± 9.0. The main sources of variation for CL, FE and TE were the season (P≤0.002) and hCGAI1 (P≤0.001) effects. For the interval of 24h of hCGIA1 the number of TE was increased by 4 compared with the TE obtained for the 41 h interval of hCGIA1. Maternal and paternal genetic effects were estimated using restricted maximum likelihood methodology applied to the univariate animal model, whereas genetic covariance components were estimated in bivariate models. Estimates of maternal heritability were 0.45 for CL, 0.32 for FE, 0.29 for TE and 0.05 for UE whereas for the paternal effect, heritabilities were very low (<0.06). Genetic correlation between CL, FE and TE variables were very high (>0.89) for the maternal effect. A breeding scheme based on CL selection in response to superovulation could thus improve the number of transferable embryos.

Correlated responses of pre- and postweaning growth and backfat thickness to six generations of selection for ovulation rate or prenatal survival in French Large White pigs

Correlated effects of selection for components of litter size on growth and backfat thickness were estimated using data from 3 pig lines derived from the same base population of Large White. Two lines were selected for 6 generations on either high ovulation rate at puberty (OR) or high prenatal survival corrected for ovulation rate in the first 2 parities (PS). The third line was an unselected control (C). Genetic parameters for individual piglet BW at birth (IWB); at 3 wk of age (IW3W); and at weaning (IWW); ADG from birth to weaning (ADGBW), from weaning to 10 wk of age (ADGPW), and from 25 to 90 kg of BW (ADGT); and age (AGET) and average backfat thickness (ABT) at 90 kg of BW were estimated using REML methodology applied to a multivariate animal model. In addition to fixed effects, the model included the common environment of birth litter, as well as direct and maternal additive genetic effects as random effects. Genetic trends were estimated by computing differences between OR or PS and C lines at each generation using both least squares (LS) and mixed model (MM) methodology. Average genetic trends for direct and maternal effects were computed by regressing line differences on generation number. Estimates of direct and maternal heritabilities were, respectively, 0.10, 0.12, 0.20, 0.24, and 0.41, and 0.17, 0.33, 0.32, 0.41, and 0.21 (SE = 0.03 to 0.04) for IWB, IW3W, IWW, ADGBW, and ADGPW. Genetic correlations between direct and maternal effects were moderately negative for IWB (-0.21 +/- 0.18), but larger for the 4 other traits (-0.59 to -0.74). Maternal effects were nonsignificant and were removed from the final analyses of ADGT, AGET, and ABT. Direct heritability estimates were 0.34, 0.46, and 0.21 (SE = 0.03 to 0.05) for ADGT, AGET, and ABT, respectively. Direct and maternal genetic correlations of OR with performance traits were nonsignificant, with the exception of maternal correlations with IWB (-0.28 +/- 0.13) and ADGPW (0.23 +/- 0.11) and direct correlation with AGET (-0.23 +/- 0.09). Prenatal survival also had low direct but moderate to strong maternal genetic correlations (-0.34 to -0.65) with performance traits. The only significant genetic trends were a negative maternal trend for IBW in the OR line and favorable direct trends for postweaning growth (ADGT and AGET) in both lines. Selection for components of litter size has limited effects on growth and backfat thickness, although it slightly reduces birth weight and improves postweaning growth.

Correlated responses for litter traits to six generations of selection for ovulation rate or prenatal survival in French Large White pigs

Effects of selection for reproductive traits were estimated using data from 3 pig lines derived from the same Large White population base. Two lines were selected for 6 generations on high ovulation rate at puberty (OR line) or high prenatal survival corrected for ovulation rate in the first 2 parities (PS line). The third line was an unselected control line. Genetic parameters for age and BW at puberty (AP and WP); number of piglets born alive, weaned, and nurtured (NBA, NW, and NN, respectively); proportions of stillbirth (PSB) and survival from birth to weaning (PSW); litter and average piglet BW at birth (LWB and AWB), at 21 d (LW21 and AW21), and at weaning (LWW and AWW) were estimated using REML methodology. Heritability estimates were 0.38 +/- 0.03, 0.46 +/- 0.03, 0.16 +/- 0.01, 0.08 +/- 0.01, 0.09 +/- 0.01, 0.04 +/- 0.01, 0.04 +/- 0.02, 0.19 +/- 0.02, 0.10 +/- 0.02, 0.10 +/- 0.02, 0.36 +/- 0.02, 0.27 +/- 0.01, and 0.24 +/- 0.01 for AP, WP, NBA, PSB, NW, NN, PSW, LWB, LW21, LWW, AWB, AW21, and AWW, respectively. The measures of litter size showed strong genetic correlations (r(a) >/= 0.95) and had antagonistic relations with PSB (r(a) = -0.59 to -0.75) and average piglet BW (r(a) = -0.19 to -0.46). They also had strong positive genetic correlations with prenatal survival (r(a) = 0.67 to 0.78) and moderate ones with ovulation rate (r(a) = 0.36 to 0.42). Correlations of litter size with PSW were negative at birth but positive at weaning. The OR and PS lines were negatively related to PSW and average piglet BW. Puberty traits had positive genetic correlations with OR and negative ones with PS. Genetic trends were estimated by computing differences between OR or PS and control lines at each generation using least squares and mixed model methodologies. Average genetic trends were computed by regressing line differences on generation number. Significant (P < 0.05) average genetic trends were obtained in OR and PS lines for AP (respectively, 2.1 +/- 0.9 and 3.2 +/- 1.0 d/generation) and WP (respectively, 2.0 +/- 0.5 and 1.8 +/- 0.5 d/generation) and in the PS line for NBA (0.22 +/- 0.10 piglet/generation). Tendencies (P < 0.10) were also observed for LWB (0.21 +/- 0.12 kg/generation) and AWW (-0.25 +/- 0.14 kg/generation) in the PS line. Selection on components of litter size can be used to improve litter size at birth, but result in undesirable trends for preweaning survival.