Relationship between Somatic Cell Score and Fat Plus Protein Yield in the First Three Lactations in Spanish Florida Goats

: The aim of this study was to estimate genetic parameters of somatic cell score (SCS) and fat plus protein yield (FPY) using repeatability (RM) and random regression (RRM) models in Florida goats. The data consisted of 340,654 test-day controls of the first three lactations, and the pedigree contained 36,144 animals. Covariance components were estimated with a bivariate RM and RRM using the REML approach. Both models included as fixed effects the combination of herd and control date, litter size, kidding number and lactation length, and as random effects, the additive genetic and permanent environmental effects. A variation in the shape of the genetic parameters along the lactation curve was observed for both traits, and h 2 oscillated between 0.272 and 0.279 for SCS and 0.099 and 0.138 for FPY. The genetic correlation between SCS and FPY was negative and medium ( − 0.304 to − 0.477), indicating that a low-SCS EBV is associated with a genetic predisposition to high FPY production. Our results showed that given the magnitude of h 2 for SCS and its r g with FPY, the SCS could be used as a selection criterion to increase resistance to mastitis, thus obtaining an improved dairy and cheese aptitude in this breed.


Introduction
According to Miller and Lu [1], the goat dairy industry is increasing its activity worldwide.This trend is due not only to the nutritional value of goat milk [2] but also to the additional attributes of its consumption behavior and habits, which give this species an essential role in maintaining the environment of the region where it is farmed by decreasing the risk of soil deterioration, reducing the possibility of fires and keeping many marginal areas of Europe productive [3].On the other hand, increases of over 7% in the demand for goat milk for the cosmetics industry are predicted up to 2025 [4], which means that the prospects for goat farming are secure.
At the European level, Spain is the second largest producer of goat milk [5], with the Andalusia region being where the largest number of goats are concentrated and where around 40% of the milk of this species is produced, mostly for the manufacture of different types of cheese.Among the indigenous goat breeds in this country, the Florida breed is the one with the greatest productive potential, and it has been subjected to an improvement program for over 20 years, with highly satisfactory results [6].A recent publication indicates that cumulative milk production in the period up to 240 days of lactation in the Florida breed has increased at a rate of 1.7% per year.However, this may, as a collateral effect, have increased potential infection problems, which is reflected in aspects such as somatic cell count (SCC), due to the constant interaction between the udder and the milking equipment [7].Although a certain level of SCC may be tolerated as an indicator of resistance to infection, immunity and health, it may also represent risks to product quality and could result in economic losses.
High SCC or total bacterial count (TBC) represents a significant problem in the dairy industry, leading to important economic losses, not to mention the hygienic and legal consequences or the major impact it has on food safety and public health [8].One of the factors leading to an increase in these parameters in milk is subclinical mastitis (SCM).
Bulk milk bacterial cells can enter from the environment during milking or from bacterial growth because of insufficient cleaning and sanitation of the system, but in some cases, they can also enter from the goat's mammary gland.Therefore, the TBC of bulk milk can sometimes be used to monitor the level of mastitis on a farm [9].Indeed, the identification of risk factors for high SCC is of great importance in regions where SCC regulations are already applied in the goat industry, as well as in regions undergoing changes towards the implementation of programs and regulations for goat milk quality.However, there is limited information currently available about the on-farm risk factors associated with SCC levels in goat bulk milk [10,11].
The Florida is a native Spanish breed of dairy goat distributed mostly in the center and south of Spain, raised under a wide variety of systems of production, ranging from semi-extensive to semi-intensive systems [12].The main selection criterion in the breeding program of this breed is the fat plus protein yield per lactation [13].
In this breed, Jiménez-Granado et al. [7] showed how the levels of this parameter have increased alarmingly over the last decade, which is consistent with the levels of mastitis perceived by the farmer.It is therefore considered necessary to establish a plan to control and prevent mastitis, which combines sanitary and management measures with the search for greater resistance to this disease.Since the year 2004, the collection of records on SCC in goat milk has become mandatory throughout the EU.So, with the availability of this information and the well-known relationship between mastitis and SCC, the genetic study of this variable has become possible, with the aim of improving resistance to mastitis.In fact, there is already evidence of positive results of a breeding program to decrease the level of SCC [14,15], with low genetic relationships with milk production and its components.In order to determine the optimal strategy that takes into account not only somatic cell levels but also their interrelationships with the other selection criteria, it is essential to estimate the genetic relationship between this parameter and the main selection criteria in the breed.
Accordingly, the aim of this study was to estimate genetic parameters of somatic cell score (SCS) and fat plus protein yield across lactations and kidding number, using repeatability and random regression models.

Phenotypic Data and Pedigree
For this study, a total of 1,031,143 test-day (TD) records collected between 2005 and 2019 were accessed from the National Association of Florida Goat Breeders (ACRIFLOR).These TD records belonged to 90 herds and included a total of 58,606 females, bred on farms in the regions of Andalusia and Extremadura, which share similar climatic conditions, with an average temperature of 16.38 • C (between 14 • C and 24 • C) and an average relative humidity of 65.7% (between 37% and 78%; see supplementary materials Map S1).All the data were subjected to a data-editing process, and all lactations longer than 305 days, parity numbers greater than 6 and records with somatic cell count (SCC) outside the usual range (<40 × 10 3 to >10,000 × 10 3 ) were excluded.In the end, we used for the variance component analysis a total of 340,654 TDs recorded in the first 43 weeks of lactation of the first three kiddings from January 2006 to November 2019, from 27,479 daughters of 941 sires and 16,243 dams, of which 8788 were in the data vector.The pedigree contained all the known ancestors of the phenotyped animals, with a total of 36,144 animals.The variables analyzed were SCC transformed into SCS = log2 (SCC) + 3 [16] and the total daily amount of fat plus protein (FPY) expressed in grams.

Statistical Analysis
The traits studied were recorded over a fixed time scale (over the first three lactations), so their statistical processing was applied in stages.First, SCS and FPY were analyzed using a linear fixed effects model in order to obtain initial information on the causes of variation and also to be able to represent the forms of responses of both dependent variables at each kidding along the lactation trajectory.In this analysis, we considered the fixed effects of a combination of herd and control date (HCD i 8618 levels), litter size (LS j 4 levels), kidding number (NP k 3 classes) and lactation length expressed in weeks (DIM wl 43 weeks).
In the initial analysis, the (co)variance components were studied using bivariate models, according to a repeatability model (RM), which assumed no variation across lactation, and each dependent variable was expressed as an average for each kidding, as follows: where y 1 and y 2 represent the phenotypic values for SCS and FPY, respectively.The elements b i represent the same fixed effects mentioned previously; µ i are the vectors for the additive genetic random effect and p i is the permanent environmental random effect due to repetitions of the same observations in the animal.The matrices X i , Z i and W i are incidence matrices connecting the fixed and random effects to the vector of dependent variables.Finally, e i represents the random vector of the error.Subsequently, a second analysis was performed with a random regression model (RRM) analyzing the variables SCS and FPY recorded in each TD throughout lactation (DIM wI ).The representation of this model was similar to the RM, but in the RRM, the elements of a Legendre polynomial of order r were included as covariates.(Φr) of order r, as a fixed covariate in X i and a random covariate in Z i , was included to estimate the evolution of the (co)variance components between both dependent variables along the lactation trajectory expressed as DIM wi .
Both models were run for each of the first three lactations and the total number of lactations, using the statistical program Asreml 3 [17].The expected (co)variance components for the RM were: where the elements σ 2 µ 1 ; σ 2 pi and σ 2 ei stand for the genetic, permanent environment and residual variances for lscc (X 1 ) and fpy (X 2 ), respectively, σ µ 12 is their covariance, I W and l n are identity matrices and A is the denominator of the kinship relationship.The genetic heritability parameters (h 2 xi ) and the genetic correlations between both dependent variables (r g12 ) were estimated according to classical formulas [18].
Unlike the RM, in RRM, the genetic (co)variance components for both dependent variables were not considered to be the same for each point on the lactation scale, which means that there will be parameters of h 2 xi ; r g12 and EGV i for each i th point of DIM wi .Traits which are expressed longitudinally require an additional procedure presented by Jamrozik and Schaeffer [19], although basically the structure of expected (co)variances and estimates was similar: In this representation, I only changes the point expression of µ i by α i , which is a linear function for SCS and FPY throughout the lactation expressed in terms of a Legendre polynomial (Φ d ).The term ⊗ is a symbol of the Kronecker operator and C α is a matrix containing the eigen elements of the polynomial used, while the rest of the terms are the same as previously presented.The matrix C α t has a complex structure consisting of four , with, on the diagonal (c α1 and c α2 ), the genetic (co)variance components for each trait, with the elements corresponding to a polynomial of order r = 2, which best fits the data, and, outside the diagonal (c α12 = c α21 ), the covariances between all the terms of each variable.The structure of each of these submatrices was: iα ; σ 2 sα and σ 2 qα are the variances of the intercept, slope and quadratic term, respectively, while σ isα ; σ iqα and σ sqα are the corresponding covariances of the elements of the Legendre polynomial Φ d .This longitudinal procedure is more complex but offers multiple advantages, as shown below.
With the results of C αi and the coefficients of Φ d , the heritability values (h 2 α i ) and the genetic correlations (r g α ), both within and between both dependent variables at each point of the dimw scale including Φ d , were calculated as follows: , the same formula was used for i = SCS and FPY.
The intrapartum genetic correlations used for each variable throughout lactation were: while the intrakidding genetic correlations for each pair of variables throughout lactation were: Finally, the expected genetic values for SCS and FPY (EGV X i ) were estimated by direct solution in the RM, while RRM at each of the three births and the total at any point along the DIM w trajectory were estimated by: where α i is the solution of the RRM corresponding to animal i and consists of the corresponding elements of the Legendre polynomial (Φ) of order 2, in this case a quadratic equation.By manipulating the terms of Φ, the EGV for cumulative production up to any point in the lactation can be estimated.The procedure was similar for both variables.

Phenotypic Parameters
Table 1 gives an overview of the database used.As can be seen, we followed the performance of 27,749 does, with an average number of 12 TD and 2.06 lactations per animal.These were daughters of 16,243 goats and 941 sires.The results of the preliminary fixed effects model determined that all effects included in the model using both dependent variables were highly significant.Figure 1 shows the main phenotypic trends for FPY and SCS throughout lactation length and kidding number.

Phenotypic Parameters
Table 1 gives an overview of the database used.As can be seen, we followed the performance of 27,749 does, with an average number of 12 TD and 2.06 lactations per animal.These were daughters of 16,243 goats and 941 sires.The results of the preliminary fixed effects model determined that all effects included in the model using both dependent variables were highly significant.Figure 1 shows the main phenotypic trends for FPY and SCS throughout lactation length and kidding number.For both variables, there was a clear tendency to increase their level as the number of lactations increased, although this was more marked between the first and the other lactations for FPY.If the different lactations are compared, similar behavior can be observed in the three lactations for the SCC, so that as the lactation progresses, the SCC increases.The number of somatic cells increased rapidly until week 15, with a few weeks of delay compared to the productive peak of FPY.After that, it continued to grow more slowly until reaching the maximum peak at the end of lactation at almost week 40 (approx.10 months in lactation).In the case of FPY, there was an inverse, similar trend in all three lactations, but with a clearly lower magnitude in the case of the first lactation.This difference was probably due to the fact that the first kidding usually takes place when physical development is not yet complete (the animal has not yet reached adult size), so that during For both variables, there was a clear tendency to increase their level as the number of lactations increased, although this was more marked between the first and the other lactations for FPY.If the different lactations are compared, similar behavior can be observed in the three lactations for the SCC, so that as the lactation progresses, the SCC increases.The number of somatic cells increased rapidly until week 15, with a few weeks of delay compared to the productive peak of FPY.After that, it continued to grow more slowly until reaching the maximum peak at the end of lactation at almost week 40 (approx.10 months in lactation).In the case of FPY, there was an inverse, similar trend in all three lactations, but with a clearly lower magnitude in the case of the first lactation.This difference was probably due to the fact that the first kidding usually takes place when physical development is not yet complete (the animal has not yet reached adult size), so that during the first lactation, many nutrients still have to be used for the animal's growth as well as for milk production.

Repeatability Model
The estimated results obtained according to this model are presented in Table 2.The variation between lactations of the variance components can be observed, for the first three lactations analyzed, both for the SCS trait and the FPY.Thus, heritability ranged between 0.272 and 0.279 for the SCS trait and between 0.099 and 0.138 for FPY.Table 3 presents the genetic correlations between the EGVs for each kidding for each variable and between both variables within each lactation.When using the RM, the correlations between the EGVs for each variable in the three lactations showed a positive pattern of relationships within each variable.The EGV Xi of each variable in the three kiddings also showed a pattern of positive relationships within the variable but of moderate magnitude, while there was a slight antagonism between SCS and FPY (ranging from −0.307 for the first kidding to −0.592 for the third one).
Table 3. Correlations between the expected genetic values estimated according to the repeatability model in the Florida goat breed: above the diagonal, correlations between kiddings for SCS; below the diagonal, correlations between kiddings for FPY and on the diagonal, correlations for each kidding between SCS and FPY.This correlation showed that animals with a higher predisposition to present mastitis (higher genetic values for SCS) have a lower potential for milk production (FPY, in our case).It therefore follows that genetic selection for higher resistance to mastitis should result in a better response for milk production.

Random Regression Model
The evolution of the (co)variance components, heritabilities and within-trait genetic correlations of both variables along the DIM wi scale and each kidding is represented in Figure 2. Genetic variances and h 2 estimates showed similar shapes, with a decreasing trend as DIM w increases and a slight increase in the middle of lactation.In the same way, the within-trait genetic correlations across the DIM wi scale and number of kiddings were different from unity.The genetic correlations between FPY and SCS throughout the lactation and in the three different parities are presented in Figure 3. Their estimated values were negative in all cases and oscillated between −0.53; −0.17, −0.44; −0.24 and −0.40; −0.27, for the first, second and third kidding, respectively.These results indicate that an expected low genetic somatic cell count is associated with a high genetic production potential.The correlations were higher in the first phase, when the animal initiates milk production, as well as in the final phase, when after a lactation, the udder epithelium has suffered the stress of milk production for a whole lactation.In the middle of lactation another peak occurs, coinciding with the productive stress after the lactation peak (although it is less evident in the case of the first lactation).The genetic correlations between FPY and SCS throughout the lactation and in the three different parities are presented in Figure 3. Their estimated values were negative in all cases and oscillated between −0.53; −0.17, −0.44; −0.24 and −0.40; −0.27, for the first, second and third kidding, respectively.These results indicate that an expected low genetic somatic cell count is associated with a high genetic production potential.The correlations were higher in the first phase, when the animal initiates milk production, as well as in the final phase, when after a lactation, the udder epithelium has suffered the stress of milk production for a whole lactation.In the middle of lactation another peak occurs, coinciding with the productive stress after the lactation peak (although it is less evident in the case of the first lactation).
The results shown in Figure 2 indicate that the variables studied should not be considered as expressing the same traits throughout the lactation: in other words, there is genetic variation in the type of response.To illustrate this variation, it was considered that, in practical terms, the objective can be interpreted as higher cumulative FPY volume at 240 days of lactation and lower average SCS.To obtain these values, the corresponding EGVs were estimated following the suggestion by Jamrozik and Schaeffer [19], whose general formulation is presented above in the Materials and Methods section.Using this method, the best 500 animals were selected and the results are presented in Figure 4.
It should be noted that, among the best selected animals, there were contrasting forms of response in both variables, with animals that can be described as robust maintaining the same level of EGVs throughout the lactation, while others, which can be described as plastic, showing a tendency to decrease or increase in the same trajectory.This illustrates one of the key advantages of using RRM, in that it allows us to identify the animals that present special attributes in their response pattern throughout lactation.Our work demonstrates that there is genetic variation in the shape of the lactation curve for the traits studied, so that the results of the RM may be biased by assuming the opposite.The correlations between the EGVs estimated by both models corroborate this view, as can be seen in Figure 5.The results shown in Figure 2 indicate that the variables studied should not be considered as expressing the same traits throughout the lactation: in other words, there is genetic variation in the type of response.To illustrate this variation, it was considered that, in practical terms, the objective can be interpreted as higher cumulative FPY volume at 240 days of lactation and lower average SCS.To obtain these values, the corresponding EGVs were estimated following the suggestion by Jamrozik and Schaeffer [19], whose general formulation is presented above in the Materials and Methods section.Using this method, the best 500 animals were selected and the results are presented in Figure 4.It should be noted that, among the best selected animals, there were contrasting forms of response in both variables, with animals that can be described as robust maintaining the same level of EGVs throughout the lactation, while others, which can be described as plastic, showing a tendency to decrease or increase in the same trajectory.This  For SCS, there was a slight increase until the 10th week of lactation for each of the kiddings, followed by a decreasing response until the 34th week of lactation.This correlation of less than 0.8 in all cases implies that selection for increased resistance to mastitis (lower somatic cell count) has a different direct response (much lower in the case of the RM) depending on the assessment method used (RM or RRM), especially in the early stages of lactation.In the same way, the indirect correlated response for resistance to mastitis, if the selection criterion is FPY, will be lower when using RM than when using RRM, since the genetic correlations between both traits is always lower in RM (Tables 2 and 3).
For FPY, the pattern was similar, but with correlations below 0.53 (Figure 5), so that selection for milk production appears to also have an indirect response for resistance to mastitis, although in this case of a lesser magnitude than in the reverse selection.Also, in this case, the indirect response is higher when selecting RRM.

Discussion
The issue of udder health, and particularly subclinical and clinical mastitis and its simplest and most economical indicator, somatic cell count (SCC), has prompted numerous studies to incorporate these indicators into selection programs.The positive results published in Holstein cattle in the USA and Canada [20], with reductions ranging from around 10% to 40% in SCC levels, have provided very encouraging support for work in small ruminants.Results in dairy goats (Scholtens et al. [15] in New Zealand breeds; Arnal et al. [21] in French Alpine goats) and dairy sheep (Serrano et al. [22] in the Manchega breed; Riggio et al. [23]; Tolone et al. [24] in the Valle del Belice breed) also point in the same direction.The usefulness of recording and analyzing SCC has the major advantage that it is a very simple character which can be incorporated in bulk into the periodic dairy control systems in these species.However, the level of SCC might not distinguish the threshold from which infection could be identified, so, there has been a certain degree of controversy about the periodicity in which SCC should be evaluated: following the average lactation, when using the well-known RM model or on each test-day (TD), when RRM is used.A detailed study published by Windig et al. [25] comparing both models has suggested that more benefits could be obtained with the use of TD.For SCS, there was a slight increase until the 10th week of lactation for each of the kiddings, followed by a decreasing response until the 34th week of lactation.This correlation of less than 0.8 in all cases implies that selection for increased resistance to mastitis (lower somatic cell count) has a different direct response (much lower in the case of the RM) depending on the assessment method used (RM or RRM), especially in the early stages of lactation.In the same way, the indirect correlated response for resistance to mastitis, if the selection criterion is FPY, will be lower when using RM than when using RRM, since the genetic correlations between both traits is always lower in RM (Tables 2 and 3).
For FPY, the pattern was similar, but with correlations below 0.53 (Figure 5), so that selection for milk production appears to also have an indirect response for resistance to mastitis, although in this case of a lesser magnitude than in the reverse selection.Also, in this case, the indirect response is higher when selecting RRM.

Discussion
The issue of udder health, and particularly subclinical and clinical mastitis and its simplest and most economical indicator, somatic cell count (SCC), has prompted numerous studies to incorporate these indicators into selection programs.The positive results published in Holstein cattle in the USA and Canada [20], with reductions ranging from around 10% to 40% in SCC levels, have provided very encouraging support for work in small ruminants.Results in dairy goats (Scholtens et al. [15] in New Zealand breeds; Arnal et al. [21] in French Alpine goats) and dairy sheep (Serrano et al. [22] in the Manchega breed; Riggio et al. [23]; Tolone et al. [24] in the Valle del Belice breed) also point in the same direction.The usefulness of recording and analyzing SCC has the major advantage that it is a very simple character which can be incorporated in bulk into the periodic dairy control systems in these species.However, the level of SCC might not distinguish the threshold from which infection could be identified, so, there has been a certain degree of controversy about the periodicity in which SCC should be evaluated: following the average lactation, when using the well-known RM model or on each test-day (TD), when RRM is used.A detailed study published by Windig et al. [25] comparing both models has suggested that more benefits could be obtained with the use of TD.
It can be seen how the values of h 2 lie within the published range of the literature available for goat species, with estimates of h 2 for SCS of between 0.18 and 0.32 (Rupp et al. [26] and Arnal et al. [27] in Alpine and Saanen breeds; Scholtens et al. [15] in New Zealand breeds).It should be noted, however, that previous references were not homogeneous with regard to how to quantify the level of SCC, and a wide range of models were also used.
Estimates of h 2 for total fat plus protein production have shown ranges of h 2 between 0.04 and 0.40, as published by numerous authors, in goat species (Rupp et al. [26] in Alpine and Saanen breeds; García-Peniche et al. [28] and Castañeda-Bustos et al. [29] in US breeds; Scholtens et al 2019 [15] in New Zealand breeds; Arnal et al 2019 [27] in Alpine and Saanen breeds).The results of this study for both variables were within the lower levels of the above references, while the same trends can be seen in the wide variability of genetic origin that could be exploited in a selection program and the negative but weak relationships between SCS and milk components.However, this negative correlation does not coincide with the farmer's viewpoint, at least on those farms with higher production levels, which could lead us to think that production levels in this breed are not yet high enough for the production stress suffered by the animals to compromise their immune response capacity (i.e., cause a lower resistance to udder infection).However, Figure 1 shows that the inverse relationship between the two variables at the phenotypic level follows the same tendency as at the genetic level.An increase in the number of somatic cells might be related to natural desquamation by simple repletion of the udder in highly productive animals, without there being a defensive inflammatory reaction caused by bacterial growth and without being related to a higher propensity of genetic origin.
In accordance with Stinchcombe and Kirkpatrick [30], traits which are expressed across a temporal trajectory (or various environmental levels) are classified as 'function value traits'.These traits should be studied using RRM, such as the one employed in this study.In line with the findings of Zumbach et al. [31], Apodaca-Sarabia et al. [32] and Arnal et al. [21], our study concurs that a Legendre polynomial of order 2 provides the best fit to data on goat milk production under diverse conditions.
The few references available on SCS present values of h 2 from 0.12 to 0.25 [32], while Arnal et al. [21] found an h 2 of between 0.10 and 0.155; both studies showed slight increases in the middle of lactation.Applying the same statistical model, Arnal et al. [21] presented values of h 2 from 0.14 to 0.23 for total fat and protein production in the three kiddings of Alpine goats in France, with both variables showing higher values of h 2 in the middle of lactation.Similar trends were indicated by Zumbach et al. [31], but with higher levels of h 2 , from 0.28 to 0.47, while Oliveira et al. [33] showed results of h 2 from 0.40 to 0.50 for fat and protein percentages, respectively, recorded along the lactation scale.
Regarding the genetic correlations (r g ) between different points in the lactation, the trends for each variable were very similar for each kidding, with a decreasing response with distance between lactation stages.This same response was presented by Zumbach et al. [31] (r g ≈ 0.48 to 0.80 for total protein and r g ≈ 0.37 to 0.73 for total fat) in different goat breeds in Germany; Arnal et al. [21], with (r g ≈ 0.50 to 0.85) for fat and total protein in Alpine goats from France, while Oliveira et al. [33] presented the same pattern (r g ≈ 0.50) for percentages of fat and protein in Saanen and Alpine goats from Brazil.The problem is even more complex if we consider the r g between SCS and FPY, which were median but negative throughout lactation.References of this type are scarce, however: Rupp et al. [26] presented values of r g ≈ −0.13 to − 0.20, which are very close to those in Figure 2, following the same trend as in dairy sheep [23].In general, the relationships between milk components in goat species tend to be negative, according to several references presented in tables in the article by Scholtens et al. [15].
In general, references to r g between milk production traits and SCS are scarce and have only been reported for fat yield and protein yield separately, using a repeatability model.Rupp et al. [26] observed negative r g between SCS and FY and PY (−0.13 and −0.04, respectively) in the French Saanen breed.Conversely, other studies have reported positive correlations between these traits, varying between 0.06 and 0.23 (Rupp et al. [26] in Alpine goats; Valencia-Posadas et al. [34] in Nigerian Dwarf goats).

Figure 1 .
Figure 1.Least-squares means of somatic cell score (SCS) and fat plus protein yield (FPY) across lactation length (in weeks) and number of kiddings in the Florida breed.

Figure 1 .
Figure 1.Least-squares means of somatic cell score (SCS) and fat plus protein yield (FPY) across lactation length (in weeks) and number of kiddings in the Florida breed.

Figure 2 .
Figure 2. Evolution of variance components, heritabilities and genetic correlations across lactations and each parity for somatic cell score (SCS) and fat plus protein yield (FPY) using a random regression model in the Florida breed.

Figure 2 .
Figure 2. Evolution of variance components, heritabilities and genetic correlations across lactations and each parity for somatic cell score (SCS) and fat plus protein yield (FPY) using a random regression model in the Florida breed.

Figure 3 .
Figure 3. Evolution of the genetic correlations between somatic cell score (SCS) and fat plus protein yield (FPY) throughout lactation and in each parity using a random regression model in the Florida breed.

Figure 3 .
Figure 3. Evolution of the genetic correlations between somatic cell score (SCS) and fat plus protein yield (FPY) throughout lactation and in each parity using a random regression model in the Florida breed.Dairy 2024, 5, x FOR PEER REVIEW 9 of 13

Figure 4 .
Figure 4. Frequency distribution of estimated genetic values (EGVs) for somatic cell score (SCS) and fat plus protein yield (FPY) up to 240 days of lactation and variation in shapes of the lactation curves of the 500 best animals, using a random regression model in the Florida breed.

Figure 4 .
Figure 4. Frequency distribution of estimated genetic values (EGVs) for somatic cell score (SCS) and fat plus protein yield (FPY) up to 240 days of lactation and variation in shapes of the lactation curves of the 500 best animals, using a random regression model in the Florida breed.

Figure 5 .
Figure 5. Correlations between the estimated genetic values (EGVs) of somatic cell score (SCS) and fat plus protein yield (FPY) estimated by the repeatability and random regression models throughout lactation for each kidding in the Florida breed.

Figure 5 .
Figure 5. Correlations between the estimated genetic values (EGVs) of somatic cell score (SCS) and fat plus protein yield (FPY) estimated by the repeatability and random regression models throughout lactation for each kidding in the Florida breed.

Table 1 .
Descriptive statistics of the data analyzed (standard deviation in parenthesis) in the Florida goat breed.

Table 1 .
Descriptive statistics of the data analyzed (standard deviation in parenthesis) in the Florida goat breed.

Table 2 .
Estimation of genetic parameter components for somatic cell score (SCS) and total daily fat plus protein yield (FPY) in the Florida goat breed using a bivariate repeatability model.