Postnatal development of collagen structure in ovine articular cartilage
© van Turnhout et al; licensee BioMed Central Ltd. 2010
Received: 29 January 2010
Accepted: 7 June 2010
Published: 7 June 2010
Articular cartilage (AC) is the layer of tissue that covers the articulating ends of the bones in diarthrodial joints. Across species, adult AC shows an arcade-like structure with collagen predominantly perpendicular to the subchondral bone near the bone, and collagen predominantly parallel to the articular surface near the articular surface. Recent studies into collagen fibre orientation in stillborn and juvenile animals showed that this structure is absent at birth. Since the collagen structure is an important factor for AC mechanics, the absence of the adult Benninghoff structure has implications for perinatal AC mechanobiology. The current objective is to quantify the dynamics of collagen network development in a model animal from birth to maturity. We further aim to show the presence or absence of zonal differentiation at birth, and to assess differences in collagen network development between different anatomical sites of a single joint surface. We use quantitative polarised light microscopy to investigate properties of the collagen network and we use the sheep (Ovis aries) as our model animal.
Predominant collagen orientation is parallel to the articular surface throughout the tissue depth for perinatal cartilage. This remodels to the Benninghoff structure before the sheep reach sexual maturity. Remodelling of predominant collagen orientation starts at a depth just below the future transitional zone. Tissue retardance shows a minimum near the articular surface at all ages, which indicates the presence of zonal differentiation at all ages. The absolute position of this minimum does change between birth and maturity. Between different anatomical sites, we find differences in the dynamics of collagen remodelling, but no differences in adult collagen structure.
The collagen network in articular cartilage remodels between birth and sexual maturity from a network with predominant orientation parallel to the articular surface to a Benninghoff network. The retardance minimum near, but not at, the articular surface at all ages shows that a zonal differentiation is already present in the perinatal animals. In these animals, the zonal differentiation can not be correlated to the collagen network orientation. We find no difference in adult collagen structure in the nearly congruent metacarpophalangeal joint, but we do find differences in the dynamics of collagen network remodelling.
The articulating ends of bones in the diarthrodial joints are covered with a thin layer of articular cartilage (AC). During early development AC functions as a surface growth plate for the underlying bone [1–3]. In adult life, AC functions as a load bearing surface that transmits loads and provides a low friction environment. Between birth and skeletal maturity, AC has to accommodate both developmental and load bearing demands.
AC has been shown to develop from a fairly homogeneous tissue to a tissue with site specific composition [1, 4, 5] and site specific mechanical properties . It is accepted that this development is affected by mechanical loads [7–9] that vary over the joint surfaces [10–12]. It is reasonable to assume that mechanical loads also affect the development of the collagen network, and that the differences in mechanical properties in adult life partly reflect differences in collagen network.
It is well known that the composition of the collagen network in AC changes during postnatal development: e.g. collagen type I is replaced by collagen type II [13–15], and the amount of collagen increases [5, 16, 17]. Structural remodelling of the collagen network in terms of the three zones has also been described for postnatal development [1, 18]. However, the remodelling capacity of the collagen network in adult AC is limited [8, 19–23]. Therefore, differences in collagen structure in adult AC are most likely already present at birth, or developed during postnatal cartilage maturation.
In adults, collagen fibre orientation and organisation are important parameters for the mechanical functions of AC [7, 18, 24, 25]. The collagen network partly determines AC stiffness [26–29] and Poisson ratio  and in transient loads it also influences interstitial fluid flow [31–34]. Across species and anatomical sites, adult AC shows what is known as 'a Benninghoff structure' [1, 17, 18, 25, 35–37]. This structure is characterised by three layers in AC: from articular surface to tidemark there is first a thin layer with collagen fibres mainly oriented parallel to the articular surface, second there is a thicker transitional zone where the collagen fibres appear to lack a predominant orientation, and third there is the deep zone where collagen fibres are mainly oriented perpendicular to the tidemark.
Recently, (semi) quantitative studies into collagen fibre orientation in stillborn and juvenile animals versus adults have confirmed that the Benninghoff structure is absent at birth and that it is developed during early postnatal life [17, 36, 38]. Van Turnhout et al.  looked at equine AC in stillborn and mature animals. They found that collagen fibre orientation in the stillborn animals was predominantly parallel to the articular surface throughout the complete cartilage layer. In the adult animals, they confirmed the Benninghoff structure. They also noticed that some zonal differentiation was already present in the stillborn AC, although this was not reflected in the depth-dependent collagen orientation. Rieppo et al.  looked at juvenile and adult porcine AC. For the juvenile animals, they also found that a major part of the collagen fibrils was predominantly oriented parallel to the articular surface throughout the complete cartilage layer. In the adult animals, they confirmed the Benninghoff structure. Julkunen et al.  looked at rabbit AC from age 4 weeks to 18 months at six time points. Their results confirm the findings by Van Turnhout et al.  and Rieppo et al. .
Since the collagen network is an important factor for AC mechanics, these 'remarkable structural alterations during development and growth'  in part explain the stiffening of the AC layer during development . Investigations into how the collagen network affects cartilage mechanics during postnatal development and the mechanical implications of the absence of a Benninghoff structure, can be performed with fibril reinforced finite element models, e.g. [40–42], provided that we have (quantitative) data on the collagen network at the time points of interest. Knowledge on the dynamics of collagen network development is also important for the validation of collagen remodelling algorithms, e.g. . And finite element modelling of the time line of collagen structure development in AC will provide insight into the mechanical environment that goes with this development. Finally, better insight into the mechanobiology of cartilage development will help to mimic this process in vitro, which is important in the field of cartilage tissue engineering.
Last year's studies [17, 36, 38] showed evidence of postnatal collagen network remodelling in AC, but they had their limitations in showing its dynamics. Van Turnhout et al.  used two sample points (stillborn and adult) and had a limited number of animals. Rieppo et al.  used three sample points (4 months, 11 months and 21 months) that did not include perinatal animals. The study by Julkunen et al.  does contain more and earlier sample points in a time resolved manner. However, perinatal animals were not included in this study and the number of early sample points is still too limited to convincingly show the dynamics of collagen structure remodelling. These authors did show differences in collagen structure development between different joint surfaces (femur and tibia).
The objective of the present research is therefore to quantify the development of the collagen network in a model animal from birth to maturity. We further aim to shed more light on the zonal differentiation at birth as noted by Van Turnhout et al. , and to assess differences in (the development of) the collagen network between different anatomical sites of a single joint surface. We will look at anatomical sites that are expected to have different loading histories. We use quantitative polarised light microscopy (qPLM), sometimes called 'the gold standard of histology' [44, 45], to investigate properties of the collagen network and we use the sheep (Ovis aries) as our model animal.
Results for the statistical model for cartilage depth
For the factor hind leg/fore leg, there is a main effect (table 1: B k , p < 0.001) and an interaction with age (table 1: (AB) jk , p = 0.021). For the main effect, the model predicts that cartilage is 6.40% ± 2.15% thicker in the fore legs than in the hind legs. For the interaction with age, cartilage in the fore legs is thicker than the cartilage in the hind legs at all ages but 0 weeks and 36 weeks (hind legs thicker, figure 2b). The interaction with age results in significant differences at 52 weeks (fore legs 14.3% ± 7.29% thicker) and 72 weeks (fore legs 13.7% ± 7.25% thicker).
For the factor lateral site/medial site, there is no main effect (table 1: C l , p = 0.053), but there is an interaction with age (table 1: (AC) jl , p < 0.001). This interaction does not show a trend (figure 2c) and is only significant at age 4 weeks (lateral site 15.8% ± 7.40% thicker).
For the factor caudal site/distal site/rostral site, there is a main effect (table 1: D m , p < 0.001) and an interaction with age (table 1: (AD) jm , p < 0.001). For the main effect, the model predicts that cartilage is 4.50% ± 2.12% thicker at the distal site than at the caudal site, and cartilage is 7.36% ± 2.18% thicker at the rostral site than at the caudal site. The interaction with age does not show a trend (figure 2d), but there are significant differences at ages 0 weeks (rostral 14.1% ± 7.31% thicker than caudal), 2 weeks (distal 25.0% ± 8.00% thicker than caudal, and distal 14.6% ± 7.33% thicker than rostral), 20 weeks (rostral 14.3% ± 7.32% thicker than distal), 52 weeks (distal 16.3% ± 7.45% thicker than caudal) and 72 weeks (rostral 16.3% ± 7.43% thicker than distal).
with D f the cartilage depth in μm and t the age in weeks.
Results for the statistical model for orientation index
There is a main effect for lateral sites versus medial sites (table 2: C l , p = 0.013). The model shows that is 0.22 higher at the lateral site. Inspection of the average orientation patterns for these sites per age (additional file 2) indicates two differences: 1) the medial site appears to be further in the development of the Benninghoff structure at ages 4 weeks, 12 weeks and 52 weeks; and 2) the depth where the predominant collagen fibre orientation is halfway in the transition from superficial zone to deep zone, is smaller for the medial site than for the lateral site at ages 28 weeks, 36 weeks and 52 weeks. Both effects result in a lower value of .
For the factor caudal site/distal site/rostral site, there is a main effect (table 2: D m , p < 0.001) and an interaction with age (table 2: (AD) jm , p = 0.011). For the main effect, the model predicts that is 0.06 higher at the distal site than at the caudal site, and is 0.08 higher at the rostral site than at the caudal site. The interaction with age shows significant differences between 12 weeks and 52 weeks (figure 4b). is higher at the rostral site than at the caudal site at 12 weeks (0.16), 20 weeks (0.09), 28 weeks (0.20) and 36 weeks (0.10). is higher at the distal site than at the caudal site at 28 weeks (0.21) and 52 weeks (0.17). Inspection of the average orientation patterns for these sites per age (additional file 3) shows that the caudal site develops its Benninghoff structure earlier in life than the rostral and distal site: predominant collagen orientation at a depth of 100 μm at 28 weeks is 104° at the caudal site against 127° at the rostral site and 125° at the distal site. The average orientation patterns do not show a trend for differences between development for the rostral site and distal site. At 72 weeks of age, all site-dependent variation in collagen orientation has vanished (figure 4b).
Results for the statistical model for the maximum retardance
There are no other main effects, but there are two significant interactions with the factor caudal site/distal site/rostral site (table 3: (AD) jm , p = 0.034 and (BD) km , p = 0.005). The interaction with age (AD) jm shows significant differences at ages 0 weeks (distal 30.3% ± 12.8% higher than rostral), 2 weeks (caudal 25.0% ± 12.2% higher than rostral, and distal 41.3% ± 13.8% higher than rostral) and 20 weeks (rostral 28.4% ± 12.6% higher than distal), but no trend (figure 6b). The interaction of caudal site/distal site/rostral site with hind legs/fore legs is due to the caudal site: Δmax is 18.7% ± 4.92%, higher in the hind legs than in the fore legs for the caudal site. The model shows no differences between hind legs and fore legs for the rostral and distal site.
Results for the statistical model for the position of the retardance minimum near the articular surface
The interaction of caudal site/distal site/rostral site with age (table 4: (AD) jm , p < 0.001) has significant differences at ages 0 weeks (rostral 17.2% ± 8.59% higher than caudal), 2 weeks (rostral 19.1% ± 8.72% higher than caudal, and rostral 21.7% ± 8.91% higher than distal), 4 weeks (rostral 28.4% ± 9.40% higher than distal), 28 weeks (rostral 17.4% ± 8.59% higher than caudal, and distal 16.4% ± 8.53% higher than caudal), 36 weeks (rostral 17.7% ± 8.62% higher than distal) and 52 weeks (distal 18.7% ± 8.69% higher than caudal, and distal 30.1% ± 9.52% higher than rostral), but these effect do not show a trend (figure 6d). There is also a significant main effect for lateral site/medial site (table 4: C l , p = 0.002): d|Δmin is 7.14% ± 2.31% higher at the lateral site.
The results on cartilage thickness are mostly in line with previous work. Brommer et al.  for instance, analysed site-dependent development of cartilage thickness in the equine proximal phalanx. As we did in our study, they found that cartilage thickness decreases during development, and they found no site-dependent differences in cartilage thickness in mature animals. In our study, cartilage at the caudal site is thinner than at the rostral and distal site and cartilage in the fore legs is thicker than in the hind legs, although the differences are small (< 8%). Cartilage thickness has been shown to respond to mechanical loads during development [7, 46–48] and the results for the development of the collagen network suggest influences of different loading regimes at different anatomical sites. However, cartilage thickness in the mature animal is thought to relate to joint congruency, i.e. the similarity of the two articulating cartilage surfaces [49, 50]. This explains why we find only small site-dependent differences in cartilage thickness in the nearly congruent joint that we investigated. This leaves us with the difference in cartilage thickness between hind legs and fore legs. Contrary to the remark by Simon  that cartilage is in general thicker in the hind limbs, we find that cartilage is thicker in the fore limbs, and that the difference increases slightly between birth and maturity.
As in previous studies on collagen structure in horses , pigs  and rabbits , we find that the perinatal predominant collagen fibre orientation is parallel to the articular surface throughout most of the cartilage depth and that this arrangement changes to an arcade-like Benninghoff structure as the cartilage matures. We also confirmed the presence of a retardance minimum in the perinatal animals as noted earlier by Van Turnhout et al. . The collagen remodelling process appears finished at 36 weeks of age in our sheep in this joint (e.g. figures 3b and 5b): the differences in qPLM results are small after this age and three of the four measured parameters (D, Δmax, and d|Δmin) do not change significantly between 36 weeks and 72 weeks of age. Hunziker et al.  found no noteworthy changes after puberty in the structural organization of rabbit AC and 36 weeks is also approximately when the sheep are sexually mature: the first ovulation or puberty in the sheep occurs at ≈ 6 months of age or at a mass of 30 kg to 35 kg, and this transition lasts ≈ 5 weeks . It is well before skeletal maturity, which occurs around 12 months of age [53–55]. Also, the slaughter mass results (additional file 4) show that the animals still developed after the age of 36 weeks.
The fact that collagen orientation appears to vary in the deep zone (figure 3a) is related to our decision to average sample profiles up to the mean depth of the sample pool. There are samples that are shorter than this mean depth and the number of samples that we (can) use to average therefore decreases from the minimum depth in the sample pool to the mean depth. The benefit of this procedure is that it better captures the superficial and transitional zones: e.g., the observation that the absolute depth of the retardance minimum appears invariant to age is evident when analysed in this manner. The downside is that from the minimum depth to the mean depth in the sample pool, we have an increasing proportion of data that is very close to the cartilage/bone interface. The azimuth data in particular may deviate from the deep zone data at the cartilage/bone interface, where the collagen ensures the integration of these two tissues. Inspection of the individual orientation patterns has shown that the azimuth is in fact fairly constant in the deep zone (additional file 5), and that the deviations towards the cartilage/bone interface in figure 3a are largely due to the increasing influence of cartilage/bone interfaces as described above. We find site dependent differences in collagen network remodelling (additional files 2 and 3). Both and d|Δmin are higher at the lateral site than at the medial site (tables 2 and 4, C l , p = 0.013 and p = 0.002). is higher for less advanced Benninghoff structures and for Benninghoff structures with higher d|Δmin. The results for d|Δmin suggest that its smaller value for the medial sites is correlated to the smaller values of . Thus, the transitional zone is at a larger depth at the lateral side, i.e. the superficial zone is thicker at the lateral side.
Also, the caudal site develops different from the distal and rostral sites: it develops its Benninghoff structure earlier in life (figure 4b, additional file 2) and maximum retardance is higher in the hind legs for the caudal site, but not for the distal and rostral site (table 3, (BD) km , p = 0.005). These observations suggests a loading regime during development that is different for the caudal sites compared to the rostral and distal sites. In the last sample point however, we find no site dependent differences in predominant collagen orientation patterns (additional file 2) or the orientation index (figure 4b). Gründer  reports a faster structural development for sheep tibial AC compared to femural AC, Julkunen et al.  showed that collagen structure development and adult collagen structure both differ between femur and tibia in rabbits, and Brama et al.  showed changes in collagen structure in juvenile foals as a result of exercise training. However, Egger et al.  did not find differences in collagen structure in adult elephant knee cartilage and they propose that this is due to the congruency in the elephant knee. This hypothesis is supported by our data: in the examined nearly congruent joint that was subjected to physiological loads, we also find no differences in collagen orientation between sites on the same cartilage surface for the oldest (adult) animals.
Our data shows that the predominant collagen orientation first changes at a depth that marks the proximal boundary of the transitional zone in the mature animal (figure 3a). The resulting intermediate pattern for predominant collagen orientation that we find at 12 weeks, is similar to the intermediate pattern found in a 4.5 month year old horse , in 11 month old pigs , and in rabbits of 4 weeks to 6 weeks . This intermediate pattern is present for all anatomical sites in the current study. This suggests that the depth-dependent pattern of collagen remodelling is similar across species, joints and anatomical sites, except for the timing (e.g. earlier development for the caudal site in this study, figure 4b and additional file 2). The difference in species (porcine versus ovine) and anatomical location (metacarpus versus knee) can explain why Rieppo et al.  found that collagen remodelling was not finished when their pigs reached sexual maturity.
The absence of a Benninghoff collagen structure at birth may have implications for postnatal morphogenesis. Results of simulations with finite element models [58, 59] and in vitro experimental results  show correlations between hydrostatic pressure and inhibition of ossification and between shear strains and promotion of ossification. An effect of the Benninghoff collagen structure is that it limits cartilage swelling, and it thereby effectively increases the hydrostatic (osmotic) pressure in the deep zone [61–63]. The development of the Benninghoff structure might thus serve to prevent ossification and maintain a functional cartilage layer in the adult animal. It would be interesting to see if cartilage ossification is promoted when the Benninghoff collagen structure cannot develop.
The retardance minimum near the articular surface is traditionally associated with the transitional zone, a layer with a weak anisotropic collagen arrangement, in adult animals [18, 35]. Our special interest in this retardance minimum appears to have been justified. Not only is it present at birth without evidence of differentiation in collagen orientation, which is confirmed by earlier work on equine  and porcine cartilage , it is also the only parameter that is similar in the youngest and oldest animals (figure 6c). During cartilage maturation, its thickness decreases to about one third of the perinatal thickness; collagen fibre orientation changes from a network predominantly parallel to the articular surface over the cartilage depth to a Benninghoff structure; and maximum retardance increases to about twice that of neonatal cartilage. Compared to these developmental changes, the depth of the retardance minimum (≈ 50 μm) appears relatively indifferent to the remodelling process (figure 6c). Thus, a layer with a weak anisotropic collagen arrangement, or less collagen, or both, is already formed before birth. This can for instance be a result of in utero movement and muscular loads on the joint [64, 65]. If so, then apparently the mechanical state in this layer does not change enough after birth to result in a much stronger anisotropic collagen fibre arrangement. If in utero joint loading is responsible for the emergence of this layer, we would predict its absence in animal models that e.g. lack skeletal muscle activity . Finite element models may help to explain what is so special about this layer, and its location.
With this paper, we provide essential information for analysis of the role of the collagen fibre network during development, e.g. by fibril reinforced finite element models [40–42]. It is the first time that the dynamics of post natal remodelling of collagen fibre organisation have been time and space resolved from birth to maturity. Therefore, we were able to present the first solid observation of a retardance minimum in perinatal animals, more or less at a fixed distance (≈ 50 μm) from the articular surface. This is an important step towards a better understanding of the mechanobiology of articular cartilage development. With additional information on collagen densities during development, we will be able to interpret the retardance minimum as a region with less collagen or a region with more anisotropic collagen and see how much this interpretation changes during development. With additional information on the development of glycosaminoglycan concentrations and fixed charge densities in the AC, it becomes possible to estimate the mechanical environment that drives the depth-dependent AC development in general, and depth-dependent collagen orientation remodelling in particular.
The collagen network in articular cartilage remodels between birth and sexual maturity from a network with predominant orientation parallel to the articular surface to a Benninghoff network. Changes in predominant collagen orientation occur first at a depth that marks the proximal boundary of the transitional zone in adult animals. This is the first solid observation of the dynamics of collagen structure remodelling, and provides important information for the further investigation of articular cartilage mechanobiology.
The fact that the retardance minimum near the articular surface is present at all ages shows that a zonal differentiation is already present in the perinatal animals. However, the zonal differentiation can not be correlated to changes in collagen network orientation in these animals. Furthermore, we were able to show that the depth of this retardance minimum is fairly constant during development. Without information on the development of collagen densities, we cannot tell if this retardance minimum is the result of a region with less collagen or a region with collagen with less anisotropic predominant orientations; and we cannot tell whether the relative contributions (less collagen and less anisotropic arrangement) change during development. We intent to address this in a future publication.
We observe differences in the dynamics of collagen remodelling between different anatomical sites, but in the adult animals we see no difference in collagen structure (i.e. they are all similar Benninghoff structures) or cartilage thickness. We expect that this is due to the congruency of the joint that we investigated. Comparison of the data in the current paper with other data on collagen remodelling in literature, suggests that the general pattern of collagen remodelling (i.e. where the first changes in collagen orientation occur) is similar between joints and species.
Animals were obtained from a local sheep farm. We obtained 50 female sheep for 9 sample points (five sheep for each sample point, and five spare sheep), at ages 2 weeks, 4 weeks, 8 weeks, 12 weeks, 20 weeks, 28 weeks, 36 weeks, 52 weeks and 72 weeks. An additional four stillborn lambs were used (labelled age = 0 weeks). Animals were kept at the farm with their mother until sacrifice or the age of 12 weeks. Animals that were 12 weeks old were collected at the farm and housed at the universities laboratory animal facility 'Ossekampen' until sacrifice.
Animals were kept together at pasture with ample space. They had ad lib access to water and food. The animals were born between 24th of February 2007 and 21st of March 2007. In September and October the disease blue tongue spread over the Netherlands and we had to house the animals temporarily inside a barn. Due to the disease, three infected animals had to be excluded from the experiment. All animals were closely monitored, and the other animals showed no clinical signs of blue tongue infection. A further three animals were excluded from the experiment due to accidents with a dog and a male sheep. Therefore, the total number of animals at the end of the experiment was 48, and the number of animals for the first sample point (0 weeks, stillborn) and the last sample point (72 weeks) was four. Animals were weighed prior to sacrifice. The experiment was approved was by the the Wageningen University Animal Experiments Committee.
Animals for the sample points from 2 weeks to 12 weeks were collected at the farm and sacrificed with an overdose of T61. Animals for the remaining sample points were brought to an abattoir to be sacrificed. The animals' legs were collected immediately following sacrifice. For each animal we randomly selected either the two left legs or the two right legs for qPLM analysis. At each sample point, the ratio left:right was either 3:2, 2:3 or 2:2.
We used a dental saw to take the medial and lateral hemispheres from the distal end of each cannon bone. These hemispheres were fixed with formalin and decalcified with EDTA (10% EDTA, pH 7.4) until the hemispheres could be cut with a razor blade. The hemispheres were then divided into a rostral, a distal and a caudal sample (figure 7). Of these, the distal site is expected to be subject to a more static load and the rostral and caudal sites are expected to be subject to a more intermittent load during locomotion . These samples were washed and infiltrated with sucrose (25% sucrose on PBS) overnight, snap frozen in liquid nitrogen and stored at -80°C until further processing, and finally cut parallel to the superficial split lines (as observed in two of the five animals) to 7 μm thick histological slices with a cryostat (Reichert 2800N).
Quantitative polarised light microscopy
The birefringent material in cartilage are the fibrils in the collagen network and qPLM measures two parameters of birefringent structures: azimuth and retardance. The azimuth φ does not measure individual fibril orientations, but is the predominant orientation of the collagen fibres in the network in the pixel. The retardance Δ is a measure of the amount of birefringent material (e.g. collagen fibrils) that is associated with azimuth φ : low retardance indicates one of three things: that 1) there is fewer collagen in this pixel; or that 2) the measured orientation belongs to a collagen network with a low level of anisotropy; or 3) a combination of these two.
Macroscopically normal histological samples were mounted with water and analysed with the LC-PolScope system for qPLM [68, 69]. Images were obtained with a Zeiss Axiovert 200M microscope at a 5×/1.6 magnification, equipped with a Q-imaging monochrome HR Retiga EX 1350 camera. Recorded intensity images had a resolution of 1.59 μm2/pixel and were stored in 8 bit TIFF format. We used the five frame setting with background correction as described by . The recorded images were analysed for predominant collagen fibril orientation and tissue retardance with custom written scripts implemented in Matlab (version 7.8.0 R2009a, The MathWorks, Inc., 1984-2005).
Azimuth and Retardance patterns
for W = 101 pixels over the width of the ROI, and for on the interval 0 < ≤ π . The summation over the absolute values of the inner dot product found the line with the smallest difference in angles compared to all W lines described by the azimuth values φ(w). The multiplication with the retardance Δ (w) assigned a smaller importance to the azimuth in pixels with low retardance, i.e. to pixels with less collagen associated with the predominant orientation φ(w) .
with D the number of points (pixels) over the depth of the sample (figure 8).
Average azimuth patterns between samples
Data were analysed with generalized linear mixed models because some of the variables analysed are not normally distributed. Also, measurements on the same animal and position within an animal are dependent. This excludes conventional analyses such analysis of variance or regression that are intended for normally distributed and independent data. We therefore used the penalized quasi-likelihood methodology described by Schall , Breslow & Clayton  and Engel & Keen . Calculations were performed with GenStat . The models comprised random effects with associated components of variance, that allowed for dependence between observations of the same animals and the same anatomical sites. Thus, we used a nested structure within animal for hind leg/fore leg, lateral/medial and caudal/distal/rostral sites. In particular, this allowed for additional dependence within animals between duplicate observations on the same site. We are interested in the development of differences between the different anatomical sites with age. Therefore, fixed effects (systematic effects) comprised main effects and all second order interactions for factors age, hind leg/fore leg, lateral/medial site and caudal/distal/rostral site in the initial models. Models were fitted separately to four response variables: cartilage thickness D, the orientation index , the maximum retardance value near the articular surface Δmax, and the position of the retardance minimum near the articular surface d|Δmin. For the variables D, Δmax and d|Δmin, we used a log link and gamma variance function, with a multiplicative dispersion parameter. For the variable , we used a logit link (logit p = log p/(1 - p)) and binomial variance function with an additional multiplicative dispersion parameter. Random effects on the link scale were assumed to follow normal distributions. Tests were based on an approximate F-test  applied to the adjusted dependent variate from the last iteration step of the iterative re-weighted restricted maximum likelihood algorithm  that we used. The link functions provide the relationship between the linear predictor and the mean of the distribution function and the chosen link and variance functions were needed to achieve satisfactory (normally distributed) residuals for the models. Non-significant (p > 0.05) higher order interactions were dropped from the initial models. We used the following symbols in the models: μ: intercept; A j , j = 0, 2, 4, 8, 12, 20, 28, 36, 52, 72: age in weeks; B k , k = 1, 2: fixed factor hind leg/fore leg; C l , l = 1, 2: fixed factor lateral site/medial site; D m , m = 1, 2 3: fixed factor caudal site/distal site/rostral site; L i : random factor individual lamb; and (LB) ik , (LBC) ikl and (LBCD) iklm nested random factors within lamb. The final model that we fitted for each covariate, is presented in the results section. In the text, we quantify significant differences as mean ± standard error as predicted by the model. For the model with the logit link (for ), we can only present the mean, and not the standard errors. In the figures, we use raw means and associated standard deviations, and not model predictions, to illustrate the results.
We kindly thank Wim van de Laan for providing the sheep. From Wageningen University, we thank the people at De Haar/Ossekampen laboratory animal facility for taking care of the sheep; and Dr. Norbert de Ruijter at the Plant Cell Biology group for support with the qPLM equipment. In our own group we thank Adwin van den Brink for help with tissue processing; and Dr. Martin Lankheet for his critical review of the manuscript.
- Hunziker E, Kapfinger E, Geiss J: The structural architecture of adult mammalian articular cartilage evolves by a synchronized process of tissue resorption and neoformation during postnatal development. Osteoarthritis & Cartilage. 2007, 15 (4): 403-413. 10.1016/j.joca.2006.09.010. [Corrigendum in Osteoarthritis & Cartilage 15 (2007) 403-413.]View ArticleGoogle Scholar
- Archer CW, Morrison H, Pitsillides AA: Cellular aspects of the development of diarthrodial joints and articular cartilage. Journal of Anatomy. 1994, 184 (Pt 3): 447-456.PubMed CentralPubMedGoogle Scholar
- Hayes A, MacPherson S, Morrison H, Dowthwaite G, Archer C: The development of articular cartilage: evidence for an appositional growth mechanism. Anatomy and Embryology. 2001, 203 (6): 469-479. 10.1007/s004290100178.View ArticlePubMedGoogle Scholar
- Brama P, Tekoppele J, Bank R, Karssenberg D, Barneveld A, Van Weeren P: Topographical mapping of biochemical properties of articular cartilage in the equine fetlock joint. Equine Veterinary Journal. 2000, 32: 19-26. 10.2746/042516400777612062.View ArticlePubMedGoogle Scholar
- Brama P, TeKoppele J, Bank R, Barneveld A, Van Weeren P: Functional adaptation of equine articular cartilage: the formation of regional biochemical characteristics up to age one year. Equine Veterinary Journal. 2000, 32 (3): 217-221. 10.2746/042516400776563626.View ArticlePubMedGoogle Scholar
- Brommer H, Brama P, Laasanen M, Helminen H, Van Weeren J PRand Jurvelin: Functional adaptation of articular cartilage from birth to maturity under the influence of loading: a biomechanical analysis. Equine Veterinary Journal. 2005, 37 (2): 148-154. 10.2746/0425164054223769.View ArticlePubMedGoogle Scholar
- Helminen HJ, Hyttinen MM, Lammi MJ, Arokoski JP, Lapveteläinen T, Jurvelin J, Kiviranta I, Tammi MI: Regular joint loading in youth assists in the establishment and strengthening of the collagen network of articular cartilage and contributes to the prevention of osteoarthrosis later in life: a hypothesis. Journal of Bone and Mineral Metabolism. 2000, 18 (5): 245-257. 10.1007/PL00010638.View ArticlePubMedGoogle Scholar
- Murray R, Birch H, Lakhani K, Goodship A: Biochemical composition of equine carpal articular cartilage is influenced by short-term exercise in a site-specific manner. Osteoarthritis & Cartilage. 2001, 9 (7): 625-632. 10.1053/joca.2001.0462.View ArticleGoogle Scholar
- Brama P, TeKoppele J, Bank R, Van Weeren P, Barneveld A: Influence of different exercise levels and age on the biochemical characteristics of immature equine articular cartilage. Equine Veterinary Journal Supplement. 1999, 55-61. 31Google Scholar
- Brama P, Karssenberg D, Barneveld A, Van Weeren P: Contact areas and pressure distribution on the proximal articular surface of the proximal phalanx under sagittal plane loading. Equine Veterinary Journal. 2001, 33: 26-32. 10.2746/042516401776767377.View ArticlePubMedGoogle Scholar
- Hodge W, Fijan R, Carlson K, Burgess R, Harris W, Mann R: Contact pressures in the human hip joint measured in vivo. Proceedings of the National Academy of Sciences of the United States of America. 1986, 83 (9): 2879-2883. 10.1073/pnas.83.9.2879.PubMed CentralView ArticlePubMedGoogle Scholar
- Palmer JL, Bertone AL, Litsky A: Contact area and pressure distribution changes of the equine third carpal bone during loading. Equine Veterinary Journal. 1994, 26 (3): 197-202. 10.1111/j.2042-3306.1994.tb04369.x.View ArticlePubMedGoogle Scholar
- Bland YS, Ashhurst DE: Development and ageing of the articular cartilage of the rabbit knee joint: distribution of the fibrillar collagens. Anatomy and Embryology. 1996, 194 (6): 607-619. 10.1007/BF00187473.View ArticlePubMedGoogle Scholar
- Bland YS, Ashhurst DE: Changes in the content of the fibrillar collagens and the expression of their mRNAs in the menisci of the rabbit knee joint during development and ageing. Histochemical Journal. 1996, 28 (4): 265-274. 10.1007/BF02409014.View ArticlePubMedGoogle Scholar
- Morrison E, Ferguson M, Bayliss M, Archer C: The development of articular cartilage: I. The spatial and temporal patterns of collagen types. Journal of Anatomy. 1996, 189 (Pt 1): 9-22.PubMed CentralPubMedGoogle Scholar
- Williamson AK, Chen AC, Masuda K, Thonar EJM, Sah RL: Tensile mechanical properties of bovine articular cartilage: variations with growth and relationships to collagen network components. Journal of Orthopaedic Research. 2003, 21 (5): 872-880. 10.1016/S0736-0266(03)00030-5.View ArticlePubMedGoogle Scholar
- Rieppo J, Hyttinen M, Halmesmaki E, Ruotsalainen H, Vasara A, Kiviranta I, Jurvelin J, Helminen H: Changes in spatial collagen content and collagen network architecture in porcine articular cartilage during growth and maturation. Osteoarthritis & Cartilage. 2009, 17 (4): 448-455. 10.1016/j.joca.2008.09.004.View ArticleGoogle Scholar
- Hughes L, Archer C, Gwynn Ia: The ultrastructure of mouse articular cartilage: collagen orientation and implications for tissue functionality. A polarised light and scanning electron microscope study and review. European Cells & Materials Journal. 2005, 9: 68-84.Google Scholar
- Hyttinen M, Arokoski J, Parkkinen J, Lammi M, Lapveteläinen T, Mauranen K, Király K, Tammi M, Helminen H: Age matters: collagen birefringence of superficial articular cartilage is increased in young guinea-pigs but decreased in older animals after identical physiological type of joint loading. Osteoarthritis & Cartilage. 2001, 9 (8): 694-701. 10.1053/joca.2001.0466.View ArticleGoogle Scholar
- Brama P, Tekoppele J, Bank R, Barneveld A, Firth E, Van Weeren P: The influence of strenuous exercise on collagen characteristics of articular cartilage in Thoroughbreds age 2 years. Equine Veterinary Journal. 2000, 32 (6): 551-554. 10.2746/042516400777584613.View ArticlePubMedGoogle Scholar
- Saadat E, Lan H, Majumdar S, Rempel DM, King KB: Long-term cyclical in vivo loading increases cartilage proteoglycan content in a spatially specific manner: an infrared microspectroscopic imaging and polarized light microscopy study. Arthritis Research & Therapy. 2006, 8 (5): R147-10.1186/ar2040.View ArticleGoogle Scholar
- Tran-Khanh N, Hoemann CD, McKee MD, Henderson JE, Buschmann MD: Aged bovine chondrocytes display a diminished capacity to produce a collagen-rich, mechanically functional cartilage extracellular matrix. Journal of Orthopaedic Research. 2005, 23 (6): 1354-1362. 10.1016/j.orthres.2005.05.009.1100230617.View ArticlePubMedGoogle Scholar
- Vaughan-Thomas A, Dudhia J, Bayliss MT, Kadler KE, Duance VC: Modification of the composition of articular cartilage collagen fibrils with increasing age. Connective Tissue Research. 2008, 49 (5): 374-382. 10.1080/03008200802325417.View ArticlePubMedGoogle Scholar
- Rieppo J, Töyräs J, Nieminen MT, Kovanen V, Hyttinen MM, Korhonen RK, Jurvelin JS, Helminen HJ: Structure-function relationships in enzymatically modified articular cartilage. Cells Tissues Organs. 2003, 175 (3): 121-132. 10.1159/000074628.View ArticlePubMedGoogle Scholar
- Egger GF, Witter K, Weissengruber G, Forstenpointner G: Articular cartilage in the knee joint of the African elephant, Loxodonta africana, Blumenbach 1797. Journal of Morphology. 2008, 269: 118-127. 10.1002/jmor.10600.View ArticlePubMedGoogle Scholar
- Wu J, Kirk T: Study of altered mechanical properties of articular cartilage in relation to the collagen network. Advanced Materials Research. 2008, 41-42: 9-14. 10.4028/www.scientific.net/AMR.41-42.9.View ArticleGoogle Scholar
- Shirazi R, Shirazi-Adl A, Hurtig M: Role of cartilage collagen fibrils networks in knee joint biomechanics under compression. Journal of Biomechanics. 2008, 41 (16): 3340-3348. 10.1016/j.jbiomech.2008.09.033.View ArticlePubMedGoogle Scholar
- Park S, Nicoll SB, Mauck RL, Ateshian GA: Cartilage mechanical response under dynamic compression at physiological stress levels following collagenase digestion. Annals of Biomedical Engineering. 2008, 36 (3): 425-434. 10.1007/s10439-007-9431-6.View ArticlePubMedGoogle Scholar
- Korhonen RK, Julkunen P, Wilson W, Herzog W: Importance of collagen orientation and depth-dependent fixed charge densities of cartilage om mechanical behavior of chondrocytes. Journal of Biomedical Engineering. 2008, 130 (2): 21003-1-21003-11.Google Scholar
- Kiviranta P, Rieppo J, Korhonen RK, Julkunen P, Töyräs J, Jurvelin JS: Collagen network primarily controls Poisson's ratio of bovine articular cartilage in compression. Journal of Orthopaedic Research. 2006, 24 (4): 690-699. 10.1002/jor.20107.View ArticlePubMedGoogle Scholar
- Visser Sd, Bowden J, Wentrup-Byrne E, Rintoul L, Bostrom T, Pope J, Momot K: Anisotropy of collagen fibre alignment in bovine cartilage: comparison of polarised light microscopy and spatially resolved diffusion-tensor measurements. Osteoarthritis & Cartilage. 2008, 16 (6): 689-697. 10.1016/j.joca.2007.09.015.View ArticleGoogle Scholar
- Lötjönen P, Julkunen P, Töyröas J, Lammi MJ, Jurvelin JS, Nieminen HJ: Strain-dependent modulation of ultrasound speed in articular cartilage under dynamic compression. Ultrasound in Medicine and Biology. 2009Google Scholar
- Li L, Korhonen R, Iivarinen J, Jurvelin J, Herzog W: Fluid pressure driven fibril reinforcement in creep and relaxation tests of articular cartilage. Medical Engineering and Physics. 2008, 30 (2): 182-189. 10.1016/j.medengphy.2007.03.001.View ArticlePubMedGoogle Scholar
- Federico S, Herzog W: On the anisotropy and inhomogeneity of permeability in articular cartilage. Biomechanics and Modeling in Mechanobiology. 2008, 7 (5): 367-378. 10.1007/s10237-007-0091-0.View ArticlePubMedGoogle Scholar
- Benninghoff A: Form und Bau der Gelenkknorpel in ihren Beziehungen zur Funktion. Zweiter Teil: Der Aufbau des Gelenkknorpels in seinen Beziehungen zur Funktion. Zeitschrift für Zellforschung und Mikroskopische Anatomie. 1925, 2: 783-862. 10.1007/BF00583443.View ArticleGoogle Scholar
- Van Turnhout MC, Haazelager MB, Gijsen MA, Schipper H, Kranenbarg S, Van Leeuwen JL: Quantitative description of collagen structure in the articular cartilage of the young and adult equine distal metacarpus. Animal Biology. 2008, 58 (4): 353-370. 10.1163/157075608X383674.View ArticleGoogle Scholar
- Zambrano NZ, Montes GS, Shigihara KM, Sanchez EM, Junqueira LC: Collagen arrangement in cartilages. Acta Anatomica. 1982, 113: 26-38.View ArticlePubMedGoogle Scholar
- Julkunen P, Iivarinen J, Brama P, Arokoski J, Jurvelin J, Helminen H: Maturation of collagen fibril network structure in tibial and femoral cartilage of rabbits. Osteoarthritis & Cartilage. 2010, 18 (3): 406-415. 10.1016/j.joca.2009.11.007.View ArticleGoogle Scholar
- Julkunen P, Harjula T, Iivarinen J, Marjanen J, Seppänen K, Närhi T, Arokoski J, Lammi MJ, Brama P, Jurvelin J, Helminen H: Biomechanical, biochemical and structural correlations in immature and mature rabbit articular cartilage. Osteoarthritis & Cartilage. 2009, 17 (12): 1628-1638. 10.1016/j.joca.2009.07.002.View ArticleGoogle Scholar
- Korhonen RK, Laasanen MS, Töyräs J, Lappalainen R, Helminen HJ, Jurvelin JS: Fibril reinforced poroelastic model predicts specifically mechanical behavior of normal, proteoglycan depleted and collagen degraded articular cartilage. Journal of Biomechanics. 2003, 36 (9): 1373-1379. 10.1016/S0021-9290(03)00069-1.View ArticlePubMedGoogle Scholar
- Wilson W, Huyghe J, Van Donkelaar C: Depth-dependent compressive equilibrium properties of articular cartilage explained by its composition. Biomechanics and Modeling in Mechanobiology. 2007, 6 (1-2): 43-53. 10.1007/s10237-006-0044-z.View ArticlePubMedGoogle Scholar
- Li L, Cheung J, Herzog W: Three-dimensional fibril-reinforced finite element model of articular cartilage. Medical and Biological Engineering and Computing. 2009, 47 (6): 607-615. 10.1007/s11517-009-0469-5.View ArticlePubMedGoogle Scholar
- Wilson W, Driessen N, Van Donkelaar C, Ito K: Prediction of collagen orientation in articular cartilage by a collagen remodeling algorithm. Osteoarthritis & Cartilage. 2006, 14 (11): 1196-1202. 10.1016/j.joca.2006.05.006.View ArticleGoogle Scholar
- Xia Y, Ramakrishnan N, Bidthanapally A: The depth-dependent anisotropy of articular cartilage by Fourier-transform infrared imaging (FTIRI). Osteoarthritis & Cartilage. 2007, 15 (7): 780-788. 10.1016/j.joca.2007.01.007.View ArticleGoogle Scholar
- Alhadlaq H, Xia Y, Moody J, Matyas J: Detecting structural changes in early experimental osteoarthritis of tibial cartilage by microscopic magnetic resonance imaging and polarised light microscopy. Annals of the Rheumatic Diseases. 2004, 63 (6): 709-717. 10.1136/ard.2003.011783.PubMed CentralView ArticlePubMedGoogle Scholar
- Kiviranta I, Tammi M, Jurvelin J, Arokoski J, Säämänen A, Helminen H: Articular cartilage thickness and glycosaminoglycan distribution in the young canine knee joint after remobilization of the immobilized limb. Journal of Orthopaedic Research. 1994, 12 (2): 161-167. 10.1002/jor.1100120203.View ArticlePubMedGoogle Scholar
- Kiviranta I, Tammi M, Jurvelin J, Arokoski J, Säämänen AM, Helminen HJ: Articular cartilage thickness and glycosaminoglycan distribution in the canine knee joint after strenuous running exercise. Clinical Orthopaedics and Related Research. 1992, 302-308. 283Google Scholar
- Oettmeier R, Arokoski J, Roth AJ, Helminen HJ, Tammi M, Abendroth K: Quantitative study of articular cartilage and subchondral bone remodeling in the knee joint of dogs after strenuous running training. Journal of Bone and Mineral Research. 1992, 7 (Suppl 2): S419-S424. 10.1002/jbmr.5650071410.View ArticlePubMedGoogle Scholar
- Simon WH, Friedenberg S, Richardson S: Joint congruence. A correlation of joint congruence and thickness of articular cartilage in dogs. Journal of Bone and Joint Surgery. American Volume. 1973, 55 (8): 1614-1620.Google Scholar
- Braune W, Fischer O: Bewegungen des Kniegelenks nach einer neuen Methode am lebenden Menschen gemessen. Abhandlungen der Mathematisch-Physischen Klasse der Königlich-Sächsischen Gesellschaft der Wissenschaften, Leipzig. Edited by: Hirzel S. 1891, 17: 75-150.Google Scholar
- Simon W: Scale effects in animal joints. I. Articular cartilage thickness and compressive stress. Arthritis and Rheumatism. 1970, 13 (3): 244-256. 10.1002/art.1780130305.View ArticlePubMedGoogle Scholar
- Foster DL, Jackson LM: Puberty in the Sheep. Knobil and Neill's Physiology of Reproduction. Edited by: Neill JD. 2006, Elsevier Inc, 2: 2127-2176. full_text. 3View ArticleGoogle Scholar
- Mastrangelo AN, Magarian EM, Palmer MP, Vavken P, Murray MM: The effect of skeletal maturity on the regenerative function of intrinsic ACL cells. Journal of Orthopaedic Research. 2010, 28 (5): 644-651.PubMed CentralPubMedGoogle Scholar
- Hopkins D, Stanley D, Martin L, Gilmour A: Genotype and age effects on sheep meat production 1. Production and growth. Australian Journal of Experimental Agriculture. 2007, 47 (10): 1119-1127. 10.1071/EA06236.View ArticleGoogle Scholar
- Cake M, Boyce M, Gardner G, Hopkins D, Pethick D: Genotype and gender effects on sheep limb bone growth and maturation: selection for loin depth causes bone hypotrophy. Australian Journal of Experimental Agriculture. 2007, 47 (10): 1128-1136. 10.1071/EA07058.View ArticleGoogle Scholar
- Gründer W: MRI assessment of cartilage ultrastructure. NMR in Biomedicine. 2006, 19 (7): 855-876. 10.1002/nbm.1092.View ArticlePubMedGoogle Scholar
- Brama P, Holopainen J, Van Weeren P, Firth E, Helminen H, Hyttinen M: Effect of loading on the organization of the collagen fibril network in juvenile equine articular cartilage. Journal of Orthopaedic Research. 2009, 27 (9): 1226-1234. 10.1002/jor.20866.View ArticlePubMedGoogle Scholar
- Beaupré GS, Stevens SS, Carter DR: Mechanobiology in the development, maintenance, and degeneration of articular cartilage. Journal of Rehabilitation Research and Development. 2000, 37 (2): 145-151.PubMedGoogle Scholar
- Carter DR, Wong M: Modelling cartilage mechanobiology. Philosophical Transactions of the Royal Society of London. Series B: Biological Sciences. 2003, 358 (1437): 1461-1471. 10.1098/rstb.2003.1346.PubMed CentralView ArticlePubMedGoogle Scholar
- Smith RL, Carter DR, Schurman DJ: Pressure and shear differentially alter human articular chondrocyte metabolism: a review. Clinical Orthopaedics and Related Research. 2004, S89-S95. 427 SupplGoogle Scholar
- Maroudas A: Balance between swelling pressure and collagen tension in normal and degenerate cartilage. Nature. 1976, 260 (5554): 808-809. 10.1038/260808a0.View ArticlePubMedGoogle Scholar
- Bank RA, Soudry M, Maroudas A, Mizrahi J, TeKoppele JM: The increased swelling and instantaneous deformation of osteoarthritic cartilage is highly correlated with collagen degradation. Arthritis and Rheumatism. 2000, 43 (10): 2202-2210. 10.1002/1529-0131(200010)43:10<2202::AID-ANR7>3.0.CO;2-E.View ArticlePubMedGoogle Scholar
- Ateshian GA, Soltz MA, Mauck RL, Basalo IM, Hung CT, Lai WM: The role of osmotic pressure and tension-compression nonlinearity in the frictional response of articular cartilage. Transport in Porous Media. 2003, 50 (1-2): 5-33. 10.1023/A:1020618514874.View ArticleGoogle Scholar
- Hamrick MW: A chondral modeling theory revisited. Journal of Theoretical Biology. 1999, 201 (3): 201-208. 10.1006/jtbi.1999.1025.View ArticlePubMedGoogle Scholar
- Onyekwelu I, Goldring MB, Hidaka C: Chondrogenesis, joint formation, and articular cartilage regeneration. Journal of Cellular Biochemistry. 2009, 107 (3): 383-392. 10.1002/jcb.22149.View ArticlePubMedGoogle Scholar
- Gomez C, David V, Peet NM, Vico L, Chenu C, Malaval L, Skerry TM: Absence of mechanical loading in utero influences bone mass and architecture but not innervation in MyoD-Myf5-deficient mice. Journal of Anatomy. 2007, 210 (3): 259-271. 10.1111/j.1469-7580.2007.00698.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Van Turnhout MC, Kranenbarg S, Van Leeuwen JL: Modelling optical behaviour of birefringent biological tissues for evaluation of quantitative polarised light microscopy. Journal of Biomedical Optics. 2009, 14 (5): 10.1117/1.3241986. 054018-1-11Google Scholar
- Oldenbourg R, Mei G: New polarized light microscope with precision universal compensator. Journal of Microscopy. 1995, 180 (Pt 2): 140-147.View ArticlePubMedGoogle Scholar
- Oldenbourg R: Polarization microscopy with the LC-microscope. Live Cell Imaging: A Laboratory Manual. Edited by: Goldman RD, Spector DL. 2004, Cold Spring Harbor Laboratory Press, 205-238.Google Scholar
- Shribak M, Oldenbourg R: Techniques for fast and sensitive measurements of two-dimensional birefringence distributions. Applied Optics. 2003, 42 (16): 3009-3017. 10.1364/AO.42.003009.View ArticlePubMedGoogle Scholar
- Upton GJ, Fingleton B: Spatial data analysis by example - Volume 2: Categorical and directional data. 1989, Wiley series in probability and mathematical statistics, John Wiley and SonsGoogle Scholar
- Schall R: Estimation in generalized linear models with random effects. Biometrika. 1991, 78 (4): 719-727. 10.1093/biomet/78.4.719.View ArticleGoogle Scholar
- Breslow N, Clayton D: Approximate inference in generalized linear mixed models. Journal of the American Statistical Association. 1993, 88 (421): 9-25. 10.2307/2290687.Google Scholar
- Engel B, Keen A: A simple approach for the analysis of generalized linear mixed models. Statistica Neerlandica. 1994, 48: 1-22. 10.1111/j.1467-9574.1994.tb01428.x.View ArticleGoogle Scholar
- Payne R, Harding S, Murray D, Soutar D, Baird D, Glaser A, Channing I, Welham S, Gilmour A, Thompson R, Webster R: The Guide to GenStat Release 12. VSN International, Hemel Hempstead. 2009Google Scholar
- Kenward MG, Roger JH: Small sample inference for fixed effects from restricted maximum likelihood. Biometrics. 1997, 53 (3): 983-997. 10.2307/2533558.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.