Short communication: Genetic parameter estimates for caprine arthritis encephalitis in dairy goats.

Caprine arthritis encephalitis (CAE) is a chronic disease caused by a retrovirus from the Lentivirus genus. No effective vaccines or treatments exist, and therefore genetic selection for CAE resistance might be a feasible alternative. To our best knowledge, no other studies have investigated the genetic architecture of CAE resistance in dairy goats. In this context, this study was designed to estimate genetic parameters for CAE infection in Alpine and Saanen goats using a Bayesian threshold model. A total of 542 adult goats (and >3-generation pedigree), which were group-housed in a population with high CAE prevalence, were tested based on a serological infection assessment test (negative = 1 or positive = 2) and used for this study. Genetic parameters were estimated using the BLUPF90 family programs. There was considerable genetic variability for CAE resistance, and pedigree-based heritability was significantly different from zero (0.026 < heritability < 0.128). Our findings indicate that the prevalence of CAE in goat herds can be reduced or eliminated through direct genetic selection for CAE resistance in addition to proper management strategies.

Triggering the activation of Activin A type II receptor in human adipose stem cells towards tenogenic commitment using mechanomagnetic stimulation.

Stem cell therapies hold potential to stimulate tendon regeneration and homeostasis, which is maintained in response to the native mechanical environment. Activins are members of the mechano-responsive TGF-ß superfamily that participates in the regulation of several downstream biological processes. Mechanosensitive membrane receptors such as activin can be activated in different types of stem cells via magnetic nanoparticles (MNPs) through remote magnetic actuation resulting in cell differentiation. In this work, we target the Activin receptor type IIA (ActRIIA) in human adipose stem cells (hASCs), using anti-ActRIIA functionalized MNPs, externally activated through a oscillating magnetic bioreactor. Upon activation, the phosphorylation of Smad2/3 is induced allowing translocation of the complex to the nucleus, regulating tenogenic transcriptional responses. Our study demonstrates the potential remote activation of MNPs tagged hASCs to trigger the Activin receptor leading to tenogenic differentiation. These results may provide insights toward tendon regeneration therapies.

Exploring Stem Cells and Inflammation in Tendon Repair and Regeneration.

Tendon injuries are frequent and are responsible for substantial morbidity both in sports and in the workplace. Despite the endogenous mechanisms of tendon repair and regeneration, tendon healing upon injury is slow and often insufficient to restore complete biomechanics functionality.Inflammation has a pivotal role in tendon healing and failed healing responses contribute to the progression of tendinopathies. However, the molecular and cellular mechanisms involved are poorly understood requiring further insights.During inflammation, bioactive molecules such as cytokines secreted locally at the injury site, influence resident stem cells that contribute as modulatory agents over the niche towards homeostasis, holding great promise as therapeutic agents for tendon pathological conditions associated to unresolved inflammation and failed healing.This review overviews the role of cytokines and resident cells, focusing on the participation of tendon stem cell population in inflammation and tendon healing upon injury and their potential action in resolution of pathological conditions.

Using the Small Ruminant Nutrition System to develop and evaluate an alternative approach to estimating the dry matter intake of goats when accounting for ruminal fiber stratification.

The first objective of this research was to assess the ability of the Small Ruminant Nutrition System (SRNS) mechanistic model to predict metabolizable energy intake (MEI) and milk yield (MY) when using a heterogeneous fiber pool scenario (GnG1), compared with a traditional, homogeneous scenario (G1). The second objective was to evaluate an alternative approach to estimating the dry matter intake (DMI) of goats to be used in the SRNS model. The GnG1 scenario considers an age-dependent fractional transference rate for fiber particles from the first ruminal fiber pool (raft) to an escapable pool (λr), and that this second ruminal fiber pool (i.e., escapable pool) follows an age-independent fractional escape rate for fiber particles (ke). Scenario G1 adopted only a single fractional passage rate (kp). All parameters were estimated individually by using equations published in the literature, except for 2 passage rate equations in the G1 scenario: 1 developed with sheep data (G1-S) and another developed with goat data (G1-G). The alternative approach to estimating DMI was based on an optimization process using a series of dietary constraints. The DMI, MEI, and MY estimated for the GnG1 and G1 scenarios were compared with the results of an independent dataset (n=327) that contained information regarding DMI, MEI, MY, and milk and dietary compositions. The evaluation of the scenarios was performed using the coefficient of determination (R(2)) between the observed and predicted values, mean bias (MB), bias correction factor (Cb), and concordance correlation coefficient. The MEI estimated by the GnG1 scenario yielded precise and accurate values (R(2) = 082; MB = 0.21 Mcal/d; Cb = 0.98) similar to those of the G1-S (R(2) = 0.85; MB = 0.10 Mcal/d; Cb=0.99) and G1-G (R(2) = 0.84; MB = 0.18 Mcal/d; Cb = 0.98) scenarios. The results were also similar for the MY, but a substantial MB was found as follows: GnG1 (R(2) = 0.74; MB = 0.70 kg/d; Cb = 0.79), G1-S (R(2) = 0.71; MB = 0.58 kg/d(1); Cb = 0.85) and G1-G (R(2) = 0.71; MB = 0.65 kg/d; Cb = 0.82). The alternative approach for DMI prediction provided better results with the G1-G scenario (R(2)=0.88; MB = -71.67 g/d; Cb = 0.98). We concluded that the GnG1 scenario is valid within mechanistic models such as the SRNS and that the alternative approach for estimating DMI is reasonable and can be used in diet formulations for goats.

Dynamics of chromosomal evolution in the genus Hypsiboas (Anura: Hylidae).

Hylidae is one of the most species-rich families of anurans, and 40% of representatives in this group occur in Brazil. In spite of such remarkable diversity, little is known about this family and its taxonomical and systematic features. Most hylids have 2n = 24, even though most of the cytogenetic data are mainly obtained based on the conventional chromosomal staining and are available for only 16% of Hypsiboas species, a genus accounting for about 10% of the hylid diversity. In this study, cytogenetic data of distinct species and populations of Hypsiboas were analyzed, and the evolutionary dynamics of chromosomal macro- and microstructure of these amphibians were discussed. Contrary to the conservativeness of 2n = 24, this genus is characterized by a high variation of chromosomal morphology with as much as 8 karyotype patterns. Differences in the number and location of nucleolus organizer regions and C-bands allowed the identification of geographical variants within nominal species and cytotaxonomical chromosomal markers. Comparative analyses revealed a strong phylogeographic relationship between chromosomal patterns in this group.

Genetic evaluation of milk yield in Alpine goats for the first four lactations using random regression models.

Random regression models have been used in evaluating test-day milk yield, providing accurate estimates of genetic values in animals. However, herd evaluation with only information from the first lactation may not be the best option from an economic perspective. Other factors should be taken into account, particularly other lactations. Our objective in this study was to analyze the genetic divergence between the first four lactations of Alpine goats. The RENPED software was used to perform descriptive statistics, check for errors in pedigree, recode the data, and for Pearson's and Spearman's correlations. The WOMBAT software was used to estimate the variance components and predict the breeding values. The CALC software was adopted to calculate the percentage of coincidence between the ranking of the animals and the animals kept in common at each lactation evaluation. The results show that selection using only the first lactation in small herds with a low degree of technology can be employed as a palliative measure, in view of the difficulty in evaluating all lactations. However, the selection of breeding goats and the production of catalogues should not be based only on the first lactation, because the results demonstrate inversions in the classification of the best breeders when other lactations are analyzed.

Factors that influence the test day milk yield and composition.

The objective of this research was to study the factors that influence the test day milk yield (TDMY) and percentages of protein, fat, lactose, and total dry extract obtained on test day. We evaluated 110,732 weekly milk production records from 1496 goats and 19,271 monthly milk constituent records from 1245 Alpine and Saanen goats, which were collected from 1997 to 2010 in the goat sector at Universidade Federal de Viçosa. To ensure greater record reliability, only lactation data with kidding order between 1 to 6, type of kidding data including 0 to 3 kids, milk control years after 1997, and genetic groupings other than types 7 or 9 were considered, due to the relative lack of information recorded for some classes of these factors. Data in which the reported milk days were less than 7 or greater than 315 were also eliminated. Goats aged greater than 300 days at calving and those aged less than 6 years at control were considered in this study. Milk production was higher in the dry season in comparison to the rainy season. Genetic grouping did not influence all traits in both breedings. The TDMY tended to increase along with increasing age of the goats at kidding, while the opposite trend was observed relative to kidding order. Factors that significantly influenced all of the studied traits varied, and the factors that significantly influenced each trait were altered between the relationship of Alpine and Saanen breeds. Thus, the analysis of factors that influence traits to be evaluated in the herd under study is critical for defining the best evaluation model.

Genetic evaluation of Alpine goats using different milk control intervals.

The objective of this study was to compare the results of genetic evaluations by using different milk control intervals to reduce the cost of milk yield controls without harming the quality of genetic evaluation of the animals. We analyzed test day milk yield data from the Goat Sector of Universidade Federal de Viçosa. After editing and checking for errors in the database, there were 20,710 records of test day milk yield for the 667 first lactations of Alpine goats, constituting the complete file, with 7-day control intervals. Information on specific weeks was excluded from the complete file to create files with data on control intervals of 15, 21, and 28 days. The RENPED program was used to recode the pedigree and data files and correct pedigree errors; the WOMBAT program was used for genetic evaluations of the 4 files. The following comparison criteria of analysis results were used: logarithm of the function of the restricted maximum likelihood, length of the analyses in seconds, Pearson and Spearman correlations, and common elimination percentage among the areas below the regression curve of the genetic values of the animals. Overall, it is recommended that a 7-day interval among milk controls should be used in breeding programs and farms with a high technical level. Intervals of 14 and 21 days can achieve satisfactory results combined with a lower data collection cost for farms with an average-to-low technical level, less effective size, and genetic variability that depend on external genetic material for genetic improvement.

Random regression models using Legendre orthogonal polynomials to evaluate the milk production of Alpine goats.

The objective of this study was to identify the best random regression model using Legendre orthogonal polynomials to evaluate Alpine goats genetically and to estimate the parameters for test day milk yield. On the test day, we analyzed 20,710 records of milk yield of 667 goats from the Goat Sector of the Universidade Federal de Viçosa. The evaluated models had combinations of distinct fitting orders for polynomials (2-5), random genetic (1-7), and permanent environmental (1-7) fixed curves and a number of classes for residual variance (2, 4, 5, and 6). WOMBAT software was used for all genetic analyses. A random regression model using the best Legendre orthogonal polynomial for genetic evaluation of milk yield on the test day of Alpine goats considered a fixed curve of order 4, curve of genetic additive effects of order 2, curve of permanent environmental effects of order 7, and a minimum of 5 classes of residual variance because it was the most economical model among those that were equivalent to the complete model by the likelihood ratio test. Phenotypic variance and heritability were higher at the end of the lactation period, indicating that the length of lactation has more genetic components in relation to the production peak and persistence. It is very important that the evaluation utilizes the best combination of fixed, genetic additive and permanent environmental regressions, and number of classes of heterogeneous residual variance for genetic evaluation using random regression models, thereby enhancing the precision and accuracy of the estimates of parameters and prediction of genetic values.

Phylogeographic structure and karyotypic diversity of the Brazilian shrew mouse ( Blarinomys breviceps, Sigmodontinae) in the Atlantic Forest.

Blarinomys breviceps possesses cryptic and burrowing habits with poorly documented genetics and life history traits. Due to its rarity, only a few specimens and DNA sequences have been deposited in collections worldwide. Here, we present the most comprehensive cytogenetic and molecular characterization of this rare genus. Phylogenetic analyses based on partial cytochrome b sequences were performed, attempting to establish the relationships among individuals with distinct karyotypes along the geographic distribution of the genus in the Atlantic Forest. Classical and molecular cytogenetics, using banding patterns and FISH of telomeric and whole chromosome X-specific painting probes (obtained from the Akodontini Akodon cursor) were used to characterize and compare the chromosomal complements. Molecular phylogenetic analyses recovered 2 main geographically structured clades, northeastern and southeastern with pairwise sequence divergences among specimens varying between 4.9 and 8.4%. Eight distinct karyomorphs are described: (A) 2n = 52 (50A, XX), (B) 2n = 52 (48A, XY+2Bs), (C) 2n = 45 (42A, XY+1B), (D) 2n = 43 (37A, XX+4Bs), (E) 2n = 37 (34A, XY+1B), (F) 2n = 34 (32A, XX), (G) 2n = 31 (27A, XX+2Bs), (H) 2n = 28 (26A, XY), all with the same number of autosomal arms (FN(A) = 50). Variation of 0-4 supernumerary chromosomes (Bs) presenting heterogeneity in morphology and distribution of interstitial telomeric sequences (ITSs) is reported. ITSs are also found in some metacentric autosomes. The phylogeographic separation between 2 major lineages with high levels of genetic divergence, and the wide karyotypic diversity indicate that B. breviceps is a diverse group that warrants taxonomic re-evaluation.

Genetic and environmental factors that influence production and quality of milk of Alpine and Saanen goats.

Data from 1279 lactations of 783 Alpine and Saanen goats of the herd of our university in Minas Gerais, Brazil, were used to study environmental effects on and to estimate genetic parameters for milk production until 270 days of lactation (MP270) and for production and percentages of fat (PFAT and %FAT), protein (PPROT and %PROT), lactose (PLACT and %LACT), and total dry extract (PEXTR and %EXTR). Environmental effects were estimated by a statistical model that included contemporary group effect, type of kidding, genetic grouping, and kidding order. A multi-trait animal model with animal and permanent environment random effects was used to estimate genetic parameters and the significant environmental effects (fixed). Contemporary group influenced all traits; genetic grouping did not influence %LACT; type of kidding did not influence PFAT, %PROT or %LACT, and kidding order did not influence %FAT or %EXTR. Heritability and repeatability estimates were, respectively, 0.19 and 0.37 (MP270); 0.10 and 0.20 (PFAT); 0.12 and 0.24 (PPROT); 0.15 and 0.27 (PLACT); 0.13 and 0.24 (PEXTR); 0.21 and 0.34 (%FAT); 0.39 and 0.44 (%PROT); 0.17 and 0.29 (%LACT); 0.31 and 0.47 (%EXTR). Estimates of genetic correlations among MP270 and production of milk constituents were positive and high, but correlations between MP270 and %FAT, MP270 and %PROT, MP270 and %ESTR were moderate and negative. These heritability estimates show that satisfactory genetic gains can be obtained by selection, especially for milk constituents.

Energy requirements and efficiency of Alpine goats in early lactation.

Dairy goats may rely heavily on body fat and protein reserves in early lactation. Therefore, we aimed to determine the energy requirement and estimate the efficiency of utilization the nutrients of tissues mobilized in the first 8 weeks of lactation for milk production using the comparative slaughter technique. The average initial body mass of 51 multiparous goats was 57.19 ±â€¯8.38 kg and a body condition score of 3.0 ±â€¯0.5. Three goats were slaughtered at the beginning of the experiment to serve as baseline animals to estimate initial empty BW and initial body composition. We used a complete randomized design in which the factor was the day of lactation for slaughtering (the 7th, 14th, 21st, 28th, 35th, 42nd, 49th and 56th day), with six repeats, totalling 48 goats. No fasting before slaughtering. All animals received a single experimental diet. The efficiency of transferring energy from body reserves to milk was estimated using a multiple linear regression equation yielding a value of 0.76. The total energy stored in the empty body decreased over the eight lactation weeks, from 726.47 ±â€¯26.19 to 316.18 ±â€¯49.21 MJ, a 56.47% reduction, mainly because of a reduction in the energy from internal fat of 3.96 ±â€¯1.98 MJ/day. In conclusion, the net energy required for maintenance is 60 ±â€¯30 kJ/BW0.75 per day, and the net energy required for lactation decreases 70 ±â€¯30 kJ/day during the first eight lactation weeks.

Magnetic responsive materials modulate the inflammatory profile of IL-1ß conditioned tendon cells.

Tendinopathies represent half of all musculoskeletal injuries worldwide. Inflammatory events contribute to both tendon healing and to tendinopathy conditions but the cellular triggers leading to one or the other are unknown. In previous studies, we showed that magnetic field actuation modulates human tendon cells (hTDCs) behavior in pro-inflammatory environments, and that magnetic responsive membranes could positively influence inflammation responses in a rat ectopic model. Herein, we propose to investigate the potential synergistic action of the magnetic responsive membranes, made of a polymer blend of starch with polycaprolactone incorporating magnetic nanoparticles (magSPCL), and the actuation of pulsed electromagnetic field (PEMF): 5 Hz, 4mT of intensity and 50% of duty cycle, in IL-1ß-treated-hTDCs, and in the immunomodulatory response of macrophages. It was found that the expression of pro-inflammatory (TNFα, IL-6, IL-8, COX-2) and ECM remodeling (MMP-1,-2,-3) markers tend to decrease in cells cultured onto magSPCL membranes under PEMF, while the expression of TIMP-1 and anti-inflammatory genes (IL-4, IL-10) increases. Also, CD16++ and CD206+ macrophages were only found on magSPCL membranes with PEMF application. Magnetic responsive membranes show a modulatory effect on the inflammatory profile of hTDCs favoring anti-inflammatory cues which is also supported by the anti-inflammatory/repair markers expressed in macrophages. These results suggest that magnetic responsive magSPCL membranes can contribute for inflammation resolution acting on both resident cell populations and inflammatory cells, and thus significantly contribute to tendon regenerative strategies. Statement of significance Magnetically-assisted strategies have received great attention in recent years to remotely trigger and guide cell responses. Inflammation plays a key role in tendon healing but persistent pro-inflammatory molecules can contribute to tendon disorders, and therefore provide a therapeutic target for advanced treatments. We have previously reported that magnetic fields modulate the response of human tendon cells (hTDCs) conditioned to pro-inflammatory environments (IL-1ß-treated-hTDCs), and that magnetic responsive membranes positively influence immune responses. In the present work, we combined pulsed electromagnetic field (PEMF) and magnetic responsive membranes to guide the inflammatory profile of IL-1ß-treated-hTDCs and of macrophages. The results showed that the synergistic action of PEMF and magnetic membranes supports the applicability of magnetically actuated systems to regulate inflammatory events and stimulate tendon regeneration.

Chromosomal evolution in the Brazilian geckos of the genus Gymnodactylus (Squamata, Phyllodactylidae) from the biomes of Cerrado, Caatinga and Atlantic rain forest: evidence of Robertsonian fusion events and supernumerary chromosomes.

Chromosomes of the South American geckos Gymnodactylus amarali and G. geckoides from open and dry areas of the Cerrado and Caatinga biomes in Brazil, respectively, were studied for the first time, after conventional and AgNOR staining, CBG- and RBG-banding, and FISH with telomeric sequences. Comparative analyses between the karyotypes of open areas and the previously studied Atlantic forest species G. darwinii were also performed. The chromosomal polymorphisms detected in populations of G. amarali from the states of Goiás and Tocantins is the result of centric fusions (2n = 38, 39 and 40), suggesting a differentiation from a 2n = 40 ancestral karyotype and the presence of supernumerary chromosomes. The CBG- and RBG-banding patterns of the Bs are described. G. geckoides has 40 chromosomes with gradually decreasing sizes, but it is distinct from the 2n = 40 karyotypes of G. amarali and G. darwinii due to occurrence of pericentric inversions or centromere repositioning. NOR location seems to be a marker for Gymnodactylus, as G. amarali and G. geckoides share a medium-sized subtelocentric NOR-bearing pair, while G. darwinii has NORs at the secondary constriction of the long arm of pair 1. The comparative analyses indicate a non-random nature of the Robertsonian rearrangements in the genus Gymnodactylus.

Follicle and hormone dynamics in single versus double ovulating heifers.

Follicles > or = 5 mm were ablated at 4 day post-ovulation in heifers to induce a follicular wave, and prostaglandin F(2alpha) was given at day 6 to increase the incidence of double ovulations. Follicle diameters and plasma hormone concentrations were compared between single ovulators (n=12) and double ovulators (n=8). In double ovulators, the interval from follicle deviation to the peak of the pre-ovulatory LH surge was shorter (1.9+/-0.2 vs 2.5+/-0.2 days; P<0.02) and diameter of the largest pre-ovulatory follicle was smaller (12.2+/-0.5 vs 13.3+/-0.3 mm; P<0.02). The LH concentrations of the pre-ovulatory surge did not differ between single and double ovulators for 24 h on each side of the peak. When data were normalised to LH peak, the peak of the pre-ovulatory FSH and oestradiol (E(2)) surges occurred in synchrony with the peak of LH surge for both groups. Concentration of FSH for 24 h on each side of the peak showed a group effect (P<0.0001) from lower concentration in the double ovulators. A group-by-hour interaction (P<0.008) for E(2) reflected greater concentration in the double ovulators before and at the peak. Results indicated that two pre-ovulatory follicles resulted in an earlier and greater E(2) increase, leading to lower FSH concentration, an earlier LH surge, and ovulation at a smaller diameter. In conclusion, the difference in hormone concentrations during the pre-ovulatory period was an effect rather than a cause of double ovulations.

Human adipose tissue-derived tenomodulin positive subpopulation of stem cells: A promising source of tendon progenitor cells.

Cell-based therapies are of particular interest for tendon and ligament regeneration given the low regenerative potential of these tissues. Adipose tissue is an abundant source of stem cells, which may be employed for the healing of tendon lesions. However, human adult multipotent adipose-derived stem cells (hASCs) isolated from the stromal vascular fraction of adipose tissue originate highly heterogeneous cell populations that hinder their use in specific tissue-oriented applications. In this study, distinct subpopulations of hASCs were immunomagnetic separated and their tenogenic differentiation capacity evaluated in the presence of several growth factors (GFs), namely endothelial GF, basic-fibroblast GF, transforming GF-ß1 and platelet-derived GF-BB, which are well-known regulators of tendon development, growth and healing. Among the screened hASCs subpopulations, tenomodulin-positive cells were shown to be more promising for tenogenic applications and therefore this subpopulation was further studied, assessing tendon-related markers (scleraxis, tenomodulin, tenascin C and decorin) both at gene and protein level. Additionally, the ability for depositing collagen type I and III forming extracellular matrix structures were weekly assessed up to 28 days. The results obtained indicated that tenomodulin-positive cells exhibit phenotypical features of tendon progenitor cells and can be biochemically induced towards tenogenic lineage, demonstrating that this subset of hASCs can provide a reliable source of progenitor cells for therapies targeting tendon regeneration.

The ontogenetic allometry of body morphology and chemical composition in dairy goat wethers.

We studied the ontogenetic growth of goat wethers (castrated male goats) of the Saanen and Swiss Alpine breeds based on a large range of intraspecific body mass (BM). The body parts and the chemical constituents of the empty body were described by the allometric function by using BM and the empty body mass (EBM) as the predictors for morphological traits and chemical composition, respectively. We fitted the allometric scaling function by applying the SAS NLMIXED procedure, but to evaluate assumptions regarding variances in morphological and compositional traits, we combined the scaling function with homoscedastic (MOD1), and the heteroscedastic exponential (MOD2) and power-of-the-mean (MOD3) variance functions. We also predicted the ontogenetic growth by using the traditional log-log transformation and back-transformed results into the arithmetic scale (MOD4). We obtained predictions from MOD4 in the arithmetic scale by a two-step process, and evaluated MOD1, MOD2 and MOD3 by a model selection framework, and compared MOD4 with MOD1, MOD2 and MOD3 based on goodness-of-fit measures. Based on information criteria for model selection, heterogeneous variance functions were more likely to describe 10 over 36 traits with a low level of model selection uncertainty. One trait was predicted by averaging the MOD1 and MOD2 variance functions; and nine traits were better described by averaging the MOD2 and MOD3 variance functions. The predictions for other 16 traits were averaged from MOD1, MOD2 and MOD3. However, MOD4 better described 11 traits according to the goodness-of-fit measures. Depending on the variable being analyzed, the body parts and the chemical amounts exhibited the three types of allometric behavior with respect to BM and EBM, that is, positive, negative and isometric ontogenetic growth. Reference BMs, that is, 20, 27, 35 and 45 kg, were used to compute the net protein and energy requirements based on the first derivative of the scaling function, and the results were presented in reference to the EBM and EBM0.75. Both the net protein and energy requirements scaled to EBM0.75 increased from 20 to 45 kg of BM.

Combining different functions to describe milk, fat, and protein yield in goats using Bayesian multiple-trait random regression models.

We proposed multiple-trait random regression models (MTRRM) combining different functions to describe milk yield (MY) and fat (FP) and protein (PP) percentage in dairy goat genetic evaluation by using Bayesian inference. A total of 3,856 MY, FP, and PP test-day records, measured between 2000 and 2014, from 535 first lactations of Saanen and Alpine goats, including their cross, were used in this study. The initial analyses were performed using the following single-trait random regression models (STRRM): third- and fifth-order Legendre polynomials (Leg3 and Leg5), linear B-splines with 3 and 5 knots, the Ali and Schaeffer function (Ali), and Wilmink function. Heterogeneity of residual variances was modeled considering 3 classes. After the selection of the best STRRM to describe each trait on the basis of the deviance information criterion (DIC) and posterior model probabilities (PMP), the functions were combined to compose the MTRRM. All combined MTRRM presented lower DIC values and higher PMP, showing the superiority of these models when compared to other MTRRM based only on the same function assumed for all traits. Among the combined MTRRM, those considering Ali to describe MY and PP and Leg5 to describe FP (Ali_Leg5_Ali model) presented the best fit. From the Ali_Leg5_Ali model, heritability estimates over time for MY, FP. and PP ranged from 0.25 to 0.54, 0.27 to 0.48, and 0.35 to 0.51, respectively. Genetic correlation between MY and FP, MY and PP, and FP and PP ranged from -0.58 to 0.03, -0.46 to 0.12, and 0.37 to 0.64, respectively. We concluded that combining different functions under a MTRRM approach can be a plausible alternative for joint genetic evaluation of milk yield and milk constituents in goats.