|
|
||||||||
Department of Dairy and Animal Science, The Pennsylvania State University, University Park 16802
Corresponding author: A. J. Heinrichs; e-mail: ajh{at}psu.edu.
| ABSTRACT |
|---|
|
|
|---|
Key Words: rumen tissue sampling papillae correlation power analysis
Abbreviation key: A = caudal portion of caudal ventral blind sac, LSN = least significant number, LB = left side caudal dorsal sac, LC = left side cranial dorsal sac, LD = left side cranial ventral sac, LE = left side ventral portion of caudal ventral blind sac, PL = papillae length, PW = papillae width, PC = papillae per square centimeter, RB = right side caudal dorsal sac, RC = right side cranial dorsal sac, RD = right side cranial ventral sac, RE = right side ventral portion of caudal ventral blind sac, RWT = rumen wall thickness
| INTRODUCTION |
|---|
|
|
|---|
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
![]() |
where
ij | = | the correlation between PL, PW, RWT, or PC of the ith and jth areas, samples, or measurements;
| ij | = | the covariance of the ith and jth areas, samples, or measurements for PL, PW, RWT, or PC;
| i, ( j) | = | the standard deviation of the ith (jth) area, sample, or measurement for PL, PW, RWT, or PC ( i j); and
| i, j | = | A, RB, RC, RD, RE, LB, LC, LD, or LE for area, 1 to 5 for sample, or 1 to 4 for measurement (i j).
|
Comparisons.
Comparisons between experiments, areas, samples, and measurements for PL, PW, RWT, and PC were conducted across all calves and experiments using the MIXED procedure of SAS (1999) with a repeated measures statement. A separate model was utilized for PL and PW analysis than for RWT and PC analysis. The model for PL and PW analysis was
![]() |
where
| ytasm | = | an observed value for PL or PW for the mth measurement, taken from the sth sample, collected from the ath area from a calf in the tth experiment;
| µ | = | the overall mean of the population;
| t | = | the fixed effect of the tth experiment where t = 1 to 4;
| ßa | = | the random effect of the ath area where a = A, RB, RC, RD, RE, LB, LC, LD, or LE;
| s | = | the random effect of the sth sample where s = 1 to 5;
| m | = | the random effect of the mth measurement where m = 1 to 4; and
| etasm | = | the error associated with the mth measurement, taken from the sth sample, collected from the ath area from a calf in the tth experiment; ).
|
Calf nested within experiment and sample nested within area were included in the RANDOM statement of the model, and measurement was utilized as the repeated variable. The model for RWT and PC analysis was
![]() |
where
| ytas | = | an observed value for RWT or PC for the sth sample, collected from the ath area from a calf in the tth experiment;
| µ | = | the overall mean of the population;
| t | = | the fixed effect of the tth experiment where t = 1 to 4;
| ßa | = | the random effect of the ath area where a = A, RB, RC, RD, RE, LB, LC, LD, or LE;
| s | = | the random effect of the sth sample where s = 1 to 5; and
| etas | = | the error associated with the sth sample, collected from the ath area from a calf in the tth experiment; ).
|
Calf nested within experiment was included in the RANDOM statement, and sample was utilized as the repeated variable. Treatment influence was accounted for by including treatment in all comparison models. Age variation is accounted for by the inclusion of an experiment effect in the model. Differences between experiments were considered significant at P < 0.10 and between areas, samples, or measurements at P < 0.001. A strict P value was utilized for area, sample, and measurement comparisons to insure that significant differences were also physiologically measurable differences.
Regression analysis.
The REGRESSION procedure of SAS (1999) was utilized to determine relationships between PL, PW, RWT, and PC across all calves and experiments. Linear models were initially fit, subsequently followed by higher order models if applicable.
Power analysis.
To determine the ability of the procedure to distinguish differences between rumen developmental levels, the necessary number of calves per treatment, samples per calf, measurements per calf, and the resultant power of each test, a power analysis was conducted for each described variable using a power macro (version 1.2) obtained from SAS (Latour, 2003). The power macro enabled the description of the desired P value (
), root mean square error (
), and treatment difference or effect size (
) calculated as
= [SS(Hyp)/n]0.5 where
| SS | = | sum of squares for the comparison of interest;
| Hyp | = | hypotheses stating that the effect size between 2 observed measurements for PL, PW, RWT, or PC equals 0 (null hypothesis) or that the effect size is not equal to 0 and is some nonzero number (alternate hypothesis); and
| n | = | number of observations within the comparison of interest.
|
Effect size is determined by the difference between the observed effect size and the null hypothesis effect size. Levels for
were set at 0.01, 0.05, and 0.10. Levels for
and
were not defined by the researchers, but were calculated by the power macro through incorporation of the GLM procedure of SAS (1999) using the following model:
![]() |
where
| ypr | = | an observed value for PL, PW, RWT, or PC from the rth calf, sample, or measurement in the pth rumen development group;
| µ | = | the overall mean of the population;
| p | = | the fixed effect of the pth rumen development group where p = low or high;
| ßr | = | the random effect of the rth calf, sample, or measurement where r = 1 to 42 for calf; 1 to 5 for sample; or 1 to 4 for measurement; and
| epr | = | the error associated with the rth calf, sample, or measurement in the pth rumen development group; ).
|
The power of the test at n, the least significant number (LSN), and the power of the test when n = LSN were provided for each
level using the estimated
and
values. The LSN is defined as the sample size required to produce a significant test for a sample having an
,
, and
values equal to those of the data set used. Average values for PL, PW, RWT, and PC were calculated for each calf. The average values of PL, PW, and RWT were sorted from lowest to highest, with low average measurements considered as indicative of a low level of rumen development and high average measurements representative of a high level of rumen development. Average values for PC were sorted from highest to lowest with high PC values indicative of low rumen development (Anderson et al., 1982; Klein et al., 1987; Zitnan et al., 1999). The calves were then divided into 2 treatment groups, a low or high rumen development group, based on average variable measurements. Consideration was given to simply selecting the 21 lowest and highest averages and their corresponding calves, as this method would possibly decrease the observed variance for each rumen development variable. However, with the desire to make the results of this paper as realistic and applicable as possible, a selection criterion was developed. Division into treatment groups was as follows:
Papillae length was the primary grouping factor due to analyses indicating a high level of correlation between PL throughout the rumen and stronger relationships between PL and the other 3 variables than between PW, RWT, or PC and the other 3 variables. The importance order in step 3 was developed using the same reasoning and statistics, with higher correlation values and stronger relationships for PW, followed by RWT, with PC having the lowest correlation values and weakest relationships. The power analysis step was conducted last to avoid any influence on previous analyses due to the selection criteria used. The power analyses utilized all 42 calves and included observations from all areas. Calves from all 4 experiments were present in both the low and high group accounting for experiment and age differences. All
and
values presented were calculated by the power analyses and are unique to the data set utilized. Each calf, in the initial data set, had 180 observations for PL and PW (labeled 1 to 180) and 45 observations for RWT and PC (labeled 1 to 45) taken from 45 tissue samples (labeled 1 to 45) within the rumen. Average PL, PW, RWT, and PC were calculated for each calf to create the CALF power analysis data set. Average PL, PW, RWT, and PC were calculated for each sample (1 to 45) across calves within the given treatment group to create the SAMPLE power analysis dataset. For creation of the MEASUREMENT dataset, average PL and PW for each measurement (1 to 180) were calculated across calves within the treatment group. Data sets for SAMPLE and MEASUREMENT power analyses essentially represented 2 rumens, one with a high level of rumen development and one with a low level of rumen measurement, and contained 45 and 180 observations per treatment, respectively. Values averaged across calves within treatments were utilized due to a separate power analysis being conducted for each effect of interest (calf, sample, and measurement); therefore each effect analyzed was considered a main effect within their respective power analysis (Dawson and Lagakos, 1993). Calves within treatment and observations within calves were the same for all data sets. Variation attributable to each of the dependent variables of interest (i.e., calves, samples, measurements) was captured by the manner in which averages were calculated for each distinct data set, and the number of observations within each data set.
| RESULTS |
|---|
|
|
|---|
|
Comparisons.
Values for least square means, standard error, and comparisons of PL, PW, RWT, and PC for areas across calves and experiments are presented in Table 2
. Significant differences (P < 0.001) were observed between some areas for PL, PW, RWT, and PC. Occurrences of difference between areas may be primarily explained by intrarumen variation. In addition, standard errors for PL and PW made up a larger portion of the mean than standard errors for RWT and PC, indicating greater variation within the former rumen development parameters. Similar values and minimal significant differences between areas corresponding to the caudal portion of the caudal ventral blind sac and the caudal dorsal sac (A and B) were observed, indicating that samples from either area could represent the caudal rumen. Numerous differences between right and left rumen areas were observed for PL, with only one difference observed for PW and RWT, and no right and left rumen area differences seen for PC. More differences between right and left rumen areas for PL may be explained by higher variation within PL when compared with PW, RWT, and PC.
|
Similar correlation results were observed between measurements taken across samples, areas, calves, and experiments. Measurement correlation values for variables PL and PW were all significantly different (P < 0.001) from 0, and ranged from 0.90 to 0.92 for PL and from 0.72 to 0.76 for PW, indicating a strong correlation between measurements for PL and PW.
Comparisons.
Values for least square means of PL (average 1.11 mm), PW (average 0.70 mm), RWT (average 1.38 mm), and PC (average 75.86) were similar for different samples across areas, calves, and experiments with no significant differences between samples. Standard errors as a percentage of means were highest for variable PL (0.09) followed by PW (0.05), PC (3.30), then RWT (0.04), as observed in the area analysis.
No significant differences were observed between least square means of PL (average 1.11 mm) and PW (average 0.70 mm) from measurements taken across samples, areas, calves, and experiments. Standard errors as a percentage of means were once again highest for variable PL (0.09) followed by PW (0.04), as seen in the area and sample analyses.
Rumen Variable Relationships
All calves and experiments were included in the regression analyses. Relationships between most rumen variable pairs were minimal, except for PL and PW. There was a significant linear and quadratic relationship between PL and PW. The linear regression was:
![]() |
The quadratic relationships between PL and PW were:
![]() |
with
![]() |
with
![]() |
Power Analysis
Calves, samples, and measurements required.
Table 3
presents the results of power analyses conducted to determine the number of calves, samples, and measurements necessary to find significant differences between treatments with an
of 0.01, 0.05, or 0.10 for PL, PW, RWT, and PC. Results of the power analysis for PL, PW, and RWT indicate the ability of the procedure to detect treatment differences, at an acceptably high power, with a data set containing 21 calves per treatment, 45 samples, and 180 measurements per calf. Results of the power analysis for PC indicate an inability of the procedure to detect significant differences between treatments with the number of samples utilized.
|
,
, and LSN values are unique to this data set, LSN values may not be universally applicable, but do provide a starting point for research design and planning. With the observed
,
, and
values, the results from Table 3
of 0.01, 0.05, and 0.10, respectively. With the same format, the results suggest 6, 4, or 3 calves per treatment, 3, 3, or 2 samples per calf, and 5, 3, or 3 measurements per calf to detect treatment difference for PW. Results for RWT suggest 18, 11, or 8 calves per treatment and 13, 8, or 6 samples/measurements per calf to detect treatment differences with an
of 0.01, 0.05, or 0.10, respectively. Values of LSN for RWT are higher than for PL and PW, possibly due to low
values (observed difference between treatments) for this variable. Due to low
values and low relationship between
and
, values of LSN for PC are high, suggesting that PC may not be a feasible variable in rumen development research. | DISCUSSION |
|---|
|
|
|---|
Rumen Variable Relationships
The relationship among and explained variation of the rumen variables is greatest for PL, followed by PW, then RWT, with PC having little to no relationship with the other variables. In addition, the regression of PL, PW, and RWT on PC indicated a slight negative relationship. This suggests that as PC increases, values for PL, PW, and RWT decrease, a trait reported in previous research (Anderson et al., 1982; Klein et al., 1987; Zitnan et al., 1999). This relationship is not surprising as increased values for PL, PW, and RWT are indicative of increased rumen and papillary growth and development. An increase in the size of papillae in a fixed 1-cm2 section of the rumen wall must result in a decrease in the number of papillae in that finite area. In addition, an increased rumen volume, possibly represented by an increased RWT, enlarges the area covered by 1-cm2 at birth over a larger area at a later age, thereby decreasing the number of papillae in the finite section once again.
Comparison and Power Analyses
Papillae length.
Values for PL indicate a greater occurrence of treatment and area differences and high variability within PL. The high variation and differences may explain the increased power of the procedure to detect treatment differences in PL, as presented in Table 3
. Therefore, PL may be the most important variable for rumen development research, and may represent the greatest influence of treatment on rumen development. However, high PL variability may also increase the number of areas that require sampling to sufficiently represent development of the entire rumen. Previous research has detected treatment differences for PL when weaning date or the dry portion of the ration was chemically or physically altered (Nocek et al., 1984; Greenwood et al., 1997; Zitnan et al., 1999). These previous findings are expected as rumen development is greatly influenced by age and the presence of butyrate from microbial and protozoal degradation of readily fermentable carbohydrate sources in the rumen (Brownlee, 1956; Warner et al., 1956; Klein et al., 1987). Other research has not detected treatment differences in PL, but has detected age differences (Anderson et al., 1982; Klein et al., 1987).
Papillae width.
Values for PW indicate a lower occurrence of treatment and area differences and decreased variability within PW than PL, but higher than RWT and PC. These values possibly explain why the power of PW falls between the power for PL and RWT (Table 3
). However, the power of the procedure to detect treatment differences for PW is also very high. Therefore, these results suggest that PW is likely a secondarily important rumen variable for rumen development research. Detection of differences in PW has been limited in previous research, and weaning age appears to have a greater effect than chemical or physical alteration of the dry ration (Klein et al., 1987; Greenwood et al., 1997; Zitnan et al., 1999).
Rumen wall thickness.
It appears that RWT is tertiary in importance as a rumen development variable. Statistical treatment differences can be obtained with the procedure, but at a decreased ability than indicated for PL and PW, possibly due to insufficient observations per area or the inability of treatment differences to influence this variable to the same extent as PL and PW. Previous research has not detected significant differences between treatments for rumen wall and/or epithelial thickness when measurements were obtained in areas similar to those described for this procedure (Anderson et al., 1982; Greenwood et al., 1997). In addition, it has been suggested that an increase in rumen muscularization may occur independently of rumen epithelial growth (Brownlee, 1956; Harrison et al., 1960). However, as no attempts were made to separate the rumen muscle from the rumen epithelium in this procedure, the results in these analyses may be confounded.
Papillae per square centimeter.
Results for PC indicate a low power, or possible inability, of the procedure to detect treatment differences in PC, suggesting inapplicability of PC as a rumen development variable. The inability of the procedure to detect treatment differences is likely a result of high
values and low
values. It is also possible that the described procedure is incapable of detecting treatment differences for PC due to sampling error. In addition, the possibility of high genetic control over PC, subsequently limiting environmental influence, should also be considered. Age of the calf, age at weaning, and subsequent length of time that concentrates make up an appreciable portion of the daily ration are likely the environmental factors having the greatest influence on detectable differences in PC and therefore have the highest possibility of confounding results (Klein et al., 1987; Zitnan et al., 1998; Zitnan et al., 1999). Furthermore, age of the calf and related rumen volume has a tremendous effect on the finite area of the rumen sampled, as described in the rumen variable relationships section. However, some researchers have reported significant differences between treatments for this variable (Anderson et al., 1982; Nocek et al., 1984).
Calf Numbers, Samples, and Measurements
It is apparent from the results that the areas sampled, the number of tissue samples taken, and the number of measurements per area can be reduced from that originally suggested by this procedure. Papillae length is the only variable appearing to require sampling from the right and left side of the rumen. However, a high power of the test for PL may overcome the need for right and left side representation. In addition, the occurrence of only one right and left side difference for PW and RWT, and no differences for PC, coupled with relatively high correlations between right and left rumen areas, indicates a limited need for sampling both the right and left side of the rumen. Therefore, the rumen dissection procedure reported by McGavin and Morrill (1976a) may be effective for dissecting and sampling the reticulo-rumen. However, the visible guides utilized for dissection may not be readily apparent in the young calf, increasing the difficulty of the McGavin and Morrill (1976a) dissecting procedure. Due to the results indicated, it is recommended that tissue samples be obtained from the cranial and caudal sacs of the ventral and dorsal rumen (area n = 4), 3 random tissue samples collected from each area (sample n = 12), and 2 measurements per sample for PL, PW, and RWT recorded (measurement n = 24). For purposes of representing physiological growth throughout the entire rumen, multiple sampling sites are suggested. However, for statistical purposes, samples collected from a specific rumen area may also be acceptable, provided that rumen sampling site is identical across calves and treatments. Twenty-four measurements per calf reduces the originally suggested sample numbers by 87% for PL and PW and by 47% for RWT, and should result in statistically acceptable power for PL and PW, but may be limited in statistical power for RWT. Therefore, more RWT measurements per calf may be desirable but will require additional tissue samples. However, the possible inapplicability of RWT may not warrant additional time spent on this variable, and PL or PW may sufficiently address rumen epithelial changes. It is not suggested to record PC measurements due to the high numbers required to find statistical differences for this variable. The required calves per treatment to find differences range from 3 to 18; however, these values result in detectable treatment difference with a possibly unacceptable power. In addition, it is not suggested to use less than 3 calves per treatment to avoid outlier influence and maintain degrees of freedom within treatment. It appears that 3 calves per treatment should result in an analysis of sufficient power for PL and PW, valid indicators of rumen development, and minimally limited power for RWT.
| CONCLUSION |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication May 29, 2003. Accepted for publication October 23, 2003.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. A. Khan, H. J. Lee, W. S. Lee, H. S. Kim, S. B. Kim, S. B. Park, K. S. Baek, J. K. Ha, and Y. J. Choi Starch Source Evaluation in Calf Starter: II. Ruminal Parameters, Rumen Development, Nutrient Digestibilities, and Nitrogen Utilization in Holstein Calves J Dairy Sci, March 1, 2008; 91(3): 1140 - 1149. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. A. Khan, H. J. Lee, W. S. Lee, H. S. Kim, S. B. Kim, K. S. Ki, S. J. Park, J. K. Ha, and Y. J. Choi Starch Source Evaluation in Calf Starter: I. Feed Consumption, Body Weight Gain, Structural Growth, and Blood Metabolites in Holstein Calves J Dairy Sci, November 1, 2007; 90(11): 5259 - 5268. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. A. Khan, H. J. Lee, W. S. Lee, H. S. Kim, K. S. Ki, T. Y. Hur, G. H. Suh, S. J. Kang, and Y. J. Choi Structural Growth, Rumen Development, and Metabolic and Immune Responses of Holstein Male Calves Fed Milk Through Step-Down and Conventional Methods J Dairy Sci, July 1, 2007; 90(7): 3376 - 3387. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. E. Lesmeister and A. J. Heinrichs Effects of Adding Extra Molasses to a Texturized Calf Starter on Rumen Development, Growth Characteristics, and Blood Parameters in Neonatal Dairy Calves J Dairy Sci, January 1, 2005; 88(1): 411 - 418. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. E. Lesmeister and A. J. Heinrichs Effects of Corn Processing on Growth Characteristics, Rumen Development, and Rumen Parameters in Neonatal Dairy Calves J Dairy Sci, October 1, 2004; 87(10): 3439 - 3450. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |