|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
ANIMAL GENETICS |
USDA, ARS, US Meat Animal Research Center, Clay Center, NE 68933-0166
| Abstract |
|---|
|
|
|---|
Key Words: fetal survival placental growth pig uterine capacity
| INTRODUCTION |
|---|
|
|
|---|
Despite clear evidence of a genetic basis for litter size (Johnson et al., 1999
), chromosomal regions containing QTL for the component traits of litter size have been difficult to detect (Rohrer et al., 1999
; Cassady et al., 2001
), indicating the genetic architecture is composed of many loci with small effects, which influence complex interactions. Gene expression data can provide complementary information to genomic scans for detection of the underlying genetic variation that influence messenger RNA profiles (Caetano et al., 2004
; Blomberg et al., 2005
). However, a better understanding of the temporal and spatial expression of line differences for embryo or fetal survival is needed to facilitate future gene discovery experiments within the OR and UC lines.
Our specific experimental objective was to characterize timing and magnitude of line differences in prenatal mortality in UHO gilts that were serially slaughtered at 25, 45, 65, 85, and 105 d of gestation. Additionally, the experiment generated data to characterize these selection lines for the proportionality of growth of fetal and placental components of developing conceptuses.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Experimental Design and Traits
A 4-breed composite with equal contributions from Chester White, Landrace, Large White, and Yorkshire breeds was formed to estimate breed, heterosis, and recombination effects (Cassady et al., 2002
). From a common base generation of this composite produced in 1986, selection was initiated during 1988 within 2 replicated seasons for ovulation rate estimated from laparoscopic examination of gilts in the OR line and increased uterine capacity of UHO surgically-altered gilts in the UC line, and an unselected control (CO) line was maintained. The selection experiment proceeded for 11 generations, and then all 3 lines were evaluated for responses in the component traits (Leymaster and Christenson, 2000
). Relative to the unselected control, selection for ovulation rate increased ovulation rate by 3.2 ova, decreased uterine capacity by 0.97 pigs, decreased prenatal survival by 10.3%, and increased litter size by 0.30 pigs. Selection for uterine capacity increased ovulation rate by 0.13 ova, increased uterine capacity by 2.15 pigs, increased prenatal survival by 3.5%, and increased litter size by 0.62 pigs. Subsequent to this terminal evaluation, each line has been maintained under no intentional selection pressure, in 2 replicated seasons farrowing annually in March or September. Genetic diversity within each selection line and season was maintained by an intended target of 10 boars and 40 gilts farrowed.
Unilateral-hysterectomy-ovariectomy surgery was performed at approximately 160 d of age on 50 gilts of each line in each of 4 farrowing seasons, for a total of 593 gilts. Removal of the left or right uterine horn and ovary was done on an alternating basis. Justification and details of UHO as a procedure to estimate uterine capacity have been described previously (Christenson et al., 1987
). All UHO-treated gilts that expressed estrus were single-sire mated at approximately 280 d of age to boars within the same selection line during a 21-d breeding season, avoiding paternal half-sib matings. A total of 502 gilts were mated and assigned to 1 of 5 gestational time points (d 25, 45, 65, 85, or 105) for slaughter. Gilts were assigned to a specific slaughter date to match one of the intended time points as closely as possible (± 2 d) and also to distribute paternal half-sibs across as many of the 5 time points as possible.
Gilts were slaughtered at the US Meat Animal Research Center abattoir, where electrical stunning rendered them unconscious and they were exsanguinated. The reproductive tracts were removed and immediately evaluated. Tracts that displayed evidence of incomplete healing from UHO were discarded, whether or not the gilt was pregnant. The number of observations is presented by line and gestational stage (Table 1
). Data collected from 402 gilts were analyzed. The ovary was excised and weighed before dissection and recorded for number of corpora lutea, representing ovulation rate. The cervix and broad ligament were closely trimmed, and gravid uterine weight was recorded. Length from the cervical end of the uterine body to the oviductal end of the uterine horn was recorded. Each tract was then opened along the antimesometrial border. Each fetus was removed by detachment from the umbilical cord, and status (alive, dead, or mummy) was recorded. Observation of a heartbeat was required for live status to be recorded. Weight of each fetus was also recorded. Placental tissue matching each fetus was separated from the endometrium and weighed. At 85 and 105 d of gestation, the tracts were allowed to sit at room temperature for 1 to 2 h after removal of the fetuses and before separation of the placental tissue. Traits calculated on a per gilt basis from individual fetal data included average and litter fetal and placental weights from live fetuses, total number of live fetuses, and total number of dead fetuses including mummified fetuses.
|
Sampling issues when evaluating reproductive traits like litter size can generate least squares means that do not make sense from a biological point of view. For example, raw means and subsequently the least squares means estimated for number of live fetuses at d 105 were slightly greater than the estimates obtained at d 85 for each selection line. Obviously, litter size cannot increase from d 85 to 105, and this discrepancy was due to the random sample of gilts selected for each stage. A second statistical analysis was utilized to assess the dynamic relationship between ovulation rate and litter size (and thus prenatal survival) as a nonlinear decay function that can only decrease with time. To quantify prenatal survival curve characteristics, data for number of live fetuses were fitted to a nonlinear equation with 3 parameters:
![]() | [1] |
where f (d) = the number of live fetuses at a specific day of gestation, A = the asymptotic estimate of number live fetuses, B = the estimate of ovum wastage, k = the shape parameter of the curve, and d = the day of gestation. The nonlinear solution procedure (PROC NLIN) supported by SAS was used to provide estimates of the parameters using the Gauss-Newton method. All gilts provided data to estimate the number of potentially live fetuses at d = 0 (based on the ovulation rate) and also 1 time point in gestation (25, 45, 65, 85, 105 d); thus, parameters were estimated from a data set with 2 observations per gilt. Parameters A, B, and k were estimated for each line and also for data pooled across 2 or 3 lines. An F-ratio was calculated to test whether the estimation of the parameters specific to each line significantly improved the fit of the data relative to estimation of the parameters from pooled data sets. The test statistic calculated was as follows:
![]() |
where RSS represents the residual sums of squares, Rdf represents the residual degrees of freedom, and the subscripts P and L represent the pooled and contributing line-specific models, respectively. Comparison of a single equation pooling all 3 lines relative to 3 line-specific equations had 6 df in the numerator, whereas both 2-way comparisons of lines, each selection line (OR or UC) pooled with the CO, compared with contributing line-specific equations were 3-df tests. Large values of each F-statistic generated evidence that the pooled, single model was inappropriate, and the relationships among ovulation rate and number of live fetuses differed between lines contributing to the pooled function.
A similar statistical approach was used to characterize and compare lines for average and total litter fetal and placental growth of live fetuses during gestation. A general exponential form of fetal growth was evaluated using the equation:
![]() | [2] |
where f(d) = the weight of the live fetus or placental tissue at a specific day (d) of gestation. Parameters estimated included an intercept W and 2 terms to describe the instantaneous growth rate X of the tissue (fetus or placenta) that decreased by the value Y each day of gestation. Parameters were again estimated for line-specific and combinations of pooled functions. An F-ratio (as described previously) was calculated to test whether estimation of the parameters specific to each line significantly improved fit of the data relative to estimation of parameters from pooled data sets. Large values of each F-statistic generated evidence that a pooled, single model was inappropriate and that the average or total fetal and placental growth rates differed between lines. As stated before for survival curves, comparison of a single equation for all lines relative to 3 line-specific equations was a 6-df test, whereas the 2 comparisons of the selection lines (OR and UC) each pooled with CO were 3-df tests.
Statistical Analyses of Individual Conceptus Growth Traits
Allometric growth rate of fetal and placental components of conceptus tissue were evaluated using the equations and concepts for part-whole relationships developed by Huxley (1932)
. The allometric equation is represented in the following form:
![]() |
where Y = a measure of the component, X = a measure of the whole, a = a constant (value of Y when X = 1), and b = the relative growth coefficient. This equation can be linearized to its logarithmic form, ln(Y) = ln(a) + b* ln(X), to allow easily computed and interpreted solutions for the b parameter.
Conceptus tissue mass was calculated as the sum of fetal and placental weights (n = 3,371 live conceptuses). Fixed effects of farrowing season (1, 2, 3, and 4), line (CO, OR, and UC), gestational age (25, 45, 65, 85, and 105 d), and 2-way interactions were fitted. Sire (n = 119) of the fetus nested within farrowing season and line and also gilt (n = 402) nested within sire of the fetus were fitted as random effects. The Kenward-Rogers option was used to approximate the denominator degrees of freedom associated with the random effects. Growth impetuses (relative growth rates) of the fetal and placental components of the conceptus were determined. Natural log-transformed values of the fetus or the placental weight were regressed on the natural log of the entire conceptus weight. A slope of 1.0 implies that the component tissue was growing at the same rate as the entire conceptus. A slope of less than 1.0 indicates the component tissue was growing at a slower rate than the entire conceptus, whereas a slope of greater than 1.0 indicates a faster rate of growth by the component tissue. Heterogeneity of these slopes throughout gestation and between selection lines was evaluated by fitting the subclass regressions within line and gestational age as well as within the interaction of line x gestational age.
| RESULTS |
|---|
|
|
|---|
Effects of farrowing season, farrowing season x gestational stage, and farrowing season x line were not tabulated. Ovulation rate, and subsequently ovarian weight, was influenced by environmental effects associated with these farrowing seasons. Gilts born in September and mated during June had greater ovulation rates (1 ovum) than gilts born in March and mated during December. This seasonal influence on ovulation rate has been consistently observed within these lines (Ford et al., 2004
).
Characterization of Traits
Least-squares means and residual SD are presented by stage of gestation after accounting for selection line and farrowing season fixed effects (Table 2
) to give an indication of the experimental variation. Inspection of the residual SD for average and litter fetal and placental weights provided justification for adjusting these raw data by log transformation to account for the scaling effect of the rapidly increasing means with increased stage of gestation.
|
|
Line Effects on Fetal Survival
Estimates of parameters used to describe the number of live fetuses as a function of gestational age are presented in Table 4
. Fetal losses occurred throughout gestation as measured by the nonlinear decay function described earlier (equation [1]). The asymptotic estimate of number of live fetuses per uterine horn, when ovulation rate is not limiting, is represented by parameter A: A is an estimate of uterine capacity with the caveat that gestation time is limited to 114 d. The B parameter is the estimated difference between the asymptote and the value of the function at d = 0 (ovulation rate), and thus in biological terms represents total ovum wastage. Estimates of the k parameter dictate the shape of the decay pattern with gestational age. Three line-specific functions fit the data better than a single pooled function (Table 4
; Figure 1
). The F-tests also rejected (P < 0.01) the null hypotheses that OR and CO were equivalent, as well as UC and CO. Analysis of residual values did not reveal any bias associated with the line-specific equations. In general, the curve for OR differed from the other 2 lines by having the greatest value of B (ovum wastage), and the lowest value of A. The curve for UC exhibited the lowest value of B and the greatest value of A. Numbers of live fetuses at d 114 were predicted for each line with values of 7.0, 5.9, and 7.8 live pigs per uterine horn for CO, OR, and UC, respectively. When expressed as deviations from CO, these results imply that single trait selection for uterine capacity increased the intended trait by 0.8 pigs per uterine horn, and single trait selection for ovulation rate decreased uterine capacity by 1.1 pigs per uterine horn. These results would support a negative genetic correlation between ovulation rate and uterine capacity.
|
|
|
|
Analysis of average placental weight (Figure 3
) rejected (P < 0.01) the hypothesis that all lines were similar. The F-tests indicated that UC and CO were similar to each other and that functions for CO and OR were different (P < 0.01). This difference was evident by d 45, and the magnitude increased until d 85 when all 3 lines reached a plateau phase of placental growth.
|
Relative Growth Rates of Fetal and Placental Tissues
Results from tests of homogeneity of growth coefficients among selection lines and gestational stages relating fetal and placental growth relative to the entire conceptus are presented in Table 6
. Differences (P < 0.05) were detected among combinations of lines and gestational stages for relative growth of the placenta and approached statistical significance (P < 0.10) for relative growth of the fetus. Based on these results, estimates of relative growth coefficients for the interaction of selection line with gestational stage are presented in Table 7
.
|
|
| DISCUSSION |
|---|
|
|
|---|
Lines of pigs that differed greatly in conceptus survival were created (Leymaster and Christenson, 2000
) by selection for ovulation rate or uterine capacity, the number of fetuses carried to term in UHO-treated gilts. The intent was not to select for litter size directly, but rather to select independently for the component traits and provide data to validate a litter size model (Bennett and Leymaster, 1989
). The OR line had no intended selection pressure for uterine capacity, and similarly no selection pressure for ovulation rate was intended in the UC line. Early embryo survival before d 25 of gestation (Youngs et al., 1994
) and uterine capacity effects from d 30 to parturition (Leymaster et al., 1986
; Johnson et al., 1999
) have both been shown to contribute significantly to prenatal losses. It was not known when increased uterine capacity provided a more favorable uterine environment to increase survival rates. Selection in the UC line was based on number born at farrowing, allowing for expression of accumulated effects throughout gestation. The current experiment was specifically designed to address the timing of prenatal mortality in these unique lines. Line differences in timing of prenatal mortality have implications on approaches to investigate genetic variation underlying the loci affected by selection. Physiological and metabolic pathway targets considered for evaluation may differ depending on timing of mortality.
Differences between lines in ovulation rate were similar to those estimated in the terminal evaluation of the selection phase (Leymaster and Christenson, 2000
), with OR exceeding the other 2 lines by 3 to 4 ova. Similar rates of embryonic loss to d 25 for all 3 lines indicate that losses due to reduced ovum quality or maturation, fertilization failures, and genetic abnormalities were not affected by OR selection. Also, direct selection for uterine capacity did not result in a significant increase of live embryos before d 25. Numbers of dead embryos observed at this early stage accounted for less than 30% of the ova lost to that point. If there is a genetic basis for potential embryo viability as defined in Bennett and Leymaster (1989)
, associated with selection for ovulation rate or uterine capacity, it was not evident in this study. Least-squares means for ovulation rate and number of live fetuses indicated 67, 74, and 72% survival to d 25 for CO, OR, and UC lines, respectively. These estimates of embryonic viability at this time point were slightly lower than the modeled value of 0.82 estimated from a summary of literature at that time (Bennett and Leymaster, 1990
).
Dramatic changes between lines in fetal losses between d 25 and 45 were evident (Figure 1
) and provided the majority of the differences between lines. Gilts from the OR line had less endometrial space available per embryo at implantation and this severely crowded environment resulted in lower average placental weights and greater fetal losses during the critical period for placental development. Compared with the CO line, the UC line had a similar initial density of embryos per unit of endometrial space, yet provided a superior uterine environment to increase number of live fetuses and maintain similar average fetal and placental weights. A greater number of UC embryos developed a sufficient placenta and survived to and beyond d 45, producing greater litter placental and fetal weights. Despite the smallest number of live fetuses, the OR line had the lowest average fetal and placental weights. Apparently, the OR line was not able to capitalize on uterine space associated with early-stage embryos that died, suggesting that the space allotted to each placenta is fixed during early gestation. This conclusion is consistent with a report by Lamberson and Eckardt (1996)
that also indicated surviving fetuses do not make use of uterine space vacated by nonsurviving littermates.
The period of rapid embryo elongation and communication among potential littermates to allocate available uterine space and the subsequent initiation of attachment occurs by d 13, with complete attachment over the entire surface of the conceptus by d 26 (Foxcroft et al., 2000
). A rapid increase of allantoic fluid occurs between d 20 and 30 that is associated with expansion of the chorioallantoic membranes and direct intimate contact between the placenta and endometrial surface (Knight et al., 1977
). Development of the placental membrane length during this period is critical to the overall allocation of endometrial space occupied by an individual embryo. Placental length increases rapidly from d 20 to 30 and then continues to increase at a slower rate to d 60 when it reaches a plateau, with placental weight reaching a plateau approximately 10 d later (Pomeroy, 1960
; Knight et al., 1977
). This rapid placental growth period requires a high level of coordination for appropriate amino acid metabolism to optimize the DNA and protein synthesis machinery to support this level of cell proliferation. During this specific time of gestation (from d 20 to 40), placental branched chain amino acid transport, transamination, and activities of related enzymes increase markedly (Self et al., 2004
). It has been demonstrated that d 40 of gestation is also coincident with maximum placental glutamine synthesis and when glutamine is most abundant in fetal allantoic fluid (Self et al., 2004
). Uptake of glutamine by the embryo becomes critical in its role as a substrate for fetal synthesis of arginine because uterine uptake of arginine is not sufficient to support fetal growth during late pregnancy (Wu et al., 1999
). Pathways and enzymes associated with amino acid metabolism changes in this critical time period would make reasonable candidate genes for intensive evaluation of altered allele frequencies and functional analysis in these selected lines.
Development of fetal erythropoiesis also coincides with the critical interval from d 20 to 40. It has been suggested that intrauterine crowding results in impaired fetal erythropoiesis (Pearson et al., 1998
) and breed differences exist in the rate of blood cell development (Vallet et al., 2003
). Allelic variation at the erythropoietin receptor gene (EPOR) was discovered in 2 populations and the fetal genotype was associated with a functional increase in the expression of the EPOR gene and improved fetal survival (Vallet et al., 2005
; Vallet and Freking, 2006
). Although this functional polymorphism was encountered within the 3 selection lines from the current experiment, the frequency of the favorable allele was very low (less than 5%) and only marginally greater in the UC line. In spite of a rather strong indication of this single locus affecting embryo survival and litter size (Vallet et al., 2005
), the frequency of this polymorphism was not increased by direct selection for uterine capacity. This evidence is consistent with the concept that genetic architecture of uterine capacity is polygenic with many loci each contributing small effects on the phenotype.
After the critical d 25 to 45 interval, line differences in average placental weights became evident (decreased OR placental weights at d 45) and preceded differences in fetal weights between d 65 and 85. This result is consistent with concepts of maximum placental growth occurring earlier than fetal growth (Pomeroy, 1960
) and that restriction of placental growth leads to a subsequent restriction of fetal growth. At what point restriction of placental growth becomes a critical factor is not clear. Conceptus estrogen production may play a critical role in promoting growth of the conceptus trophectoderm by stimulating endometrial production of growth factors (Wilson and Ford, 2000
; Jaeger et al., 2001
; Blomberg et al., 2006
). Variation in mitotic activity of trophectoderm between and among conceptuses may play a critical role in determining ultimate placental development and function among embryos competing for limited maternal resources. Changes in embryonic estrogen secretion or production in these lines have been investigated with mRNA from d 11 to 12, and preliminary reports indicate decreased levels of steroidogenic transcripts by the OR line (Miles et al., 2006
). However, it is not yet clear whether this is a cause of reduced elongation in the OR line or a consequence of increased intrauterine crowding experienced by these blastocysts. It seems reasonable that greater synchrony of early embryo development can lead to more uniform placental development and fewer embryos at risk for loss. However, typical embryonic growth can also occur up to approximately d 30 in the pig even in the absence of a functional placenta, as in the case of gynogenetic embryos (Surani et al., 1990
; J. Piedrahita, North Carolina State University, Raleigh, NC, personal communication). Vallet et al. (2002)
reported that small conceptus diameter on d 11 was not correlated with uterine capacity or with placental weights during later stages of pregnancy. Fetal growth, in terms of mass, typically begins to exceed placental growth by d 70 (Pomeroy, 1960
). This is supported by our analysis of data for relative growth coefficients (Table 7
). Relative growth coefficients were greater for placental tissue compared with fetal tissue to d 65 and then decreased to a similar rate as that of fetal tissues near term. Relative growth coefficients for the UC line (Table 7
) indicated a later maturing placental tissue, and it appears that the steepest portion of the growth curve (Figure 3
) for placenta coincides with greatest relative growth coefficients.
Wilson et al. (1999)
hypothesized that the increased ratio of fetal mass to placental mass could be used as a criterion to increase litter size based on observations of greater ratios in Meishan pigs compared with occidental breeds and a single generation of divergent selection yielding differences in litter size. This concept was further examined in a line of pigs selected for components of litter size, and results indicated that uterine capacity increases in the selected line were achieved through a different mechanism than altering the ratio of fetal weight to placental weight, as the control line exceeded the selected line for this ratio (Mesa et al., 2003
). Similar to Mesa et al. (2003)
, changes in uterine capacity from the current experiment were not achieved through altered placental efficiency measured as fetal weight relative to placental weight. Functions for number of live fetuses were significantly different between UC and CO (Table 4
), but UC and CO differences were not detected for either average fetal weight (Table 5
) or average placental weight (Table 6
). It is likely that placental function rather than weight is most critical to support fetal development. Divergent selection for 4 generations for this ratio generated little evidence to suggest that selection for improved placental efficiency would increase litter size (Mesa et al., 2005
). All of these studies are based on the concept that placental and fetal mass components at term are predictive of litter size, when in reality it is likely that the critical time period in terms of placental support of larger litters occurs much earlier in gestation. In the current experiment, the UC line had greater relative growth coefficients for the placenta at d 65 and 105, and litter size differences were established by d 45.
Anecdotal observations while separating placental and uterine tissues suggested that placental tissue from the UC line was physically stronger and could withstand greater pressure before tearing. Perhaps altered collagen or other structural proteins were the result of direct selection for uterine capacity, improving the function of individual placenta to support fetal growth under restricted conditions. Understanding improved placental function is critically important because ovulation rate has increased within swine industry populations and a concomitant increase in uterine capacity must also occur to maximize the opportunities in litter size. Even if the value for litter size has approached the industry optimum, postnatal survival is a problem. The importance of prenatal programming to postnatal performance has been documented (Foxcroft et al., 2006
). The current experiment is the largest one to report values for number of live fetuses throughout the gestation interval and has shed light on the temporal expression of line differences for fetal survival. Results indicated the importance of placental tissues in expression of these differences. Our intention is to pursue placental transcriptome differences that resulted from selection for components of litter size and resulted in line differences for fetal survival.
Increased ovulation rates have resulted in less than optimal piglet survival rates during prenatal and early postnatal periods and have hindered the desired outcome of increased numbers of pigs weaned. Understanding the biological basis for improved uterine capacity as a consequence of genetic selection would allow exploitation of greater levels of reproductive potential. Unique fetal survival curves resulted from selection of pigs for ovulation rate or uterine capacity. Embryonic losses from ovulation to d 25 were proportional to ova shed in all 3 lines. Direct selection for uterine capacity has increased survival of fetuses to term, primarily by increasing survival between d 25 and 45 of gestation without altering fetal or placental weights. Losses continued to occur from d 45 to 105, but line differences in number of fetuses at d 45 were essentially maintained to d 105. The critical period of time identified coincides with rapid development of the placenta and focuses attention on metabolic pathways that can create and develop a placental environment capable of maintaining more fetuses of a similar birth weight to term.
| Footnotes |
|---|
2 Corresponding author: Brad.Freking{at}ars.usda.gov
Received for publication November 21, 2006. Accepted for publication May 8, 2007.
| LITERATURE CITED |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. F Trott, K. C Horigan, J. M Gloviczki, K. M Costa, B. A Freking, C. Farmer, K. Hayashi, T. Spencer, J. E Morabito, and R. C Hovey Tissue-specific regulation of porcine prolactin receptor expression by estrogen, progesterone, and prolactin J. Endocrinol., July 1, 2009; 202(1): 153 - 166. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. R. Miles, B. A. Freking, L. A. Blomberg, J. L. Vallet, and K. A. Zuelke Conceptus development during blastocyst elongation in lines of pigs selected for increased uterine capacity or ovulation rate J Anim Sci, September 1, 2008; 86(9): 2126 - 2134. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. L. Vallet and B. A. Freking Differences in placental structure during gestation associated with large and small pig fetuses J Anim Sci, December 1, 2007; 85(12): 3267 - 3275. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |